降水小波分析matlab程序,小波分析MATLAB程序

%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('小波方差');