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

總結

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

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。