你可以给我提供GMRES求解二维泊松偏微分方程边值问题的matlab代码吗

当然可以，以下是一个使用MATLAB实现GMRES求解二维泊松偏微分方程边值问题的示例代码：

% 设置偏微分方程的参数
N = 100; % 网格数
h = 1/N; % 步长
x = linspace(0,1,N+1); % 离散化的网格点
[X,Y] = meshgrid(x,x); % 二维网格点
f = ones(N-1,N-1); % 右端项
A = gallery('poisson',N-1); % 离散化后的系数矩阵

% 将系数矩阵转化为稀疏矩阵格式
A = sparse(A);

% 将右端项展开为向量形式
f = reshape(f,[],1);

% 定义GMRES算法的参数
maxit = 100; % 最大迭代次数
tol = 1e-6; % 收敛精度
restart = 10; % 重启次数

% 使用GMRES求解线性方程组
x0 = zeros(size(f)); % 初始解
[x,flag,relres,iter,resvec] = gmres(A,f,restart,tol,maxit);

% 将解展开为二维矩阵形式
u = reshape(x,[N-1,N-1]);

% 将边界值加上
u = [zeros(N-1,1),u,zeros(N-1,1)];
u = [zeros(1,N+1);u;zeros(1,N+1)];

% 绘制解
surf(X,Y,u)

该代码使用GMRES求解了一个简单的二维泊松方程，其中$f=1$，边界条件为$u(x,y)=0$。在实际使用中，需要根据具体问题进行适当的修改。