偏微分方程-数值解).docVIP

学 生 实 验 报 告 实验课程名称 偏微分方程数值解 开课实验室 数统104机房 学 院 数统学院 年级 xx专业班 xx班 学 生 姓 名 xxxxxx 学 号 xxxxx 开 课 时 间 2013 至 2014 学年第 2 学期 总 成 绩 教师签名 数学与统计学院制 开课学院、实验室： 实验时间 ： 2014 年5 月 7 日 实验项目 名 称 扩散问题 实验项目类型 验证 演示 综合 设计 其他 指导教师 成 绩 一．实验目的 用显示和隐示方法求解热传导方程，比较其不同 二．实验内容 分别用显示和隐示方法求解，并画出t=0.02,0.04,0.1的图像。 三．实验原理、方法（算法）、步骤 程序1: function U=explicit(M) N=20;U=zeros(6,22); g=sin(2*pi*x); h=1/(N+1);k=0.1/M;X=0:h:1;w=k/h^2; A=zeros(N,N); u=zeros(M+1,20); for i=1:20 x=h*i; u(1,i)=eval(g); end for k=2:N-1 A(k,k-1)=w;A(k,k)=1-2*w;A(k,k+1)=w; end for j=2:M+1 u(j,:)=A*u(j-1,:); end for j=1:M+1 for i=2:21 U(j,i)=u(j,i-1); end end end function Time2 X=0:1/21:1;U1=explicit(5);U2=explicit(10);U3=explicit(20); fun=exp(-4*(pi^2)*0.02)*sin(2*pi*x); plot(X,U1(2,:),-o,X,U2(3,:),-s,X,U3(5,:),--); axis([0,1,-0.5,0.5]); hold on plot(fun,X); legend(M=5,M=10,M=20,Exact,3); function Time4 X=0:1/21:1; U1=explicit(5); U2=explicit(10); U3=explicit(20); fun=exp(-4*(pi^2)*0.04)*sin(2*pi*x); plot(X,U1(3,:),-o,X,U2(5,:),-s,X,U3(9,:),--); axis([0,1,-0.5,0.5]); hold on plot(X,fun); legend(M=5,M=10,M=20,Exact,3); function Time1 X=0:1/21:1; U1=explicit(5); U2=explicit(10); U3=explicit(20); fun=exp(-4*(pi^2)*1)*sin(2*pi*x); plot(X,U1(6,:),-o,X,U2(11,:),-s,X,U3(21,:),--); axis([0,1,-60,60]); hold on plot(X,fun); legend(M=5,M=10,M=20,Exact,3); 分别运行的结果为: 方法2: function U=implicit(M) N=20;U=zeros(M+1,22); g=sin(2*pi*x); h=1/(N+1);k=0.1/M;X=0:h:1; w=k/h^2; B=zeros(N,N); B(1,1)=1+2*w; B(1,2)=-w; B(N,N-1)=-w; B(N,N)=1+2*w; u=zeros(M+1,20); for i=1:20 x=h*i; u(1,i)=eval(g); end for k=2:N-1 B(k,k-1)=-w;B(k,k)=1+2*w; B(k,k+1)=-w; end for j=2:M+1 u(j,:)=inv(B)*u(j-1,:); end for j=1:M+1 for i=2:21 U(j,i)=u(j,i-1); end end end function Time22 X=0:1/21:1; U1=implicit(5); U2=impli