Ingyenes tananyagok

Analízis 3

Fourier-sorok

Egy periodikus függvény felbontása szinuszok és koszinuszok összegére: honnan jönnek az együtthatók, mit egyszerűsít a szimmetria, és hova konvergál a sor.

15 perc olvasás

1. Az alapötlet

A Taylor-sornál egy függvényt hatványfüggvények összegére bontottunk, és a közelítés egyetlen pont körül volt jó. Periodikus jeleknél ez rossz üzlet: egy négyszögjelet vagy egy fűrészfogat polinomokkal reménytelen leírni. Kézenfekvő olyan építőkockákat választani, amik maguk is periodikusak — szinuszokat és koszinuszokat. Ebben a leckében végigvezetlek azon, hogyan találjuk meg ezeknek az építőkockáknak a súlyát, és hol lehet a számolás felét megspórolni.

Fourier-sor (2π szerint periodikus esetben)

f(x)∼a02+∑n=1∞(ancos⁡nx+bnsin⁡nx)f(x)\sim\frac{a_{0}}{2}+\sum_{n=1}^{\infty}\Big(a_{n}\cos nx+b_{n}\sin nx\Big)

Egy egyszerű példa: ami már eleve ilyen alakú

Az f(x)=3+2cos⁡x−sin⁡2xf(x)=3+2\cos x-\sin 2x függvény maga is egy (véges) Fourier-sor. Olvasd le az együtthatóit: a02=3\frac{a_0}{2}=3 (tehát a0=6a_0=6), a1=2a_1=2, b2=−1b_2=-1, az összes többi nulla. A Fourier-sorfejtés célja pontosan ez: egy tetszőleges periodikus függvényről kideríteni, mekkora súllyal szerepel benne az egyes szinuszok és koszinuszok. Itt ránézésre látszik, a következő fejezetben megtanuljuk kiszámolni.

A a02\frac{a_0}{2} tag a jel átlagos értéke (az „egyenkomponens”), az n=1n=1 tag az alaphang, a többi pedig a felharmonikusok. Ez a fizikai olvasat: a Fourier-sor megmondja, milyen frekvenciákból és milyen arányban áll össze a jel.

Miért működik: ortogonalitás

A sin⁡nx\sin nx és cos⁡nx\cos nx függvények ortogonális rendszert alkotnak: bármely kettőnek a szorzatát integrálva egy teljes perióduson nullát kapunk, kivéve ha ugyanarról a függvényről van szó. Például ∫−ππsin⁡nxcos⁡mx dx=0\int_{-\pi}^{\pi}\sin nx\cos mx\,dx=0 minden n,mn,m-re. Épp ez teszi lehetővé, hogy az együtthatókat egyenként „kihalásszuk”: beszorzunk azzal a taggal, amelyiknek az együtthatóját keressük, és integrálunk — a többi tag magától eltűnik.

Elakadtál ennél a résznél?

Egy órán megnézzük, pontosan hol csúszik el, és addig gyakoroljuk, amíg magabiztos nem leszel benne.

2. A Fourier-együtthatók

Az előző bekezdés gondolatát végigvíve adódnak az Euler–Fourier-képletek. Ezt a hármat kell tudni:

Együtthatók 2π periódusra

an=1π∫−ππf(x)cos⁡nx dx,bn=1π∫−ππf(x)sin⁡nx dxa_{n}=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\cos nx\,dx,\qquad b_{n}=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin nx\,dx

Az a0a_0 ugyanennek a képletnek az n=0n=0 esete, hiszen cos⁡0=1\cos 0=1:

a0=1π∫−ππf(x) dx⟹a02=12π∫−ππf(x) dxa_{0}=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\,dx\qquad\Longrightarrow\qquad \frac{a_{0}}{2}=\frac{1}{2\pi}\int_{-\pi}^{\pi}f(x)\,dx

Egy egyszerű példa: a konstans függvény

Legyen f(x)=1f(x)=1. Ekkor a0=1π∫−ππ1 dx=2ππ=2a_0=\frac1\pi\int_{-\pi}^{\pi}1\,dx=\frac{2\pi}{\pi}=2, és n≥1n\ge1-re an=1π∫−ππcos⁡nx dx=0a_n=\frac1\pi\int_{-\pi}^{\pi}\cos nx\,dx=0, bn=1π∫−ππsin⁡nx dx=0b_n=\frac1\pi\int_{-\pi}^{\pi}\sin nx\,dx=0, mert egy teljes periódusra integrált szinusz vagy koszinusz nulla. A sor tehát a02=1\frac{a_0}{2}=1: pontosan a függvény, ahogy kell. Figyeld meg, hogy a felezés nélkül 2-t kaptunk volna, ami nyilván rossz: ez a legegyszerűbb ellenőrzés az a02\frac{a_0}{2}-re.

Az a₀ és a felezés

Sok tankönyv azért írja a sor elejére a02\frac{a_0}{2}-t, hogy az ana_n képlete minden nn-re ugyanaz lehessen, a nullát is beleértve. Van olyan jegyzet is, ahol a sor elején sima a0a_0 áll, és akkor annak a képletében szerepel a 12π\frac{1}{2\pi}. A kettő ugyanaz — de a vizsgán azt a konvenciót kövesd, amit az előadáson használtak, és a felezést egyszer se hagyd le.

Az együtthatók számolása kivétel nélkül parciális integrálás — a legtöbb feladatban ∫xsin⁡nx dx\int x\sin nx\,dx típusú integrálok jönnek ki. Érdemes fejben tartani két gyakran előkerülő értéket:

cos⁡nπ=(−1)n,sin⁡nπ=0(n∈Z)\cos n\pi=(-1)^{n},\qquad \sin n\pi=0\quad (n\in\mathbb{Z})

3. Páros és páratlan függvények

Ez az a fejezet, ami a dolgozatban a legtöbb időt spórolja. Ha a függvénynek van paritása, az együtthatók fele automatikusan nulla, és nem kell kiszámolni őket — elég egy mondattal indokolni.

  • Páros függvény (f(−x)=f(x)f(-x)=f(x)): minden bn=0b_n=0, csak koszinuszok maradnak. Ezt hívják koszinuszsornak.
  • Páratlan függvény (f(−x)=−f(x)f(-x)=-f(x)): minden an=0a_n=0, beleértve az a0a_0-t is; csak szinuszok maradnak. Szinuszsor.

Az indoklás egyszerű: páratlan függvényt szimmetrikus tartományon integrálva nullát kapunk. Ha ff páros, akkor f(x)sin⁡nxf(x)\sin nx páratlan (páros szor páratlan), így az integrálja nulla. Ha megvan a paritás, a megmaradó együtthatókat is elég a fél tartományon számolni, kettővel szorozva:

Páros f esetén

an=2π∫0πf(x)cos⁡nx dx,bn=0a_{n}=\frac{2}{\pi}\int_{0}^{\pi}f(x)\cos nx\,dx,\qquad b_{n}=0

Páratlan f esetén

an=0,bn=2π∫0πf(x)sin⁡nx dxa_{n}=0,\qquad b_{n}=\frac{2}{\pi}\int_{0}^{\pi}f(x)\sin nx\,dx

Egy egyszerű példa: a paritás leolvasása

Az f(x)=x2f(x)=x^2 páros, tehát minden bn=0b_n=0, és csak koszinuszsor lesz. Az f(x)=x3f(x)=x^3 páratlan: minden an=0a_n=0, az a0a_0-lal együtt, csak szinuszok maradnak. Az f(x)=x+1f(x)=x+1 viszont egyik sem, ott mind a két fajta együtthatót ki kell számolni. Ezt a mondatot a dolgozatban le is kell írni, mert ez indokolja, miért maradt ki a számolás fele.

4. Kidolgozott példa: a fűrészfog

Fejtsük Fourier-sorba az f(x) = x függvényt a (−π, π) intervallumon

  1. Paritás. Az f(x)=xf(x)=x páratlan, tehát minden an=0a_n=0. Csak a bnb_n együtthatókat kell számolni.
  2. Az integrál. Parciális integrálás u=xu=x, v′=sin⁡nxv'=\sin nx választással:
    ∫0πxsin⁡nx dx=[−xcos⁡nxn]0π+1n∫0πcos⁡nx dx\int_{0}^{\pi}x\sin nx\,dx=\left[-\frac{x\cos nx}{n}\right]_{0}^{\pi}+\frac{1}{n}\int_{0}^{\pi}\cos nx\,dx
    A második tag nulla, mert [sin⁡nx]0π=0\big[\sin nx\big]_{0}^{\pi}=0, marad:
    ∫0πxsin⁡nx dx=−πcos⁡nπn=π(−1)n+1n\int_{0}^{\pi}x\sin nx\,dx=-\frac{\pi\cos n\pi}{n}=\frac{\pi(-1)^{n+1}}{n}
  3. Az együttható.
    bn=2π⋅π(−1)n+1n=2(−1)n+1nb_{n}=\frac{2}{\pi}\cdot\frac{\pi(-1)^{n+1}}{n}=\frac{2(-1)^{n+1}}{n}
  4. A sor.
    x∼2∑n=1∞(−1)n+1nsin⁡nx=2(sin⁡x−sin⁡2x2+sin⁡3x3−… )x\sim 2\sum_{n=1}^{\infty}\frac{(-1)^{n+1}}{n}\sin nx=2\left(\sin x-\frac{\sin 2x}{2}+\frac{\sin 3x}{3}-\dots\right)

Ajándék: a Leibniz-sor

Helyettesíts be x=π2x=\frac{\pi}{2}-t. A bal oldal π2\frac{\pi}{2}, a jobb oldalon a páros indexű tagok kiesnek (sin⁡nπ=0\sin n\pi=0), és épp a Leibniz-sor marad:

π4=1−13+15−17+…\frac{\pi}{4}=1-\frac13+\frac15-\frac17+\dots

A Fourier-sorok egyik szép mellékterméke, hogy nevezetes numerikus sorok összegét lehet velük megkapni — olyanokét is, amiket másképp nagyon nehéz.

5. Tetszőleges periódus

Ha a periódus nem 2π2\pi, hanem 2L2L, akkor csak át kell skálázni a változót: az nxnx helyére nπxL\frac{n\pi x}{L} kerül, az 1π\frac1\pi szorzó helyére 1L\frac1L:

Együtthatók 2L periódusra

an=1L∫−LLf(x)cos⁡nπxL dx,bn=1L∫−LLf(x)sin⁡nπxL dxa_{n}=\frac{1}{L}\int_{-L}^{L}f(x)\cos\frac{n\pi x}{L}\,dx,\qquad b_{n}=\frac{1}{L}\int_{-L}^{L}f(x)\sin\frac{n\pi x}{L}\,dx
f(x)∼a02+∑n=1∞(ancos⁡nπxL+bnsin⁡nπxL)f(x)\sim\frac{a_{0}}{2}+\sum_{n=1}^{\infty}\left(a_{n}\cos\frac{n\pi x}{L}+b_{n}\sin\frac{n\pi x}{L}\right)

A 2π2\pi-s képletek ennek az L=πL=\pi esetei, tehát elég ezt az egy alakot megjegyezni. A számolás menete szóról szóra ugyanaz, a paritásos rövidítések is változatlanul érvényesek.

Egy egyszerű példa: a fűrészfog 2 periódussal

Legyen f(x)=xf(x)=x a (−1;1)(-1;1) intervallumon, tehát L=1L=1. Páratlan, ezért csak bnb_n kell: bn=2∫01xsin⁡(nπx) dxb_n=2\int_0^1 x\sin(n\pi x)\,dx. Ugyanaz a parciális integrálás, mint az előbb, csak a nevezőben nπn\pi áll nn helyett: bn=2(−1)n+1nπb_n=\dfrac{2(-1)^{n+1}}{n\pi}. A sor: x∼2π(sin⁡πx−sin⁡2πx2+… )x\sim\dfrac{2}{\pi}\left(\sin\pi x-\dfrac{\sin2\pi x}{2}+\dots\right). Hasonlítsd össze a 2π2\pi-s változattal: a szinuszok argumentumába π\pi került, és az együtthatók π\pi-vel osztódtak.

6. Hova konvergál a sor?

A sor és a függvény között eddig ∼\sim jel állt, nem egyenlőség. Hogy mikor lehet egyenlőséget írni, arról a Dirichlet-féle konvergenciatétel szól.

Dirichlet-tétel

a sor o¨sszege=f(x−)+f(x+)2\text{a sor összege}=\frac{f(x^{-})+f(x^{+})}{2}

Ha ff szakaszonként folytonos és szakaszonként monoton (ez a szokásos feltétel, és minden vizsgapéldában teljesül), akkor a Fourier-sor minden pontban konvergens, és az összege:

  • ott, ahol a függvény folytonos: maga az f(x)f(x) érték;
  • a szakadási helyeken: a bal és jobb oldali határérték számtani közepe — vagyis a sor „megfelezi” az ugrást, akkor is, ha a függvény ott mást vesz fel.

A fenti példában x=πx=\pi-nél a periodikus kiterjesztés π\pi-ről −π-\pi-re ugrik. A sor összege itt π+(−π)2=0\frac{\pi+(-\pi)}{2}=0 — és tényleg: minden tag sin⁡nπ=0\sin n\pi=0.

A Gibbs-jelenség

Egy ugrás közelében a részletösszegek mindig „túllőnek”: a grafikonon egy kis púp jelenik meg, ami a tagszám növelésével nem tűnik el, csak egyre közelebb kúszik a szakadáshoz. Ez a Gibbs-jelenség; nem számolási hiba, hanem a szakaszonkénti konvergencia természetes velejárója.

Parseval-formula

A jel „energiája” is kiolvasható az együtthatókból:

1π∫−ππf2(x) dx=a022+∑n=1∞(an2+bn2)\frac{1}{\pi}\int_{-\pi}^{\pi}f^{2}(x)\,dx=\frac{a_{0}^{2}}{2}+\sum_{n=1}^{\infty}\big(a_{n}^{2}+b_{n}^{2}\big)

A fűrészfogas példára alkalmazva a bal oldal 2π23\frac{2\pi^{2}}{3}, a jobb oldal pedig ∑4n2\sum\frac{4}{n^{2}}. Innen egy sorral adódik a híres összeg:

∑n=1∞1n2=π26\sum_{n=1}^{\infty}\frac{1}{n^{2}}=\frac{\pi^{2}}{6}

7. Tipikus hibák

  • Az a0a_0 felezésének elhagyása. A sor első tagja a02\frac{a_0}{2}, nem a0a_0. Ez a leggyakoribb pontvesztés, és a végeredmény ránézésre helyesnek tűnik.
  • A paritás kihasználásának elmulasztása. Ha a függvény páros vagy páratlan, az együtthatók fele nulla. Ezt nem kiszámolni kell, hanem leírni egy mondatban — és a maradék munka a felére csökken.
  • A cos⁡nπ\cos n\pi kiírásának elhagyása. A végeredményben ne maradjon cos⁡nπ\cos n\pi, írd át (−1)n(-1)^{n} alakba; enélkül a sor jellege nem látszik.
  • Egyenlőség a szakadási pontokban. A sor ott az ugrás felezőpontjához konvergál, nem a függvényértékhez. A jelölésben ezért szerepel ∼\sim és nem ==.
  • Rossz periódusú képlet. Ha az intervallum nem [−π,π][-\pi,\pi], ne a 1π\frac1\pi-s képletet használd. Írd fel először, mennyi LL, és onnan indulj.

Elakadtál menet közben?

Egy tananyag megmutatja, hogyan működik a dolog. Azt viszont, hogy pontosan hol csúszik el nálad, egy óra alatt derítjük ki. Foglalj egy alkalmat, vagy beszéljük meg egy ingyenes konzultáción, mire van szükséged.