| 查看: 4537 | 回复: 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朵）
董大瓜» 猜你喜欢
时间戳又变了
已经有15人回复
-
科研孤儿太难了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
2楼2013-04-12 09:26:14
6楼2013-04-12 09:33:50
7楼2013-04-12 09:51:38
11楼2013-04-12 20:30:26
17楼2013-04-19 13:56:10
18楼2013-04-23 10:35:51
20楼2013-05-13 16:54:02
简单回复
qyyos3楼
2013-04-12 09:33
回复
五星好评 顶一下，感谢分享！
sunguov4楼
2013-04-12 09:33
回复
五星好评 
qyyos5楼
2013-04-12 09:33
回复
2013-04-12 09:57
回复
五星好评 谢谢分享
awmc20089楼
2013-04-12 11:53
回复
五星好评 
chuzhaonan10楼
2013-04-12 19:10
回复
五星好评
beekingxxx6912楼
2013-04-14 11:35
回复
五星好评 顶一下，感谢分享！
han2008200813楼
2013-04-14 21:17
回复
五星好评 顶一下，感谢分享！
han2008200814楼
2013-04-14 21:17
回复
顶一下，感谢分享！
xylz618815楼
2013-04-17 17:08
回复
五星好评 顶一下，感谢分享！
xylz618816楼
2013-04-17 20:46
回复
顶一下，感谢分享！
董大瓜19楼
2013-05-03 07:24
回复
小小冰竹21楼
2013-05-14 11:07
回复
五星好评 顶一下，感谢分享！
董大瓜22楼
2013-07-25 08:22
回复
顶
嗯嗯浅悟23楼
2013-12-27 19:35
回复
五星好评 顶一下，感谢分享！
戴世杰24楼
2013-12-28 23:30
回复
五星好评 好 [ 发自小木虫客户端 ]
boshin_liu25楼
2014-01-13 09:27
回复
五星好评 顶一下，感谢分享！
szzpanda26楼
2014-03-31 15:02
回复
五星好评 顶一下，感谢分享！
suibianba27楼
2014-04-15 18:15
回复
五星好评 顶一下，感谢分享！
sam_lucy28楼
2015-04-22 10:16
回复
五星好评 顶一下，感谢分享！
xx140829楼
2019-11-17 16:24
回复
五星好评 顶一下，感谢分享！
40