時間序列分析是地學研究中能夠經(jīng)常遇到的問題张咳。在時間序列的研究中鼎俘，時域和頻域是常用的兩種基本形式。
- 時域具有時間定位功能纲酗，但無法得到關(guān)于時間序列變化的更多信息材鹦。
- 頻域（Fourier）分析雖然具有準確的頻率定位功能逝淹，但僅僅適合平穩(wěn)時間序列分析。
但是地學的許多特征（如徑流桶唐、降雨栅葡、地震波等）隨時間的辯護往往受到多種因素的綜合影響，大都屬于非平穩(wěn)時間序列尤泽。他們不但具有趨勢性欣簇、周期性等特征，還存在著隨機性坯约、突變性以及“多時間尺度”結(jié)構(gòu)熊咽，具有多層次演變的規(guī)律。
小波分析步驟：
- 小波函數(shù)的選取
- 數(shù)據(jù)格式轉(zhuǎn)換
- 邊界效應的消除或者減小
- 計算小波系數(shù)
- 計算復小波的實部闹丐、模横殴、模方、方差
- 繪制小波系數(shù)實部卿拴、模衫仑、模方等值線圖
- 繪制小波方差圖
- 繪制主周期趨勢圖
實例應用：
- import excel -> file->import data-> save as .mat
- wavemenu
- signal extesion -> direction to extend-> both || extension mode -> symmetic(half-point) 均衡-> save transformed signal -> e_data.mat
- (core) One-dimensional Complex continuous wavelet analysis(CWT)-> Load a signal-> choose the wavelet -> choose the sampling period-> choose the scale settings-> Analyze->save cofficients-> coff.mat
- load the coff.mat to workplace
% 求實部值
shibu = real(coefs)
% 求模、模方
mo = abs(coefs)
mofang = (mo).^2
% 求方差
fangcha = sum(abs(coefs).^2,2)
----------繪制實部等值線圖----------------
contourf(shibu, 10, '-') % 實部等值線圖
colorbar;
time = starttime: step: endtime
xlabel('Year');
ylabel('Time scales');
set(gca, 'xticklabel', time); % 確定橫坐標的范圍
print(gcf, '-r300', '-dpng', 'name.png');
小波分析matlab 操作總結(jié)：
- 導入數(shù)據(jù)
- 保存 * .mat file
- 保存 *.png file
- matlab 輸出圖片
- print(gcf, '-r300', '-dpng', 'name.png'); 這種方法最好用堕花。圖片質(zhì)量高文狱。
- save as 相當于截圖 圖片質(zhì)量不高
