给我1个有限差分法求解二维非线性Klein-Gordon方程的matlab代码

以下是一个使用有限差分法求解二维非线性Klein-Gordon方程的matlab代码：

% 定义常数和初始条件
N = 100; % 空间点数
M = 100; % 时间点数
L = 10; % 空间长度
T = 5; % 时间长度
h = L/N; % 空间步长
k = T/M; % 时间步长
a = 1; b = 1; c = 1; d = 1; % 常数
x = linspace(-L/2,L/2,N); y = linspace(-L/2,L/2,N);
[X,Y] = meshgrid(x,y);
u = exp(-X.^2-Y.^2); % 初始条件

% 定义差分算子
Dxx = toeplitz(sparse([1,1],[1,2],[-2,1]/h^2,1,N));
Dyy = toeplitz(sparse([1,1],[1,2],[-2,1]/h^2,1,N));
Dx = toeplitz(sparse([1,1],[1,2],[-1,1]/(2*h),1,N));
Dy = toeplitz(sparse([1,1],[1,2],[-1,1]/(2*h),1,N));

% 使用显式欧拉方法求解
for n=1:M
    lap_u = Dxx*u + Dyy*u;
    u_xx = Dxx*u;
    u_yy = Dyy*u;
    u_t = a*u - b*u.^3 - c*u_xx - d*u_yy;
    u = u + k*u_t + k^2/2*lap_u;
end

% 绘制图形
figure;
surf(X,Y,u);
xlabel('x'); ylabel('y'); zlabel('u');
title('Solution of Klein-Gordon Equation')

其中，使用了显式欧拉方法进行求解。需要注意的是，此代码仅供参考，可能需要针对具体问题做出一些调整。