%1.xiaozao函數(shù)，是需要對(duì)標(biāo)準(zhǔn)化的序列進(jìn)行消除數(shù)據(jù)噪音分析；

%2.Db3函數(shù)，是對(duì)數(shù)列進(jìn)行Db3趨勢分析；

%3.period函數(shù)，是求得時(shí)間序列的實(shí)部和模的平方。

%其中周期變化圖是實(shí)部的等值線圖

%而小波方差是模的平方的算數(shù)平均。

s=load('11.txt'); % input SST time series

s=zscore(s);

scales=[1:1:32];

%進(jìn)行連續(xù)小波變換得到小波系數(shù)矩陣,選擇復(fù)morlet小波函數(shù)

wf=cwt(s,scales,'cmor1-1.5','plot');

% 求得系數(shù)的實(shí)部

shibu=real(wf);

subplot(221);

contourf(shibu,10,'-');

colormap('HSV');

colorbar;

time=1960:5:2010

Xlabel('年份');

Ylabel('時(shí)間尺度/年');

set(gca,'XTickLabel', time) %更新XTickLabel

% 小波方差是模的平方的算數(shù)平均

mo=abs(wf);

mofang=mo.^2;

fangcha=mean(mofang,2);

subplot(222);

plot(fangcha,'k-','linewidth',1.5);

Xlabel('時(shí)間尺度/年');

Ylabel('小波方差');

總結(jié)

以上是生活随笔為你收集整理的降水小波分析matlab程序,小波分析MATLAB程序的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。