Furijeov red je matematička operacija kojom se periodična funkcija razlaže na svoje „spektralne komponente“ radi jednostavnije analize. Nekoliko prvih članova takvog razvoja se u tehnici često uzimaju kao veoma korisna vrsta aproksimacije.
Diskretna furijeova transformacija pretvara diskretne vrednosti (vektor) u Furijeove koeficijente. Neprekidna furijeova transformacija radi to isto sa funkcijom. Naziv je dobila po francuskom matematičaru Žozefu Furijeu (1768—1830).
Matematička osnova
Uzmimo neku periodičnu funkciju
sa periodom T, za koju važi
. Zbog periodičnosti možemo da je razdelimo na N sinus i kosinus funkcija:
,
, gde je
osnovna frekvencija, odnosno harmonik.
Treba imati na umu da je sinus samo kosinus sa faznim pomerajem:
![{\displaystyle f(t)=\sum _{n=0}^{N}A_{n}\cos(n\omega t+\varphi _{n})=A_{0}+\sum _{n=1}^{N}(A_{n}\cos \varphi _{n}\cdot \cos(n\omega t)-A_{n}\sin \varphi _{n}\cdot \sin(n\omega t))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/711867aac84e4c656023db870e802224f0d5f07d)
Kada definišemo
, a potom
i
dobijamo isti izraz, ovog puta bez faze:
![{\displaystyle f(t)=a_{0}+\sum _{n=1}^{N}(a_{n}\cos(n\omega t)-b_{n}\sin(n\omega t)).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/268cb86fa6f2da16b94fda2620a664d612bb39e2)
Zašto se ne uzima tan ili recimo cosh? Zašto baš cos i sin? Razlog je ortogonalnost sin i cos funkcija.
Ideja iza furijeove transformacije je sledeća: ceo prostor koji ima „normalne“ ose transformišemo u prostor u kome su nove ortogonalne ose kosinus i sinus talasi i njihovi viši harmonici. Signal koji transformišemo je samo jedna tačka (mesni vektor), a vrednosti na svakoj osi su amplitude svakog harmonika pojedinačno (
).
Sada se uključuje Ojlerov identitet uz pomoć koga ove trigonometrijske funkcije možemo da zamenimo kompleksnim pandanima:
i ![{\displaystyle \sin(x)={\frac {1}{2\mathrm {i} }}\left(e^{\mathrm {i} x}-e^{-\mathrm {i} x}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e83cc2e0d09ad21d7755d331bb9fc3dd7329edc)
Iz toga dalje sledi
![{\displaystyle f(t)=a_{0}+\sum _{n=1}^{N}{\frac {1}{2}}\left(a_{n}(e^{\mathrm {i} n\omega t}+e^{-\mathrm {i} n\omega t})-{1 \over \mathrm {i} }b_{n}(e^{\mathrm {i} n\omega t}-e^{-\mathrm {i} n\omega t})\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c1294eb5e20457830dfcde2ba90f1de9b945ead)
![{\displaystyle =a_{0}+\sum _{n=1}^{N}{\frac {1}{2}}\left(a_{n}(e^{\mathrm {i} n\omega t}+e^{-\mathrm {i} n\omega t})+\mathrm {i} b_{n}(e^{\mathrm {i} n\omega t}-e^{-\mathrm {i} n\omega t})\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/272cc89015ab8472e5df74e717836ee6e7394296)
![{\displaystyle =a_{0}+\sum _{n=1}^{N}{\frac {1}{2}}\left((a_{n}+\mathrm {i} b_{n})e^{\mathrm {i} n\omega t}+(a_{n}-\mathrm {i} b_{n})e^{-\mathrm {i} n\omega t}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c4dea0ed583d173006db67071293ec023a2e065a)
Zamenimo realne koeficijente kompleksnim:
,
i ![{\displaystyle c_{-n}:={\frac {1}{2}}(a_{n}-\mathrm {i} b_{n})={\overline {c_{n}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f40ed7f148a63891410274358b0690be656ce41)
dobijamo sumu sa negativnim indeksima:
![{\displaystyle f(t)=\sum _{k=-N}^{N}c_{k}e^{\mathrm {i} k\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0fbae59e739a7091c27f2bcf945bb08c9094cfa8)
Takođe, ne treba gubiti iz vida da su
funkcije isto ortonormalne baze (svaki vektor koji predstavlja osu ima dužinu 1 i normalan je u odnosu na sve ostale vektore):
U slučaju
![{\displaystyle (e^{\mathrm {i} jt},e^{\mathrm {i} jt})={\frac {1}{2\pi }}\int _{0}^{2\pi }e^{\mathrm {i} jt}{\overline {e^{\mathrm {i} jt}}}dt={\frac {1}{2\pi }}\int _{0}^{2\pi }e^{\mathrm {i} jt}e^{-\mathrm {i} jt}dt={\frac {1}{2\pi }}\int _{0}^{2\pi }1dt=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6444d4f5cdbaff534f035cf39c0517ee7a3af611)
A za
važi:
![{\displaystyle (e^{\mathrm {i} jt},e^{\mathrm {i} kt})={\frac {1}{2\pi }}\int _{0}^{2\pi }e^{\mathrm {i} jt}{\overline {e^{\mathrm {i} kt}}}dt={\frac {1}{2\pi }}\int _{0}^{2\pi }e^{\mathrm {i} jt}e^{-\mathrm {i} kt}dt={\frac {1}{2\pi }}\int _{0}^{2\pi }e^{\mathrm {i} (j-k)t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94239daaa450bcc9328745bec61d4fa98d8f230c)
![{\displaystyle ={\frac {1}{2\pi \mathrm {i} (j-k)}}\left[e^{\mathrm {i} (j-k)t}\right]_{0}^{2\pi }={\frac {1}{2\pi \mathrm {i} (j-k)}}\left(e^{\mathrm {i} (j-k)2\pi }-1\right)=}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e1a31212b6347a8915e289b603978131e74f6309)
![{\displaystyle ={\frac {1}{2\pi \mathrm {i} (j-k)}}\cdot 0=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d6d705e96e22a1eefbbe3cf9ecef03297e4edf82)
Furijeovi redovi
No, želimo sada da neku periodičnu i neprekidnu funkciju približno izračunamo uz pomoć sume trigonometrijskih funkcija (konkretno: kosinusa i sinusa). Videli smo kako možemo da dođemo do
; gornju jednačinu množimo sa
i naposletku integrišemo sa obe strane po intervalu [0,T] odnosno u trajanju jedne periode:
![{\displaystyle e^{-\mathrm {i} m\omega t}f(t)=\sum _{n=-N}^{N}c_{n}\left(e^{\mathrm {i} (n\omega t)}e^{-\mathrm {i} m\omega t}\right)=\sum _{n=-N-m}^{N-m}c_{n+m}e^{\mathrm {i} (n+m)\omega t-\mathrm {i} m\omega t}=\sum _{n=-N-m}^{N-m}c_{n+m}e^{\mathrm {i} n\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7d336bf6775da1222c602a932efe97855c9f4e9c)
![{\displaystyle \Leftrightarrow \int _{0}^{T}e^{-\mathrm {i} m\omega t}f(t)dt=\sum _{n=-N-m}^{N-m}c_{n+m}\int _{0}^{T}e^{\mathrm {i} n\omega t}dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f025631e3496eaafada0ba6150eb3e956cce61a)
Za integrale sa desne strane važi:
- kada je n=0:
![{\displaystyle \int _{0}^{T}e^{\mathrm {i} 0\omega t}dt=\int _{0}^{T}e^{0}=\left[1\right]_{0}^{T}=T}](https://wikimedia.org/api/rest_v1/media/math/render/svg/32663cf201e76d795deb90f5406895c8b3f689c3)
- a kada je n≠0:
![{\displaystyle ={\frac {1}{\mathrm {i} n\omega }}(e^{\mathrm {i} n\omega T}-1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e76a96023bf724a04ad026c7d19b4df4353aaeba)
Iz
sledi
, a to dalje možemo da primenimo na gore navedeni integral:
![{\displaystyle \int _{0}^{T}e^{\mathrm {i} n\omega t}dt=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5489cf53435db3f8167eec8573eda9270d60a10e)
Na kraju se cela računica uprošćava:
![{\displaystyle \int _{0}^{T}f(t)e^{-\mathrm {i} m\omega t}dt=\sum _{n=-N-m}^{N-m}c_{n+m}\int _{0}^{T}e^{\mathrm {i} n\omega t}dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6869dae79c9ba57cd900d4c2b11829c0cf3225a5)
![{\displaystyle =\sum _{n=-N-m}^{-1}c_{n+m}\int _{0}^{T}e^{\mathrm {i} n\omega t}dt+c_{m}\cdot \int _{0}^{T}e^{\mathrm {i} 0\omega t}dt+\sum _{n=1}^{N-m}c_{n+m}\int _{0}^{T}e^{\mathrm {i} n\omega t}dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/699595e32e8a3b64da20954a26d59630f0d2b779)
![{\displaystyle =0+c_{m}T+0=c_{m}T=\int _{0}^{T}f(t)e^{-\mathrm {i} m\omega t}dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e121585b575374901293645d0597bf2b9dfbbd5c)
U celom računu neka nas ne zbunjuje korišćenje promenljive
, njena svrha je puko uprošćavanje jednačine. Sve je stoga samo dosetljivost, odnosno umetnost kako napisati jedno te isto na drugačiji način.
Na kraju, Furijeov red definišemo:
![{\displaystyle f_{N}(t)=\sum _{n=-N}^{N}c_{n}e^{in\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f3af30edd19caefdcc32f989b94b20637a86f3c3)
Konvergentnost Furijeovog reda
Furijeov red konvergira ka mnogim funkcijama; tu spadaju pored ostalih sve funkcije koje imaju izvod ili su kvadratno integrabilne (L2 prostor).
Pretpostavimo da je
jedna takva funkcija. Kada namestimo
, onda ona takođe može da se napiše i ovako:
![{\displaystyle f(t)=\sum _{n=-N}^{N}c_{n}e^{\mathrm {i} n\omega t}=\sum _{n=-N}^{N}{\frac {1}{T}}(\int _{0}^{T}f(t)e^{-\mathrm {i} n\omega t}dt)\cdot e^{\mathrm {i} n\omega t}=\sum _{n=-N}^{N}{\frac {e^{\mathrm {i} n\omega t}}{T}}\cdot \int _{0}^{T}f(t)e^{-\mathrm {i} n\omega t}dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b27e501bca993440e21d4943b229d878e64ced61)
![{\displaystyle f(t)=\sum _{n=-\infty }^{\infty }{\frac {e^{\mathrm {i} n\omega t}}{T}}\cdot \int _{0}^{T}f(t)e^{-\mathrm {i} n\omega t}dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9beff70d43addbeb20d0d66165ad357c59cc7513)
Povezano
Literatura
- William E. Boyce and Richard C. DiPrima (2005). Elementary Differential Equations and Boundary Value Problems (8th izd.). New Jersey: John Wiley & Sons, Inc.. ISBN 978-0-471-43338-5.
- Joseph Fourier, translated by Alexander Freeman (published 1822, translated 1878, re-released 2003). The Analytical Theory of Heat. Dover Publications. ISBN 978-0-486-49531-6. 2003 unabridged republication of the 1878 English translation by Alexander Freeman of Fourier's work Théorie Analytique de la Chaleur, originally published in 1822.
- Enrique A. Gonzalez-Velasco (1992). „Connections in Mathematical Analysis: The Case of Fourier Series”. American Mathematical Monthly 99 (5): 427-441. DOI:10.2307/2325087.
- Katznelson, Yitzhak (1976). An introduction to harmonic analysis (Second corrected izd.). New York: Dover Publications, Inc. ISBN 978-0-486-63331-2.
- Felix Klein, Development of mathematics in the 19th century. Mathsci Press Brookline, Mass, 1979. Translated by M. Ackerman from Vorlesungen über die Entwicklung der Mathematik im 19 Jahrhundert, Springer, Berlin, 1928.
- Walter Rudin (1976). Principles of mathematical analysis (3rd izd.). New York: McGraw-Hill, Inc.. ISBN 978-0-07-054235-8.
- A. Zygmund (2002). Trigonometric series (third izd.). Cambridge: Cambridge University Press. ISBN 978-0-521-89053-3. The first edition was published in 1935.