应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » MATLAB 重绘MS导出的*.csv格式的能带图，可以自动分开a画lpha和beta。
291/1返回列表
查看: 3780  |  回复: 28
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价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
    • 2013-04-12 09:11:46, 85.32 K
  • 附件 2 : plotCsvBandStructure.m
    • 2013-04-13 15:58:11, 2.08 K

» 收录本帖的淘帖专辑推荐

第一性原理 科研技巧 我的第一原理学习贴 MS & VASP
兴趣

» 本帖已获得的红花（最新10朵）

董大瓜

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-04-12 09:12:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

董大瓜

新虫 (小有名气)
笑脸的地方时这样的。一个冒号和右半块括号
x(L/Num_Band,: )=[];
y(L/Num_Band,: )=[];
赞一下(1人) 回复此楼
2楼2013-04-12 09:26:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei520jing

铁虫 (初入文坛)

★★★★★ 五星级,优秀推荐

师弟，继续努力……前途无量！
赞一下 回复此楼
6楼2013-04-12 09:33:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

董大瓜

新虫 (小有名气)
引用回帖:
6楼: Originally posted by wei520jing at 2013-04-12 09:33:50
师弟，继续努力……前途无量！

刚开始学MS。。。新手，向大师学习。。。
赞一下 回复此楼
7楼2013-04-12 09:51:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)

★★★★★ 五星级,优秀推荐

继续努力……前途无量！
赞一下 回复此楼
11楼2013-04-12 20:30:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

董大瓜

新虫 (小有名气)
顶起来
回复此楼
17楼2013-04-19 13:56:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

董大瓜

新虫 (小有名气)
顶起来
回复此楼
18楼2013-04-23 10:35:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

董大瓜

新虫 (小有名气)
送红花一朵
送花一朵
回复此楼
20楼2013-05-13 16:54:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
qyyos3楼
2013-04-12 09:33   回复  
五星好评  顶一下，感谢分享！
sunguov4楼
2013-04-12 09:33   回复  
五星好评  
qyyos5楼
2013-04-12 09:33   回复  
wangjb4018楼
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   回复  
五星好评  顶一下，感谢分享！
相关版块跳转 我要订阅楼主 董大瓜 的主题更新
291/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 感觉在瑞典不适合做科研 +14 毕生所学 2024-09-19 19/950 2024-09-22 15:08 by 毕生所学
[有机交流] 核磁解谱 50+4 luffyyyyy 2024-09-16 17/850 2024-09-22 13:33 by luffyyyyy
[论文投稿] 作者署名：nature大子刊共一还是共同通讯？ 10+8 Eriknpeng 2024-09-16 21/1050 2024-09-22 12:29 by gdd2018
[有机交流] 跑板分纯 15+3 DY12@ 2024-09-21 3/150 2024-09-22 10:59 by shouwe
[论文投稿] 有没有便宜一点的开源期刊 100+5 zhuzhuyaya 2024-09-19 8/400 2024-09-22 10:58 by shouwe
[考研] 大家慎重选择南科大李文甲 +3 小红豆 2024-09-21 7/350 2024-09-22 10:37 by chen5805
[论文投稿] 论文投稿 5+4 love319 2024-09-20 6/300 2024-09-22 10:14 by bobvan
[职场人生] 直系领导让我做他的博后 +3 木之角 2024-09-21 3/150 2024-09-22 09:37 by ciompman
[论文投稿] 请问有composites science and technology投稿模板吗？ 5+3 固定资产 2024-09-21 3/150 2024-09-22 09:29 by nono2009
[教师之家] 大专在职教师是否需要读一个在职博士？ +19 伍乘柒35 2024-09-19 35/1750 2024-09-22 08:36 by Quakerbird
[有机交流] 如何除去高沸点有机溶剂 10+3 lyfsl 2024-09-19 4/200 2024-09-22 06:15 by xiaomei1031
[论文投稿] PCCP投稿 +5 wuxiaohuang1 2024-09-20 5/250 2024-09-21 23:38 by 超爱烤冷面
[硕博家园] 关于博士选择 +4 翰墨古风 2024-09-20 9/450 2024-09-21 20:54 by 翰墨古风
[硕博家园] 现博士在读，硕士学校的硕导又让干活 +8 也就这样 2024-09-18 9/450 2024-09-21 18:45 by 芒果与橙子
[论文投稿] 提交proof后，SI还可以修改吗？ 20+4 力行成士 2024-09-20 6/300 2024-09-21 17:32 by 莱茵润色
[论文投稿] 顶级水刊，几乎不拒稿的SCI! +5 a433587770 2024-09-18 5/250 2024-09-21 12:39 by 88817753
[基金申请] 今年中了国家课题，明年聘讲师，国家级课题可以用于进副教授吗 20+10 春天播种，秋天� 2024-09-16 23/1150 2024-09-21 11:28 by gwmgyp
[论文投稿] 文章改了几轮，该投稿了，选了几篇杂志发给导师 +5 evan0489 2024-09-19 6/300 2024-09-20 20:13 by lyd782
[基金申请] 山东省自然科学基金有查到专家意见么 +12 乱码哈哈哈 2024-09-18 16/800 2024-09-20 10:30 by andywei1028
[基金申请] 结项4年了，论文还能挂基金号吗 +10 xiandaiqi 2024-09-15 12/600 2024-09-19 22:00 by yuch
信息提示
关闭
请填处理意见
关闭