使用MATLAB绘制原子轨道和电子云图形.doc

下载文档 降价啦

234
0
约2.39千字
约 5页
2018-12-10 发布于江西
举报
版权申诉
保障服务

使用MATLAB绘制原子轨道和电子云图形.doc

1、本文档共5页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

使用MATLAB绘制原子轨道和电子云图形

使用MATLAB绘制原子轨道和电子云图形截图6.jpg (34.91 KB) 2009-2-1 20:13 截图7.jpg (80.81 KB) 2009-2-1 20:13 复制内容到剪贴板代码: %%%氢原子电子云图%by dynamic%see also %2009.1.20function matlabskyclc;close all;figure(numbertitle,off,name,波函数图——by Matlabsky)r=0:80;subplot(221)plot(r,R53(r));xlabel(r);ylabel(R(r));title(5f);grid onsubplot(222)plot(r,R53(r).^2)xlabel(r);ylabel(R(r)^2);title(5f);grid onsubplot(223)D=r.^2.*(R53(r).^2);plot(r,D)xlabel(r);ylabel(D(r));title(5f);grid onsubplot(224)D=r.^2.*R63(r).^2;plot(r,D)xlabel(r);ylabel(D(r));title(6f);grid on 以下内容需要回复才能看到复制内容到剪贴板代码: %%figure(numbertitle,off,name,原子轨道图——by Matlabsky)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;subplot(221)R=Ypx(theta,phi);[X,Y,Z]=sph2cart(T,P,R);mesh(X,Y,Z)axis equalxlabel(x);ylabel(y);zlabel(z);title(Ypx)subplot(222)R=Ypx(theta,phi);R=R.^2;[X,Y,Z]=sph2cart(T,P,R);mesh(X,Y,Z)xlabel(x);ylabel(y);zlabel(z);title(|Ypx|^2)subplot(223)R=Ydxy(theta,phi);R=R.^2;[X,Y,Z]=sph2cart(T,P,R);mesh(X,Y,Z)xlabel(x);ylabel(y);zlabel(z);title(|Ydxy|^2)subplot(224)R=Yfzxy(theta,phi);R=R.^2;[X,Y,Z]=sph2cart(T,P,R);mesh(X,Y,Z)xlabel(x);ylabel(y);zlabel(z);title(|Yfzxy|^2)%%figure(numbertitle,off,name,电子云图——by Matlabsky)x=-15:15;y=x;[X,Y]=meshgrid(x,y);Z=zeros(size(X));psi=psi_3dxy(X,Y,Z);psi2=psi.^2;subplot(221)mesh(X,Y,psi)xlabel(x);ylabel(y);zlabel(3dxy)subplot(222)mesh(X,Y,psi2)xlabel(x);ylabel(y);zlabel(|3dxy|^2)subplot(223)contour(psi2,10);xlabel(x);ylabel(y);zlabel(3dxy)z=x;[X,Y,Z]=meshgrid(x,y,z);psi=psi_3dxy(X,Y,Z);psi2=psi.^2;subplot(224)slice(X,Y,Z,psi2,[],[],0);xlabel(x);ylabel(y);zlabel(|3dxy|^2)view(0,90);colorbar%%function y=R53(r)p=2*r/5;y=1/300/sqrt(70)*(8-p).*p.^3.*exp(-p/2);%%function y=R63(r)p=r/3;y=1/2592/sqrt(35)*(72-18*p+p.^2).*p.^3.*exp(-p/2);%%function y=Ypx(theta,phi)y=sqrt(3/4/pi)*abs(sin(theta).*cos(phi));%%function y=Ydxy(theta,phi)y=sqrt(15/16/pi)*abs(sin(theta).^2.*sin(2*phi)