数值求解二维扩散方程的初边值问题.docVIP

数值求解二维扩散方程的初边值问题古典显式格式:将原格式化为:附源程序:%-------------------------------------------运用古典显式差分格式求解二维扩散方程的初边值问题;function gdxs(ti,h,t)%-------------------------------------------ti:时间步长;%-------------------------------------------h:空间步长;k=t/ti;m=1/h+1;r=ti/h^2; %------------------------------ r为网格比;w=ones(m,m);u=ones(m,m);for i=2:m-1for j=2:m-1u(i,j)=sin(pi*(i-1)*h)*sin(2*pi*h*(j-1));endendticfor l=1:kfor i=2:m-1for j=2:m-1w(i,j)=r*u(i-1,j)+r*u(i,j-1)+r*u(i+1,j)+r*u(i,j+1)+(1-4*r)*u(i,j);endendu=w;endtoct=tocumesh(u)交替方向隐式格式(P-R格式):将原差分格式化为:代入边界条件,转化为三对角矩阵附追赶法源程序:%--