matlab脚本画投影能带没有横坐标和网格线

Qiang0508

投影能带

您好，我用vaspkit处理投影能带，用群里提供的matlab脚本画图，但是画出来的图片为什么和给的示例不一样，没有横坐标，也没有图片中的白色网格线，麻烦您帮忙看一下，谢谢。我把我的文件都上传到附件了，麻烦帮忙看一下，谢谢。

Dirac-learn

1. % First run vaspkit (version => 1.2.3) with task 216, 257 or 285 to get MOMENTUM.grd, ENERGY.grd and WEIGHT.grd
   
2. 
   
3. clc;clear;clear all;
   
4. 
   
5. momentum=load('MOMENTUM.grd');
   
6. energy=load('ENERGY.grd');
   
7. weight=load('WEIGHT.grd');
   
8. ymin=-5.0;
   
9. ymax=5.0;
   
10. fig_title =' Nb-d'    ;
    
11. 
    
12. % Read content from KLABELS, assign values to kpoint_value and k_value
    
13. fid = fopen('KLABELS','r');%打开文件句柄
    
14. FormatString='%s';
    
15. High_Symmetry_Kpoints=textscan(fid,FormatString,'HeaderLines',1,'Delimiter',''); %跳过第1行，以空格为数据的分隔符
    
16. High_Symmetry_Kpoints{1}(end)=[];  %删除最后一行
    
17. %初始化
    
18. kpoint_value = cell(1,length(High_Symmetry_Kpoints{1}));
    
19. k_value = length(High_Symmetry_Kpoints{1});
    
20. fclose(fid);%关闭文件句柄
    
21. for i = 1:length(High_Symmetry_Kpoints{1})
    
22. a  = strsplit(High_Symmetry_Kpoints{1,1}{i,1},' ');  
    
23. kpoint_value(i)  =  a(1);
    
24. k_value(i) = str2double(cell2mat(a(2)));
    
25. end            
    
26. momentum_min = min(momentum(:));
    
27. momentum_max = max(momentum(:));
    
28. k_value(length(k_value))=momentum_max;
    
29. k_value(1)=momentum_min;
    
30. energy_min   = min(energy(:));
    
31. energy_max   = max(energy(:));
    
32. colormap('jet')
    
33. set(gca,'Fontsize',20)
    
34. 
    
35. %subplot(1, 5,subfig_number)
    
36. fig = pcolor(momentum, energy, weight);
    
37. shading interp;
    
38. set(fig, 'LineStyle', 'none')
    
39. title(fig_title)
    
40. ylabel('E-E_F (eV)')
    
41. ylim([ymin ymax]);
    
42. set(gca,'ytick',ymin:1:ymax)
    
43. set(gca, 'XTickLabel', kpoint_value , 'XTick', k_value,'Fontname', 'Times newman','Fontsize',20)
    
44. 
    
45. line([momentum_min momentum_max], [0 0], 'LineStyle', '--', 'Color', 'w', 'LineWidth',1)                                
    
46. for a=2:1:length(k_value)-1
    
47. line([1 1]*k_value(a), [energy_min energy_max], 'LineStyle', '--', 'Color', 'w', 'LineWidth',1)
    
48. end
    

复制代码

Dirac-learn

简单修改了一下，你用这个代码算一下看看

Qiang0508

Dirac-learn 发表于 2023-5-24 23:08
简单修改了一下，你用这个代码算一下看看

太感谢您了，我去试一下，多谢