偏微分方程数值解法题解教程.docVIP

PAGE  PAGE - 34 - 偏微分方程数值解法（带程序） 例1 求解初边值问题要求采用树脂格式 ，，完成下列计算： 取分别计算时刻的数值解。 取分别计算时刻的数值解。 取分别计算时刻的数值解。 并与解析解进行比较。 解：程序 function A=zhongxinchafen(x,y,la) U=zeros(length(x),length(y)); for i=1:size(x,2) if x(i)0x(i)=0.5 U(i,1)=2*x(i); elseif x(i)0.5x(i)1 U(i,1)=2*(1-x(i)); end end for j=1:length(y)-1 for i=1:length(x)-2 U(i+1,j+1)=U(i+1,j)+la*(U(i+2,j)-2*U(i+1,j)+U(i,j)); end end A=U(:,size(U,2)) function u=jiexijie1(x,t) for i=1:size(x,2) k=3; a1=(1/(1^2)*sin(1*pi/2)*sin(1*pi*x(i))*exp(-1^2*pi^2*t)); a2=a1+(1/(2^2)*sin(2*pi/2)*sin(2*