关于小波分析中的lag1

雪颜

lag1 = 0.72;  % 红噪声背景谱的自相关检验参数

程序默认是0.72，可是我看torrence的文献，里面不同时间序列的lag1是不同的。

我觉得直接用0.72很不妥。

雪颜

如图，都是95%，但是lag1的值是不同的（也就是α）

The 95% confidence level was derived using (31) and assuming red-noise spectra
(16) with α = 0.72 for Nino3 SST and α = 0.61 for the SOI.

原文

雪颜

原文：
where α is the assumed lag-1 autocorrelation, x0 = 0, and zn is taken from Gaussian white noise.
Thus, by choosing an appropriate lag-1 autocorrelation, one can use (16) to model a red-noise spectrum.
This red-noise was estimated from (α1+ √α2)/2, where α1 and α2 are the lag-1 and lag-2 autocorrelations of the Nino3 SST.

也就是说，所有用lag-1=0.72的小波分析都是错的，因为0.72只是对于nino3 sst的红噪声自相关系数，是要自己算的。

雪颜

我把论坛上的小波程序看完了，和原程序进行比对，我发现很多都是瞎改的，画出来的图也是错的。
大家还是好好学习小波，多看点文献，不要只做伸手党，画完图就算了，错了没都不知道。

gegbt_Sky

非常认可楼主的观点，我也有相同的思考。那么请问楼主，红噪声检验背景谱的自相关检验参数
lag1是不就是待处理数据的一阶自相关系数？如果是的话，在修改Torrence程序的时候，lag1改成自己数据的一阶自相关系数就可以了吧？

雪颜

gegbt_Sky 发表于 2017-2-14 00:04
非常认可楼主的观点，我也有相同的思考。那么请问楼主，红噪声检验背景谱的自相关检验参数
lag1是不就是待 ...

是的，lag1必须改。

吧啦吧啦

The red-noise background spectrum for α = 0.72 is shown by the lower dashed curve in Fig. 3. This red-noise was estimated from ( α1+ √α2)/2, where α1 and α2 are the lag-1 and lag-2 autocorrelations of the Niño3 SST.
我很好奇为什么α 这么取，不是只是取lag-1就好了吗？

看云淡风清

非常认同楼主说的，我用了几个小波程序，结果都不一样。还有那篇文章很经典，只是我还看不懂额。。。。

雪颜

吧啦吧啦 发表于 2017-3-22 20:52
The red-noise background spectrum for α = 0.72 is shown by the lower dashed curve in Fig. 3. This r ...

α 就是lag-1

想飞的三叶草

我也发现了这个问题，求问楼主，lag-1公式是什么？怎么算