| 查看: 4943 | 回复: 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» 猜你喜欢
确定了，国自然21号放榜
已经有4人回复
-
静等基金结果
已经有16人回复
-
filecode
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- origin画ms电荷密度图
已经有21人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- 新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度
已经有15人回复
- 【讨论】请教用origin画一种图
已经有15人回复
- 一般用什么软件计算离子的键能和电子云密度
已经有11人回复
- 【求助】band.dat用origin画能带图
已经有5人回复
- 【求助】求电子云密度大小，满意追加金币
已经有10人回复
- 用origin画了个图，似乎里面有Type 3 字体，生成pdf以后，提交的时候通不过。
已经有22人回复
- 【求助】如何采用Matlab绘制电子云密度图
已经有8人回复
YUN杰2070
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 2376.3
- 红花: 3
- 帖子: 119
- 在线: 95.6小时
- 虫号: 1279373
- 注册: 2011-04-27
- 性别: GG
- 专业: 无机纳米化学
3楼2011-10-28 21:20:32
★
小木虫(金币+0.5):给个红包，谢谢回帖
小木虫(金币+0.5):给个红包，谢谢回帖
|
本帖内容被屏蔽 |
2楼2011-10-11 15:40:54
风之梦
木虫 (小有名气)
- 应助: 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
