Matlab处理气象数据（十二）DTR及其异常的计算

DTR（diurnal temperature range）是日较差，即最高温和最低温度的差值。一段时期平均DTR的增加或降低可以反映出一个区域内的气候变化状况。

load('Temmax1.mat');
load('Temmax2.mat');
load('Temmin1.mat');
load('Temmin2.mat');
DTR1=Temmax1-Temmin1;
DTR2=Temmax2-Temmin2;

%画两套数据的DTR
plot(DTR1,'r.-','linewidth',2);
hold on
plot(DTR2,'b.-','linewidth',2);
legend({'NCEP','Observed'},'Location','Northeast');
set(gca,'xtick',[2 7 12 17 22 27 32],'xticklabel',{'1980','1985','1990','1995','2000','2005','2010'});%在x轴特定位置上添加标注
set(gca, 'FontSize',10,'FontWeight','Bold','tickdir','out') %设置标注为10号字、加粗、标记线向外
h=xlabel('Year'); %设置x轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
h=ylabel('Temperarure(\circC)'); %设置y轴名称
set(h, 'FontSize',10,'FontWeight','Bold') 
xlim([1 33])%x轴范围锁定为1~33
box off %去掉外框

两套数据的DTR

%画两套数据DTR的异常及其趋势线
aDTR1=DTR1-mean(DTR1);
aDTR2=DTR2-mean(DTR2);
plot(aDTR1,'r.-','linewidth',2);%画NCEP数据最高气温折线图，红色实线实心点
hold on
plot(aDTR2,'b.-','linewidth',2);
axis([ -inf inf -1 1]);
legend({'NCEP','Observed'},'Location','Northwest');
set(gca,'xtick',[2 7 12 17 22 27 32],'xticklabel',{'1980','1985','1990','1995','2000','2005','2010'});%在x轴特定位置上添加标注
set(gca, 'FontSize',10,'FontWeight','Bold','tickdir','out') %设置标注为10号字、加粗、标记线向外
h=xlabel('Year'); %设置x轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
h=ylabel('Temperarure(\circC)'); %设置y轴名称
set(h, 'FontSize',10,'FontWeight','Bold') 
xlim([1 33])%x轴范围锁定为1~33
box off %去掉外框
plot(aDTR2-detrend(aDTR2),'b--','linewidth',2);
plot(aDTR1-detrend(aDTR1),'r--','linewidth',2);
hold off

两套数据的DTR距平及其趋势线

---

相关链接：
Matlab处理气象数据——目录