Fourierova transformace

Fourierova transformace je integrální transformace převádějící signál mezi časově a frekvenčně závislým vyjádřením pomocí harmonických signálů, tj. funkcí \sin a \cos, obecně tedy funkcí komplexní exponenciály. Slouží pro převod signálů z časové oblasti do oblasti frekvenční. Signál může být buď ve spojitém či diskrétním čase.

Spojitý čas

Definice

Fourierova transformace S(\omega) funkce s(t) je definována integrálním vztahem :S(\omega)=\int\limits_{-\infty}^{\infty} s(t){\mathrm e}^{-{\mathrm i}\omega t}\,{\mathrm d}t

Funkci s(t) vypočteme z S(\omega) inverzní Fourierovou transformací :s(t)=\frac{1}{2\pi}\int\limits_{-\infty}^{\infty} S(\omega){\mathrm e}^{{\mathrm i}\omega t}\,{\mathrm d}\omega

Nevlastní integrály chápeme ve smyslu Cauchyovy hlavní hodnoty, tj. :\int\limits_{-\infty}^{\infty} [.]\,{\mathrm d}=\lim_{T \to \infty}\int\limits_{-T}^{T} [.]\, {\mathrm d}

Dvojice ve Fourierově transformaci se nazývají originál (zde s(t)) a obraz (zde S(\omega)). Vztah mezi originálem a obrazem vyjadřujeme zápisem :S(\omega)=\mathcal{F}[s(t)] a s(t)=\mathcal{F}^{-1}[S(\omega)].

V technické oblasti je \omega úhlová frekvence, S(\omega) představuje spektrum signálu s(t).

Spektrum je komplexní veličina a lze vyjádřit ve tvaru S(\omega)=\left|S(\omega)\right|{\mathrm e}^{{\mathrm i} \mathrm{arg}\,S(\omega)}. Velikost \left|S(\omega)\right| nazýváme amplitudové spektrum a úhel \mbox{arg}\,S(\omega) fázové spektrum signálu.

Vlastnosti Fourierovy transformace

Věta o linearitě

Lineární kombinaci signálů odpovídá lineární kombinace jejich spekter :\mathcal{F}\left[\sum_{i} c_is_i(k)\right]=\sum_{i} c_i\mathcal{F}[s_i(k)]

Věta o změně měřítka (Podobnost)

Má-li signál s(t) spektrum S(\omega), má signál s(at), a\neq0 spektrum :\frac{1}{\left|a\right|}S\left(\frac{\omega}{a}\right).

Tedy rozšíření signálu v časové oblasti odpovídá zúžení ve spektrální oblasti a naopak.

Posun signálu v čase (Posunutí)

Má-li signál s(t) spektrum S(\omega), má signál zpožděný o veličinu a spektrum :\mathcal{F}[s(t-a)]={\mathrm e}^{-{\mathrm i}\omega a}S(\omega)

Amplitudové spektrum posunutého signálu se nemění, mění se jen fázové spektrum a to přímo úměrně zpoždění a kmitočtu. Na rozdíl od věty o translaci v Laplaceově transformaci platí věta pro libovolné a, tedy i pro a.

Spektrum reálného signálu

Je-li signál reálný, pak pro jeho spektrum platí: * amplitudové spektrum je sudou funkcíplyne to z F(\omega) = F^*(-\omega); viz [url=https://web. archive. +moreorg/web/20070331195935/http://www. exos. sk/GDFS/Infokomunikacia/fourtrans. pdf],]strana 6 * fázové spektrum je [[Sudé a liché funkce#Lichá funkce|lichou[/url]] funkcí * spektrum sudého signálu je sudou reálnou funkcí * spektrum lichého signálu je lichou ryze imaginární funkcí.

Diskrétní čas

Definice

Fourierova transformace S(\Omega) posloupnosti s(k) je definována vztahem :S(\Omega)=\sum_{k=-\infty}^{\infty} s(k){\mathrm e}^{-{\mathrm i}\Omega k}

Posloupnost s(k) vypočteme z S(\Omega) inverzní Fourierovou transformací :s(k)=\frac{1}{2\pi}\int\limits_{0}^{2\pi} S(\Omega){\mathrm e}^{{\mathrm i}\Omega k}\,{\mathrm d}\Omega

Někteří autoři označují tuto transformaci DtFT (discrete-time Fourier transformation), aby ji odlišili od Fourierovy transformace spojitého signálu. Zde nebudeme značením nijak odlišovat Fourierovu transformaci spojitého a diskrétního signálu. +more Vztah mezi signálem a jeho spektrem budeme tedy značit :S(\Omega)=\mathcal{F}[s(k)] a :s(k)=\mathcal{F}^{-1}[S(\Omega)].

Spektrum diskrétního signálu se od spektra spojitého signálu liší tím, že je periodické s periodou 2\pi.

Diskrétní Fourierova transformace

Definiční vztahy Fourierovy transformace vyžadují znalost matematického vyjádření signálu či spektra. Pokud však zpracováváme naměřené hodnoty, tj. +more známe vzorky signálu či spektra z konečného intervalu, stojíme před problémem, jak určit spektrum z vzorků signálu či signál ze vzorků spektra. K tomu účelu používáme numerické metody, která je známa jako diskrétní Fourierova transformace (DFT).

Diskrétní Fourierova transformace mezi posloupnostmi \{ d(k) \}_{k=0}^{N-1}, \{ D(n) \}_{n=0}^{N-1}, je definována vztahy:

* přímá diskrétní Fourierova transformace :D(n)=\sum_{k=0}^{N-1} d(k){\mathrm e}^{-{\mathrm i} nk2\pi/N}, n=0,...,N-1

* a zpětná (inverzní) diskrétní Fourierova transformace :d(k)=\frac{1}{N}\sum_{n=0}^{N-1} D(n){\mathrm e}^{{\mathrm i} nk2\pi/N}, k=0,...,N-1

Diskrétní Fourierova transformace našla velké uplatnění zejména s rozvojem výpočetní techniky. Součástí řady přístrojů jsou jednoúčelové procesory realizující tuto transformaci. +more Výpočet DFT podle definičního vztahu vyžaduje N^2 komplexních součinů a N^2 komplexních součtů. Toto množství operací výrazně snižovalo možnost aplikace DFT na výpočty v reálném čase.

Situace se změnila po roce 1965, kdy J. W. +more Cooley a J. W. Tukey popsali velmi efektivní algoritmus výpočtu DFT, tzv. rychlou Fourierovu transformaci (FFT - Fast Fourier Transform), který vyžaduje jen N\log_2 N komplexních součinů a N\log_2 N komplexních součtů. Díky tomuto algoritmu se stala diskrétní Fourierova transformace nejrozšířenějším prostředkem pro numerický výpočet Fourierovy transformace. Algoritmus FFT je také implementován ve všech nejrozšířenějších matematických programech jako je např. GNU Octave, Mathcad, Mathematica, Maple, Matlab atd.

Diskrétní Fourierova transformace se dá vyjádřit jako násobení vektoru d_n Vandermondovou maticí, která která obsahuje komplexní odmocniny z 1.

Zpětná Fourierova transformace

:f(t) = \frac{1}{2 \pi} \int\limits_{-\infty}^{\infty} S (\omega){\mathrm e}^{{\mathrm i}\omega t}\,{\mathrm d}\omega

Integrál je nutno chápat ve smyslu hlavní hodnoty. Popisuje význam F. +more t. - rozklad obecného (i neperiodického) signálu f(t) na směs harmonických kmitů, jejichž amplituda je daná amplitudovým spektrem a fázový posun fázovým spektrem.

Související články

Fourierova řada * Diskrétní kosinová transformace * Rychlá Fourierova transformace

Reference

Externí odkazy

Václav Hlaváč: [url=https://web. archive. +moreorg/web/20070221224422/http://cmp. felk. cvut. cz/~hlavac/Public/TeachingLectures/p2LinIntTxFourier. pdf]LINEÁRNÍ INTEGRÁLNÍ TRANSFORMACE - Fourierova transformace[/url], Fakulta elektrotechnická ČVUT v Praze, katedra kybernetiky, Centrum strojového vnímání * [url=http://matematika. cuni. cz/dl/analyza/37-fou/lekce37-fou-pmax. pdf]FOURIEROVA TRANSFORMACE[/url] M. Hušek, P. Pyrih et al. , Matematicko-fyzikální fakulta, Univerzita Karlova v Praze * [url=http://ivankuckir. blogspot. cz/2013/11/fourierova-transformace-srozumitelne. html]Fourierova transformace - názorné vysvětlení[/url] * [url=https://tomasboril. cz/fourierseries3d/cz/]Fourier Series 3D - interaktivní demonstrace principu Fourierových řad[/url] HTML5 a JavaScript: Unikátní interaktivní 3D zobrazení propojující časovou, frekvenční, amplitudovou a fázovou osu.

Kategorie:Matematická analýza Kategorie:Zpracování signálu Kategorie:Integrální transformace