《抛物型方程差分格式.doc》由会员分享，可在线阅读，更多相关《抛物型方程差分格式.doc（2页珍藏版）》请在三一文库上搜索。
1、;实验：抛物型方程差分格式 (1).cleara=0;b=1;J=10;h=(b-a)/J;T=1;lamda=0.25;taol=lamda*h2;x=a:h:b;t=0:taol:T;N=length(t);u=zeros(N,J+1);u(1,:)=sin(pi*x);for n=2:N for j=2:J u(n,j)=lamda*u(n-1,j+1)+(1-2*lamda)*u(n-1,j)+lamda*u(n-1,j-1); endendxx,tt=meshgrid(x,t);zz=exp(-pi2*tt).*sin(pi*xx);subplot(1,2,1)mesh(xx,tt,
2、zz)subplot(1,2,2)mesh(xx,tt,u)(2).cleara=0;b=1;J=10;h=(b-a)/J;T=1;lamda=8;taol=lamda*h2;x=a:h:b;t=0:taol:T;N=length(t);u=zeros(N,J+1);u(1,:)=sin(pi*x);S=zeros(J-1,J-1);for i=1:J-1 for j=1:J-1 if abs(i-j)=1 S(i,j)=1; end endendA=(1+2*lamda)*eye(J-1)-lamda*S;B=inv(A);for n=2:N u(n,2:J)=(B*(u(n-1,2:J)&
3、#39;)'endxx,tt=meshgrid(x,t);zz=exp(-pi2*tt).*sin(pi*xx);subplot(1,2,1)mesh(xx,tt,zz)subplot(1,2,2)mesh(xx,tt,u)(3).cleara=0;b=1;J=10;h=(b-a)/J;T=1;lamda=0.5;taol=lamda*h2;x=a:h:b;t=0:taol:T;N=length(t);u=zeros(N,J+1);u(1,:)=sin(pi*x);S=zeros(J-1,J-1);for i=1:J-1 for j=1:J-1 if abs(i-j)=1 S(i,j)=1; end endendA=(1-lamda)*eye(J-1)+0.5*lamda*S;B=(1+lamda)*eye(J-1)-0.5*lamda*S;for n=2:N u(n,2:J)=(inv(B)*A*(u(n-1,2:J)')'endxx,tt=meshgrid(x,t);zz=exp(-pi2*tt).*sin(pi*xx);subplot(1,2,1)mesh(xx,tt,zz)subplot(1,2,2)mesh(xx,tt,u).'