Transformată Fourier în timp discret - Discrete-time Fourier transform

În matematică , transformata Fourier în timp discret ( DTFT ) este o formă de analiză Fourier aplicabilă unei secvențe de valori.

DTFT este adesea folosit pentru a analiza eșantioane ale unei funcții continue. Termenul timp discret se referă la faptul că transformarea operează pe date discrete, de multe ori eșantioane al căror interval are unități de timp. Din probe uniform distanțate, produce o funcție de frecvență care este o însumare periodică a transformatei Fourier continue a funcției continue originale. În anumite condiții teoretice, descrise de teorema eșantionării , funcția continuă originală poate fi recuperată perfect din DTFT și astfel din eșantioanele discrete originale. DTFT în sine este o funcție continuă a frecvenței, dar eșantioane discrete ale acestuia pot fi ușor calculate prin intermediul transformatei Fourier discrete (DFT) (a se vedea § Eșantionarea DTFT ), care este de departe cea mai comună metodă de analiză Fourier modernă.

Ambele transformări sunt inversabile. DTFT invers este secvența de date eșantionată originală. DFT invers este o însumare periodică a secvenței originale. Transformare rapidă Fourier (FFT) este un algoritm de calcul un ciclu de DFT și inversa ei produce un ciclu inversului DFT.

Definiție

Transformata Fourier în timp discret a unei secvențe discrete de numere reale sau complexe x [ n ] , pentru toate numerele întregi n , este o serie Fourier , care produce o funcție periodică a unei variabile de frecvență. Când variabila de frecvență, ω, are unități normalizate de radiani / probă , periodicitatea este , iar seria Fourier este :

 

 

 

 

( Ec. 1 )

Utilitatea acestei funcții de domeniu de frecvență este înrădăcinată în formula de însumare Poisson . Fie X ( f ) transformata Fourier a oricărei funcții, x ( t ) , ale cărei probe la un anumit interval T ( secunde ) sunt egale (sau proporționale) cu secvența x [ n ] , adică Tx ( nT ) = x [ n ] . Atunci funcția periodică reprezentată de seria Fourier este o însumare periodică a lui X ( f ) în termeni de frecvență f în hertz ( cicluri / sec ) :

 

 

 

 

( Ec. 2 )

Fig 1. Reprezentarea unei transformate Fourier (stânga sus) și a însumării sale periodice (DTFT) în colțul din stânga jos. Colțul din dreapta jos prezintă mostre ale DTFT care sunt calculate de o transformată Fourier discretă (DFT).

Numărul întreg k are unități de cicluri / eșantion , iar 1 / T este rata de eșantionare, f s ( eșantioane / sec ). Deci X 1 / T ( f ) cuprinde copii exacte ale lui X ( f ) care sunt deplasate cu multipli de f s hertz și combinate prin adunare. Pentru f s suficient de mare , termenul k = 0 poate fi observat în regiunea [- f s / 2, f s / 2 ] cu distorsiuni mici sau deloc ( aliasing ) față de ceilalți termeni. În Fig.1, extremitățile distribuției în colțul din stânga sus sunt mascate prin aliasing în sumarea periodică (stânga jos).

De asemenea, observăm că e - i2πfTn este transformata Fourier a lui δ ( t - nT ) . Prin urmare, o definiție alternativă a DTFT este:

 

 

 

 

( Ec. 3 )

Funcția de pieptene Dirac modulat este o abstracție matematică uneori denumită eșantionare prin impuls .

Transformarea inversă

O operație care recuperează secvența de date discretă din funcția DTFT se numește DTFT invers . De exemplu, transformata Fourier continuă inversă a ambelor părți ale ecuației 3 produce secvența sub forma unei funcții de pieptene Dirac modulate:

Cu toate acestea, menționând că X 1 / T ( f ) este periodică, toate informațiile necesare sunt conținute în orice interval de lungime 1 / T . In ambele Eq.1 si Eq.2 , a însumările peste n sunt o serie Fourier , cu coeficienți x [ n ] . Formulele standard pentru coeficienții Fourier sunt, de asemenea, transformatele inverse:

 

 

 

 

( Ec. 4 )

Date periodice

Când secvența de date de intrare x [ n ] este N -periodică , Eq.2 poate fi redusă prin calcul la o transformată Fourier discretă (DFT), deoarece :

  • Toate informațiile disponibile sunt conținute în N eșantioane.
  • X 1 / T ( f ) converge la zero peste tot, cu excepția multiplilor întregi de 1 / ( NT ) , cunoscuți cafrecvențe armonice . La aceste frecvențe, DTFT diverg la diferite rate de frecvență. Și aceste rate sunt date de DFT-ul unui ciclu alsecvenței x [ n ] .
  • DTFT este periodic, deci numărul maxim de amplitudini armonice unice este (1 / T ) / (1 / ( NT )) = N

Coeficienții DFT sunt dați de :

    iar DTFT este :
     

Înlocuirea acestei expresii în formula de transformare inversă confirmă :

( toate numerele întregi )

cum era de așteptat. DFT invers în linia de mai sus este uneori denumită o serie de Fourier discret (DFS).

Eșantionarea DTFT

Când DTFT este continuu, o practică obișnuită este de a calcula un număr arbitrar de eșantioane ( N ) dintr-un ciclu al funcției periodice X 1 / T : 

unde este o însumare periodică :

    (vezi seria Fourier discretă )

Secvența este invers DFT. Astfel, eșantionarea noastră a DTFT face ca transformarea inversă să devină periodică. Matricea de | X k | 2 valori sunt cunoscute ca periodogramă , iar parametrul N se numește NFFT în funcția Matlab cu același nume.

Pentru a evalua un ciclu numeric, avem nevoie de o secvență de lungime finită x [ n ] . De exemplu, o secvență lungă ar putea fi trunchiată de o funcție de fereastră de lungime L, rezultând trei cazuri demne de menționat special. Pentru simplitatea noțională, luați în considerare valorile x [ n ] de mai jos pentru a reprezenta valorile modificate de funcția fereastră.

Caz: Decimarea frecvenței. L = NI , pentru un număr întreg I (de obicei 6 sau 8)

Un ciclu de reduce la o însumare a I segmente de lungime N . DFT are apoi diferite nume, cum ar fi :

  • fereastră-prezum FFT
  • Greutate, suprapunere, adăugare (WOLA)
  • polifază DFT
  • banc filtru polifazic
  • fereastră cu blocuri multiple și aliasing în timp .

Amintiți-vă că decimarea datelor eșantionate într-un domeniu (timp sau frecvență) produce suprapuneri (uneori cunoscute sub numele de aliasing ) în celălalt și invers. Comparativ cu un DFT cu lungime L , însumarea / suprapunerea determină decimarea frecvenței, lăsând doar probele DTFT cel mai puțin afectate de scurgeri spectrale . Aceasta este de obicei o prioritate atunci când implementați un filtru FFT (canalizator). Cu o funcție convențională de fereastră de lungime L , pierderea prin scoici ar fi inacceptabilă. Deci, ferestrele multi-bloc sunt create folosind instrumentele de proiectare a filtrelor FIR . Profilul lor de frecvență este plat în cel mai înalt punct și cade rapid în punctul mediu dintre restul de probe DTFT. Cu cât valoarea parametrului I este mai mare, cu atât este mai bună performanța potențială.

Caz: L = N +1 .

Când o funcție de fereastră simetrică, cu lungimea L ( ) este trunchiată cu 1 coeficient, se numește periodică sau uniformă DFT . Trunchierea afectează DTFT. O DFT a probelor de secvențe trunchiate DTFT la intervale de frecvență de 1 / N . Pentru a preleva probe la aceleași frecvențe, pentru comparație, DFT este calculat pentru un ciclu al însumării periodice,

Fig 2. DFT de e i2πn / 8 pentru L = 64 și N = 256
Fig 3. DFT de e i2πn / 8 pentru L = 64 și N = 64

Caz: interpolare de frecvență. LN

În acest caz, DFT se simplifică într-o formă mai familiară:

Pentru a profita de un algoritm rapid de transformare Fourier pentru calcularea DFT, însumarea se efectuează de obicei pe toți N termenii, chiar dacă N - L dintre ei sunt zerouri. Prin urmare, cazul L < N este adesea denumit zero-padding .

Scurgerea spectrală, care crește pe măsură ce L scade, este în detrimentul anumitor valori importante ale performanței, cum ar fi rezoluția componentelor de frecvență multiple și cantitatea de zgomot măsurată de fiecare probă DTFT. Dar aceste lucruri nu contează întotdeauna, de exemplu atunci când secvența x [ n ] este o sinusoidă fără zgomot (sau o constantă), modelată de o funcție de fereastră. Apoi, este o practică obișnuită să folosiți zero-padding pentru a afișa grafic și a compara modelele detaliate de scurgere ale funcțiilor ferestrei. Pentru a ilustra acest lucru pentru o fereastră dreptunghiulară, luați în considerare secvența:

și

Figurile 2 și 3 sunt graficele cu mărimea a două DFT de dimensiuni diferite, așa cum este indicat în etichetele lor. În ambele cazuri, componenta dominantă este la frecvența semnalului: f = 1/8 = 0,125 . De asemenea, vizibil în Fig 2 este modelul de scurgere spectrală a ferestrei dreptunghiulare L = 64 . Iluzia din Fig 3 este un rezultat al eșantionării DTFT doar la trecerea la zero. Mai degrabă decât DTFT-ul unei secvențe de lungime finită, dă impresia unei secvențe sinusoidale infinit de lungi. Factorii care contribuie la iluzie sunt utilizarea unei ferestre dreptunghiulare și alegerea unei frecvențe (1/8 = 8/64) cu exact 8 cicluri (un număr întreg) la 64 de probe. O fereastră Hann ar produce un rezultat similar, cu excepția faptului că vârful ar fi lărgit la 3 eșantioane (a se vedea fereastra DFT-chiar Hann ).

Convoluţie

Teorema convoluție pentru secvențe este :

Un caz special important este convoluția circulară a secvențelor x și y definite de unde este o însumare periodică. Natura cu frecvență discretă înseamnă că produsul cu funcția continuă este, de asemenea, discret, ceea ce duce la o simplificare considerabilă a transformatei inverse :

Pentru secvențele x și y a căror durată diferită de zero este mai mică sau egală cu N , o simplificare finală este :

Semnificația acestui rezultat este explicată la algoritmii de convoluție circulară și de convoluție rapidă .

Proprietăți de simetrie

Când părțile reale și imaginare ale unei funcții complexe sunt descompuse în părțile lor pare și impare , există patru componente, notate mai jos prin indicii RE, RO, IE și IO. Și există o mapare unu-la-unu între cele patru componente ale unei funcții de timp complexe și cele patru componente ale transformării sale frecvente complexe :

Din aceasta, apar diferite relații, de exemplu :

  • Transformarea unei funcții cu valoare reală ( x RE + x RO ) este funcția simetrică uniformă X RE + i X IO . În schimb, o transformare simetrică uniformă implică un domeniu de timp cu valoare reală.
  • Transformarea unei funcții cu valoare imaginară ( i x IE + i x IO ) este funcția simetrică impar X RO + i X IE , iar inversul este adevărat.
  • Transformarea unei funcții simetrice pare ( x RE + i x IO ) este funcția cu valoare reală X RE + X RO , iar inversul este adevărat.
  • Transformarea unei funcții impar-simetrice ( x RO + i x IE ) este funcția cu valoare imaginară i X IE + i X IO , iar inversul este adevărat.

Relația cu transformata Z.

este o serie Fourier care poate fi exprimată și în termeni ai transformatei Z bilaterale . Adică:

unde notația distinge transformata Z de transformata Fourier. Prin urmare, putem exprima și o porțiune a transformatei Z în termenii transformatei Fourier:

Rețineți că atunci când parametrul T se modifică, termenii de rămân o separare constantă , iar lățimea lor se ridică în sus sau în jos. Termenii lui X 1 / T ( f ) rămân o lățime constantă, iar separarea lor 1 / T se ridică în sus sau în jos.

Tabelul transformatelor Fourier în timp discret

Unele perechi comune de transformare sunt prezentate în tabelul de mai jos. Se aplică următoarea notație:

  • este un număr real care reprezintă frecvența unghiulară continuă (în radiani pe eșantion). ( este în cicluri / sec și este în sec / eșantion.) În toate cazurile din tabel, DTFT este 2π-periodic (în ).
  • desemnează o funcție definită pe .
  • desemnează o funcție definită pe și zero în altă parte. Atunci:
  • este funcția delta Dirac
  • este funcția sinc normalizată
  • este funcția dreptunghi
  • este funcția triunghi
  • n este un număr întreg care reprezintă domeniul în timp discret (în eșantioane)
  • este funcția pasului de unitate de timp discret
  • este delta Kronecker
Domeniul de timp
x [ n ]
Domeniul de frecvență
X (ω)
Observații Referinţă
întreg

    nui M     chiar M

întreg

Termenul trebuie să fie interpretată ca o distribuție în sensul unei valoare principale Cauchy în jurul său poli la .
    -π <a <π

numar real

număr real cu
număr real cu
numere întregi și
numere reale cu
număr real ,
funcționează ca un filtru diferențiator
numere reale cu
Transformarea lui Hilbert
Semnal trapezoidal.svg numere reale complexe

Proprietăți

Acest tabel prezintă câteva operații matematice în domeniul timpului și efectele corespunzătoare în domeniul frecvenței.

Proprietate Domeniul de timp
x [ n ]
Domeniul de frecventa
Observații Referinţă
Linearitatea numere complexe
Inversarea timpului / Inversarea frecvenței
Conjugarea în timp
Inversarea timpului și conjugarea
Partea reală în timp
O parte imaginară în timp
Parte reală în frecvență
Partea imaginară în frecvență
Schimbare în timp / Modulare în frecvență număr întreg k
Schimbare în frecvență / Modulare în timp numar real
Decimare   întreg
Expansiunea timpului întreg
Derivată în frecvență
Integrarea în frecvență
Diferențierea în timp
Suma în timp
Convoluție în timp / Multiplicare în frecvență
Multiplicarea în timp / Convoluția în frecvență Convoluție periodică
Corelarea încrucișată
Teorema lui Parseval

Vezi si

Note

Citări de pagină

Referințe

Lecturi suplimentare

  • Porat, Boaz (1996). Un curs de procesare a semnalului digital . John Wiley și Sons. pp. 27-29 și 104-105. ISBN 0-471-14961-6.
  • Siebert, William M. (1986). Circuite, semnale și sisteme . Seria MIT de inginerie electrică și informatică. Cambridge, MA: MIT Press. ISBN 0262690950.
  • Lyon, Richard G. (2010). Înțelegerea procesării digitale a semnalului (ediția a 3-a). Prentice Hall. ISBN 978-0137027415.