Matlab 科研绘图汇总

 条形图

%% 绘制条形图
figure
x=1:3;
y=indicator_ztd_all;
h=bar(x,y);
ylabel('Accurasy(mm)');
xlabel('Map Function');
set(gca,'XTickLabel',type); %type=[string('GMF'),string('NMF'),string('VMF1')];
ylim([-2 10])  %y轴范围

%柱形图数字注记
for i=1:length(x)
text(x(i)+0.1,-1.3,num2str(y(i,1),'%.2f'),'HorizontalAlignment','center','VerticalAlignment','bottom')
text(x(i)-0.05,y(i,2)+0.1,num2str(y(i,2),'%.2f'),'HorizontalAlignment','center','VerticalAlignment','bottom')
text(x(i)+0.25,y(i,3)+0.1,num2str(y(i,3),'%.2f'),'HorizontalAlignment','center','VerticalAlignment','bottom')
end
%设置figure图像大小
set(gcf,'unit','centimeters','position',[6,4,12,6.18])
%设置图例
lgd=legend('bias','std','rms');
%设置图例大小和位置
 set(lgd,'unit','centimeters','position',[7.85,3,0.5,0.5]); 

 

折线图

 %绘制折线图
    plot(ztdSeris{1,1}(:,1),ztdSeris{1,1}(:,2),'LineWidth',1)
    hold on
    plot(ztdSeris{2,1}(:,1),ztdSeris{2,1}(:,2),'LineWidth',1)
    hold on
    plot(ztdSeris{3,1}(:,1),ztdSeris{3,1}(:,2),'LineWidth',1)
    
    set(gca,'FontSize',10);
    xlabel('Doy');
    ylabel('Ztd Bias(mm)');
    xlim([143 150])   %x轴范围
    grid on;
    set(gcf,'unit','centimeters','position',[6,4,8,4.944]) 
    
    lgd=legend('GMF','NMF','VMF1');
    set(lgd,'unit','centimeters','position',[5.5,3.7,0.5,0.5],'FontSize',8);
    %     set(lgd,'Orientation','horizon');  横排
    picff=sprintf('%s%s_map',savedir,list{i,1})
    print(gcf,picff,'-r300','-dtiff');

 线性拟合相关系数图

GNSS_PWV=diff(:,4);
Radio_PWV=diff(:,5);
rms=sqrt(sumsqr(diff(:,6))/(length(diff(:,6))));
%相关系数
corrArray=corrcoef(GNSS_PWV,Radio_PWV);
corr=corrArray(1,2);

%绘制线性拟合图
t=polyfit(GNSS_PWV,Radio_PWV,1);
plot(GNSS_PWV,Radio_PWV,'.',GNSS_PWV,polyval(t,GNSS_PWV),'MarkerSize',24,'LineWidth',3)
% Label axes
xlabel('GNSS/PWV (mm)');
ylabel('Radio/PWV (mm)');
set(gca,  'FontName','times new roman','FontSize',10) ;
set(gcf,'unit','centimeters','position',[6,4,14,7])
%绘制文字框
RMS='RMS:';
RMS=[RMS,num2str(roundn(rms,-2)),' mm'];
Correlation='Correlation:';
Correlation=[Correlation,num2str(roundn(corr,-4))];
strings={RMS;Correlation};
annotation(  ...
    'textbox',[0.15,0.75,0.25,0.15],...
    'String',strings,...
    'LineWidth',1,...
    'FontSize',10,...
    'FontName','times new roman',...
    'FitBoxToText','off');
grid on
%限制x轴范围
min=min(GNSS_PWV);
max=max(GNSS_PWV);
xlim([min max]);
picfn=sprintf('seris_%d%3.3d_%d%3.3d%',sy,sdoy,ey,edoy);
picff=fullfile(workDir,picfn);
print(gcf,picff,'-r300','-dpng');
close(gcf)