4

我有以下代码可用于绘制氢 2pz 轨道的 xy 切片:

%probability densities
pd_psi_210 = @(r,theta,phi) exp(-r).*(r.^2).*(cos(theta).^2)/(32*pi);

%configuring the range
[x y z] = meshgrid(-10:.1:10,-10:.1:10,-2:.1:2);
[THETA,PHI,R] = cart2sph(x,y,z);

%create array of probability density magnitudes
psi_210_vals = pd_psi_210(R,THETA,PHI);

%plotting
imagesc(psi_210_vals(:,:,1)); %x-y plane

我想绘制轨道的 3d 等高线图。我试过这个(它似乎没有得到我想要的东西):

isosurface(psi_210_vals(:,:,:)); %3D contour

我怎样才能让它工作?

4

1 回答 1

6

您只需指定底层网格和所需的级别。例如:

>> isosurface(-10:.1:10, -10:.1:10, -2:.1:2, psi_210_vals, 0.001);
>> axis equal

在此处输入图像描述

于 2011-11-07T18:43:36.720 回答