记一个绘制态密度与能带的matlab脚本

题记

记一个能带与态密度绘制的脚本，用于处理p4VASP提取的数据后的数据（需要p4vasp的可以在我博客里找）

用法

能带数据需保存在名为band.dat的文件里
态密度数据保存在名为dos.dat的文件里

输入

脚本当中有些数据需要先手动输入（共3处）
比如排列方式
能带的xy轴的取值范围
能带大高对称点

%%%%%%%%%%%%%%%%%%%%能带%%%%%%%%%%%%%%%
dat= importdata('band.dat');%输入能带数据
x=dat(:,1);
y=dat(:,2);
[m,n]=size(x);
mun1=1;
subplot(1,2,1)%能带图的放置位置，这里1，2，1分别表示把图按1行2列的形式排列，其中这幅图排1的位置
hold on
for i=2:m
    if x(i)== 0
        plot(munx,muny,'k','linewidth',1)%绘制能带图默认颜色为黑色‘k’线的粗细为1
        hold on
        mun1=1;
        x_long=x(i-1);
        clear munx muny
        i=i+1;
    else
        munx(mun1)=x(i);
        muny(mun1)=y(i);
        mun1=mun1+1;
        i=i+1;
    end
end
%%%%%%%%%%%%%%%%%%%%%输入%%%%%%%