AO和NAO指数

wlzhongouc

闲着没事，画画图，AO和NAO是两种不同的现象，但是它们又有着很大的关联，下面用Matlab处理了一下（数据来自：http://www.cpc.ncep.noaa.gov/products/precip/CWlink/daily_ao_index/ao.shtml#current）
程序如下：
%-------------AO and NAO Index data----------------
%----------------Created by 猛牛, Apr 22, 2012------------------------
clear;clf;clc;
load monthly.ao.index.b50.current.ascii
load norm.nao.monthly.b5001.current.ascii
AO=monthly_ao_index_b50_current;
NAO=norm_nao_monthly_b5001_current;
AO_p=AO(637:744,3);                   %-------Time: 2003 to 2011------
NAO_p=NAO(637:744,3);
A_index1=find(AO_p>0);
AO_positive=AO_p(AO_p>0);
A_index2=find(AO_p<0);
AO_negtive=AO_p(AO_p<0);
up=0;
for i=1:length(AO_p)
    up=up+(AO_p(i)-mean(AO_p))*(NAO_p(i)-mean(NAO_p));
end
down_1=0;down_2=0;
for i=1:length(AO_p)
    down_1=down_1+(AO_p(i)-mean(AO_p))^2;
    down_2=down_2+(NAO_p(i)-mean(NAO_p))^2;
end
r=up/(sqrt(down_1*down_2));           %-----Caculate the correlation coefficient-----
t=abs(r*sqrt(length(AO_p)-2)/sqrt(1-r^2));  %-----t-test------
bb=bar(A_index1,AO_positive,'b');
hold on
cc=bar(A_index2,AO_negtive,'r');
plot(NAO_p,'k','linewidth',2);
set(gca,'xtick',[1,13,25,37,49,61,73,85,97],'xticklabel',...
        {'2003-Jan','2004-Jan','2005-Jan','2006-Jan','2007-Jan','2008-Jan','2009-Jan','2010-Jan','2011-Jan'},'xlim',[1 110],'fontsize',16);
set(gca,'ylim',[-4.5 3],'fontsize',16);
xlabel('Year(Monthly mean)','FontSize',21);
ylabel('AO and NAO Index','FontSize',21);
set(gca,'XColor','k','YColor','k','linewidth',2);
text(3,-3,strcat('Correlation Coefficient R:',num2str(r)),'color','g','fontsize',16,'fontweight','bold');
legend('AO-Positive','AO-Negtive','NAO',3);
legend('boxoff');
title('AO and NAO Index','FontSize',23);
saveas(gcf,strcat('AO and NAO Index(Monthly mean)','.emf'));

成图如下：



可以看出AO和NAO的密切关系，相关系数达到0.68888，通过了显著水平0.05的检验。

PS：不是做气象的，欢迎讨论与批评指正 ：）

言深深

哈哈，不错嘛，matlab的绘图蛮好看的啊

majianyong

好羡慕会matlab的...

mofangbao

majianyong 发表于 2012-4-22 11:35

                               
登录/注册后可看大图

好羡慕会matlab的...

俺也还不会呢。。。

兰溪之水

言深深 发表于 2012-4-22 11:25

                               
登录/注册后可看大图

哈哈，不错嘛，matlab的绘图蛮好看的啊

兰溪哥哥教你~

言深深

兰溪之水 发表于 2012-4-22 12:32

                               
登录/注册后可看大图

兰溪哥哥教你~

切···

wode_q_x

感谢楼主分享自己的小发现！

wlzhongouc

wode_q_x 发表于 2012-6-13 10:10

                               
登录/注册后可看大图

感谢楼主分享自己的小发现！

共同努力 呵呵

butterflying

学习一下

/xin村儿/

时间序列正负值双色画法