| 查看: 4333 | 回复: 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朵）
» 猜你喜欢
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有18人回复
-
垃圾破二本职称评审标准
已经有12人回复
-
职称评审没过，求安慰
已经有17人回复
-
EST投稿状态问题
已经有7人回复
-
谈谈两天一夜的“延安行”
已经有15人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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