| 查看: 4399 | 回复: 8 | |||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||
guolianshun金虫 (著名写手)
|
[交流]
用Origin画MS中电子云密度的技术！！！超漂亮 已有8人参与
|
||||||||||
|
用Origin画MS中电子云密度的技术！！！超漂亮 有心人可以研究一下！！！ 版主请加精啊， 费了好大力气的！ 呵呵！ 欢迎大家学习！ MATLAB程序用来处理MS电荷密度的等高线做法，分四步： 1.找到电荷密度图的存储文件，一般是在自己所建的计算project中，隐藏的，castep模块是.charg_frm格式的，可以用txt打开后另存为txt格式。（这个地方注意把文件单独拷贝出来，以免破坏原始文件） 2.读取数据，在matlab命令里输入一下命令： [path,fn]=uigetfile('*.txt','Open'); fp=fopen([fn,path],'r'); head=fscanf(fp,'%s',4); data=fscanf(fp,'%f',[4,18*18*144]); fclose(fp); 这是读取刚才保存的txt数据的，其中18*18*144是可以改的，分别对应a,b,c的重复单元。 3.作矩阵。 a=reshape(data(1,,[18,18,144]); b=reshape(data(2,,[18,18,144]); c=reshape(data(3,,[18,18,144]); d=reshape(data(4,,[18,18,144]); 其中18*18*144是可以改的，分别对应a,b,c的重复单元，即长宽高的值。 4.作等高线图 layer=1; contour(reshape(a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8); layer=1（2，3，4....） 可以修改作图的层数，分别是从前面到后面。最后的8表示显示线条数目，可以自己修改成9，10或其它。 hold on 表示累积作图，hold off表示清空前面的 作图 举例；做表面层和第九层的等高线图 >> layer=9; contour(reshape(a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8); >> hold on >> layer=1; contour(reshape(a(:,layer,  ,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8);这是沿侧面做等高线（即xz平面），其它方向作图类似。只要修改其中的参数就可以了，如把a(:,layer, ,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8);改为a(layer,:, ,[18,144]),reshape(c(layer,l:,,[18,144]),reshape(d(layer,:,,[18,144]),8); 一个牛人的总结，转自分子论坛。 |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 | 物理前沿 | 第一性原理计算 | @个人收集@ |
| Photochemistry | 计算化学相关 | 精华帖子 | 第一性 |
» 本帖已获得的红花（最新10朵）
awmc2008» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有4人回复
-
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- origin画ms电荷密度图
已经有21人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- 新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度
已经有15人回复
- 【讨论】请教用origin画一种图
已经有15人回复
- 一般用什么软件计算离子的键能和电子云密度
已经有11人回复
- 【求助】band.dat用origin画能带图
已经有5人回复
- 【求助】求电子云密度大小，满意追加金币
已经有10人回复
- 用origin画了个图，似乎里面有Type 3 字体，生成pdf以后，提交的时候通不过。
已经有22人回复
- 【求助】如何采用Matlab绘制电子云密度图
已经有8人回复