豆粑粑 matlab 画电子云
别人的程序，我就是转一下。虽然没完全看明白。
%%
%氢原子电子云图
%by PengjuZhang
%2012.03.02
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/60:pi/2;
t=0:2*pi/100:2*pi;
=meshgrid(p,t);
theta=pi/2-P;
phi=T;
R=Ypx(theta,phi);
=sph2cart(T,P,R);%球坐标转化为笛卡尔坐标的XYZ
mesh(X,Y,Z);
axis equal;
hidden off
xlabel('x'),ylabel('y'),zlabel('z');
title('Ypx');%角度波函数Y
subplot(2,2,2);
R=Ypx(theta,phi);%调用函数
R=R.*R;
=sph2cart(T,P,R);
mesh(X,Y,Z);
axis equal;
hidden off
xlabel('x'),ylabel('y'),zlabel('z');
title('1Ypx1^2');
subplot(2,2,3);
R=Ydxy(theta,phi);
R=R.*R;
=sph2cart(T,P,R);
mesh(X,Y,Z);
axis equal;
hidden off
xlabel('x'),ylabel('y'),zlabel('z');
title('1Ydxy1^2');
subplot(2,2,4);
R=Yfzxy(theta,phi);
R=R.*R;
=sph2cart(T,P,R);
mesh(X,Y,Z);
axis ('equal');
hidden off
xlabel('x'),ylabel('y'),zlabel('z');
title('1Yfzxy1^2');
h=gcf;
fig_na='./fig_ele_cloud';
fun_work_li_035_myfig_out(h,fig_na,3)
%%
% figure('numbertitle','off','name','电子云图');
% x=-15:15;
% y=x;
% =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;
% =meshgrid(x,y,z);
% psi=psi_3dxy(X,Y,Z);
% psi2=psi.*psi;
% subplot(2,2,4);
% slice(X,Y,Z,psi2,[],[],);
% axis tight;
% xlabel('x'),ylabel('y'),title('13dxy1^2');
% view(0,90);
% colorbar;
%
%
页:
[1]