Mathematica中文论坛-非官方

标题: 豆粑粑　matlab 画势能面 [打印本页]

---

作者: meatball1982    时间: 2018-11-21 09:33
标题: 豆粑粑　matlab 画势能面
别人的问题，记录一下。

这是个记录的帖子：
等我不在这个世界上的时候，希望有人能通过这个帖子的内容，少走些弯路，多学些。


用surf画势能面。分别调整透明程度。
用surf画原子。
用text标注。
设置其它。

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. %% outline
   
6. % plot diff layer potential surfs
   
7. 
   
8. %% main
   
9. % the x,y
   
10. x=[1:10];
    
11. y=[0:8];
    
12. [X,Y]=meshgrid(x,y);
    
13. 
    
14. % your potential values,
    
15. Z1  = sin(0.2*X+0.7*Y );
    
16. Z2  = exp(0.1*X.*cos(0.5*Y.^1.1))+3;
    
17. Z3  = 0.1*((X-5).^2 + (Y-6).^2)+7;
    
18. 
    
19. % atom x,y,z
    
20. [s_x,s_y,s_z]=sphere(50);
    
21. x_shift = 3;
    
22. y_shift = 2;
    
23. z_shift = 12;
    
24. r=0.6;
    
25. 
    
26. 
    
27. %% plot
    
28. hold on
    
29. % surf 3 layers
    
30. h1=surf(X,Y,Z1);
    
31. h2=surf(X,Y,Z2,Z2+10); % you can modify 10 to change the color of the layer
    
32. h3=surf(X,Y,Z3,Z3+20); % you can modify 20 to change the color of the layer
    
33. 
    
34. alpha(h1,0.95);
    
35. alpha(h2,0.9);
    
36. alpha(h3,0.8);
    
37. 
    
38. % plot atom
    
39. surf(r*s_x+x_shift,r*s_y+y_shift,r*s_z+z_shift,'edgecolor','none','markerfacecolor','w')
    
40. 
    
41. % Text
    
42. text(3,4,10,'S_2','fontsize',16)
    
43. 
    
44. 
    
45. % plot point on surf
    
46. plot3(X(6,5),Y(6,5),Z3(6,5),'yo','markerfacecolor','y','markersize',10);
    
47. 
    
48. % other settings
    
49. view(150,30);
    
50. axis off
    
51. axis equal
    
52. camlight
    
53. 
    
54. 
    
55. h=gcf;
    
56. fi_na = 'fig_3layers_pot'
    
57. fun_work_li_035_myfig_out(h,fi_na,3)
    
58. %% logs
    
59. % mod : 2018年 11月 21日 星期三 08:52:09 CST
    
60. % mod by mm
    
61. % contact : meatball1982@163.com

复制代码

fig_3layers_pot.png (145.11 KB, 下载次数: 2174)

下载附件

2018-11-21 09:31 上传

---

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

Powered by Discuz! X3.2