版块导航: 正在加载中...

客户端APP下载

登录注册

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

24小时热门版块排行榜

【奖励】本帖被评价18次，作者董大瓜增加金币 14.4 个

董大瓜

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 3194.1
帖子: 158
在线: 57.9小时
虫号: 2326901

[资源] 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朵）

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响已经有17人回复

1楼2013-04-12 09:12:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

董大瓜

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 3194.1
帖子: 158
在线: 57.9小时
虫号: 2326901

笑脸的地方时这样的。一个冒号和右半块括号
x(L/Num_Band,: )=[];
y(L/Num_Band,: )=[];

赞一下(1人)

2楼2013-04-12 09:26:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wei520jing

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 373.6
帖子: 31
在线: 9.3小时
虫号: 1833940

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

师弟，继续努力……前途无量！

6楼2013-04-12 09:33:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

董大瓜

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 3194.1
帖子: 158
在线: 57.9小时
虫号: 2326901

引用回帖:

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

铁杆木虫 (著名写手)

应助: 41 (小学生)
金币: 6014.2
帖子: 1875
在线: 940小时
虫号: 321641

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

继续努力……前途无量！

11楼2013-04-12 20:30:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

董大瓜

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 3194.1
帖子: 158
在线: 57.9小时
虫号: 2326901

顶起来

17楼2013-04-19 13:56:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

董大瓜

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 3194.1
帖子: 158
在线: 57.9小时
虫号: 2326901

顶起来

18楼2013-04-23 10:35:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

董大瓜

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 3194.1
帖子: 158
在线: 57.9小时
虫号: 2326901

送红花一朵

送花一朵

20楼2013-05-13 16:54:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

简单回复

qyyos3楼

2013-04-12 09:33 回复

五星好评顶一下，感谢分享！

2013-04-12 09:33 回复

五星好评

qyyos5楼

2013-04-12 09:33 回复

2013-04-12 09:57 回复

五星好评谢谢分享

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 回复

顶一下，感谢分享！

2013-04-17 17:08 回复

五星好评顶一下，感谢分享！

2013-04-17 20:46 回复

顶一下，感谢分享！

2013-05-03 07:24 回复

小小冰竹21楼

2013-05-14 11:07 回复

五星好评顶一下，感谢分享！

2013-07-25 08:22 回复

顶

嗯嗯浅悟23楼

2013-12-27 19:35 回复

五星好评顶一下，感谢分享！

2013-12-28 23:30 回复

五星好评好 [ 发自小木虫客户端 ]

boshin_liu25楼

2014-01-13 09:27 回复

五星好评顶一下，感谢分享！

2014-03-31 15:02 回复

五星好评顶一下，感谢分享！

2014-04-15 18:15 回复

五星好评顶一下，感谢分享！

2015-04-22 10:16 回复

五星好评顶一下，感谢分享！

xx140829楼

2019-11-17 16:24 回复

五星好评顶一下，感谢分享！

相关版块跳转我要订阅楼主董大瓜的主题更新