matlab画径向波函数,关于使用MATLAB画氢原子波函数的问题

参考下面的程序命令，修改一下，希望可以自己定义一些命令，参数，最后得到图片的效果，需要相应的函数程序，可以发邮件联系yxzsyxzy

@163.com

%氢原子电子云图

clc;

close all;

figure('numbertitle','off','name','波函数图');

subplot(2,2,1);

r=0:80;

plot(r,R53(r));

xlabel('r');

ylabel('R(r)');

title('5f');

grid on; % 径向波函数R(r)

subplot(2,2,2);

plot(r,R53(r).^2);

xlabel('r');

ylabel('R(r)^2');

title('5f');

grid on; %径向波函数密度R(r)*R(r)

subplot(2,2,3);

D=(r.^2).*(R53(r).^2);

plot(r,D);

xlabel('r');

ylabel('D(r)');

title('5f');

grid on;

subplot(2,2,4);

D=(r.^2).*(R63(r).^2);

plot(r,D);

xlabel('r');

ylabel('D(r)');

title('6f');

grid on;%径向分布函数D(r)

%%

figure('numbertitle','off','name','原子轨道图');

subplot(2,2,1);

p=-pi/2:pi/30:pi/2;

t=0:2*pi/50:2*pi;

[P,T]=meshgrid(p,t);

theta=pi/2-P;

phi=T;

R=Ypx(theta,phi);

[X,Y,Z]=sph2cart(T,P,R);%球坐标转化为笛卡尔坐标的XYZ

mesh(X,Y,Z);

axis equal;

xlabel('x'),ylabel('y'),zlabel('z');

title('Ypx');%角度波函数Y

subplot(2,2,2);

R=Ypx(theta,phi);%调用函数

R=R.*R;

[X,Y,Z]=sph2cart(T,P,R);

mesh(X,Y,Z);

axis equal;

xlabel('x'),ylabel('y'),zlabel('z');

title('1Ypx1^2');

subplot(2,2,3);

R=Ydxy(theta,phi);

R=R.*R;

[X,Y,Z]=sph2cart(T,P,R);

mesh(X,Y,Z);

axis equal;

xlabel('x'),ylabel('y'),zlabel('z');

title('1Ydxy1^2');

subplot(2,2,4);

R=Yfzxy(theta,phi);

R=R.*R;

[X,Y,Z]=sph2cart(T,P,R);

mesh(X,Y,Z);

axis ('equal');

xlabel('x'),ylabel('y'),zlabel('z');

title('1Yfzxy1^2');

figure('numbertitle','off','name','电子云图');

x=-15:15;

y=x;

[X,Y]=meshgrid(x,y);

Z=zeros(size(X));

psi=psi_3dxy(X,Y,Z);

psi2=psi.*psi;

subplot(2,2,1);

mesh(X,Y,psi);

axis tight;

xlabel('x'),ylabel('y'),zlabel('3dxy');

subplot(2,2,2);

mesh(X,Y,psi2);

axis tight;

xlabel('x'),ylabel('y'),zlabel('13dxy1^2');

subplot(2,2,3);

contour(psi2,10);

axis tight;

xlabel('x'),ylabel('y'),title('13dxy1^2');

z=x;

[X,Y,Z]=meshgrid(x,y,z);

psi=psi_3dxy(X,Y,Z);

psi2=psi.*psi;

subplot(2,2,4);

slice(X,Y,Z,psi2,[],[],[0]);

axis tight;

xlabel('x'),ylabel('y'),title('13dxy1^2');

view(0,90);

colorbar;

2016-5-7 17:45 上传

效果图