写一个3D拉普拉斯算子的MATLAB代码

好的，以下是3D拉普拉斯算子的MATLAB代码：

function [L] = laplacian_3D(x,y,z)
% x, y, z分别为三个方向上的坐标数组

hx = x(2) - x(1); % 步长
hy = y(2) - y(1);
hz = z(2) - z(1);

Dx = spdiags([ones(size(x)),-2*ones(size(x)),ones(size(x))],[-1,0,1],length(x),length(x))/hx^2; % 生成x方向上的二阶差分矩阵
Dy = spdiags([ones(size(y)),-2*ones(size(y)),ones(size(y))],[-1,0,1],length(y),length(y))/hy^2; % 生成y方向上的二阶差分矩阵
Dz = spdiags([ones(size(z)),-2*ones(size(z)),ones(size(z))],[-1,0,1],length(z),length(z))/hz^2; % 生成z方向上的二阶差分矩阵

L = kron(kron(Dz,speye(length(y))),speye(length(x))) + kron(kron(speye(length(z)),Dy),speye(length(x))) + kron(kron(speye(length(z)),speye(length(y))),Dx); % 生成拉普拉斯算子L

end

笑话时间：为什么程序员总是用黑色的键盘？因为黑色的键盘才有 CTRL。（CTRL是电脑键盘上重要的一组功能键）