前些日子科研上突然遇到了需要對(duì)序列進(jìn)行頻譜分析的問(wèn)題，查閱了大量的python資料，雖然諸多科學(xué)計(jì)算庫(kù)中也提供了多種頻譜分析的方法和函數(shù)，但是計(jì)算結(jié)果總是不太對(duì)勁，沒(méi)有算法是針對(duì)氣象氣候序列的分析，不得已又掏出老本行去NCL上運(yùn)行，搞來(lái)搞去十分麻煩，于是一咬牙一狠心，沒(méi)有算法自己開(kāi)發(fā)好了，翻閱了魏鳳英老師的統(tǒng)計(jì)方法一書(shū)，結(jié)合氣象家園現(xiàn)有的fortran和matlab代碼，終于是成功將適用于氣候序列頻譜分析的功率譜和小波代碼移植了出來(lái)。
具體還需要再完善，目前只是適配了測(cè)試數(shù)據(jù)，還需要再進(jìn)行封裝和注釋。
上一下結(jié)果圖進(jìn)行對(duì)比吧。
功率譜
該部分參考氣象家園原貼
移植python后出圖：
小波變換
原代碼參考見(jiàn)：點(diǎn)擊這里
(Written January 1998 by C. Torrence)
氣象家園也提供了代碼和測(cè)試數(shù)據(jù)下載
http://bbs.06climate.com/forum.php?mod=viewthread&tid=31901&highlight=%D0%A1%B2%A8
matlab出圖效果：
python移植結(jié)果：
逐步debug調(diào)整算法，終于得到了正確的結(jié)果嗎，抑制不住自己激動(dòng)的內(nèi)心，至于tick那個(gè)密集的東西，現(xiàn)在實(shí)在是沒(méi)有心情去調(diào)整了，先發(fā)文記錄，封裝好后應(yīng)該會(huì)上傳至github或者我自己的網(wǎng)站。
