0 svar
25 visningar
naytte Online 7909 – Moderator
Postad: Igår 21:06 Redigerad: Igår 21:18

Hur ska man bestämma konstanten C_n, PDE?

Hej!

Jag sitter med uppgiften nedan:

Jag började med att konstatera att det finns radiell symmetri, så vi söker en lösning u(r,t)u(r,t). Jag ansatte

u(r,t)=R(r)T(t)\displaystyle u(r,t)=R(r)T(t) och fick då

R''(r)+r-1R(r)R(r)=T'(t)0.05T(t)=:λ\displaystyle \frac{R''(r)+r^{-1}R(r)}{R(r)}=\frac{T'(t)}{0.05T(t)}=:\lambda

Vi börjar med RR:

r2R''(r)+rR'(r)-r2λR(r)=0\displaystyle r^2R''(r)+rR'(r)-r^2\lambda R(r)=0

Här kommer min första fråga. Jag ansatte här att λ<0\lambda < 0 eftersom jag såg min föreläsare göra detta i en uppgift, men jag förstår inte hur man kommer fram till att man ska göra detta.

Hur som helst ansatte jag x:=r|λ|x:=r\sqrt{|\lambda |} och definierade en funktion FF enligt R(r)=:F(x(r))R(r)=: F(x(r)). Då transformerades ekvationen till

x2F''(x)+xF'(x)+x2F(x)\displaystyle x^2F''(x)+xF'(x)+x^2F(x)

Eftersom Y0Y_0 beter sig dåligt då r0r \to 0 är måste lösningen bestå av J0J_0. Vi har alltså, upp till en skalär,

F(x)=J0(x)R(r)=J0(r-λ)F(x)=J_0(x) \implies R(r) = J_0(r\sqrt{-\lambda})

Randvillkor ger

Rn(r)=J0(r-λn)R_n(r)=J_0(r\sqrt{-\lambda_n}), λn=-(πn5)2\lambda_n = -(\frac{\pi_n}{5})^2

Vi har då

Tn(t)=Cne-0.05(πn5)2t\displaystyle T_n(t)=C_ne^{-0.05(\frac{\pi_n}{5})^2t}

Med superposition:

u(r,t)=n1Rn(r)Tn(t)=n1J0(rπn5)Cne-0.05(πn5)2t\displaystyle u(r,t)=\sum_{n\ge 1}R_n(r)T_n(t)=\sum_{n\ge 1}J_0(\frac{r\pi_n}{5})C_ne^{-0.05(\frac{\pi_n}{5})^2t}

Nu vet vi ju också att vid t=0t=0 har vi

n1J0(rπn5)Cn=30\displaystyle \sum_{n\ge 1}J_0(\frac{r\pi_n}{5})C_n=30

Jag tänkte att man skulle utveckla 3030 i basen J0(rπn5)J_0(\frac{r\pi_n}{5}) för att bestämma CnC_n genom att jämföra koefficienter, men hur gör man det? Jag vet att bara J0(πnr)J_0(\pi_n r) är en bas för r  2(0,1)\mathcal{L}^{\;\; 2}_r(0,1), men vi har ju inte dessa funktioner i summan... Blir våra funktioner en bas för L2(0,5) eftersom vi delar med 5 inuti?

Svara
Close