| 查看: 1544 | 回复: 17 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【讨论】电子结构图的制作
|
|||
| 我看了一篇文献,该文章也是用castep做的,文章的电子结构图很漂亮,请高手指点该图是否是用castep做的.如何做? |
回复此楼» 猜你喜欢
面上项目申报
已经有3人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有5人回复
-
酰胺脱乙酰基
已经有9人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有6人回复
-
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
meimeidyh
铁杆木虫 (著名写手)
进城务工人员
- 应助: 116 (高中生)
- 金币: 11474.2
- 散金: 4112
- 红花: 19
- 帖子: 2702
- 在线: 702.4小时
- 虫号: 133720
- 注册: 2005-12-14
- 性别: GG
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuchenwf(金币+5,VIP+0):谢谢 7-2 22:28
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuchenwf(金币+5,VIP+0):谢谢 7-2 22:28
|
转一个好帖给你看看（用matlab做电荷密度等高线）： 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);[ Last edited by meimeidyh on 2009-7-2 at 13:17 ] |
14楼2009-07-02 13:14:26
liuxiaocunde
铁杆木虫 (著名写手)
古德里安将军
- 应助: 83 (初中生)
- 贵宾: 2.195
- 金币: 9333.7
- 散金: 3014
- 红花: 20
- 帖子: 1927
- 在线: 716小时
- 虫号: 536180
- 注册: 2008-03-30
- 性别: GG
- 专业: 固体无机化学
2楼2009-06-25 21:24:17
3楼2009-06-25 23:50:22
tcwzb168
新虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 2776.6
- 散金: 200
- 红花: 12
- 帖子: 663
- 在线: 76.8小时
- 虫号: 479526
- 注册: 2007-12-14
- 性别: GG
- 专业: 质谱分析
4楼2009-06-26 00:01:44