| 查看: 4309 | 回复: 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朵）
» 猜你喜欢
读博
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有13人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
origin画ms电荷密度图
已经有21人回复- 【活动】Show 脚本，赢奖励
已经有23人回复
- 新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度
已经有15人回复
- 【讨论】请教用origin画一种图
已经有15人回复
- 一般用什么软件计算离子的键能和电子云密度
已经有11人回复
- 【求助】band.dat用origin画能带图
已经有5人回复
- 【求助】求电子云密度大小，满意追加金币
已经有10人回复
- 用origin画了个图，似乎里面有Type 3 字体，生成pdf以后，提交的时候通不过。
已经有22人回复
- 【求助】如何采用Matlab绘制电子云密度图
已经有8人回复
★
小木虫(金币+0.5):给个红包，谢谢回帖
小木虫(金币+0.5):给个红包，谢谢回帖
|
本帖内容被屏蔽 |
2楼2011-10-11 15:40:54
YUN杰2070
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 2376.3
- 红花: 3
- 帖子: 119
- 在线: 95.6小时
- 虫号: 1279373
- 注册: 2011-04-27
- 性别: GG
- 专业: 无机纳米化学
3楼2011-10-28 21:20:32
ustc
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 9239.4
- 散金: 207
- 红花: 7
- 帖子: 1218
- 在线: 480.4小时
- 虫号: 28859
- 注册: 2003-11-13
- 性别: GG
- 专业: 催化化学
4楼2011-10-28 21:27:57
送鲜花一朵
 |
5楼2011-10-29 23:07:11
风之梦
木虫 (小有名气)
- 应助: 37 (小学生)
- 金币: 5863.2
- 红花: 6
- 帖子: 278
- 在线: 352.4小时
- 虫号: 211914
- 注册: 2006-03-08
- 性别: GG
- 专业: 无机非金属材料
6楼2011-10-31 09:27:44
fanny0007
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 972.6
- 散金: 12
- 帖子: 220
- 在线: 120.5小时
- 虫号: 533141
- 注册: 2008-03-26
- 性别: MM
- 专业: 电路与系统
7楼2011-10-31 10:24:55
余昭福
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 0.9
- 散金: 280
- 帖子: 69
- 在线: 27.6小时
- 虫号: 3079522
- 注册: 2014-03-22
- 性别: GG
- 专业: 金属材料的力学行为
8楼2015-11-09 11:12:55
9楼2019-01-03 13:14:01