2. Fourier-raðir

2.1. Inngangur

Rifjum upp að ef við höfum grunn af vigrum x1,x2,,xn í Rn má rita sérhvern vigur y sem

y=a1x1++anxn

þar sem stuðlarnir a1,,an eru ótvírætt ákvarðaðir. Við munum nú spyrja okkur spurningarinnar, er hægt að gera eitthvað sambærilegt þegar vigurrúmið er óendanlega vítt, t.d. þegar það samanstendur af föllum. Við þekkjum dæmi um slíkt, þegar rita má óendanlega oft diffranlegt fall f með Taylor-röð þess

f(x)=a01+a1x+a2x2+,þar seman=f(n)(0)n!.

Í þessum kafla munum við skilgreina svokallaðar Fourier-raðir sem líta svipað út en í stað fallanna 1,x,x2, munum við liða lotubundin föll f í grunn sem samanstendur af hornaföllum (eða jafngilt, veldisvísisföllum), finna formúlur fyrir stuðlunum í framsetningunni og loks skoða hvernig má nota raðirnar við lausn hlutafleiðujafna.

2.1.1. Skilgreining

Fall f:RR er sagt vera T-lotubundið ef f(x+T)=f(x) fyrir öll xR.

2.2. Fourier-raðir

Við munum skoða tilfellið þegar föllin eru lotubundin með lotu 2π og sjá hvaða skilyrði tryggja að hægt sé að liða slík föll í grunn sem samanstendur af hornaföllunum sin(nx) og cos(nx) annars vegar eða enix hins vegar, þar sem n0.

2.2.1. Skilgreining

Látum IR vera bil.

  1. Rúmið L1(I) er mengi þeirra falla f:IC þannig að
I|f(x)|dx<.
  1. Rúmið L2(I) er mengi þeirra falla f:IC þannig að
I|f(x)|2dx<.

Ef f og g eru föll í L2(I) kallast

f,g=1|I|If(x)¯g(x)dx

innfeldi þeirra (misjafnt er hvort deilt er með |I|, lengdinni á I, í skilgreiningunni). Ef f,g=0 segjum við að f og g séu hornrétt.

Athugið

Lj(I), j=1,2 eru vigurrúm, af því að

  1. Ef fLj(I) og gLj(I) þá er fallið f+gLj(I)
  2. Ef fLj(I) þá er αfLj(I), þar sem αR

2.2.2. Skilgreining

Ef fL1([π,π]) er 2π-lotubundið þá skilgreinum við Fourier-stuðla þess með

cn=cn(f)=12πππeinxf(x)dx,n=,2,1,0,1,2,,

Fourier-kósínus-stuðla f með

an=an(f)=1πππf(x)cos(nx)dx,n=0,1,2,,

og Fourier-sínus-stuðla f með

bn=bn(f)=1πππf(x)sin(nx)dx,n=1,2,.

Raðirnar

a02+n1(ancos(nx)+bnsin(nx))ogn=cneinx

kallast Fourier-raðir f og til aðgreiningar er sú fyrri oft nefnd hornafallaröð f.



2π-lotubundna fallið er skilgreint með því að gefa formúlu fyrir því á bilinu [0,2π].



Athugið

Þegar T-lotubundið fall er heildað yfir eina lotu skiptir ekki máli hvar upphafspunktur heildisins er valinn, þ.e.

T/2T/2f(x)dx=T0f(x)dx=α+Tαf(x)dx,fyrir öll αR.

2.2.3. Setning - Reiknireglur

Látum f,gL1([π,π]) vera 2π-lotubundin föll.

  1. Fourier-stuðlarnir eru línulegar varpanir á L1([π,π]),
an(αf+βg)=αan(f)+βan(g)bn(αf+βg)=αbn(f)+βbn(g)cn(αf+βg)=αcn(f)+βcn(g)
  1. Eftirfarandi samband gildir
a0=2c0,an=cn+cn,bn=i(cncn),c0=a02,cn=12(anibn),cn=12(an+ibn),n1.
  1. Ef g(x)=f(x+α) , þar sem αR þá er cn(g)=einαcn(f) fyrir öll n=0,±1,±2,.
  2. Ef f er raungilt fall þá eru an(f) og bn(f) rauntölur og cn(f)=¯cn(f).
  3. Ef f er jafnstætt fall þá er bn(f)=0 fyrir öll n=1,2,3, og
an(f)=2ππ0f(x)cos(nx)dx.

6 Ef f er oddstætt fall þá er an(f)=0 fyrir öll n=0,1,2, og

bn(f)=2ππ0f(x)sin(nx)dx.
  1. Ef f,f,,f(m) eru í L1([π,π]) þá er
cn(f(k))=(in)kcn(f),0km,nZ.

2.2.4. Skilgreining

Ef fL1([T/2,T/2]) er T-lotubundið þá setjum við ω=2π/T og skilgreinum Fourier-stuðla þess með

cn=cn(f)=1TT/2T/2einωxf(x)dx,n=,2,1,0,1,2,,

Fourier-kósínus-stuðla f með

an=an(f)=2TT/2T/2f(x)cos(nωx)dx,n=0,1,2,,

og Fourier-sínus-stuðla f með

bn=bn(f)=2TT/2T/2f(x)sin(nωx)dx,n=1,2,.

Raðirnar

a02+n1(ancos(nωx)+bnsin(nωx))ogn=cneinωx

kallast Fourier-raðir f og til aðgreiningar er sú fyrri oft nefnd hornafallaröð f.

Athugið

Sambærilegar reiknireglur fyrir T-lotubundin fást út frá reglunum fyrir 2π-lotubundin föll, með því að “skipta 2π út fyrir T ” á viðeigandi stöðum.

2.3. Samleitni Fourier-raða

Í þessari grein fjöllum við um skilyrði sem tryggja samleitni Fourier-raða falls og hvenær og í hvaða skilningi fallið er jafnt Fourier-röð sinni. Við munum notast talsvert við innfeldið sem skilgreint er á L2([π,π]) og setjum því fram nokkrar reiknireglur um innfeldi

2.3.1. Reiknireglur um innfeldi

Ef u,v,wL2([π,π]) og α,βC þá gilda eftirfarandi reiknireglur

αu+βv,w=αu,w+βv,w,u,αv+βw=¯αu,v+¯βu,w,u,v=¯v,u,u,u0.

Síðasta reglan leyfir okkur að skilgreina lengd eða staðal fallsins u sem

u=u,u.

Ein mikilvægasta ójafna stærðfræðinnar er Cauchy-Schwarz ójafnan

2.3.2. Cauchy-Schwarz ójafna

Fyrir u,vL2([π,π]) gildir

|u,v|12πππ|u(x)v(x)|dxuv.

Athugum nú að föllin einx og eimx eru hornrétt ef nm því þá gildir

einx,eimx=12πππe(nm)ixdx=[e(nm)ixi(nm)]ππ=0.

Ef n=m gildir hins vegar að einx,eimx=1.

Ef rita má 2π-lotubundið fall f með röð á forminu

f(x)=n=cneinx

og ef víxla má á heildi og óendanlegri summu í eftirfarandi reikningum þá fæst

12πππf(x)eimxdx=f,eimx=n=cneinx,eimx=cm.

Þar með eru stuðlarnir cn ótvírætt ákvarðaðir og jafnir Fourier-stuðlum fallsins f og f er jafnt Fourier-röð sinni. Í framhaldinu munum við fjalla betur um þessa reikninga og undir hvaða skilyrðum þeir eru rættlætanlegir.

2.3.3. Regla Pýþagórasar

Ef u,vL2[π,π] eru hornrétt, þá er

u+v2=u2+v2.

Nokkuð einfalt er að sanna eftirfarandi ójöfnu.

2.3.4. Bessel-ójafnan

Ef fL2([π,π]) er 2π–lotubundið og hefur Fourier-stuðla cn=cn(f), þá er

+n=|cn|212πππ|f(x)|2dx.

Losaralegir reikningar leyfa okkur að færa rök fyrir því að sterkari niðurstaða gildir, ójafnan er í raun jafna:

Ef rita má

f(x)=n=cneinx

og að því gefnu að víxla megi á óendanlegum summum og heildum í eftirfarandi reikningum þá er

12πππ|f(x)|2dx=f,f=n=cneinx,m=cmeimx=n=m=cn¯cmeinx,eimx=n=m=cn¯cmδnm=n=|cn|2.

Táknið δnm sem kallast Kronecker-δ og uppfyllir δmn=1 ef m=n en δmn=0 annars. Það er talsvert flóknara að réttlæta þessa niðurstöðu með fullnægjandi hætti en það er hægt og við ræðum niðurstöðuna aftur þegar við fjöllum um Parseval-jöfnuna.

2.3.5. Skilgreining

Fall f á R er sagt vera samfellt deildanlegt á köflum ef skipta má R í endanlega mörg bil með skiptipunktum x1,x2,,xk þannig að fallið er samfellt diffranlegt á opnu bilunum ]xj,xj+1[ og afleiðan hefur markgildi frá hægri í vinstri endapunkti bils og markgildi frá vinstri í hægri endapunkti bils. Mengi falla sem eru samfellt deildanleg á köflum er táknað með PC1(R).

Við munum skoða föll sem eru 2π-lotubundin og tilheyra menginu PC1(R)C(R), þ.e. eru samfellt diffranleg á köflum og samfelld. Dæmi um slíkt fall er 2π-lotubundna fallið sem er skilgreint með formúlunni f(x)=x2 á [π,π].

2.3.6. Setning

Ef fPC1(R)C(R) er 2π–lotubundið, þá er cn(f)=incn(f),

+n=|cn(f)|<+,

og þar með er Fourier–röðin +cn(f)einx samleitin í jöfnum mæli á R.

Meginniðurstaða þessarar greinar er eftirfarandi setning sem sýnir undir hvaða skilyrðum og í hvaða skilningi fall er jafnt Fourier-röð sinni. Rifjum upp ritháttinn

f(x+)=limyx+f(y)ogf(x)=limyxf(y)

2.3.7. Setning - Andhverfuformúla Fouriers

Ef fPC1(R) er 2π–lotubundið fall með Fourier–stuðla cn=cn(f), Fourier-kósínus–stuðla an=an(f) og Fourier–sínus–stuðla bn=bn(f), þá gildir

12(f(x+)+f(x))=+n=cneinx=limN+Nn=Ncneinx=12a0+n=1(ancosnx+bnsinnx).

Í punktum x þar sem f er samfellt gildir f(x)=12(f(x+)+f(x)) og þar með er

f(x)=+n=cneinx=12a0+n=1(ancosnx+bnsinnx).

Ef fPC1(R)C(R), þá eru raðirnar samleitnar í jöfnum mæli á R.

Þegar 2π-lotubundið fall fL2([π,π]) er ósamfellt gildir almennt ekki að það sé jafnt Fourier-röð sinni í ósamfelldnipunktunum. Við getum samt spurt okkur hvort hægt sé að tala um að fallið sé jafnt Fourier-röð sinni í einhverjum öðrum skilningi. Eftirfarandi setning segir okkur að hlutsumman

sN=Nn=Ncn(f)einx

stefnir á fallið f í staðlinum á L2([π,π]).

2.3.8. Setning - Parseval-jafnan

Ef fL2[π,π] er 2π–lotubundið, þá gildir

fsN2=12πππ|f(x)Nn=Ncn(f)einx|2dx0,N+

og af þessu leiðir jafna Parseval

+n=|cn(f)|2=12πππ|f(x)|2dx,

Athugið

Mismunurinn fNn=Ncneinx2 nefnist ferskekkja nálgunar f með Nn=Ncneinx. Hægt er að sýna að lágmarks ferskekkja fæst með því að velja stuðlana cn=cn(f).

2.4. Úrlausn á hlutafleiðujöfnum

Í þessari grein munum við líta á dæmi þar sem hagnýta má Fourier-raðir við lausn jaðargildisverkefna. Byrjum á tveimur mikilvægum skilgreiningum.

Þegar fengist er við ákveðnar tegundir jaðargildisverkefna getur verið gagnlegt að skilgreina lotubundna framlengingu af falli á bili sem annað hvort er oddstæð eða jafnstæð. Með þeim hætti má skilgreina raðir sem uppfylla sjálfkrafa jaðarskilyrðin sem gefin eru.

2.4.1. Jafnstæð framlenging og kósínus-röð

Ef L>0 og f:[0,L]C er fall á endanlegu bili skilgreinum við jafnstæða 2L-lotubundna framlengingu á f með því að setja

fJ(x)={f(x),x[0,L],f(x),x[L,0],

og framlengja fJ í 2L-lotubundið fall.

_images/jafnstaett.png

Jafnstæð framlenging falls f:[0,L]C í 2L -lotubundið fall fJ.

Fourier-stuðlar fJ eru gefnir með

an(fJ)=1LLLfJ(x)cosnπLxdx=2LL0fJ(x)cosnπLxdx=2LL0f(x)cosnπLxdx,n=0,1,2,,bn(fJ)=0n=1,2,3,.

Stuðlarnir an nefnast Fourier–kósínus–stuðlar fallsins f og röðin

12a0+n=1ancosnπLx

kallast Fourier–kósínus–röð fallsins f.

2.4.2. Oddstæð framlenging og sínus-röð

Ef L>0 og f:[0,L]C er fall á endanlegu bili skilgreinum við oddstæða 2L-lotubundna framlengingu á f með því að setja

fO(x)={f(x),x[0,L],f(x),x[L,0],

og framlengja fO í 2L-lotubundið fall.

_images/oddstaett.png

Oddstæð framlenging falls f:[0,L]C í 2L -lotubundið fall fO.

Fourier-stuðlar fO eru gefnir með

an(fO)=0n=0,1,2,,bn(fO)=1LLLfO(x)sinnπLxdx=2LL0fO(x)sinnπLxdx=2LL0f(x)sinnπLxdx,n=1,2,.

Stuðlarnir bn nefnast Fourier–sínus–stuðlar fallsins f og röðin

n=1bnsinnπLx

kallast Fourier–sínus–röð fallsins f.

Athugið

Hægt er að yfirfæra allar reiknireglur og fræðilegar niðurstöður líkt og t.d. andhverfusetninguna beint á Fourier-kósínus og Fourier-sínus raðir. Vísað er í kennslubók fyrir frekari smáatriði.

2.4.3. Setning

Látum P vera margliðu af stigi m og lítum á jöfnuna

P(D)u=(amDm+am1Dm1++a1D+a0)u=f(x),

þar sem fPC1(R)C(R) er T–lotubundið fall og setjum ω=2π/T. Ef cn(f)=0 fyrir öll n þannig að P(inω)=0, þá fæst T–lotubundin lausn af gerðinni

u(x)=+n=P(inω)0cn(f)P(inω)einωx,xR.

Eftirfarandi dæmi má finna í kennslubók og þar eru reikningar framkvæmdir í smáatriðum.

2.4.4. Dæmi

Notum Fourier-raðir til að leysa jaðargildisverkefnið

u+ω2u=f(x),u(0)=u(1)=0.

Það hefur ótvírætt ákvarðaða lausn fyrir sérhvert f ef og aðeins ef ω er ekki heiltölumargfeldi af π. Prófum að liða u í Fourier-sínus-röð en þá eru jaðarskilyrðin uppfyllt.

Lausnin er

u(x)=n=1fnω2n2π2sin(nπx)

þar sem fn eru Fourier-sínus-stuðlar fallsins f.

2.4.5. Dæmi - Sveiflandi strengur

Lítum á einvíðan streng af lengd L sem festur er í báða enda. Táknum frávik hans frá jafnvægi í punkti x á tíma t með u(x,t). Fallið u(x,t) uppfyllir þá bylgjujöfnuna í einni rúmbreytu ásamt jaðarskilyrðunum

2ut2c22ux2=0,u(0,t)=u(L,t)=0.

Gerum einnig ráð fyrir því að upphafsstaðan og hraðinn séu þekkt

u(x,0)=φ(x),tu(x,0)=ψ(x),x]0,L[.

Þetta verkefni má leysa með því að liða u(x,t) í Fourier-sínus–röð með miðað við breytuna x. Þannig eru jaðarskilyrðin sjálfkrafa uppfyllt.

Lausnin verður

u(x,t)=n=1(φncos(nπct/L)+ψnLnπcsin(nπct/L))sin(nπx/L)

þar sem ϕn og ψn eru Fourier-sínus-stuðlar fallanna ϕ og ψ.

Lausnina má einnig rita

u(x,t)=n=1Cncos(nπct/Lαn)sin(nπx/L)

þar sem

Cn=φ2n+(ψnL/nπc)2

kallast sveifluvídd og αn kallast fasahliðrun og uppfyllir

cosαn=φn/Cn,sinαn=(ψnL)/(nπcCn).

2.4.6. Dæmi - Varmaleiðni

Reiknum út hitastig u(x,t), í punkti x á tíma t, í einvíðri stöng af lengd L, sem er einangruð í báðum endapunktunum. Jaðarskilyrðin eru þá að ekkert varmaflæði er í endapunktum stangarinnar, sem þýðir að afleiða hitastigsins er núll í jaðarpunktunum 0 og L. Fallið u uppfyllir varmaleiðnijöfnuna og við höfum því eftirfarandi jaðargildisverkefni

{utκ2ux2=f(x,t),0<x<L,t>0,xu(0,t)=xu(L,t)=0,t>0

með upphafsskilyrðinu

u(x,0)=φ(x),x]0,L[.

Föllin f og ϕ eru ótiltekin.

Fallið u er liðað í Fourier-kósínus-röð til þess að jaðarskilyrði séu uppfyllt. Þá má sýna að lausnin er

u(x,t)=n=0(φneκ(nπ/L)2t+t0eκ(nπ/L)2(tτ)fn(τ)dτ)cos(nπx/L)

þar sem ϕn og fn eru Fourier-kósínus-stuðlar fallanna ϕ og f.