绘制SOI曲线

1、数据下载
http://www.cpc.ncep.noaa.gov/data/indices/
2、转成Excel格式
这是20年的SOI数据：

2、MATLAB制图
%读入SOI.xls，将12个月的数据（211）合并 在SOI_DATA
%为2112----A，进行reshape(A,[],1)变成252*1
%Y为SOI标准化值
%绘图
data=xlsread(‘G:\Experiment\oceanRS\SOI\StandardData_SOI.xlsx’,‘Sheet1’,‘A2:M22’);
X=data(1:21,2:13);
XX=X’;
Y=reshape(XX,1,252);
figure;
tt=linspace(1997,2017,252);
x=[1997:1:2017];
y1=Y;
y2=Y;
y1(Y<0)=0;
y2(Y>0)=0;
area(tt,y1,‘FaceColor’,‘b’);
hold on;
area(tt,y2,‘FaceColor’,‘r’);
grid on;
axis([1996,2017,-3,3]);
set(gca,‘xtick’,x)
x=get(gca,‘xlim’);
xlabel(‘YEAR’) ;
title(‘SOI-Standard’);
line([1996,2017],[0,0],‘linestyle’,’–’,‘color’,‘k’);
3、结果图