分子谱线数据绘图——matab

最新推荐文章于 2021-11-13 21:28:28 发布

原创最新推荐文章于 2021-11-13 21:28:28 发布 · 654 阅读

5 ·

CC 4.0 BY-SA版权

code 专栏收录该内容

15 篇文章

订阅专栏

本文通过使用Matlab绘制分子谱线数据图，并与Python的matplotlib库进行比较，展示了如何导入不同文件中的数据并利用Matlab进行图表的精确定制，包括设置坐标轴范围、刻度、图例等。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

用了分子谱线数据，matlab画图与python的matplotlib库函数是一样的。

%import data
[x80,y80]=textread('Adata_cut/A8_158228-158680.txt');
[x81,y81]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_A/A8/xclass_spectrum_output.dat');
[x82,y82]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_A/A8/xclass_spectrum_output.dat');
[x83,y83]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_A/A8/xclass_spectrum_output.dat');
[x84,y84]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_A/A8/xclass_spectrum_output.dat');
[x85,y85]=textread('Bdata_cut/B8_158228-158680.txt');
[x86,y86]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__C2H5OH_B/B8/xclass_spectrum_output.dat');
[x87,y87]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3COCH3_B/B8/xclass_spectrum_output.dat');
[x88,y88]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v18=1_B/B8/xclass_spectrum_output.dat');
[x89,y89]=textread('/home/luoluo/Documents/python_scripts/job__XCLASS/job__CH3OCHO_v=0_B/B8/xclass_spectrum_output.dat');

%plot
figure(1);
subplot(2,1,1);
subplot('Position',[0.13 0.52 0.77 0.37]);
plot(x80,y80,'k-',x81,y81,'r-',x82,y82,'g-',x83,y83,'cyan',x84,y84,'b','LineWidth',1.5);
axis([158228 158680 -1 15]);
set(gca,'xlim',[158228 158680]);
set(gca,'xTick',[158228:50:158680]);
set(gca,'yLim',[0 15]);
set(gca,'yTick',[-1:2:15]);
xlabel('Freq(MHz)');
ylabel('T(K)');
legend('DATA','C2H5OH','CH3COCH3','CH3OCHO v18=1','CH3OCHO v=0','Location','northwest','Orientation','vertical');
text(158630,13,'IRAS16293A');
%fontsize



subplot(2,1,2);
plot(x85,y85,'k-',x86,y86,'r-',x87,y87,'g-',x88,y88,'cyan',x89,y89,'b','LineWidth',1.5);
axis([158228 158680 -1 15]);
set(gca,'xlim',[158228 158680]);
set(gca,'xTick',[158228:50:158680]);
set(gca,'yLim',[0 15]);
set(gca,'yTick',[-1:2:15]);
xlabel('Freq(MHz)');
ylabel('T(K)');
legend('DATA','C2H5OH','CH3COCH3','CH3OCHO v18=1','CH3OCHO v=0','Location','northwest','Orientation','vertical');
text(158630,13,'IRAS16293B');