爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
楼主: houtian

[程序设计] matlab 可以绘制泰勒图么?

[复制链接]
发表于 2016-11-10 17:21:18 | 显示全部楼层
下载下来学习一下
密码修改失败请联系微信:mofangbao
发表于 2017-8-5 16:36:58 | 显示全部楼层
密码修改失败请联系微信:mofangbao
发表于 2017-8-5 16:38:13 | 显示全部楼层
% This is a test to use the taylordiag.m ploting function.
%
% The data file taylordiag_egdata.mat is required together
% with the function "allstats" and "ptable".
% Both are available at: http://code.google.com/p/guillaumemaze/
%
%
% This function runs the following command lines:
%
% clear
% load taylordiag_egdata.mat
%
% % Get statistics from time series:
% for ii = 2:size(BUOY,1)
%     C = allstats(BUOY(1,:),BUOY(ii,:));
%     statm(ii,:) = C(:,2);
% end
% statm(1,:) = C(:,1);
%
% % Plot:
% figure
% ax = ptable([2 3],[2 2;4 6]);
% iw=2;jw=3;
% alphab = 'ABCDEFG';
%
% subplot(iw,jw,2);
% plot(BUOY');
% grid on,xlabel('time (day)');ylabel('heat fluxes (W/m^2)');
% title(sprintf('%s: These are the different time series of daily heat fluxes (W/m2)','A'),'fontweight','bold');
%
% subplot(iw,jw,5); hold on
% [pp tt axl] = taylordiag(squeeze(statm(:,2)),squeeze(statm(:,3)),squeeze(statm(:,4)),...
%             'tickRMS',[25:25:150],'titleRMS',0,'tickRMSangle',135,'showlabelsRMS',0,'widthRMS',1,...
%             'tickSTD',[25:25:250],'limSTD',250,...
%             'tickCOR',[.1:.1:.9 .95 .99],'showlabelsCOR',1,'titleCOR',1);
%
% for ii = 1 : length(tt)
%     set(tt(ii),'fontsize',9,'fontweight','bold')
%     set(pp(ii),'markersize',12)
%     if ii == 1
%         set(tt(ii),'String','Buoy');
%     else
%         set(tt(ii),'String',alphab(ii-1));
%     end
% end
% title(sprintf('%s: Taylor Diagram at CLIMODE Buoy','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');


clear
load taylordiag_egdata.mat

% Get statistics from time series:
for ii = 2:size(BUOY,1)
    C = allstats(BUOY(1,:),BUOY(ii,:));
    statm(ii,:) = C(:,2);
end
statm(1,:) = C(:,1);

% Plot:
figure
ax = ptable([2 3],[2 2;4 6]);
iw=2;jw=3;
alphab = 'ABCDEFG';

subplot(iw,jw,2);
plot(BUOY');
grid on,xlabel('time (day)');ylabel('heat fluxes (W/m^2)');
title(sprintf('%s: These are the different time series of daily heat fluxes (W/m2)','A'),'fontweight','bold');

subplot(iw,jw,5); hold on
[pp tt axl] = taylordiag(squeeze(statm(:,2)),squeeze(statm(:,3)),squeeze(statm(:,4)),...
            'tickRMS',[25:25:150],'titleRMS',0,'tickRMSangle',135,'showlabelsRMS',0,'widthRMS',1,...
            'tickSTD',[25:25:250],'limSTD',250,...
            'tickCOR',[.1:.1:.9 .95 .99],'showlabelsCOR',1,'titleCOR',1);

for ii = 1 : length(tt)
    set(tt(ii),'fontsize',9,'fontweight','bold')
    set(pp(ii),'markersize',12)
    if ii == 1
        set(tt(ii),'String','Buoy');
    else
        set(tt(ii),'String',alphab(ii-1));
    end
end
title(sprintf('%s: Taylor Diagram at CLIMODE Buoy','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');
密码修改失败请联系微信:mofangbao
发表于 2017-11-20 18:30:01 | 显示全部楼层
同问,不知道怎么画
密码修改失败请联系微信:mofangbao
发表于 2017-11-20 18:30:23 | 显示全部楼层
好像NCL是可以画的
密码修改失败请联系微信:mofangbao
发表于 2017-11-20 18:30:50 | 显示全部楼层
onetwothree 发表于 2015-11-2 22:56
论坛里面有代码~~话说option value值怎么设置呢

不知道诶。。。。
密码修改失败请联系微信:mofangbao
发表于 2017-11-20 18:31:06 | 显示全部楼层
曲梓祎 发表于 2016-4-23 22:46
http://bbs.06climate.com/forum.php?mod=viewthread&tid=35864
这里有人分享啊。。。不过没下下来呢。。

我也没下下来。。。。
密码修改失败请联系微信:mofangbao
发表于 2018-3-31 10:36:23 | 显示全部楼层
{:5_275:}{:5_275:}{:5_275:}{:5_275:}
密码修改失败请联系微信:mofangbao
发表于 2018-4-17 15:10:55 | 显示全部楼层
泰勒图是什么意思啊
密码修改失败请联系微信:mofangbao
发表于 2018-7-2 21:02:32 | 显示全部楼层
grace639 发表于 2015-7-11 20:37
同画泰勒图 握个爪

这个可以直接用吗
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表