条形图
%% 绘制条形图
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)