Mathematica中文论坛-非官方

标题: 豆粑粑 matlab 画电子云 [打印本页]

---

作者: meatball1982    时间: 2017-7-18 15:40
标题: 豆粑粑 matlab 画电子云
别人的程序，我就是转一下。
虽然没完全看明白。

1. %%
   
2. %氢原子电子云图
   
3. %by PengjuZhang
   
4. %2012.03.02
   
5. clc;
   
6. close all;
   
7. % figure('numbertitle','off','name','波函数图');
   
8. % subplot(2,2,1);
   
9. % r=0:80;
   
10. % plot(r,R53(r));
    
11. % xlabel('r');
    
12. % ylabel('R(r)');
    
13. % title('5f');
    
14. % grid on; % 径向波函数R(r)
    
15. % subplot(2,2,2);
    
16. % plot(r,R53(r).^2);
    
17. % xlabel('r');
    
18. % ylabel('R(r)^2');
    
19. % title('5f');
    
20. % grid on; %径向波函数密度R(r)*R(r)
    
21. % subplot(2,2,3);
    
22. % D=(r.^2).*(R53(r).^2);
    
23. % plot(r,D);
    
24. % xlabel('r');
    
25. % ylabel('D(r)');
    
26. % title('5f');
    
27. % grid on;
    
28. % subplot(2,2,4);
    
29. % D=(r.^2).*(R63(r).^2);
    
30. % plot(r,D);
    
31. % xlabel('r');
    
32. % ylabel('D(r)');
    
33. % title('6f');
    
34. % grid on;%径向分布函数D(r)
    
35. 
    
36. %%
    
37. figure('numbertitle','off','name','原子轨道图');
    
38. subplot(2,2,1);
    
39. p=-pi/2:pi/60:pi/2;
    
40. t=0:2*pi/100:2*pi;
    
41. [P,T]=meshgrid(p,t);
    
42. theta=pi/2-P;
    
43. phi=T;
    
44. R=Ypx(theta,phi);
    
45. [X,Y,Z]=sph2cart(T,P,R);%球坐标转化为笛卡尔坐标的XYZ
    
46. mesh(X,Y,Z);
    
47. axis equal;
    
48. hidden off
    
49. xlabel('x'),ylabel('y'),zlabel('z');
    
50. title('Ypx');%角度波函数Y
    
51. subplot(2,2,2);
    
52. R=Ypx(theta,phi);%调用函数
    
53. R=R.*R;
    
54. [X,Y,Z]=sph2cart(T,P,R);
    
55. mesh(X,Y,Z);
    
56. axis equal;
    
57. hidden off
    
58. xlabel('x'),ylabel('y'),zlabel('z');
    
59. title('1Ypx1^2');
    
60. subplot(2,2,3);
    
61. R=Ydxy(theta,phi);
    
62. R=R.*R;
    
63. [X,Y,Z]=sph2cart(T,P,R);
    
64. mesh(X,Y,Z);
    
65. axis equal;
    
66. hidden off
    
67. xlabel('x'),ylabel('y'),zlabel('z');
    
68. title('1Ydxy1^2');
    
69. subplot(2,2,4);
    
70. R=Yfzxy(theta,phi);
    
71. R=R.*R;
    
72. [X,Y,Z]=sph2cart(T,P,R);
    
73. mesh(X,Y,Z);
    
74. axis ('equal');
    
75. hidden off
    
76. xlabel('x'),ylabel('y'),zlabel('z');
    
77. title('1Yfzxy1^2');
    
78. 
    
79. h=gcf;
    
80. fig_na='./fig_ele_cloud';
    
81. fun_work_li_035_myfig_out(h,fig_na,3)
    
82. 
    
83. %%
    
84. % figure('numbertitle','off','name','电子云图');
    
85. % x=-15:15;
    
86. % y=x;
    
87. % [X,Y]=meshgrid(x,y);
    
88. % Z=zeros(size(X));
    
89. % psi=psi_3dxy(X,Y,Z);
    
90. % psi2=psi.*psi;
    
91. % subplot(2,2,1);
    
92. % mesh(X,Y,psi);
    
93. % axis tight;
    
94. % xlabel('x'),ylabel('y'),zlabel('3dxy');
    
95. % subplot(2,2,2);
    
96. % mesh(X,Y,psi2);
    
97. % axis tight;
    
98. % xlabel('x'),ylabel('y'),zlabel('13dxy1^2');
    
99. % subplot(2,2,3);
    
100. % contour(psi2,10);
     
101. % axis tight;
     
102. % xlabel('x'),ylabel('y'),title('13dxy1^2');
     
103. % z=x;
     
104. % [X,Y,Z]=meshgrid(x,y,z);
     
105. % psi=psi_3dxy(X,Y,Z);
     
106. % psi2=psi.*psi;
     
107. % subplot(2,2,4);
     
108. % slice(X,Y,Z,psi2,[],[],[0]);
     
109. % axis tight;
     
110. % xlabel('x'),ylabel('y'),title('13dxy1^2');
     
111. % view(0,90);
     
112. % colorbar;
     
113. %
     
114. %

复制代码

fig_ele_cloud.png (432.09 KB, 下载次数: 1695)

下载附件

2017-7-18 15:40 上传

---

欢迎光临 Mathematica中文论坛-非官方 (http://ilovemathematica.com/)

Powered by Discuz! X3.2