Fourier' teisenduse ja diskreetse Fourier' teisenduse vaheline suhe. Vasak veerg: Pidev funktsioon (ülemine) ja tema Fourier' teisendus (all). Kesk-vasak veerg: Originaalfunktsiooni perioodiline summeerimine (ülemine). Fourier' teisendus (alumine) on null, v.a diskreetsete punktide kohtadel. Pöördteisendus on sinusoidide summa – Fourier' rida. Kesk-parem veerg: Originaalfunktsioon on diskreeditud (ülemine). Selle Fourier' teisendus (alumine) on originaalteisenduse perioodiline summa. Parem veerg: DFT (alumine) arvutab diskreetseid proove pidevast diskreetse aja Fourier' teisendusest. Pöörd-DFT (ülemine) on originaalproovide perioodiline summeerimine. Kiire DFT arvutab ühe tsükli DFT-d ja pöörd-DFT-d ühes pöörd-DFT tsüklis.
Fourier' teisenduse kujutis (ülemine) ja selle perioodiline summa (DTFT) (alumine)
**Diskreetne Fourier' teisendus** (lühend DFT inglise keele sõnadest *discrete Fourier transform*) on pideva Fourier' teisenduse vaste digiteeritud (ajas diskreeditud ja nivoos kvanditud) funktsioonide ja signaalide jaoks. Teatud tingimustel on DFT tulemused vastavuses pideva Fourier' teisenduse tulemustega. Kuid paljudel juhtudel on tulemused siiski oluliselt erinevad, mistõttu DFT tulemuste ülekandmisel analoogsignaalidele tuleb olla ettevaatlik.

DFT on kõige olulisem diskreetne teisendus, mida kasutatakse Fourier' analüüsi teostamiseks paljudes praktilistes rakendustes. Signaalitöötluses on iga signaal funktsioon, millest võetakse proove (lugemeid signaali väärtustest) diskreetimissagedustel. Analüüsi objektideks võivad olla näiteks helilaine, raadiosignaal või muutuvad temperatuurinäidud. Pilditöötluses võivad proovid olla pikslite väärtused rasterkujutise igas reas või veerus. DFT-d kasutatakse ka selleks, et lahendada efektiivselt diferentsiaalvõrrandeid ja teha teisi toiminguid, nagu konvolutsioon ja suurte arvude korrutamine.

Kuna DFT tegeleb piiratud andmemahtudega, saab seda rakendada, arvutades numbriliste algoritmide abil või kasutades sihtotstarbelist riistvara. Need rakendused kasutavad tavaliselt Fourier' kiirteisendust ehk **kiiret Fourier' teisendust** (*Fast Fourier Transform*, FFT), mis põhineb teisendamiseks vajalike arvutuste mahu vähendamisel teatud (korduvate) tulemuste (korduva) ärakasutamise teel või mõne spetsiifilise vaheteisenduse kasutamise teel.

Mõiste
------

*N* kompleksarvude jada {\displaystyle x\_{0},x\_{1},\ldots ,x\_{N-1}} teisendatakse N-perioodiliseks kompleksarvude jadaks:

{\displaystyle X\_{k}\ {\stackrel {\text{def}}{=}}\ \sum \_{n=0}^{N-1}x\_{n}\cdot e^{-2\pi ikn/N},\quad k\in \mathbb {Z} \,} (täisarvud)   

Perioodilisuse tõttu on tegelikult arvutatud domeen k [0, N – 1]. See on alati nii, kui DFT-d on rakendatud kiire Fourier' teisenduse kaudu. Muud domeenid on [-A / 2, N / 2-1] (*N* paaris) ja [- (N – 1) / 2, (N – 1) / 2] (*N* paaritu), kui vasak ja parem pool DFT väljundis on vahetatud.

Diskreetne Fourier' teisendus on tavaliselt tähistatud tähega {\displaystyle {\mathcal {F}}}, kus {\displaystyle \mathbf {X} ={\mathcal {F}}\left\{\mathbf {x} \right\}} või {\displaystyle {\mathcal {F}}\left(\mathbf {x} \right)} või {\displaystyle {\mathcal {F}}\mathbf {x} }.

Definitsioonivõrrandit võib tõlgendada või tuletada mitmel viisil, näiteks:

* See kirjeldab täielikult diskreetse aja Fourier' N-perioodilist järjestust, mis sisaldab ainult diskreetse sageduse komponente.
* Samuti tagab see ühtlaste vahekaugusega proove sisalduva pideva DFT antud piiratud pikkusega.
* See on sisendandmete ristkorrelatsioon, {\displaystyle X\_{n}} ja kompleksne sinusoid sagedusel k/N. Seega see toimib nagu sobitatud filter sellele sagedusele.
* See on diskreetne analoog Fourier' rea koefitsientide leidmise valemile:

{\displaystyle x\_{n}={\frac {1}{N}}\sum \_{k=0}^{N-1}X\_{k}\cdot e^{2\pi ikn/N},\quad n\in \mathbb {Z} \,}

mis on ka **N**-perioodiline. {\displaystyle \scriptstyle n\ \in \ [0,\ N-1]} domeenis on see definitsiooni valemi pöördteisendus. Selle tõlgenduse järgi on iga {\displaystyle X\_{k}} kompleksne number, mis kodeerib funktsiooni {\displaystyle x\_{n}} komplekssesse siinuselisse komponenti  {\displaystyle (e^{2\pi ikn/N})} nii amplituudi kui ka faasi. Siinuse sagedus on k/N. Amplituud ja faas vastavad võrranditele:
{\displaystyle |X\_{k}|/N={\sqrt {\operatorname {Re} (X\_{k})^{2}+\operatorname {Im} (X\_{k})^{2}}}/N}
{\displaystyle \arg(X\_{k})=\operatorname {atan2} {\big (}\operatorname {Im} (X\_{k}),\operatorname {Re} (X\_{k}){\big )}=-i\operatorname {ln} \left({\frac {X\_{k}}{|X\_{k}|}}\right),}
kus atan2 on arkustangensi funktsioon kahe argumendi jaoks.
Euleri valemi abil saab DFT valemi teisendada trigonomeetriliseks vormiks, mida kasutatakse inseneride arvutlustes ja infotehnoloogias:

**Fourier teisendus:**
{\displaystyle X\_{k}=\sum \_{n=0}^{N-1}x\_{n}\cdot (\cos(-2\pi k{\frac {n}{N}})+i\sin(-2\pi k{\frac {n}{N}})),\quad k\in \mathbb {Z} \,}

**Pöördteisendus:**
{\displaystyle x\_{n}={\frac {1}{N}}\sum \_{k=0}^{N-1}X\_{k}\cdot (\cos(2\pi k{\frac {n}{N}})+i\sin(2\pi k{\frac {n}{N}})),\quad n\in \mathbb {Z} \,}

kus

N on proovide arv
{\displaystyle n} on hetkel töödeldava proovi number (0, ... {\displaystyle N-1})
{\displaystyle x\_{n}} on signaali väärtus hetkel {\displaystyle n}
{\displaystyle k} on sagedus (0 Hz kuni {\displaystyle N-1} Hz) töötlemise hetkel
{\displaystyle X\_{k}} on sagedusega {\displaystyle k} signaali (amplituudi ja faasi väljendav kompleksarv.

Omadused
--------

### Täielikkus

Diskreetne Fourier' teisendus on pööratav lineaarne teisendus

{\displaystyle {\mathcal {F}}\colon \mathbb {C} ^{N}\to \mathbb {C} ^{N}}
kus {\displaystyle \mathbb {C} } on kompleksarvude hulk. See tähendab, et iga *N* > 0, *N*-mõõtmeline kompleksvektor omab DFT-d ja PDFT-d (pöördteisendust), mis on omakorda N-mõõtmelised kompleksvektorid.

### Ortogonaalsus

Vektorid {\displaystyle u\_{k}=\left[e^{{\frac {2\pi i}{N}}kn}\;|\;n=0,1,\ldots ,N-1\right]^{T}} moodustavad ortogonaalse baasi üle N-mõõtmeliste kompleksvektorite hulga:

{\displaystyle u\_{k}^{T}u\_{k'}^{\*}=\sum \_{n=0}^{N-1}\left(e^{{\frac {2\pi i}{N}}kn}\right)\left(e^{{\frac {2\pi i}{N}}(-k')n}\right)=\sum \_{n=0}^{N-1}e^{{\frac {2\pi i}{N}}(k-k')n}=N~\delta \_{kk'}}
kus{\displaystyle ~\delta \_{kk'}} on Kronecker delta. (Viimane samm on triviaalsete {\displaystyle k=k'} summa, kus 1+1+⋅⋅⋅=*N*, ja selline geomeetriline jada, millest saab summeerides saada kokku nulli.) Seda ortogonaalsuse tingimust saab kasutada selleks, et tuletada pöördteisenduse valemit DFT definitsioonist.

### Planchereli teoreem ja Parsevali teoreem

Kui *X**k* ja *Y**k* on vastavalt *x**n* ja *y**n* teisendused, siis ütleb Parsevali teoreem, et:

{\displaystyle \sum \_{n=0}^{N-1}x\_{n}y\_{n}^{\*}={\frac {1}{N}}\sum \_{k=0}^{N-1}X\_{k}Y\_{k}^{\*}}
kus tärn tähendab kompleksset konjugatsiooni. Parsevali teoreemi selline erijuht, kus *x**n* = *y**n*, on Planchereli teoreem, mis väidab:

{\displaystyle \sum \_{n=0}^{N-1}|x\_{n}|^{2}={\frac {1}{N}}\sum \_{k=0}^{N-1}|X\_{k}|^{2}.}
### Perioodilisus

Perioodilisust võib näha otse definitsioonist:

{\displaystyle X\_{k+N}\ {\stackrel {\mathrm {def} }{=}}\ \sum \_{n=0}^{N-1}x\_{n}e^{-{\frac {2\pi i}{N}}(k+N)n}=\sum \_{n=0}^{N-1}x\_{n}e^{-{\frac {2\pi i}{N}}kn}\underbrace {e^{-2\pi in}} \_{1}=\sum \_{n=0}^{N-1}x\_{n}e^{-{\frac {2\pi i}{N}}kn}=X\_{k}.}
Samas on näha, et pöördteisenduse võrrand viib perioodilisele laiendamisele.

### Nihe

Korrutades {\displaystyle x\_{n}} lineaarse faasiga {\displaystyle e^{{\frac {2\pi i}{N}}nm}} mingi täisarvu *m* järgi vastab nihkele väljundis {\displaystyle X\_{k}}: {\displaystyle X\_{k}}, mis asendatakse {\displaystyle X\_{k-m}}, kus indeks on tõlgendatud kui absoluutväärtus *N* (perioodiliselt). Samamoodi nihe sisendis {\displaystyle x\_{n}} vastab väljundi {\displaystyle X\_{k}} korrutamisele lineaarse faasiga. Matemaatiliselt, kui {\displaystyle \{x\_{n}\}} tähistab vektorit **x,** siis

kui {\displaystyle {\mathcal {F}}(\{x\_{n}\})\_{k}=X\_{k}}
siis {\displaystyle {\mathcal {F}}(\{x\_{n}\cdot e^{{\frac {2\pi i}{N}}nm}\})\_{k}=X\_{k-m}}
ja {\displaystyle {\mathcal {F}}(\{x\_{n-m}\})\_{k}=X\_{k}\cdot e^{-{\frac {2\pi i}{N}}km}}
### Ringkonvolutsiooni teoreem ja ristkorrelatsiooni teoreem

Konvolutsiooni teoreem diskreetse Fourier' teisenduse jaoks näitab, et kahe lõpmatu jada konvolutsiooni võib saada kahe üksiku teisenduse tulemuste pöördteisendusest. Kui järjestused on piiratud pikkustega N, siis tekib oluline lihtsustamine. Kui vaadelda DTF-d ja pöörd-DTF-d, võib seda kirjutada järgmiselt:

{\displaystyle {\mathcal {F}}^{-1}\left\{\mathbf {X\cdot Y} \right\}\_{n}\ =\sum \_{l=0}^{N-1}x\_{l}\cdot (y\_{N})\_{n-l}\ \ {\stackrel {\mathrm {def} }{=}}\ \ (\mathbf {x\*y\_{N}} )\_{n}\ ,}
mis on konvolutsioon {\displaystyle \mathbf {x} } ja {\displaystyle \mathbf {y} } jadade vahel pikendatud perioodilise liitmisega:

{\displaystyle (\mathbf {y\_{N}} )\_{n}\ {\stackrel {\mathrm {def} }{=}}\ \sum \_{p=-\infty }^{\infty }y\_{(n-pN)}=y\_{n(modN)}.\,}
Ristkorrelatsioon {\displaystyle \mathbf {x} }  ja  {\displaystyle \mathbf {y\_{N}} }  vahel on selline**:**

{\displaystyle {\mathcal {F}}^{-1}\left\{\mathbf {X^{\*}\cdot Y} \right\}\_{n}=\sum \_{l=0}^{N-1}x\_{l}^{\*}\cdot (y\_{N})\_{n+l}\ \ {\stackrel {\mathrm {def} }{=}}\ \ (\mathbf {x\star y\_{N}} )\_{n}\ .}
### Konvolutsiooni teoreemi duaalsus

Võib ka tõestada, et**:**

{\displaystyle {\mathcal {F}}\left\{\mathbf {x\cdot y} \right\}\_{k}\ {\stackrel {\mathrm {def} }{=}}\sum \_{n=0}^{N-1}x\_{n}\cdot y\_{n}\cdot e^{-{\frac {2\pi i}{N}}kn}={\frac {1}{N}}(\mathbf {X\*Y\_{N}} )\_{k},\,} mis on ringikujuline konvolutsioon {\displaystyle \mathbf {X} } ja {\displaystyle \mathbf {Y} } vahel.
### Trigonomeetriline interpolatsiooni polünoom

{\displaystyle p(t)={\frac {1}{N}}\left[X\_{0}+X\_{1}e^{2\pi it}+\cdots +X\_{N/2-1}e^{(N/2-1)2\pi it}+X\_{N/2}\cos(N\pi t)+X\_{N/2+1}e^{(-N/2+1)2\pi it}+\cdots +X\_{N-1}e^{-2\pi it}\right]} paaris *N*-i jaoks,
{\displaystyle p(t)={\frac {1}{N}}\left[X\_{0}+X\_{1}e^{2\pi it}+\cdots +X\_{\lfloor N/2\rfloor }e^{\lfloor N/2\rfloor 2\pi it}+X\_{\lfloor N/2\rfloor +1}e^{-\lfloor N/2\rfloor 2\pi it}+\cdots +X\_{N-1}e^{-2\pi it}\right]} paaritu *N*-i jaoks,
kus koefitsiendid *X**k* on eespool oleva *x**n* DFT ja rahuldavad interpoleerimist {\displaystyle p(n/N)=x\_{n}} {\displaystyle n=0,\ldots ,N-1} jaoks.

Paaris *N-i* puhul tuleb tähele panna, et Nyquist komponenti {\displaystyle {\frac {X\_{N/2}}{N}}\cos(N\pi t)} käsitletakse eraldi.

### Ühtne DFT

Teise võimalusena saab DFT-d väljendada DFT maatriksi abil ehk Vandermonde maatriksina, mille võttis kasutusele Sylvester 1867. aastal.

{\displaystyle \mathbf {F} ={\begin{bmatrix}\omega \_{N}^{0\cdot 0}&\omega \_{N}^{0\cdot 1}&\ldots &\omega \_{N}^{0\cdot (N-1)}\\\omega \_{N}^{1\cdot 0}&\omega \_{N}^{1\cdot 1}&\ldots &\omega \_{N}^{1\cdot (N-1)}\\\vdots &\vdots &\ddots &\vdots \\\omega \_{N}^{(N-1)\cdot 0}&\omega \_{N}^{(N-1)\cdot 1}&\ldots &\omega \_{N}^{(N-1)\cdot (N-1)}\\\end{bmatrix}}}
kus

{\displaystyle \omega \_{N}=e^{-2\pi i/N}\,}
on primitiivne *N*-is ühe juur.

Pöördteisendus on selle maatriksi pöördmaatriks:

{\displaystyle \mathbf {F} ^{-1}={\frac {1}{N}}\mathbf {F} ^{\*}}
Ühtse normaliseerimise konstantidega {\displaystyle 1/{\sqrt {N}}} saab DFT ühtseks teisenduseks, mis on defineeritud ühtse maatriksi abil:

{\displaystyle \mathbf {U} =\mathbf {F} /{\sqrt {N}}}
{\displaystyle \mathbf {U} ^{-1}=\mathbf {U} ^{\*}}
{\displaystyle |\det(\mathbf {U} )|=1}
kus *det()* on maatriksi determinant. Determinant on omaväärtuste tulemus, mis on kas {\displaystyle \pm 1} või {\displaystyle \pm i}. Ortogonaalsust väljendab nüüd ortonormaalsuse tingimus.

{\displaystyle \sum \_{m=0}^{N-1}U\_{km}U\_{mn}^{\*}=\delta \_{kn}}
Kui **X** on ühtne vektori x DFT, siis

{\displaystyle X\_{k}=\sum \_{n=0}^{N-1}U\_{kn}x\_{n}}
ja Planchereli teoreem on väljendatud niimoodi:

{\displaystyle \sum \_{n=0}^{N-1}x\_{n}y\_{n}^{\*}=\sum \_{k=0}^{N-1}X\_{k}Y\_{k}^{\*}}
Erijuhul {\displaystyle \mathbf {x} =\mathbf {y} } on see Parsevali teoreem

{\displaystyle \sum \_{n=0}^{N-1}|x\_{n}|^{2}=\sum \_{k=0}^{N-1}|X\_{k}|^{2}}
### Pöördvõrdelise DFT avaldamine DFT kaudu

Kasulik DFT omadus on see, et pöördteisendust saab kergesti väljendada tavalise DFT kaudu, kasutades allolevaid nippe:

Esiteks saame arvutada pöördvõrdeline DFT, pöörates vastupidi kõik peale ühe sisendi.

{\displaystyle {\mathcal {F}}^{-1}(\{x\_{n}\})={\mathcal {F}}(\{x\_{N-n}\})/N}
Teiseks võib ka ümber pöörata sisendid ja väljundid:

{\displaystyle {\mathcal {F}}^{-1}(\mathbf {x} )={\mathcal {F}}(\mathbf {x} ^{\*})^{\*}/N}
Viimaseks võime vahetada omavahel reaal- ja imaginaarsignaali osad:

{\displaystyle {\mathcal {F}}^{-1}(\mathbf {x} )={\textrm {swap}}({\mathcal {F}}({\textrm {swap}}(\mathbf {x} )))/N}
Kus *swap*({\displaystyle x\_{n}}) on {\displaystyle x\_{n}}, milles on reaal- ja imaginaar-osad omavahel niimoodi vahetatud: {\displaystyle x\_{n}=a+bi} *swap*({\displaystyle x\_{n}}) on {\displaystyle b+ai}. Samas, *swap*({\displaystyle x\_{n}}) on võrdne {\displaystyle ix\_{n}^{\*}}.
### Omaväärtused ja omavektorid

DFT maatriksi omaväärtused on lihtsad ja hästi tuntud, kusjuures omavektorid on keerulised ja neid pole veel täiesti uuritud.

Kui DFT ühtne vorm {\displaystyle \mathbf {U} } pikkusega *N* on selline*:*

{\displaystyle \mathbf {U} \_{m,n}={\frac {1}{\sqrt {N}}}\omega \_{N}^{(m-1)(n-1)}={\frac {1}{\sqrt {N}}}e^{-{\frac {2\pi i}{N}}(m-1)(n-1)}.}
See maatriks vastab maatriksite polünomiaalsele võrrandile:

{\displaystyle \mathbf {U} ^{4}=\mathbf {I} .}
See tähendab, et omaväärtused \lambda võrrandile on:

{\displaystyle \lambda ^{4}=1.}
Samas on {\displaystyle \mathbf {U} } omaväärtused neljandad ühe juured: \lambda on +1, −1, +*i*, või −*i*.

Kuna on olemas ainult neli omaväärtust, siis need omavad mitmesust, mis annab omavektorite numbri igale omaväärtusele.

Paljususe probleemi lahendasid McClellan ja Parks, samas ekvivalentse probleemi lahendas varem Gauss. Mitmesus sõltub absoluutväärtusest *N* ja on näidatud järgmises tabelis:

| suurus *N* | λ = +1 | λ = −1 | λ = -*i* | λ = +*i* |
| --- | --- | --- | --- | --- |
| 4*m* | *m* + 1 | *m* | *m* | *m* − 1 |
| 4*m* + 1 | *m* + 1 | *m* | *m* | *m* |
| 4*m* + 2 | *m* + 1 | *m* + 1 | *m* | *m* |
| 4*m* + 3 | *m* + 1 | *m* + 1 | *m* + 1 | *m* |

### Määramatuse printsiip

Kui juhuslik suurus *Xk* on piiratud:

{\displaystyle \sum \_{n=0}^{N-1}|X\_{n}|^{2}=1~,}
siis 

{\displaystyle P\_{n}=|X\_{n}|^{2}}
võib esindada diskreetse tõenäosuse *n*-i massi-funktsiooni koos tõenäosuse massi-funktsiooni funktsiooniga, mis on saadud transformeeritud muutujast,

{\displaystyle Q\_{m}=N|x\_{m}|^{2}~.}
Pidevate funktsioonide {\displaystyle P(x)} ja {\displaystyle Q(k)} jaoks väidab Heisenbergi määramatuse printsiip, et

{\displaystyle D\_{0}(X)D\_{0}(x)\geq {\frac {1}{16\pi ^{2}}}}
kus {\displaystyle D\_{0}(X)} ja {\displaystyle D\_{0}(x)} on {\displaystyle |X|^{2}} ja {\displaystyle |x|^{2}} erinevused, koos võrdsusega, mis on saavutatud sobivalt normaliseeritud Gaussi jaotusest.

Hirschmanni entroopialine ebakindlus omab kasulikku analoogi ka DFT puhul. Hirschmani määramatuse printsiip on väljendatud Shannoni entroopia poolest kahe tõenäosusfunktsiooni abil.

Diskreetsel juhul on Shannoni entroopiad 

{\displaystyle H(X)=-\sum \_{n=0}^{N-1}P\_{n}\ln P\_{n}}
ja

{\displaystyle H(x)=-\sum \_{m=0}^{N-1}Q\_{m}\ln Q\_{m}~,}
ja entroopialine määramatuse printsiip muutub

{\displaystyle H(X)+H(x)\geq \ln(N)~.}
### Reaalse sisendiga DFT

Kui {\displaystyle x\_{0},\ldots ,x\_{N-1}} on reaalarvud, nagu nad on tavalises olukorras, siis DFT allub sümmeetriale:

{\displaystyle X\_{N-k}\equiv X\_{-k}=X\_{k}^{\*},}  kus {\displaystyle X^{\*}\,} tähendab kompleksset konjugatsiooni.
*X**0* ja *X**N/2* on reaalsed väärtused ja ülejäänud DFT on täielikult määratud ainult *N/2-1* kompleksarvudega.

Üldistatud DFT (nihkunud ja mittelineaarne faas)
------------------------------------------------

On võimalik nihutada teisenduse proovivõttu aja ja/või sageduspiirkonna mingi *a* ja *b* väärtuse võrra. Seda nimetatakse tavaliselt **üldistatud DFT-ks** (**GDFT inglise keeles**) või **nihutatud DFT-ks.** Üldistatud DFT omab analoogsed omadusi nagu tavaline DFT:

{\displaystyle X\_{k}=\sum \_{n=0}^{N-1}x\_{n}e^{-{\frac {2\pi i}{N}}(k+b)(n+a)}\quad \quad k=0,\dots ,N-1.}
Kõige tihedamini kasutatakse nihet {\displaystyle 1/2} proovi võrra. Tavaline DFT vastab perioodilisele signaalile nii aja kui ka sageduse domeenis, {\displaystyle a=1/2} {\displaystyle (b=1)} esitab vastu perioodilist signaali ({\displaystyle X\_{k+N}=-X\_{k}}) ja vastupidi, kui {\displaystyle b=1/2} {\displaystyle (a=1)}. On olemas ka erijuht {\displaystyle a=b=1/2}, mida on nimetatud paaritu aja ja sagedusega DFT-ks (O2 DFT inglise keeles). Selliseid teisendusi kasutatakse sümmeetriliste signaalide puhul, et esitada erinevaid piirsümmeetriad. Reaalsignaalide puhul vastavad need erinevatele siinus- ja koosinus-signaalidele.

Teine huvitav juhtum on, kui {\displaystyle a=b=-(N-1)/2}, mida nimetatakse **tsentreeritud DFT-ks** (**CDFT** inglise keeles). See DFT omab kasulikku omadust: kui *N* on nelja korrutis, siis kõik neli tema ühisväärtust omavad võrdset mitmesust.

Diskreetset Fourier' teisendust võib vaadelda z-teisenduse erijuhuna, mida hinnatakse ühikringil komplekstasandil; z-teisendused vastavad kompleksnihetele *a* ja *b*, nagu on näidatud eespool.

Mitmemõõtmeline DFT
-------------------

Tavaline DFT tegeleb ühemõõtmelisega andmete massiiviga {\displaystyle x\_{n}} kuna see on funktsioon täpselt ühe diskreetse muutujaga *n*. Mitmemõõtmeline DFT tegeleb mitmemõõtmelise massiiviga {\displaystyle x\_{n\_{1},n\_{2},\dots ,n\_{d}}} ja avaldub *d* diskreetsete muutujatega funktsiooniks {\displaystyle n\_{\ell }=0,1,\dots ,N\_{\ell }-1} iga {\displaystyle \ell } {\displaystyle 1,2,\dots ,d} sees on defineeritud:

{\displaystyle X\_{k\_{1},k\_{2},\dots ,k\_{d}}=\sum \_{n\_{1}=0}^{N\_{1}-1}\left(\omega \_{N\_{1}}^{~k\_{1}n\_{1}}\sum \_{n\_{2}=0}^{N\_{2}-1}\left(\omega \_{N\_{2}}^{~k\_{2}n\_{2}}\cdots \sum \_{n\_{d}=0}^{N\_{d}-1}\omega \_{N\_{d}}^{~k\_{d}n\_{d}}\cdot x\_{n\_{1},n\_{2},\dots ,n\_{d}}\right)\right)\,,}
kus {\displaystyle \omega \_{N\_{\ell }}=\exp(-2\pi i/N\_{\ell })} ja *d* väljundi indeksid on {\displaystyle k\_{\ell }=0,1,\dots ,N\_{\ell }-1}. Natuke kompaktsemalt on see väljendatud kui vektor, kus {\displaystyle \mathbf {n} =(n\_{1},n\_{2},\dots ,n\_{d})} ja {\displaystyle \mathbf {k} =(k\_{1},k\_{2},\dots ,k\_{d})} on *d*-mõõtmelised vektorid indeksitega nullist {\displaystyle \mathbf {N} -1}, kus {\displaystyle \mathbf {N} -1=(N\_{1}-1,N\_{2}-1,\dots ,N\_{d}-1)}:

{\displaystyle X\_{\mathbf {k} }=\sum \_{\mathbf {n} =\mathbf {0} }^{\mathbf {N} -1}e^{-2\pi i\mathbf {k} \cdot (\mathbf {n} /\mathbf {N} )}x\_{\mathbf {n} }\,,}
kus {\displaystyle \mathbf {n} /\mathbf {N} } on {\displaystyle \mathbf {n} /\mathbf {N} =(n\_{1}/N\_{1},\dots ,n\_{d}/N\_{d})} ja summa tähistab pesastatud summeerimist üleval olevas valemis.

Pöördvõrdeline mitmemõõtmeline DFT avaldub samamoodi nagu ühe mõõtme korral:

{\displaystyle x\_{\mathbf {n} }={\frac {1}{\prod \_{\ell =1}^{d}N\_{\ell }}}\sum \_{\mathbf {k} =\mathbf {0} }^{\mathbf {N} -1}e^{2\pi i\mathbf {n} \cdot (\mathbf {k} /\mathbf {N} )}X\_{\mathbf {k} }\,.}
DFT väljendab sisendit {\displaystyle x\_{n}} sinusoidide superpositsiooniga, mitmemõõtmeline DFT väljendab sisendit tasandi lainete superpositsiooniga, või mitmemõõtmeliste sinusoididega. Ostsillatsiooni suund on {\displaystyle \mathbf {k} /\mathbf {N} } ja amplituudid on {\displaystyle X\_{\mathbf {k} }}. See lagunemine omab suurt tähtsust paljudel aladel, alustades digitaalsest pilditöötlusest (kahemõõtmeline) kuni diferentsiaalvõrrandite lahendamiseni. Tulemus on jagatud tasandite lainetele.

Mitmemõõtmeline DFT võib olla arvutatud kui iga mõõtme ühemõõtmelise DFT-de kompositsioon. Kahemõõtmelise DFT puhul {\displaystyle x\_{n\_{1},n\_{2}}} on {\displaystyle N\_{1}}sõltumatu ridade DFT-d (näiteks mööda {\displaystyle n\_{2}}) arvutatud esimesena uueks massiiviks {\displaystyle y\_{n\_{1},k\_{2}}}. Järgmisena {\displaystyle N\_{2}} sõltumatu veergude DFT-d (mööda {\displaystyle n\_{1}}) arvutatakse, et saada lõpptulemust {\displaystyle X\_{k\_{1},k\_{2}}}. Võib teha ka vastupidi ja esimesena arvutada veerud ja pärast read. Järjekord pole oluline, sest lõpus arvutatakse pesastatud summa.

Ühemõõtmelise DFT arvutamise algoritm on piisav selleks, et arvutada ka mitmemõõtmelised DFT-d. Selle algoritmi nimetuseks on *rea-veeru* algoritm. On olemas ka mõned teised algoritmid mitmemõõtmeliste DFT-de arvutamiseks.

### Reaalse sisendiga mitmemõõtmeline DFT

Kui sisend koosneb reaalsetest arvudest {\displaystyle x\_{n\_{1},n\_{2},\dots ,n\_{d}}}, DFT väljund omab konjugatiivset sümmeetriat nagu ühemõõtmelisel juhul:

{\displaystyle X\_{k\_{1},k\_{2},\dots ,k\_{d}}=X\_{N\_{1}-k\_{1},N\_{2}-k\_{2},\dots ,N\_{d}-k\_{d}}^{\*},}
kus tärn tähendab kompleksset konjugatsiooni ja {\displaystyle \ell }-is indeks on absoluutväärtus {\displaystyle N\_{\ell }} (iga {\displaystyle \ell =1,2,\ldots ,d} jaoks).

Kiire DFT
---------

**Kiire Fourier' teisendus** (Fast Fourier Transform = FFT) põhineb teisendamiseks vajalike arvutuste mahu vähendamisel teatud (korduvate) tulemuste (korduva) ärakasutamise teel või mõne spetsiifilise vaheteisenduse kasutamise teel. Seetõttu on FFT reeglina kasutatav vaid diskreeditud funktsioonide ja signaalide puhul. Siiski on ka analoogtehnikas mõned signaalitöötluse võtted vaadeldavad kui FFT teatud vormid.

FFT kasutatavusse ja FFT abil saadud tulemustesse tuleb seetõttu suhtuda vajaliku kriitilisusega, kuna lisaks diskreetmisest tulenevatele iseärasustele tuleb lisaks arvesse võtta veel ka FFT enda olemusest tulenevaid iseärasusi (lühike ajaintervall ja sellest tulenevalt nn. vaatluse akna(funktsiooni) kasutamine, sünkroonsuse probleemid).

FFT kasutatavus on praktikas laialdane (heli- ja sidetehnika, jm), kuna paljudel juhtudel suurt täpsust selle mõõtetehnilises mõttes ei nõuta (heli puhul loetakse piisavaks 1dB ehk ~10% täpsust). Arvutuse täpsuse probleemid on aga üldjuhul lahendatavad.

Kiire DFT arvutamiseks kasutatakse mitu algoritmi: Cooley-Tukey FFT, Prime-factor FFT, Bruun's FFT, Rader's FFT, ja Bluestein's FFT.

Rakendused
----------

DFT on kasutatud paljudel aladel ja siin on näidatud ainult tähtsamad neist. Peaaegu kõik DFT rakendused otseselt põhinevad sellel põhimõttel, et DFT-d ja pöörd-DFT on võimalik efektiivselt arvutada Kiire Fourier' teisenduse abil .

### Spektraalanalüüs

Kui DFT-d kasutatakse spektraalanalüüsiks, siis {\displaystyle \{x\_{n}\}\,} jada tavaliselt kujutab piiratud koguse ühtlaste aja vahedega paigutatud mingi signaali proove {\displaystyle x(t)\,}, kus *t* on aeg. Üleminek pidevalt ajalt proovidele muudab {\displaystyle x(t)\,} Fourier' teisendust Diskreetse aja Fourier teisenduseks (DTFT), mis tavaliselt tähendab deformeerimist (*aliasing* inglise keeles). Deformeerimise minimeerimiseks tuleb õigesti valida proovide võitmise sagedust, arvestades Nyquist-i sagedust. Samas, tuleb proovide jada konverteerida väga suurest(või lõpmatust) sobivaks, et vältida teist deformeerimist, nimetatud lekkimine, ehk resolutsiooni kaotamine DFT-s. Kui kättesaadav andmete hulk on liiga väike, või andmete töötlemiseks (vajaliku resolutsiooni leidmiseks) pole piisavalt aega, kasutatakse mitu DFT-d spektrogrammi tegemiseks. Tulemusena on võimsuse spekter kus esineb müra ja/või juhuslikud andmed, mida saab filtreerida mitu DFT tulemuste võrdlemisega. See protseduur vähendab spektri dispersiooni ja on nimetatud perioodgrammiks. Tuntud näited selle meetodi kasutamisest on Welchi ja Bartletti meetodid. Üldine meetod võimsuse spektri hindamiseks lärmaka signaali puhul on nimetatud spektraalne hindamine.

DFT ise on lõplik moonutamise allikas, sest DFT on ainult diskreetne proovide võtmine pidevast funktsioonist. Seda on võimalik leevendada DFT resolutsiooni suurendamisega.

### Andmete pakkimine

Digitaalse signaalitöötluse valdkond põhineb sagedusdomeenis toimuvates operatsioonides (Fourier' teisenduse tulemustes). Näiteks andmekaostega piltide töötlus ja heli pakkimise meetodid kasutavad DFT-d: signaal on lõigatud lühikesteks segmentideks (proovide võtmine), iga lõik on transformeeritud, ja siis peaaegu märkamatu kõrgesageduslikud Fourier' koefitsiendid, kõrvaldatakse. Dekompressor arvutab pöördtransformatsiooni, mis põhineb uutel DFT väärtustel. Pakkimise rakendused kasutavad sageli spetsialiseeritud DFT vormi, diskreetse koosinus teisendust, või juba seda modifitseeritud versiooni. Mõned suhteliselt hiljutised kokkupakkimisalgoritmid kasutavad nn laineke teisendust, mis annab ühtlasema kompromissi aja ja sageduse domeenide vahel. JPEG2000 juhul, see väldib vigaste piltide tekkimist, mis ilmuvad, kui pildid on liiga palju tihendatud algse JPEG pakkimise tõttu.

### Osatuletistega diferentsiaalvõrrandid

Diskreetne Fourier' teisendust kasutatakse sageli osatuletistega diferentsiaalvõrrandite ahendamisel, kus DFT-d kasutatakse Fourier' rea liigikaudseks arvutamiseks. Selle lähenemise eeliseks on see, et signaali laiendatakse kompleks eksponentideks {\displaystyle e^{inx}}, mis on diferentseerimise tulemuste funktsioonide hulk: {\displaystyle {{\text{d}}{\big (}e^{inx}{\big )}}/{\text{d}}x=ine^{inx}}. Seega Fourier' esindamises, diferentseerimine on lihtne – see on korrutamine väärtusega {\displaystyle in}. Lineaarne diferentsiaalvõrrand konstantsete kordajatega muundub kergesti lahendatavaks algebraliseks võrrandiks. Tulemuse tagasi aja domeeni teisendamiseks tuleb kasutada pöörd teisendust. Seda ka nimetatakse spektraalmeetodiks ehk spektraalanalüüsiks.

### Polünoomiline korrutamine

Oletame, et me soovime arvutada polünoomi ***c*(*x*) = *a*(*x*) · *b*(*x*)**. Tavaline koefitsientide c väljend hõlmab lineaarset (atsüklilist) konvolutsiooni. Seda saab kirjutada tsüklilise konvolutsiooni kasutades, kui võtta koefitsient vektoriteks ***a*(*x*)** ja ***b*(*x*)** pideva väljendiga esimeseks, seejärel lisades nulle nii, et tulemus koefitsient vektorid **a** ja **b** on mõõtmes ***d* > deg(*a*(*x*)) + deg(*b*(*x*))**. Siis,

{\displaystyle \mathbf {c} =\mathbf {a} \*\mathbf {b} }
Kus **c** on koefitsientide vektor ***c*(*x*)**, ja konvolutsiooni operaator {\displaystyle \*\,} on

{\displaystyle c\_{n}=\sum \_{m=0}^{d-1}a\_{m}b\_{n-m\ \mathrm {mod} \ d}\qquad \qquad \qquad n=0,1\dots ,d-1}
Aga konvolutsioon muutub korrutamiseks DFT all:

{\displaystyle {\mathcal {F}}(\mathbf {c} )={\mathcal {F}}(\mathbf {a} ){\mathcal {F}}(\mathbf {b} )}
Siin on vektor võetud elemendi kaupa. Seega polünoomi koefitsiendid *c*(*x*) on vaid 0, ... deg(*a*(*x*)) + deg(*b*(*x*)) ja koefitsient vektorist on

{\displaystyle \mathbf {c} ={\mathcal {F}}^{-1}({\mathcal {F}}(\mathbf {a} ){\mathcal {F}}(\mathbf {b} )).}
Kiire DFT algoritmi kasutuses ajaline keerukus on **O (*N* log *N*)**. Teisenduse operatsiooni jaoks tavaliselt kasutatakse Cooley-Tukey FFT algoritmi, sest see on lihtne ja suhteliselt kiire. Sel juhul peab {\displaystyle d} olema väiksem täisarv, mis on suurem, kui sisend-polünoomi algtegurite summa.

#### Suurte täisarvude korrutamine

Kõige kiirem tuntud algoritmidest väga suurte täisarvude korrutamiseks on polünoomide korrutamise meetod. Täisarve võib käsitleda, kui polünoomide väärtused väärtustatud konkreetse arvu alusel, koos polünoomi koefitsientidega, mis on selle arvu numbrid. Pärast polünoomilist korrutamist lõpetab korrutamise suhteliselt väikese keerukusega nn *carry*-paljundamine.

#### Konvolutsioon

Kuna konvolutsioon on sagedus domeenis lihtsalt korrutamine, siis palju lihtsam on teisendada sisendid DFT abil, järgmisena teha korrutist ja viimasena teisendada tulemust pöördteisendusega tagasi aja domeeni.

Olulised DFT paarid
-------------------

| {\displaystyle x\_{n}={\frac {1}{N}}\sum \_{k=0}^{N-1}X\_{k}e^{i2\pi kn/N}} | {\displaystyle X\_{k}=\sum \_{n=0}^{N-1}x\_{n}e^{-i2\pi kn/N}} | Nimetus |
| --- | --- | --- |
| {\displaystyle x\_{n}e^{i2\pi n\ell /N}\,} | {\displaystyle X\_{k-\ell }\,} | Nihe |
| {\displaystyle x\_{n-\ell }\,} | {\displaystyle X\_{k}e^{-i2\pi k\ell /N}\,} |
| {\displaystyle x\_{n}\in \mathbb {R} } | {\displaystyle X\_{k}=X\_{N-k}^{\*}\,} | Reaalne DFT |
| {\displaystyle a^{n}\,} | {\displaystyle \left\{{\begin{matrix}N&{\mbox{if }}a=e^{i2\pi k/N}\\{\frac {1-a^{N}}{1-a\,e^{-i2\pi k/N}}}&{\mbox{mujal}}\end{matrix}}\right.} | Geomeetriline progressioon |
| {\displaystyle {N-1 \choose n}\,} | {\displaystyle \left(1+e^{-i2\pi k/N}\right)^{N-1}\,} | Binominaalne teoreem |
| {\displaystyle \left\{{\begin{matrix}{\frac {1}{W}}&{\mbox{if }}2n<W{\mbox{ or }}2(N-n)<W\\0&{\mbox{mujal}}\end{matrix}}\right.} | {\displaystyle \left\{{\begin{matrix}1&{\mbox{if }}k=0\\{\frac {\sin \left({\frac {\pi Wk}{N}}\right)}{W\sin \left({\frac {\pi k}{N}}\right)}}&{\mbox{mujal}}\end{matrix}}\right.} | {\displaystyle x\_{n}} on *W* punktiline ristkujuline aknafunktsioon keskpunktiga *n*=0, kus *W* on paaritu reaalarv, ja {\displaystyle X\_{k}} on sinc funktsioon. |
| {\displaystyle \sum \_{j\in \mathbb {Z} }\exp \left(-{\frac {\pi }{cN}}\cdot (n+N\cdot j)^{2}\right)} | {\displaystyle {\sqrt {cN}}\cdot \sum \_{j\in \mathbb {Z} }\exp \left(-{\frac {\pi c}{N}}\cdot (k+N\cdot j)^{2}\right)} | Skaleeritud Gaussi funktsioonide diskreetimine ja perioodiline liitmine, kus {\displaystyle c>0}. |

Üldistused
----------

### Esinduse teooria

DFT-d võib tõlgendada kui kompleksset piiratud tsüklilise rühma esinduse teooriat. Teisisõnu, *n* kompleksarvude jada võib olla tõlgendatud nägu *n*-mõõtmeline kompleksne ruum **C***n* või lõplikku tsüklilise kompleksarvude rühma funktsioon f, **Z***n* → **C**. *f* on tsükliline piiratud rühma klassi funktsioon, ning seda saab väljendada taandumatu tähemärkide lineaarse kombinatsioonina selles rühmas. Need märgid on ühe juured.

### Muud valdkonnad

Palju DFT omadustest sõltuvad ainult sellest, et {\displaystyle e^{-{\frac {2\pi i}{N}}}} on primitiivsed ühe juured. Need omadused on täielikkus, ortogonaalsus, Plancherel/Parseval, perioodilisus, nihe, konvolutsioon ja paljud kiire DFT algoritmid. Seetõttu DFT-d võib defineerida ühe juure mõiste kasutades ka muudes valdkondades. Neid üldistusi tavaliselt nimetatakse number-teoreetilised teisendused (*number-theoretic transforms* (NTT)).

### Muud piiratud rühmad

Tavaline DFT tegeleb *x*0, *x*1, …, *x**N*−1 kompleksarvude jadaga, mis võib olla vaadeldud funktsioonina {0, 1, …, *N* − 1} → **C**. Mitmemõõtmeline DFT tegeleb mitmemõõtmeliste jadadega, mis võivad olla vaadeldud funktsioonidena

{\displaystyle \{0,1,\ldots ,N\_{1}-1\}\times \cdots \times \{0,1,\ldots ,N\_{d}-1\}\to \mathbb {C} .}
Siis tavalist DFT-d vaadeldakse tsüklilise rühma Fourier' teisendusena, ja mitmemõõtmeline DFT on tsükliliste rühmade summa Fourier' teisendus.

Märkused
--------

1. ↑  Selles kontekstis {\displaystyle W\_{N}} defineeritakse nägu N-is primitiivne ühe juur, {\displaystyle W\_{N}=e^{-2\pi i/N}}, saab järgmist vormi:

{\displaystyle X\_{k}=\sum \_{n=0}^{N-1}x\_{n}\cdot W\_{N}^{kn}}
2. ↑  DFT definitsiooni saab ka kirjeldada DFT maatriksi abil nagu lineaarne transformatsioon piiratud-mõõtmelises vektorruumis. Kui DFT on skaleeritud sobivalt, siis see muutub ühtlaseks maatriksiks ja XK on × koefitsiendid ortonormaalsel alusel.

Viited
------

1. 1 2 3 4  Smith, Steven W. (1999). "Chapter 8: The Discrete Fourier Transform". *The Scientist and Engineer's Guide to Digital Signal Processing* (Second ed.). San Diego, Calif.: California Technical Publishing. ISBN 0-9660176-3-3. Vaadatud 31.10.2016.
2. ↑  Brigham, E. Oran (1988). *The fast Fourier transform and its applications*. Englewood Cliffs, N.J.: Prentice Hall. ISBN 0-13-307505-2.
3. ↑  Strang, Gilbert (1994). "Wavelets". *American Scientist*. **82** (3): 253. Vaadatud 31.10.2016. This is the most important numerical algorithm of our lifetime...
4. ↑  Sahidullah; Saha, Goutam (2012). "A Novel Windowing Technique for Efficient Computation of MFCC for Speaker Recognition". *IEEE Signal Processing Letters*. **20**: 149–152. Vaadatud 31.10.2016.
5. 1 2 3  Бахурин Сергей. "Дискретное преобразование Фурье (ДПФ)". Vaadatud 31.10.2016.
6. 1 2  J. Cooley, P. Lewis, and P. Welch (1969). "The finite Fourier transform". *IEEE Trans. Audio Electroacoustics*. **17** (2): 77–85.`{{cite journal}}`: CS1 hooldus: mitu nime: autorite loend (link)
7. ↑  Oppenheim, Alan V.; Schafer, R. W.; and Buck, J. R. (1999). *Discrete-time signal processing*. Upper Saddle River, N.J.: Prentice Hall.`{{cite book}}`: CS1 hooldus: mitu nime: autorite loend (link)
8. ↑  Ricardo, Henry J. (2016). "A Modern Introduction to Differential Equations". Lk 428
9. 1 2 3  Бахурин Сергей (2015). "Свойства дискретного преобразования Фурье (ДПФ)" (vene). Vaadatud 31.10.2016.
10. 1 2  P. Duhamel; B. Piron; J. M. Etcheto (1988). "On computing the inverse DFT". *IEEE Trans. Acoust. Speech and Sig. Processing*. **36** (2): 285–286.
11. 1 2  Juan G. Vargas-Rubio; Balu Santhanam (2005). "On the multiangle centered discrete fractional Fourier transform". *IEEE Sig. Proc. Lett*. **12** (4): 273–276.
12. 1 2  Nussbaumer, Henri J. (1982). *Fast Fourier Transform and Convolution Algorithms* (Second Corrected and Updated Edition. ed.). Springer-Verlag. Lk-d 134–137.
13. ↑  Arfken, G. (1985) "Convolution Theorem." §15.5 in *Mathematical Methods for Physicists, 3rd ed.* Orlando, FL: Academic Press, pp. 810-814, Vaadatud 2016-10-31
14. 1 2  G.B. Wright, M. Javed, H. Montanelli, and L.N. Trefethen,(2015), "Extension of Chebfun to periodic functions," *SIAM. J. Sci. Comput.*, 37 , C554-C573 Vaadatud 31-10-2016
15. 1 2  Kamisetty Ramam Rao; Patrick C. Yip (2000). *The Transform and Data Compression Handbook*. CRC Press. ISBN 978-1-4200-3738-8. Vaadatud 31.10.2016.
16. ↑  Julius O. Smith (2007). *Mathematics of the Discrete Fourier Transform (DFT): With Audio Applications*. Julius Smith. ISBN 978-0-9745607-4-8. Vaadatud 31.10.2016.
17. ↑  TE 302 DISCRETE SIGNALS AND SYSTEMS 2004, Vaadatud 31-10-2016
18. 1 2  J. H. McClellan; T. W. Parks (1972). "Eigenvalues and eigenvectors of the discrete Fourier transformation". *IEEE Trans. Audio Electroacoust*. **20** (1): 66–74.
19. ↑  Magdy Tawfik Hanna, Nabila Philip Attalla Seif, and Waleed Abd El Maguid Ahmed (2004). "Hermite-Gaussian-like eigenvectors of the discrete Fourier transform matrix based on the singular-value decomposition of its orthogonal projection matrices". *IEEE Trans. Circ. Syst. I*. **51** (11): 2245–2254.`{{cite journal}}`: CS1 hooldus: mitu nime: autorite loend (link)
20. 1 2  Massar, S.; Spindel, P. (2008). "Uncertainty Relation for the Discrete Fourier Transform". *Physical Review Letters*. **100** (19).
21. ↑  DeBrunner, Victor; Havlicek, Joseph P.; Przebinda, Tomasz; Özaydin, Murad (2005). "Entropy-Based Uncertainty Measures for {\displaystyle L^{2}(\mathbb {R} ^{n}),\ell ^{2}(\mathbb {Z} )}, and {\displaystyle \ell ^{2}(\mathbb {Z} /N\mathbb {Z} )} With a Hirschman Optimal Transform for {\displaystyle \ell ^{2}(\mathbb {Z} /N\mathbb {Z} )}" (PDF). *IEEE Transactions on Signal Processing*. **53** (8): 2690. Vaadatud 31.10.2016.
22. 1 2  Donoho, D.L.; Stark, P.B (1989). "Uncertainty principles and signal recovery". *SIAM Journal on Applied Mathematics*. **49** (3): 906–931.
23. 1 2  Santhanam, Balu; Santhanam, Thalanayar S. (1988) "*Discrete Gauss-Hermite functions and eigenvectors of the centered discrete Fourier transform*", Proceedings of the 32nd IEEE *International Conference on Acoustics, Speech, and Signal Processing* (ICASSP 2007, SPTM-P12.4), vol. III, pp. 1385-1388. Vaadatud 31-10-2016
24. 1 2  Akansu, Ali N.; Agirman-Tosun, Handan (2010) "*Generalized Discrete Fourier Transform With Nonlinear Phase*", IEEE *Transactions on Signal Processing*, vol. 58, no. 9, pp. 4547-4556, Vaadatud 31-10-2016
25. 1 2 3 4 5 6  Richard tolimieri; Myoung An; Chao Lu (2012). *Mathematics of Multidimensional Fourier Transform Algorithms*. Springer Science & Business Media. ISBN 978-1-4684-0205-6. Vaadatud 31.10.2016.
26. 1 2 3 4  C. Sidney Burrus, Ivan Selesnick, Markus Pueschel, Matteo Frigo, and Steven G. Johnson (2008). *Fast Fourier Transforms*, Connexions
27. ↑  Shamgar Gurevich; Ronny Hadani; Nir Sochen (2008). "The finite harmonic oscillator and its applications to sequences, communication and radar". *IEEE Transactions on Information Theory*. **54** (9): 4239–4253.
28. ↑  David Salomon (20. märts 2007). *Data Compression: The Complete Reference*. Springer Science & Business Media. ISBN 978-1-84628-603-2. Vaadatud 31.10.2016.
29. ↑  Nakhle H. Asmar (2016). "10". *Partial Differential Equations with Fourier Series and Boundary Value Problems: Third Edition*. Courier Dover Publications. ISBN 978-0-486-80737-9. Vaadatud 31.10.2016.
30. ↑  J.A. Storer (2001). "5". *An Introduction to Data Structures and Algorithms*. Springer Science & Business Media. ISBN 978-0-8176-4253-2. Vaadatud 31.10.2016.
31. ↑  Umberto Cherubini; Pietro Rossi; Sabrina Mulinacci (2010). *Fourier Transform Methods in Finance*. John Wiley & Sons. Lk 208. ISBN 978-0-470-68492-4. Vaadatud 31.10.2016.
32. ↑  Smith, Stephen W (1997). "13.Convolution". *The Scientist and Engineer's Guide to Digital Signal Processing* (1 ed.). California Technical Publishing. ISBN 0966017633. Vaadatud 31.10.2016.
33. ↑  Serre, Jean-Pierre: ''Linear Representations of Finite Groups.'' Springer Verlag, New York 1977, ISBN 0-387-90190-6.Jean-Pierre Serre (2012). *Linear Representations of Finite Groups*. Springer Science & Business Media. ISBN 978-1-4684-9458-7. Vaadatud 31.10.2016.
34. ↑  Vladimir P. Gerdt; Evgenii V. Vorozhtsov; Ernst W. Mayr (2013). *Computer Algebra in Scientific Computing: 15th International Workshop, CASC 2013, Berlin, Germany, September 9-13, 2013, Proceedings*. Springer. ISBN 978-3-319-02297-0. Vaadatud 31.10.2016.
35. ↑  Terras, Audrey (1999), *Fourier Analysis on Finite Groups and Applications*, Cambridge University Press, lk 251, ISBN 978-0-521-45718-7

Välislingid
-----------

* Interactive explanation of the DFT
* Matlab tutorial on the Discrete Fourier Transformation
* Interactive flash tutorial on the DFT
* Mathematics of the Discrete Fourier Transform by Julius O. Smith III
* FFTW: Fast implementation of the DFT - coded in C and under General Public License (GPL)
* General Purpose FFT Package: Yet another fast DFT implementation in C & FORTRAN, permissive license
* Explained: The Discrete Fourier Transform