登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
本帖最后由 lotina 于 2012-9-26 12:35 编辑
原来我在气象论坛里有一个求助小波分析的周期设置的问题，但貌似还是有很多人在帖子上问这个问题，可能是那个帖子有点乱，没有讲明白。另外就是在论坛里已经有很多小波分析的程序了，还是有很多人在问这个问题啊o(︶︿︶)o 。 现在，我把自己做小波分析时学到的一点东西写出来，希望对大家有点帮助吧，也能少点这样的提问了。 PS：我对小波的原理还是不懂，这只是我使用的过程中的一些心得。若有误，请各位大侠们指正。
1. 对程序里的参数的一点说明：
这里的NC=81暂且称作层次吧。其中NX是你要做小波分析的时间序列。而这里的2，4，8，16，32是画图时需要显示到图上的时间周期，这个周期的上限是根据时间序列长度来取的，一般取为NX的一般，而时间周期和程序里的层次NC的对应关系如上，以此类推：101-64，121-128，141-256。一般nx=365时，取nc的上限为141.此时对应的是一年的逐日（365天）资料的小波分析。
2. 在三个CTL文件中要修改相应的层次和时间。 在GS文件中，除了修改时间和层次外，还需要修改在2，4，8，16，32……的周期的位置，与NC层次对应。这个需要先画出一次图形后，根据1，21，41，81……的位置来修改上面需要显示的周期的正确位置。
3．举个列子：63年的逐年资料，则层次就对应为81，即nx=62,nc=81，运行FOR程序，然后运行GS文件，得到图形为：
根据周期与层次的对应关系，修改gs中的命令，与nc的刻度对齐或者将其覆盖，需要修改的地方：
结果为：
如果不想显示NC的层次（因为在分析中只研究周期），去掉GS中语句*'set ylab off'中的*就好了。 修改语句直到画出你满意的图形。 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~OVER~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 另外是我空间里的一篇有关小波分析前的资料预处理的日志，对做低频小波分析的童鞋可能会有点用吧：
MORLET小波分析的FORTRAN程序和GRADS画图程序：
| 
[复制链接]
窥视卡
雷达卡
发表于 2012-5-7 14:01:26
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
楼主
发表于 2012-5-7 14:15:23
发表于 2012-5-7 17:27:57
谢谢姐姐为我们小白做出的付出 特来撒花学习