泰勒图-Matlab_脚本
脚本：
clear
%data = xlsread('C:\Users\Liuy\Desktop\泰勒图制作\cfmd数据集验证.xlsx','⽇尺度','B2:NB3');
%data = xlsread('C:\Users\Liuy\Desktop\泰勒图制作\Acell76.xlsx','年降⾬','B2:AB24');
data = xlsread('C:\Users\Liuy\Desktop\泰勒图制作\Tcell139.xlsx','按照逐⽉（单点）','B2:LM3');
%data = xlsread('C:\Users\Liuy\Desktop\泰勒图制作\Acell139.xlsx','按照多年⽇平均（单点）','B2:NB3');
%获得统计结果
for ii = 2 : size(data,1)
S = allstats(data(1,:),data(ii,:));%allstats.m⽂件获取统计值
MYSTATS(ii,:) = S(:,2); % We get stats versus reference
end      %for iserie
MYSTATS(1,:) = S(:,1); % We assign reference stats to the first row
%  taylordiag(MYSTaTS(:,2),MYSTaTS(:,3),MYSTaTS(:,4));
stdev=MYSTATS(1,2);%此处时参考值的标准差
MYSTATS(:,2)=MYSTATS(:,2)/stdev;% 标准化处理，每个模式的标准差/参考值标准差
MYSTATS(:,3)=MYSTATS(:,3)/stdev;% 标准化处理，每个模式的均⽅根误差/参考值标准差
%作图
%折线图（subplot⼦图）
figure
iw=2;jw=1;
ax = ptable([2 3],[2 3;4 6]); %⼦图位置及⼤⼩x
alphab = 'ABCDEFG';    %添加图编号（预设） A B
subplot(iw,jw,1); %折线图位置
plot(data');%转置
grid on,xlabel('day');ylabel('precipitation (mm)');
title(sprintf('%s: These are the different time series of daily precipitation ','A'),'fontweight','bold');
hold on
subplot(iw,jw,2); %泰勒图展⽰位置
%泰勒图标签中使⽤了很多默认参数，可以参考taylordiag⽂件中定义，进⾏更改。
[pp tt axl] =taylordiag(MYSTATS(:,2),MYSTATS(:,3),MYSTATS(:,4), ...
'tickrms',[0:.2:1],'titleRMS', 0 ,'showlabelsRMS',1,'widthRMS',0.5,'colRMS','g',...%去中⼼的均⽅根误差（就是减去平均值再求RMS）%每⼀个属性（如标题                'tickSTD',[0:.25:1.25],'limSTD',1.25,'styleSTD','-',... %标准差
'tickCOR',[.1:.2:.9 .95 .99],'showlabelsCOR',1,'titleCOR',1);%⽪尔逊相关系数
%添加图注
for ii = 1 : length(tt)
set(tt(ii),'fontsize',9,'fontweight','bold') %标志字母⼤⼩
set(pp(ii),'markersize',9) %标志点⼤⼩
if ii == 1
set(tt(ii),'String','Obser');%参考值图上显⽰名称
else
set(tt(ii),'String',alphab(ii-1));
end
end
title(sprintf('%s: Taylor Diagram at Hanjiang','B'),'fontweight','bold'); %图名及图编号
tt = axl(2).handle;
for ii = 1 : length(tt)
set(tt(ii),'fontsize',10,'fontweight','normal');
end
set(axl(1).handle,'fontweight','normal');
fontweight属性bold效果图：

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。