你好给我一个用五点差分法求解泊松方程的MATLAB代码

以下是用五点差分法求解泊松方程的MATLAB代码：

%% 五点差分法求解泊松方程 clear all; clc;

%% 定义网格 Lx = 1; Ly = 1; % 区域大小 Nx = 101; Ny = 101; % 网格数 dx = Lx/(Nx-1); dy = Ly/(Ny-1); % 网格步长 x = linspace(0, Lx, Nx); y = linspace(0, Ly, Ny); [X, Y] = meshgrid(x, y);

%% 定义边界条件 u0 = sin(pi*X).sin(piY); % 边界值 u = u0; % 初始值

%% 迭代求解 for k = 1:10000 % 迭代次数 u_old = u; % 上一次迭代的值 for i = 2:Nx-1 % 内部网格点 for j = 2:Ny-1 % 五点差分法 u(i,j) = (u_old(i+1,j)+u_old(i-1,j)+u_old(i,j+1)+u_old(i,j-1))/4 - (dx^2+dy^2)/4*f(X(i,j),Y(i,j)); end end % 边界条件 u(1,:) = u0(1,:); u(Nx,:) = u0(Nx,:); u(:,1) = u0(:,1); u(:,Ny) = u0(:,Ny); end

%% 绘图 figure; subplot(1,2,1); surf(X,Y,u0); title('边界条件'); subplot(1,2,2); surf(X,Y,u); title('求解结果');

%% 定义泊松方程右侧函数 function z = f(x,y) z = 2pi^2sin(pi*x)sin(piy); en