前些日子科研上突然遇到了需要对序列进行频谱分析的问题，查阅了大量的python资料，虽然诸多科学计算库中也提供了多种频谱分析的方法和函数，但是计算结果总是不太对劲，没有算法是针对气象气候序列的分析，不得已又掏出老本行去NCL上运行，搞来搞去十分麻烦，于是一咬牙一狠心，没有算法自己开发好了，翻阅了魏凤英老师的统计方法一书，结合气象家园现有的fortran和matlab代码，终于是成功将适用于气候序列频谱分析的功率谱和小波代码移植了出来。
具体还需要再完善，目前只是适配了测试数据，还需要再进行封装和注释。
上一下结果图进行对比吧。
功率谱
该部分参考气象家园原贴
原贴作者Aries结果图
移植python后出图：
python计算结果
小波变换
原代码参考见：点击这里
(Written January 1998 by C. Torrence)
气象家园也提供了代码和测试数据下载
http://bbs.06climate.com/forum.php?mod=viewthread&tid=31901&highlight=%D0%A1%B2%A8
matlab出图效果：
matlab
python移植结果：
python
逐步debug调整算法，终于得到了正确的结果吗，抑制不住自己激动的内心，至于tick那个密集的东西，现在实在是没有心情去调整了，先发文记录，封装好后应该会上传至github或者我自己的网站。
网友评论