0 svar
59 visningar
katal 77
Postad: 28 jan 17:28 Redigerad: 28 jan 17:32

Differentialekvationer: Fourierserielösningar

Hej!
Har fått problem med Neumannproblem. Har försökt lösa frågan på ett sätt som jag tror (på ett väldigt Dirichlet-sätt?) den borde lösas på, men det blir uppenbarligen inte rätt. Tror jag kan förvirrat mig själv i processen, men här är iallafall mitt försök.


Lös värmeledningsproblemet:

ut-2ux2=0,0<x<1,t>0,u(x,0)=x,ux(0,t)=ux(1,t)=0\frac{\partial u}{\partial t}-\frac{\partial ^2 u}{\partial x^2}=0 , 0<x<1, t>0, u(x,0)=x , u_x(0,t)=u_x(1,t)=0

Ansätt u(x,t)=X(x)T(t)u(x,t)=X(x)T(t)

Använd separation av variabler (...) tills steget där:

X(x)=Asin(λx)+Bcos(λx)X(x)=Asin(\sqrt {\lambda}x )+Bcos(\sqrt {\lambda}x )

T(t)=C1e-n2π2tT(t)=C_1 e^{-n^2 \pi ^2 t}

 

Randvillkor ger:

X'(x)=Aλcos(λx)-Bλsin(λx)X'(x)=A\sqrt {\lambda}cos(\sqrt {\lambda}x )-B\sqrt {\lambda}sin(\sqrt {\lambda}x )

X'(0)=0=Aλcos(0)=Aλ=>A=0X'(0)=0=A\sqrt {\lambda}cos(0)=A\sqrt {\lambda}=> A=0

X'(1)=0=-Bλsin(λ)=>λ=n2π2,B=1=>X(x)=cos(nπx)X'(1)=0=-B\sqrt {\lambda}sin(\sqrt {\lambda}) => \lambda=n^2\pi ^2 , B=1=> X(x)=cos(n \pi x)

Detta insätts i en fourierserie:

u(x,t)=a0+X(x)T(t)=a0+cos(nπx)Cne-n2π2tu(x,t)=a_0+\sum {X(x)T(t)}=a_0+\sum cos(n \pi x)C_n e^{-n^2 \pi ^2 t}

För C_n tar vi Begynnelsevillkoret och insätter:

u(x,0)=x=a0+cos(nπx)Cnu(x,0)=x=a_0+\sum cos(n \pi x)C_n

Vi får att a_0 =1/2 , och C_n :

xcos(kπx)=cos(kπx)2+(cos(nπx)Cn)cos(kπx)xcos(k \pi x)=\frac{cos(k \pi x)}{2} +( \sum cos(n \pi x)C_n)cos(k \pi x)

Integrera, ta ut C_n:

an=(xcos(kπx)-cos(kπx)2)dx(cos(nπx)cos(kπx))dxa_n=\frac{\int (xcos(k \pi x)-\frac{cos(k \pi x)}{2})dx}{\int (cos(n \pi x)cos(k \pi x))dx}

 

Integralen löses då k=n:(Gränserna är 1<x<0)

an=2(xcos(kπx)-cos(kπx)2)dx=[nπ(2x-1)sin(nπx)+2cos(nπx)π2n2]ana_n=2\int (xcos(k \pi x)-\frac{cos(k \pi x)}{2})dx=[\frac{n \pi (2x-1)sin(n \pi x)+2cos(n \pi x)}{ \pi ^2 n^2}]a_n

=(2cos(nπ)π2n2)-(2π2n2)=(\frac{2cos(n \pi)}{\pi ^2n^2})-(\frac{2}{\pi ^2 n^2})

=2n2π2((-1)n-1)=2n2π2(-1)n+1=\frac{2}{n ^2 \pi ^2}((-1)^n-1)=\frac{2}{n^2 \pi ^2}(-1)^{n+1}

Min lösning:

u(x,t)=12+2n2π2(-1)n+1cos(nπx)e-n2π2tu(x,t)=\frac{1}{2}+ \sum \frac{2}{n^2 \pi ^2}(-1)^{n+1}cos(n \pi x) e^{-n^2 \pi ^2 t}

 

Vilket är helt fel, för att svaret ska vara:

u(x,t)=12-4(2n-1)2π2e-(2n+1)2π2tcos((2n+1)πx)u(x,t)=\frac{1}{2}-\sum \frac{4}{(2n-1)^2 \pi ^2}e^{-(2n+1)^2 \pi ^2 t}cos((2n+1)\pi x)

 

Någon som kan leda mig i rätt riktning eller förklara vart jag tänker fel?

Tack!

Svara
Close