În matematică , transformata Fourier în timp discret , adesea prescurtată ca DTFT (acronim al termenului englezesc Transformată Fourier în timp discret ), este o transformare care pornind de la un semnal discret oferă o descriere periodică a acesteia în domeniul frecvenței , similar cu transformă Fourier tradițional (definit pentru funcții continue).
Acesta este un caz particular al transformării zeta :
- {\ displaystyle X (z) = \ sum _ {n = - \ infty} ^ {\ infty} x [n] \, z ^ {- n}}
care se obține prin plasare {\ displaystyle z = e ^ {i \ omega} \,} . De cand {\ displaystyle | e ^ {i \ omega} | = 1 \,} , transformata Fourier în timp discret este evaluarea transformatei zeta pe cercul unității în planul complex .
Definiție
Dat fiind un set de numere întregi complexe {\ displaystyle x [n]} , cu {\ displaystyle n \ in \ mathbb {Z}} , transformata sa Fourier în timp discret este seria:
- {\ displaystyle X (\ omega) = \ sum _ {n = - \ infty} ^ {\ infty} x [n] \, e ^ {- i \ omega n}}
Transformarea inversă permite obținerea funcției originale pornind de la transformarea sa:
- {\ displaystyle x [n] = {\ frac {1} {2 \ pi}} \ int _ {2 \ pi} X (\ omega) \ cdot e ^ {i \ omega n} d \ omega}
Relația cu eșantionarea
Transformata Fourier în timp discret joacă un rol important atunci când studiază semnalele eșantionate, adică semnalele în timp discret obținute dintr-un semnal în timp continuu, luând în considerare valoarea asumată în momente precise de timp, de obicei separate printr-un interval de timp fix {\ displaystyle T} . Procedura care permite obținerea unui semnal discret pornind de la unul continuu se numește eșantionare și este baza conversiei analog-digitale (ADC). Transformă o funcție continuă {\ displaystyle x (t)} în semnalul discret:
- {\ displaystyle x [n] \ {\ stackrel {\ text {def}} {=}} \ x (nT) \ qquad \ forall \, n \ in \ mathbb {Z}}
cu {\ displaystyle 1 / T} rata de eșantionare . Teorema de eșantionare plasează o limită a frecvenței maxime a semnalului continuu, care nu poate fi mai mare de {\ displaystyle 1 / (2T)} dacă doriți să evitați pierderea informațiilor (fenomen aliasing ). Transformarea discretă a timpului oferă o aproximare a transformatei Fourier {\ displaystyle {\ mathcal {F}}} :
- {\ displaystyle {\ mathcal {F}} (x (t)) (f) = X (f) = \ int _ {- \ infty} ^ {\ infty} x (t) \ cdot e ^ {- i2 \ pi ft} dt}
De fapt, luând în considerare formula de însumare Poisson , care arată cum se obține o însumare periodică a unei funcții {\ displaystyle X (f)} pornind de la mostrele unei funcții în timp continuu, avem:
- {\ displaystyle X_ {1 / T} (f) \ {\ stackrel {\ mathrm {def}} {=}} \ sum _ {k = - \ infty} ^ {\ infty} X \ left (fk / T \ dreapta) \ equiv \ sum _ {n = - \ infty} ^ {\ infty} \ underbrace {T \ cdot x (nT)} _ {x [n]} \ \ underbrace {e ^ {- i2 \ pi fTn} } _ {{\ mathcal {F}} \ {\ delta (t-nT) \}} \ = \ {\ mathcal {F}} \ left \ {\ sum _ {n = - \ infty} ^ {\ infty } x [n] \ cdot \ delta (t-nT) \ right \}}
unde este {\ displaystyle \ scriptstyle X_ {1 / T} (f) \,} include copii exacte ale {\ displaystyle \ scriptstyle X (f) \,} tradus de un multiplu de {\ displaystyle f_ {s}} și combinate prin adăugare. Pentru {\ displaystyle f_ {s}} termenul este suficient de mare {\ displaystyle K = 0} poate fi observat în regiune {\ displaystyle [-f_ {s} / 2, f_ {s} / 2]} , cu distorsiuni reduse sau deloc. O altă modalitate de a verifica acest lucru este următoarea:
- {\ displaystyle {\ begin {align} {\ mathcal {F}} \ left \ {\ sum _ {n = - \ infty} ^ {\ infty} T \ cdot x (nT) \ cdot \ delta (t-nT ) \ right \} & = {\ mathcal {F}} \ left \ {x (t) \ cdot T \ sum _ {n = - \ infty} ^ {\ infty} \ delta (t-nT) \ right \ } \\ & = X (f) * \ underbrace {{\ mathcal {F}} \ left \ {T \ sum _ {n = - \ infty} ^ {\ infty} \ delta (t-nT) \ right \ }} _ {\ sum _ {k = - \ infty} ^ {\ infty} \ delta (fk / T)} = \ sum _ {k = - \ infty} ^ {\ infty} X \ left (f- { \ frac {k} {T}} \ right) \ end {align}}}
Mai mult, prin calcularea transformatei Fourier inverse a ambilor membri ai ecuației anterioare, se obține pieptenul Dirac modulat:
- {\ displaystyle \ sum _ {n = - \ infty} ^ {\ infty} x [n] \ cdot \ delta (t-nT) = {\ mathcal {F}} ^ {- 1} \ left \ {X_ { 1 / T} (f) \ right \} \ {\ stackrel {\ mathrm {def}} {=}} \ int _ {- \ infty} ^ {\ infty} X_ {1 / T} (f) \ cdot e ^ {i2 \ pi ft} df.}
cu:
- {\ displaystyle x [n] = T \ int _ {\ frac {1} {T}} X_ {1 / T} (f) \ cdot e ^ {i2 \ pi fnT} df = {\ frac {1} { 2 \ pi}} \ int _ {2 \ pi} X (\ omega) \ cdot e ^ {i \ omega n} d \ omega}
Intrare periodică
Dacă succesiunea de intrare {\ displaystyle x [n]} este periodic cu perioada N este posibil să se extindă pieptenul Dirac în seria Fourier , obținând transformata Fourier discretă (DFT):
- {\ displaystyle {\ mathcal {F}} \ left \ {\ sum _ {n = - \ infty} ^ {\ infty} x [n] \ cdot \ delta (t-nT) \ right \} = {\ mathcal {F}} \ left \ {\ underbrace {\ sum _ {k = - \ infty} ^ {\ infty} X [k] \ cdot e ^ {i2 \ pi {\ frac {k} {NT}} t} } _ {Fourier \ series} \ right \} = \ underbrace {\ sum _ {k = - \ infty} ^ {\ infty} X [k] \ \ cdot \ \ delta \ left (f - {\ frac {k } {NT}} \ right)} _ {DTFT}}
Această relație arată că periodicitatea în timp face ca transformarea Fourier în timp discret să fie discontinuă. Cu toate acestea, formula integrală poate fi redusă la o sumă de N termeni:
- {\ displaystyle X [k] {\ stackrel {\ text {def}} {=}} \ {\ frac {1} {NT}} \ int _ {NT} \ left [\ sum _ {n = - \ infty } ^ {\ infty} x [n] \ cdot \ delta (t-nT) \ right] e ^ {- i2 \ pi {\ frac {k} {NT}} t} dt = {\ frac {1} { NT}} \ sum _ {n = - \ infty} ^ {\ infty} x [n] \ cdot \ int _ {NT} \ delta (t-nT) \ cdot e ^ {- i2 \ pi {\ frac { k} {NT}} t} dt}
- {\ displaystyle = {\ frac {1} {NT}} \ underbrace {\ sum _ {N} x [n] \ cdot e ^ {- i2 \ pi {\ frac {k} {N}} n}} _ {DFT} = {\ frac {1} {N}} \ underbrace {\ sum _ {N} x (nT) \ cdot e ^ {- i2 \ pi {\ frac {k} {N}} n}} _ {DFT}}
care este periodic în {\ displaystyle k} .
Eșantionare DTFT
Dacă transformata Fourier în timp discret este o funcție continuă , este adesea utilizată pentru a lua în considerare un număr arbitrar de eșantioane ale unui ciclu al funcției periodice {\ displaystyle X_ {1 / T}} :
- {\ displaystyle \ underbrace {X_ {1 / T} \ left ({\ frac {k} {NT}} \ right)} _ {X_ {k}} = \ sum _ {n = - \ infty} ^ {\ infty} x [n] \ cdot e ^ {- i2 \ pi {\ frac {kn} {N}}} = \ underbrace {\ sum _ {N} x_ {N} [n] \ cdot e ^ {- i2 \ pi {\ frac {kn} {N}}}} _ {DFT} \ qquad \ quad k = 0, \ dots, N-1}
unde este {\ displaystyle X_ {N}} este suma periodică:
- {\ displaystyle x_ {N} [n] \ {\ stackrel {\ text {def}} {=}} \ \ sum _ {m = - \ infty} ^ {\ infty} x [n-mN]}
Succesiunea {\ displaystyle x_ {N}} este inversul transformatei Fourier discrete . În acest fel, eșantionarea efectuată în acest mod implică faptul că transformarea inversă este periodică.
Pentru a evalua numeric un ciclu de {\ displaystyle x_ {N}} este necesară o succesiune {\ displaystyle x [n]} de lungime finită. În acest scop, o secvență este adesea tăiată prin intermediul unei funcții de fereastră de lungime adecvată. [1] [2]
Transformări notabile
Lasa-i sa fie {\ displaystyle n} domeniul de timp discret, {\ displaystyle \ omega} frecvența unghiulară (un număr real în {\ displaystyle (- \ pi, \ \ pi)} măsurat în radiani / eșantion), {\ displaystyle u [n]} pasul Heaviside în timp discret,{\ displaystyle \ operatorname {sinc} (t)} funcția sinc normalizată, {\ displaystyle \ delta (\ omega)} Delta Dirac , {\ displaystyle \ delta [n]} delta Kronecker ,{\ displaystyle \ operatorname {rect} (t)} funcția dreptunghiular :
- {\ displaystyle \ mathrm {rect} (t) = \ sqcap (t) = {\ begin {cases} 0 & {\ mbox {if}} | t |> {\ frac {1} {2}} \\ [ 3pt] {\ frac {1} {2}} & {\ mbox {if}} | t | = {\ frac {1} {2}} \\ [3pt] 1 & {\ mbox {if}} | t | <{\ frac {1} {2}} \ end {cases}}}
Și {\ displaystyle \ operatorname {tri} (t)} funcția triunghi:
- {\ displaystyle \ operatorname {tri} (t) = \ land (t) = {\ begin {cases} 1 + t; & - 1 \ leq t \ leq 0 \\ 1-t; & 0 <t \ leq 1 \ \ 0 & {\ mbox {altfel}} \ end {cases}}}
Domeniul timpului {\ displaystyle x [n] \,} | Domeniul de frecventa {\ displaystyle X (\ omega) \,} | Observații |
---|
{\ displaystyle \ delta [n] \!} | {\ displaystyle 1 \!} | |
{\ displaystyle \ delta [nM] \!} | {\ displaystyle e ^ {- i \ omega M} \!} | {\ displaystyle M} întreg |
{\ displaystyle \ sum _ {m = - \ infty} ^ {\ infty} \ delta [n-Mm] \!} | {\ displaystyle \ sum _ {m = - \ infty} ^ {\ infty} e ^ {- i \ omega Mm} = {\ frac {1} {M}} \ sum _ {k = - \ infty} ^ { \ infty} \ delta \ left ({\ frac {\ omega} {2 \ pi}} - {\ frac {k} {M}} \ right) \,} | {\ displaystyle M} întreg |
{\ displaystyle u [n] \!} | {\ displaystyle {\ frac {1} {1-e ^ {- i \ omega}}} + \ sum \ limits _ {k = - \ infty} ^ {\ infty} {\ pi \ delta (\ omega -2 \ pi k)} \!} | Termenul {\ displaystyle 1 / (1-e ^ {- i \ omega})} trebuie interpretat ca o distribuție . |
{\ displaystyle a ^ {n} u [n] \!} | {\ displaystyle {\ frac {1} {1-ae ^ {- i \ omega}}} \!} | {\ displaystyle 0 <| a | <1} |
{\ displaystyle e ^ {- ian} \!} | {\ displaystyle 2 \ pi \ delta (\ omega + a) \,} | {\ displaystyle a \ in \ mathbb {R}} |
{\ displaystyle \ cos (an) \!} | {\ displaystyle \ pi \ left [\ delta (\ omega -a) + \ delta (\ omega + a) \ right]} | {\ displaystyle a \ in \ mathbb {R}} |
{\ displaystyle \ sin (an) \!} | {\ displaystyle {\ frac {\ pi} {i}} \ left [\ delta (\ omega -a) - \ delta (\ omega + a) \ right]} | {\ displaystyle a \ in \ mathbb {R}} |
{\ displaystyle \ mathrm {rect} \ left [{(nM / 2) \ over M} \ right]} | {\ displaystyle {\ sin [\ omega (M + 1) / 2] \ over \ sin (\ omega / 2)} \, e ^ {- i \ omega M / 2} \!} | {\ displaystyle M} întreg |
{\ displaystyle \ operatorname {sinc} [(a + n)]} | {\ displaystyle e ^ {ia \ omega} \!} | {\ displaystyle a \ in \ mathbb {R}} |
{\ displaystyle W \ cdot \ operatorname {sinc} ^ {2} (Wn) \,} | {\ displaystyle \ operatorname {tri} \ left ({\ omega \ over 2 \ pi W} \ right)} | {\ displaystyle 0 <W \ leq 0.5 \ quad \ in \ mathbb {R}} |
{\ displaystyle W \ cdot \ operatorname {sinc} (Wn)} | {\ displaystyle \ operatorname {rect} \ left ({\ omega \ over 2 \ pi W} \ right)} | {\ displaystyle 0 <W \ leq 1 \ quad \ in \ mathbb {R}} |
{\ displaystyle {\ begin {cases} 0 & n = 0 \\ {\ frac {(-1) ^ {n}} {n}} și {\ mbox {în altă parte}} \ end {cases}}} | {\ displaystyle j \ omega} | Filtru diferențiator |
{\ displaystyle {\ frac {W} {(n + a)}} \ left \ {\ cos [\ pi W (n + a)] - \ operatorname {sinc} [W (n + a)] \ right \ }} | {\ displaystyle j \ omega \ cdot \ operatorname {rect} \ left ({\ omega \ over \ pi W} \ right) e ^ {ja \ omega}} | {\ displaystyle 0 <W \ leq 1 \ quad \ in \ mathbb {R} \ qquad a \ in \ mathbb {R}} |
{\ displaystyle {\ frac {1} {\ pi n ^ {2}}} [(- 1) ^ {n} -1] \!} | {\ displaystyle | \ omega | \!} | |
{\ displaystyle {\ begin {cases} 0; & n {\ mbox {even}} \\ {\ frac {2} {\ pi n}}; & n {\ mbox {odd}} \ end {cases}} } | {\ displaystyle {\ begin {cases} j & \ omega <0 \\ 0 & \ omega = 0 \\ - j & \ omega> 0 \ end {cases}}} | Transformarea lui Hilbert |
{\ displaystyle {\ frac {C (A + B)} {2 \ pi}} \ cdot \ operatorname {sinc} \ left [{\ frac {AB} {2 \ pi}} n \ right] \ cdot \ operatorname {sinc} \ left [{\ frac {A + B} {2 \ pi}} n \ right]} | | {\ displaystyle A, B \ in R \ quad B \ in \ mathbb {C}} |
Proprietate
Lasa-i sa fie {\ displaystyle * \!} convoluția discretă a două secvențe și {\ displaystyle x [n] ^ {*} \!} complexul conjugat al {\ displaystyle x [n] \!} .
Proprietate | Domeniul timpului {\ displaystyle x [n] \!} | Domeniul de frecventa {\ displaystyle X (\ omega) \!} | Observații |
---|
Linearitatea | {\ displaystyle ax [n] + de [n] \!} | {\ displaystyle aX (e ^ {i \ omega}) + bY (e ^ {i \ omega}) \!} | |
Traducerea timpului | {\ displaystyle x [nk] \!} | {\ displaystyle X (e ^ {i \ omega}) e ^ {- i \ omega k} \!} | {\ displaystyle k} întreg |
Traducerea în frecvență | {\ displaystyle x [n] e ^ {ian} \!} | {\ displaystyle X (e ^ {i (\ omega -a)}) \!} | {\ displaystyle a \ in \ mathbb {R}} |
Inversarea timpului | {\ displaystyle x [-n] \!} | {\ displaystyle X (e ^ {- i \ omega}) \!} | |
Conjugarea temporală | {\ displaystyle x [n] ^ {*} \!} | {\ displaystyle X (e ^ {- i \ omega}) ^ {*} \!} | |
Inversarea timpului și conjugarea | {\ displaystyle x [-n] ^ {*} \!} | {\ displaystyle X (e ^ {i \ omega}) ^ {*} \!} | |
Derivat | {\ displaystyle {\ frac {n} {i}} x [n] \!} | {\ displaystyle {\ frac {dX (e ^ {i \ omega})} {d \ omega}} \!} | |
Grâu integral | {\ displaystyle {\ frac {i} {n}} x [n] \!} | {\ displaystyle \ int _ {- \ pi} ^ {\ omega} X (e ^ {i \ vartheta}) d \ vartheta \!} | |
Convoluţie | {\ displaystyle x [n] * y [n] \!} | {\ displaystyle X (e ^ {i \ omega}) \ cdot Y (e ^ {i \ omega}) \!} | |
Multiplicare | {\ displaystyle x [n] \ cdot y [n] \!} | {\ displaystyle {\ frac {1} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} {X (e ^ {i \ vartheta}) \ cdot Y (e ^ {i (\ omega - \ vartheta)}) d \ vartheta} \!} | Convoluție periodică |
Corelarea încrucișată | {\ displaystyle \ rho _ {xy} [n] = x [-n] ^ {*} * y [n] \!} | {\ displaystyle R_ {xy} (\ omega) = X (e ^ {i \ omega}) ^ {*} \ cdot Y (e ^ {i \ omega}) \!} | |
Teorema lui Parseval | {\ displaystyle E = \ sum _ {n = - \ infty} ^ {\ infty} {x [n] \ cdot y ^ {*} [n]} \!} | {\ displaystyle E = {\ frac {1} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} {X (e ^ {i \ omega}) \ cdot Y ^ {*} (și ^ {i \ omega}) d \ omega} \!} | |
Transformarea poate fi, de asemenea, descompusă atât în părți reale, cât și în părți imaginare și, respectiv, în două funcții, pare și impar:
- {\ displaystyle X (e ^ {i \ omega}) = X_ {R} (e ^ {i \ omega}) + iX_ {I} (e ^ {i \ omega}) \ qquad X (e ^ {i \ omega}) = X_ {E} (e ^ {i \ omega}) + X_ {O} (e ^ {i \ omega}) \!}
Domeniul timpului {\ displaystyle x [n] \!} | Domeniul de frecventa {\ displaystyle X (e ^ {i \ omega}) \!} |
---|
{\ displaystyle x ^ {*} [n] \!} | {\ displaystyle X ^ {*} (e ^ {- i \ omega}) \!} |
{\ displaystyle x ^ {*} [- n] \!} | {\ displaystyle X ^ {*} (e ^ {i \ omega}) \!} |
Notă
Bibliografie
- Alan V. Oppenheim și Ronald W. Schafer, Discrete-Time Signal Processing , Ediția a II-a, Prentice Hall Signal Processing Series, 1999, ISBN 0-13-754920-2 .
- William McC. Siebert, circuite, semnale și sisteme , seria MIT de inginerie electrică și informatică. Cambridge, MA, MIT Press, 1986.
- Boaz Porat,Un curs de procesare digitală a semnalului , John Wiley și Sons, 1941, pp. 27-29 și 104-105, ISBN 0-471-14961-6 .
Elemente conexe