| 查看: 3470 | 回复: 22 | ||||
| 【奖励】 本帖被评价6次，作者monowolf增加金币 3.75 个 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[资源]
电荷密度图等高线的作法
|
||||
|
这是我师兄编的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);欢迎大家指点。上传附件为样本数据和上述作图txt文档。 特别感谢师兄Mr.Tang. [ Last edited by monowolf on 2008-9-25 at 21:58 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 我的第一原理学习贴 | 第一原理资料贴 |
» 猜你喜欢
309分085801求调剂
已经有7人回复
-
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
-
环境工程 085701，267求调剂
已经有14人回复
-
308求调剂
已经有12人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有21人回复
-
311求调剂
已经有8人回复
-
294分080500材料科学与工程求调剂
已经有13人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有19人回复
-
288资源与环境专硕求调剂，不限专业，有学上就行
已经有23人回复
11楼2009-05-31 10:49:40
2楼2008-09-25 21:59:25
3楼2008-09-25 22:02:57
4楼2008-09-26 14:40:58
