7 svar
173 visningar
naytte behöver inte mer hjälp
naytte 7934 – Moderator
Postad: 12 feb 22:27 Redigerad: 12 feb 22:51

Vad är grejen med egenfunktioner och egenvärden då vi löser differentialekvationer+relaterad fråga?

Hej!

Jag ska tillstå att jag inte ännu har befattat mig mycket med teorin, men jag skulle ändå vilja ställa denna fråga för att få ett hum om varför man är så intresserad av detta. Jag går en kurs i Fourieranalys nu och i många PDEs och ODEs är vi väldigt intresserade av "egenvärden" och "egenfunktioner".

Är anledningen till att man är så intresserad av dem att egenfunktionerna (som ofta har egenvärdena i sina argument) spänner upp lösningsrummet till vår differentialekvationer? Gäller detta endast vissa typer av differentialekvationer?


Jag hade också en relaterad fråga gällande lösnigsförfarandet av vissa (P)DEs. Vi kan ta differentialekvationen nedan som exempel:

f''+λf=0f(0)=0f(2)=f'(2)=0\displaystyle \left\{\begin{array}{lll}f''+\lambda f=0 \f(0)=0\f(2)=f^\prime(2)=0\end{array}\right.

I fallet då λ=0\lambda =0 brydde vi oss inte om lösningarna eftersom de inte utgjorde "egenfunktioner" (alltså inte spelar någon roll för basen till lösningsrummet då antar jag?). Om λ0\lambda \ne 0 får vi snabbt en allmän lösning:

fx=Acoshx-λ+Bsinhx-λ\displaystyle f\left(x\right)=A\cosh\left(x\sqrt{-\lambda}\right)+B\sinh\left(x\sqrt{-\lambda}\right)

Initialvillkor ger A=0A=0. Nu är det så att vår övningsledare påstod att vi kan bestämma att B=1B=1 eftersom alla skalärmultiplar av ff är linjärt beroende "i vårt  2\mathcal{L}^{\; 2}-rum". Vad menar han med det? Vad har  2\mathcal{L}^{\;2} med detta att göra?

Trinity2 Online 4291
Postad: 12 feb 22:52

Det var länge sedan jag läste ODE (för Nisse Nilsson, en mästare på just detta), men jag tror man vill diagonalisera dem och på så sätt få enkla lösningar. LuMa vet säkert allt om detta.

naytte 7934 – Moderator
Postad: 13 feb 20:56

bump

D4NIEL 3518
Postad: 14 feb 16:46 Redigerad: 14 feb 17:37

En väldigt bred klass av differentialekvationer kallas Sturm-Liouvilleproblem om man kan skriva dem på formen

L(f)+λwf=0,  B1(f)=B2(f)=0L(f)+\lambda w f=0,\quad B_1(f)=B_2(f)=0

där λ\lambda är en godtycklig konstant. För de flesta problem är f(x)0f(x)\equiv 0 den enda möjliga lösningen, men för icke-triviala lösningar bildar λ\lambda ett egenvärde med korresponderande egenvärdeslösningar. Notera att vi har

L(f)=-λwf-w-1L(f)=λfL(f)=-\lambda w f \implies -w^{-1}L(f)=\lambda f

vilket betyder att egenvärdena gäller operatorn -w-1L-w^{-1}L. Operatorn L(f)L(f) definieras som (rf')'+pf(rf^\prime)^\prime+pf och är självadjungerad i w2\mathcal{L}^{\,\,\,2}_w . För ett (icke-singlulärt) Sturm-Liouvilleproblem gäller att:

  • Alla egenvärden är reella.
  • Egenfunktioner hörande till olika egenvärden är ortogonala (med avseende på ww).
  • Egenrummet för ett egenvärde λ\lambda har som mest två dimensioner. Om randvillkoren är separerade är varje egenvärdesrum endimensionellt.

Ett egenvärdesrum består av 0\mathbf{0}-vektorn och korresponderande egenvektorer. Det är ett linjärt rum i vanlig mening och varje linjärkombination av funktioner ff och gg, till exempel c1f+c2gc_1f+c_2g, är därmed också en lösning, (superpositionsprincipen).

Man kan visa att på varje intervall [a,b][a,b] finns en ortonormal bas {ϕn}1\{\phi_n\}_1^\infty för w2(a,b)\mathcal{L}^{\,\,\,2}_w(a,b) som spänns av egenfunktioner. Givet randvillkoren B1(f)=B2(f)=0B_1(f)=B_2(f)=0 samt fC2[a,b]f\in C^2[a,b] konvergerar serien <f,ϕn>ϕnf\sum <f,\phi_n>\phi_n \to f likformigt.

Det exempel du tar upp är överbestämt och endast lösbart för f0f\equiv 0, men jag misstänker att du skrivit fel och egentligen avser problemet

f''+λf=0,f^{\prime \prime}+\lambda f=0,

f(0)=0,  f(2)-f'(2)=0f(0)=0,\quad f(2)-f^\prime(2)=0

Det är ett reguljärt Sturm-Liouvilleproblem med r(x)=1r(x)=1, q(x)=0q(x)=0 och viktfunktionen w=1w=1. Fallet λ=0\lambda=0 ger ekvationen f''=0f^{\prime\prime}=0 vilket med insatta randvillkor ger den triviala lösningen f(x)0f(x)\equiv 0.

För λ<0\lambda <0 får vi egenfunktioner på formen fn(x)=sinh(μnc)f_n(x)=\sinh(\mu_n c) (λn=-μn2\lambda_n = -\mu^2_n). Ibland måste man normalisera lösningarna, dvs fw2\| f\|^2_w  med avseende på randvillkoren. I ditt exempelfall spelar normeringen ingen roll. Du får sätta BB till vad du vill, det är ju bara en normering av basvektorn, du måste inte göra den ortonormal.

Edit: Slutligen, för λ>0\lambda >0 (med λn=μn2\lambda_n=\mu^2_n) får man egenfunktioner på formen fn(x)=sin(μnx)f_n(x)=\sin(\mu_n x)

naytte 7934 – Moderator
Postad: 16 feb 16:47

Tack för det utförliga svaret!

Jag är fortfarande inte riktigt med på varför vi är så intresserade av L2 när vi löser ekvationen. Kan det inte finnas lösningar som t.ex. ligger utanför L2 på det intervall vi studerar? Här är ett utdrag ur vår kurslitteratur där funktionen beskriver en vibrerande sträng med lämpliga initialvillkor:

Så som jag förstår det nu är vi särskilt intresserade av just L2 eftersom vi kan bestämma koefficienterna ana_n och bnb_n genom att jämföra vår funktion uu med "rätt" Fourierserie i L2 givet initialvillkor u(x,0)=f(x)u(x,0)=f(x). Är det rätt uppfattat?

D4NIEL 3518
Postad: 18 feb 18:21 Redigerad: 18 feb 18:41

Vi vill väldigt gärna ha ett vektorrum som vi kan låta våra funktioner och fourierserier verka i. Bland annat vill vi kunna definiera likhet, konvergens i norm och så vidare. Koncept som varit så framgångsrika i andra sammanhang (till exempel linjär algebra och vektoranalys!).

Frågan är vilket vektorrum vi ska välja? En första tänkbar kandidat kan vara PC(a,b)\mathrm{PC}(a,b), dvs rummet av alla styckvist kontinuerliga funktioner på intervallet [a,b][a,b]. Det visar sig att vi bland annat kan använda det rummet för att definiera konvergens i norm. Men det visar sig också att det valda rummet har en fruktansvärd nackdel, det är inte komplett! Med komplett menar vi att varje Cauchysekvens av funktioner i rummet ska har ett gränsvärde.  

Men till vår förfäran visar det sig alltså  att man kan konstruera Cauchysekvenser av styckvist kontinuerliga funktioner som konvergerar mot funktioner som inte längre är styckvist kontinuerliga. Och det är ju lite av en bummer.

Som tur är, och som vi förmodligen redan känner till, kan vi använda Lebesgueintegralen för att definera L2(a,b)L^2(a,b) med skalärprodukten <f,g>=abf(x)g(x)¯dx<f,g>=\int_a^b\, f(x)\bar{g(x)}\,\mathrm{d}x och har då ett nytt, bättre funktionsrum, som är Cauchykomplett. 

Så jag skulle säga att främsta anledningen till att vi gillar L2(a,b)L^2(a,b) är att det är ett komplett rum.

Skalärprodukten <f,g>=abf(x)g(x)¯dx< f,g>=\int_a^b f(x)\overline{g(x)}\,dx gör det möjligt att bestämma Fourierkoefficienterna som projektioner, cn=f,φnc_n=\langle f,\varphi_n\rangle. Fullständigheten garanterar att projektionssumman konvergerar till en funktion i samma rum. Utan fullständighet skulle Fourierutvecklingen som sagt kunna konvergera mot något som inte längre tillhör funktionsrummet.

Men det här rummet med tillhörande standarddefinitioner för också med sig visa underligheter, t.ex. att om fL2(a,b)f\in L^2(a,b) och dess norm är noll, så innebär det att f(x)=0f(x)=0 "nästan överallt, men bara nästan". Eller min favorit "for almost every x[a,b]x\in[a,b]".  

På samma sätt blir det roligt att definiera vad man menar med att två funktioner är lika.

naytte 7934 – Moderator
Postad: 18 feb 19:09

Men hur vet vi att alla våra lösningar kommer ligga i detta rum?

D4NIEL 3518
Postad: 18 feb 20:50 Redigerad: 18 feb 21:04

Jag börjar misstänka att jag inte förstår vad du egentligen vill veta, jag trodde du vill veta varför just L2L^2-rummet är så poppis inom Fourieranalys. Men diffekvationen i tråden är en reguljär Sturm-Liouville uppställning med kontinuerliga koefficienter på det begränsade intervallet [0,2][0,2].

Det innebär att eventuella lösningar ff är kontinuerliga, faktiskt gäller fC2([0,2])f\in C^2([0,2]). Och det gäller ju såklart allmänt att C([a,b])PC(a,b)L2(a,b)C([a,b])\subset PC(a,b)\subset L^2(a,b).

Undrar du varför man vill använda L2L^2 för obegränsade områden eller när man gör Fouriertransformer?

Edit: Jag kanske ska svara så här istället:

Så länge du har kontinuerliga koefficienter och ett begränsat intervall är det ingen risk att lösningar skulle ligga utanför L2L^2. För singulära Sturm–Liouville-problem (oändliga intervall eller singulära ändpunkter) gäller detta inte i allmänhet, tvärtom kan det i sådana fall finnas lösningar som inte är L2L^2.

Svara
Close