%1.xiaozao函数，是需要对标准化的序列进行消除数据噪音分析；
%2.Db3函数，是对数列进行Db3趋势分析；
%3.period函数，是求得时间序列的实部和模的平方。
%其中周期变化图是实部的等值线图
%而小波方差是模的平方的算数平均。
s=load('11.txt'); % input SST time series
s=zscore(s);
scales=[1:1:32];
%进行连续小波变换得到小波系数矩阵,选择复morlet小波函数
wf=cwt(s,scales,'cmor1-1.5','plot');
% 求得系数的实部
shibu=real(wf);
subplot(221);
contourf(shibu,10,'-');
colormap('HSV');
colorbar;
time=1960:5:2010
Xlabel('年份');
Ylabel('时间尺度/年');
set(gca,'XTickLabel', time) %更新XTickLabel
% 小波方差是模的平方的算数平均
mo=abs(wf);
mofang=mo.^2;
fangcha=mean(mofang,2);
subplot(222);
plot(fangcha,'k-','linewidth',1.5);
Xlabel('时间尺度/年');
Ylabel('小波方差');