泰勒图-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 iserieMYSTATS(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 C D... 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 onsubplot(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）%每一个属性（如标题，标签）后数字，1存在，0不存在'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 == 1set(tt(ii),'String','Obser');%参考值图上显示名称elseset(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');

效果图：

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！