用matlab3维线性谐振子的能级数值仿真的代码，用散点图体现能级值与x，y，z的关系

以下是一个简单的Matlab代码来计算和可视化3D线性谐振子的能级：

% 定义线性谐振子参数
m = 1;
omega = 1;
h_bar = 1;

% 定义范围和步长
x = linspace(-2,2,100);
y = linspace(-2,2,100);
z = linspace(-2,2,100);
dx = x(2) - x(1);
dy = y(2) - y(1);
dz = z(2) - z(1);

% 计算能级
n_max = 10;
E = zeros(n_max,length(x),length(y),length(z));
for n = 1:n_max
    for i = 1:length(x)
        for j = 1:length(y)
            for k = 1:length(z)
                E(n,i,j,k) = h_bar*omega*(n + 1/2);
                + 1/2*m*omega^2*(x(i)^2 + y(j)^2 + z(k)^2);
            end
        end
    end
end

% 绘制能级与位置的关系
figure;
for n = 1:n_max
    scatter3(x,y,z,E(n,:,:,:),'.');
    hold on;
end
xlabel('x');
ylabel('y');
zlabel('z');
title('3D线性谐振子的能级');

这个代码定义了线性谐振子的参数，然后在给定的范围内计算了能级，并使用scatter3函数将能级与位置的关系可视化为散点图。