精选优质文档-倾情为你奉上有限差分法的Matlab程序（椭圆型方程）function FD_PDE(fun,gun,a,b,c,d) %用有限差分法求解矩形域上的Poisson方程 tol=10(-6); % 误差界 N=1000; % 最大迭代次数 n=20; % x轴方向的网格数 m=20; % y轴方向的网格数 h=(b-a)/n; %x轴方向的步长 l=(d-c)/m; %y轴方向的步长 fori=1:n-1 x(i)=a+i*h; function FD_PDE(fun,gun,a,b,c,d)% 用有限差分法求解矩形域上的Poisson方程tol=10(-6); % 误差界N=1000; % 最大迭代次数n=20; % x轴方向的网格数m=20; % y轴方向的网格数h=(b-a)/n; % x轴方向的步长l=(d-c)/m; % y轴方向的步长for i=1:n-1x(i)=a+i*h;end % 定义网格点坐标for j=1:m-1y(j)=c+j*l;end % 定义网格点坐标u=zeros(n-1,m-1);
