eeof程序输出PC1和PC2之后如何画图

夏茗悠

我用EEOF程序分析了1981五月到八月逐日的olr，u850和u200，输出了PC1和PC2，但是不知道要怎样画出下面这样的RMM指数的位相

沙颖凯

横坐标RMM1 纵坐标RMM2，圆圈是RMM1^2+RMM2^2=1。不过这种东西在GrADS里好像不好实现

夏茗悠

沙颖凯 发表于 2014-4-16 23:00
横坐标RMM1 纵坐标RMM2，圆圈是RMM1^2+RMM2^2=1。不过这种东西在GrADS里好像不好实现

是啊，所以求助各位大神，这种图要用什么软件来画呢？

cold_wq

Matlab可以实现，比较简单

turtleweigg

clear all
clc
fclose all;

figure(1)
alpha=0:pi/50:2*pi;%角度[0,2*pi]
R=1;%半径
x=R*cos(alpha);
y=R*sin(alpha);
plot(x,y,'k-','linewidth',2);
hold on
plot([-4,-sqrt(2)/2],[-4,-sqrt(2)/2],'k','linewidth',2,'linestyle','--');
plot([sqrt(2)/2,4],[sqrt(2)/2,4],'k','linewidth',2,'linestyle','--');
plot([-4,-sqrt(2)/2],[4,sqrt(2)/2],'k','linewidth',2,'linestyle','--');
plot([sqrt(2)/2,4],[-sqrt(2)/2,-4],'k','linewidth',2,'linestyle','--');
plot([-4,-1],[0,0],'k','linewidth',2,'linestyle','--');
plot([1,4],[0,0],'k','linewidth',2,'linestyle','--');
plot([0,0],[1,4],'k','linewidth',2,'linestyle','--');
plot([0,0],[-1,-4],'k','linewidth',2,'linestyle','--');
set(gcf,'Position',[100 100 700 700]);
set(gca,'Position',[0.1 0.1 0.8 0.8]);
set(gca,'xtick',-4:1:4);
set(gca,'ytick',-4:1:4);
set(gca,'xticklabel',-4:1:4,'fontsize',18);
set(gca,'yticklabel',-4:1:4,'fontsize',18);

text(-0.8,0,'Weak MJO','fontsize',18,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(-0.8,-3.5,'Indian Ocean','fontsize',15,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(3.2,-0.5,'Maritime','fontsize',15,'BackgroundColor',[1,1,1],'rotation',90,'fontname','Cambria');
text(3.5,-0.55,'Continent','fontsize',15,'BackgroundColor',[1,1,1],'rotation',90,'fontname','Cambria');
text(-0.8,3.5,'Western Pacific','fontsize',15,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(-3.5,-0.8,'Western Hemi','fontsize',15,'BackgroundColor',[1,1,1],'rotation',90,'fontname','Cambria');
text(-3.15,-0.5,'and Africa','fontsize',15,'BackgroundColor',[1,1,1],'rotation',90,'fontname','Cambria');

text(-3.5,-1.5,'1','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(-1.5,-3.5,'2','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(1.5,-3.5,'3','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(3.5,-1.5,'4','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(3.5,1.5,'5','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(1.5,3.5,'6','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(-1.5,3.5,'7','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');
text(-3.5,1.5,'8','fontsize',20,'BackgroundColor',[1,1,1],'fontname','Cambria');

xlabel('RMM1','fontsize',18);
ylabel('RMM2','fontsize',18);

turtleweigg

接下来打点连线这种事情会matlab的一般都会，就按照个人习惯啦~

夏茗悠

turtleweigg 发表于 2014-4-17 14:44
接下来打点连线这种事情会matlab的一般都会，就按照个人习惯啦~

恩，原来是用matlab啊，因为我之前习惯用grads，所以看到这种图就不知道怎么画了，谢谢你啦，我会去试试的

木市街8号

楼主能不能把你的EEOF程序给我分享一下，最近也急着用，麻烦了

夏茗悠

木市街8号 发表于 2014-4-22 09:32
楼主能不能把你的EEOF程序给我分享一下，最近也急着用，麻烦了

就是从这里下载的，你点进去就能看见了http://cawcr.gov.au/staff/mwheeler/maproom/RMM/

木市街8号

好的，谢谢了