| 查看: 3781 | 回复: 28 | |||||||
| 【奖励】 本帖被评价18次，作者董大瓜增加金币 14.4 个 | |||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||
[资源]
MATLAB 重绘MS导出的*.csv格式的能带图，可以自动分开a画lpha和beta。
|
|||||||
|
自己做的，提供给大家一起交流学习。 调用格式如下： %% 从*.csv文件绘制Material Studio的能带图 %********************************************************************* %plot band structure form *.csv file %plotCsvBandStructure(FileName,Num_Band,id_kpts,label_kpts,TitleName,Axis) %example:plotCsvBandStructure('ZnO_BandStr Band Structure.csv',30,[1 7 21 27 41 53 59 66],'G|A|H|K|G|M|L|H','ZnO Band Structure',[0 1 -3 10]); %FileName:文件名 %Num_Band:能带条数，在*_BandStr.castep文件里Band Structure Parameters下 % 参数number of bands / k-point的值 %id_kpts:特殊K点 %label_kpts:特殊K点名称,MS能带图横坐标上显示的字母 %TitleName:能带图名称 %Axis:坐标轴范围 %Notes:MS exported csv files of band structure have two columns % first is K points positions, and second is corresponding energy % actual points number is Num_Pts-1 % by 董大瓜 %********************************************************************* function plotCsvBandStructure(FileName,Num_Band,id_kpts,label_kpts,TitleName,Axis) M=csvread(FileName); if size(M,2)==4 Draw(M(:,1:2),Num_Band,'alpha',id_kpts,label_kpts,TitleName,Axis);%alpha Draw(M(:,3:4),Num_Band,'beta',id_kpts,label_kpts,TitleName,Axis);%beta else Draw(M,Num_Band,'',id_kpts,label_kpts,TitleName,Axis); end function Draw(BandData,Num_Band,spin,id_kpts,label_kpts,TitleName,Axis) x=BandData(:,1); y=BandData(:,2); L=size(BandData,1); x=reshape(x,L/Num_Band,Num_Band); y=reshape(y,L/Num_Band,Num_Band); x(L/Num_Band,: )=[]; y(L/Num_Band,: )=[]; h=figure; plot(x,y,'k-+','LineWidth',1); ylabel('Energy(eV)'); line([0 1],[0 0],'Color','k','LineWidth',1,'LineStyle','--'); %Fermi energy line % special K-points marker loc_kpts = x(id_kpts); set(gca,'XTick',loc_kpts,'XTickLabel',label_kpts); line([loc_kpts;loc_kpts],repmat(get(gca,'yLim')',1,size(id_kpts,2)),'Color','k'); axis(Axis); set(h,'outerposition',get(0,'screensize')); if ~isempty(spin) xlabel(spin); title(strcat(TitleName,'-',spin)); saveas(h,strcat(TitleName,'-',spin,'.jpg')); else title(TitleName); saveas(h,strcat(TitleName,'.jpg')); end display('by DongJingYu');  Ni Band Structure-alpha.jpg  Ni Band Structure-beta.jpg  ZnO Band Structure.jpg[ Last edited by 董大瓜 on 2013-4-13 at 15:58 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Ni_BandStrBandStructure.csv
- 附件 2 : plotCsvBandStructure.m
2013-04-12 09:11:46, 85.32 K
2013-04-13 15:58:11, 2.08 K
» 收录本帖的淘帖专辑推荐
 第一性原理 | 科研技巧 | 我的第一原理学习贴 | MS & VASP |
| 兴趣 |
» 本帖已获得的红花（最新10朵）
董大瓜» 猜你喜欢
还有，jjw不发布ABC的动机是什么呢？
已经有6人回复
-
像我这样的人申清北博士是不是痴心妄想
已经有14人回复
-
感觉在瑞典不适合做科研
已经有19人回复
-
套磁信没有回复
已经有9人回复
-
核磁解谱
已经有17人回复
-
博士毕业16，两中两民；
已经有10人回复
-
浙江大学生物学博导推荐
已经有6人回复
-
作者署名：nature大子刊共一还是共同通讯？
已经有21人回复
-
跑板分纯
已经有3人回复
-
有没有便宜一点的开源期刊
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 谈谈工科学生如何学习数学-邹谋炎
已经有19人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
- 理工科用神舟笔记本大家看看如何
已经有32人回复
7楼2013-04-12 09:51:38
2楼2013-04-12 09:26:14
6楼2013-04-12 09:33:50
11楼2013-04-12 20:30:26
简单回复
qyyos3楼
2013-04-12 09:33
回复
五星好评 顶一下，感谢分享！
beekingxxx6912楼
2013-04-14 11:35
回复
五星好评 顶一下，感谢分享！
han2008200813楼
2013-04-14 21:17
回复
五星好评 顶一下，感谢分享！
han2008200814楼
2013-04-14 21:17
回复
顶一下，感谢分享！
50