小波分析MATLAB程序

%1.xiaozao函数，是需要对标准化的序列进行消除数据噪音分析；%2.Db3函数，是对数列进行Db3趋势分析；%3.period函数，是求得时间序列的实部和模的平方。%其中周期变化图是实部的等值线图%而小波方差是模的平方的算数平均。s=load('11.txt');%inputSSTtimeseriess=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:2010Xlabel('年份');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('小波方差');