Fourier-muunnos

Fourier-muunnos (myös Fourier’n muunnos) on matematiikassa käytetty jatkuva integraalimuunnos. Muunnosta käytetään erityisesti analyysissä differentiaaliyhtälöiden ratkaisemisessa ja signaalinkäsittelyssä erilaisiin taajuusanalyysiä vaativiin sovelluksiin. Muunnos perustuu teoreemaan, jonka mukaan mikä tahansa riittävän säännöllinen jatkuva funktio voidaan esittää sinimuotoisten funktioiden integraalina. Fourier-muunnoksella voidaan laskea näiden sinimuotoisten komponenttien amplitudi ja vaihe.

Määritelmä

Funktion ${\displaystyle f(x)\,}$ Fourier-muunnos ${\displaystyle {\hat {f}}(\omega )}$ määritellään

${\displaystyle {\hat {f}}(\omega )=\int _{-\infty }^{\infty }f(x)e^{-i\omega x}\,dx}$,

missä ${\displaystyle \omega \,}$ on kulmataajuus. Muoto ${\displaystyle e^{-i\omega x}\ }$ liittyy trigonometrisiin funktioihin siten, että kompleksieksponentin määritelmä on ${\displaystyle e^{a+ib}=e^{a}(\cos b+i\sin b)\ }$.

Fourier-muunnokselle on olemassa käänteismuunnos, joka määritelmän mukaiselle funktiolle on

${\displaystyle f(x)={\frac {1}{2\pi }}\int _{-\infty }^{\infty }{\hat {f}}(\omega )e^{i\omega x}\,d\omega }$.

Funktion ${\displaystyle f(x)\,}$ Fourier-muunnoksesta voidaan käyttää myös vaihtoehtoista merkintätapaa ${\displaystyle {\mathcal {F}}\{f(x)\}={\hat {f}}(x)}$. Tätä merkintää käytetään esimerkiksi differentiaalianalyysissä, jolloin halutaan selventää Fourier-muunnoksen käyttäminen jonkin muun integraalimuunnoksen sijasta. Kirjallisuudessa esiintyy usein myös eri tavoin normalisoituja variaatioita muunnoksesta. Fourier-muunnettu funktio voidaan ajatella alkuperäisen funktion esityksenä taajuustasossa. Fourier-muunnetun funktion taajuuskomponenttia ${\displaystyle \omega \,}$ vastaava amplitudi on

${\displaystyle A(\omega )=|{\hat {f}}(\omega )|}$

ja vaihe

${\displaystyle \phi (\omega )=\operatorname {arg} ({\hat {f}}(\omega ))}$.

Amplitudispektri antaa vaihespektriä enemmän tietoa, nähdään kaikki taajuuskomponentit ja niitä vastaavat määrät.

Fourier-muunnoksen ominaisuuksia

Olkoon ${\displaystyle f(x)\,}$, ${\displaystyle g(x)\,}$ ja ${\displaystyle h(x)\,}$ integroituvia funktioita ja näitä vastaavat Fourier-muunnokset ${\displaystyle {\hat {f}}(\xi )}$, ${\displaystyle {\hat {g}}(\xi )}$ ja ${\displaystyle {\hat {h}}(\xi )}$. Fourier-muunnoksella on seuraavat ominaisuudet ^[1].

• Lineaarisuus

Kaikille kompleksiluvuille a ja b pätee, että jos ${\displaystyle h(x)=af(x)+bg(x)\ }$, niin ${\displaystyle {\hat {h}}(\xi )=a\cdot {\hat {f}}(\xi )+b\cdot {\hat {g}}(\xi )\ }$.

• Siirto aikatasossa

Kaikille reaaliluvuille ${\displaystyle x_{0}\,}$ pätee, että jos h(x) = ƒ(x − x₀), niin ${\displaystyle {\hat {h}}(\xi )=e^{-i2\pi \xi x_{0}}{\hat {f}}(\xi )\ }$.

• Modulointi aikatasossa

Kaikille reaaliluvuille ${\displaystyle \xi _{0}}$ pätee, että jos ${\displaystyle h(x)=e^{i2\pi \xi _{0}x}f(x)\ }$, niin ${\displaystyle {\hat {h}}(\xi )={\hat {f}}(\xi -\xi _{0})\ }$.

• Aikatason kääntö

Jos ${\displaystyle h(x)=f(-x)\ }$, niin ${\displaystyle {\hat {h}}(\xi )={\hat {f}}(-\xi )\ }$.

• Aikatason skaalaus

Kaikille ei-negatiivisille reaaliluvuille ${\displaystyle a}$ pätee, että jos ${\displaystyle h(x)=f(ax)\ }$, niin ${\displaystyle {\hat {h}}(\xi )={\frac {1}{|a|}}{\hat {f}}\left({\frac {\xi }{a}}\right)\ }$. Kun ${\displaystyle a=-1\ }$, niin ominaisuus palautuu aikatason käännöksi.

• Konjugaatio

Jos ${\displaystyle h(x)={\overline {f(x)}}\ }$, niin ${\displaystyle {\hat {h}}(\xi )={\overline {{\hat {f}}(-\xi )}}\ }$.

• Konvoluutio

Jos ${\displaystyle h(x)=\left(f*g\right)(x)\ }$, niin ${\displaystyle {\hat {h}}(\xi )={\hat {f}}(\xi )\cdot {\hat {g}}(\xi )\ }$.

• Ominaisfunktiot

Fourier-muunnoksen ominaisfunktioita ovat Hermiten funktiot. Toisin sanoen näiden funktioiden Fourier-muunnos on funktio itse, ominaisarvolla kerrottuna.

Diskreetti Fourier-muunnos

Diskreetti Fourier-muunnos eli DFT on Fourier-muunnoksen diskreettiaikainen yleistys.^[2] Siinä signaali ajatellaan jaksolliseksi, jolloin se voidaan esittää äärellisenä Fourier'n sarjana ja integraali korvautuu summalausekkeella:

${\displaystyle F_{n}=\sum _{k=0}^{N-1}f_{k}e^{-i(2\pi k/N)n}\quad ,n=0,...,N-1}$

missä ${\displaystyle f_{k}\,}$ on ${\displaystyle N\,}$:n pituinen reaali- tai kompleksiarvoinen sarja.

Vastaava käänteismuunnos on

${\displaystyle f_{k}={\frac {1}{N}}\sum _{n=0}^{N-1}F_{n}e^{i(2\pi k/N)n}}$.

FFT

FFT (engl. Fast Fourier Transform) eli nopea Fourier-muunnos tarkoittaa tehokasta algoritmia diskreetin Fourier-muunnoksen ja sen käänteismuunnoksen laskemiseksi. Yleisin FFT on Cooleyn–Tukeyn algoritmi, jonka tunsi jo C.F Gauss vuonna 1805 ja käytti sitä Pallas- ja Juno-asteroidien ratojen laskentaan.^lähde? Työ ei ollut kovin tunnettu. Erilaisia rajoitettuja versioita kehitettiin 1800-luvulla ja 1900-luvun alkupuolella. Nykyisen maineensa FFT saavutti vuonna 1965, jolloin Cooley IBM:stä ja Tukey Princetonista osoittivat työssään algoritmin soveltuvuuden tietokoneohjelmointiin.

Jos diskreetti Fourier-muunnos lasketaan suoraan DFT:n määritelmästä, tarvittavien laskentaoperaatioiden määrä on verrannollinen näytepisteiden määrän neliöön ${\displaystyle N^{2}\,}$. Sen sijaan FFT-algoritmien aikavaativuusluokka on ${\displaystyle O(N\,\operatorname {log} \,N\,)}$. FFT:n nopeusero verrattuna suoraan diskreetin muunnoksen määritelmästä laskemiseen muuttuu hyvin merkittäväksi, kun näytepisteiden määrä on suuri. Käytännön sovelluksissa diskreetti Fourier-muunnos lasketaankin lähes aina numeerisesti FFT:n avulla.

Käytännön sovelluksia

FFT:tä käytetään hyväksi tekniikan ja fysiikan sovelluksissa, jotka perustuvat ilmiöiden jaksollisuuden tai spektrin mittaamiseen. Tärkeitä FFT:n sovelluksia ovat esimerkiksi spektrianalyysi ja OFDM tietoliikennetekniikassa sekä kuvan rekonstruktio magneettikuvauksessa.

MP3-äänenpakkausmenetelmässä äänen spektri, eli äänisignaalin taajuustasoesitys, lasketaan käyttäen FFT:tä (muunnettu DCT) ja spektri pakataan (häviöllisesti) jättämällä pois spektrikomponentit, joiden energia on pieni. Spektriksi hajottamisesta on se etu, että eri taajuisille äänille voidaan käyttää eri tarkkuutta. Ihmisen kuulo erottaa tietyt taajuudet tarkemmin, jolloin ne voidaan pakata vähemmällä häviöllä. Purku tapahtuu käänteismuuntamalla eli syntetisoimalla spektri takaisin aikatasoon.

Terästeollisuudessa voidaan, tutkimalla valssatun nauhan paksuusprofiilia FFT:llä, löytää epäkeskeisesti hiotut valssit tai kuluneet laakerit. Laakereiden kunnon seurantajärjestelmät perustuvat usein FFT:hen.

Fourier-muunnoksen ominaisuutta modulaatio aikatasossa käytetään siirtämään signaalia taajuustasossa. Radiotekniikassa lähetettävä tai vastaanotettu signaali siirretään halutulle taajuusalueelle moduloimalla signaalia sekoitusasteessa paikallisoskillaattorin signaalilla. Myös moduloimalla esimerkiksi yksittäistä sinimuotoista signaalia aikatasossa, sen värähtelytaajuutta voidaan kasvattaa tai pienentää. Tätä ominaisuutta käytetään hyväksi äänenkäsittelyssä yksittäisen äänen sävelkorkeuden korjaamiseksi^lähde?. Modulointitekniikalla tuotettu sävelkorkeuden korjaus ei nopeuta tai hidasta ääniraitaa kuten alkuperäisen ääniraidan suora skaalaaminen. Toisaalta, jos se kohdistetaan samalla kertaa useampaan eritaajuiseen signaaliin, niiden taajuussuhteet muuttuvat.

Katso myös

• Diskreetti kosinimuunnos

Lähteet

1. Pinsky, Mark: Introduction to Fourier Analysis and Wavelets. Brooks/Cole, 2002. ISBN 0-534-37660-6 (englanniksi)
2. Fourier–menetelmät (.pdf) (B2. Diskreetti Fourier–muunnos) Helsingin yliopisto. Arkistoitu 13.11.2012. Viitattu 12.1.2011.

Kirjallisuutta

• Oppenheim, Alan V.; Willsky Alan S.; with Nawab, Syed Hamid: Signals and Systems, s. 1–957. Prentice-Hall Signal Processing Series, 1997 (1983). ISBN 0-13-651175-9