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.

Fourier-sor (2π szerint periodikus esetben)

f(x)a02+n=1(ancosnx+bnsinnx)f(x)\sim\frac{a_{0}}{2}+\sum_{n=1}^{\infty}\Big(a_{n}\cos nx+b_{n}\sin nx\Big)

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 sinnx\sin nx és cosnx\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 ππsinnxcosmxdx=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.

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)cosnxdx,bn=1πππf(x)sinnxdxa_{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 cos0=1\cos 0=1:

a0=1πππf(x)dxa02=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

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 xsinnxdx\int x\sin nx\,dx típusú integrálok jönnek ki. Érdemes fejben tartani két gyakran előkerülő értéket:

cosnπ=(1)n,sinnπ=0(nZ)\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)sinnxf(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)cosnxdx,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)sinnxdxa_{n}=0,\qquad b_{n}=\frac{2}{\pi}\int_{0}^{\pi}f(x)\sin nx\,dx

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=sinnxv'=\sin nx választással:
    0πxsinnxdx=[xcosnxn]0π+1n0πcosnxdx\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 [sinnx]0π=0\big[\sin nx\big]_{0}^{\pi}=0, marad:
    0πxsinnxdx=πcosnπ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.
    x2n=1(1)n+1nsinnx=2(sinxsin2x2+sin3x3)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 (sinnπ=0\sin n\pi=0), és épp a Leibniz-sor marad:

π4=113+1517+\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=1LLLf(x)cosnπxLdx,bn=1LLLf(x)sinnπxLdxa_{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(ancosnπxL+bnsinnπ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.

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 sinnπ=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=11n2=π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 cosnπ\cos n\pi kiírásának elhagyása. A végeredményben ne maradjon cosnπ\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.