墨西哥帽小波变换的影响域和计算方案新探讨
New Research on Cone of Influence and Computing Scheme of Mexican Hat Wavelet Transform
-
摘要: 该文分析了Torrence等人研究的墨西哥帽 (Mexican Hat) 小波变换的影响域和高频失真问题并探讨了解决方案。平移点b处、小波尺度为a的墨西哥帽小波函数的有效定义域为[ b-2.12a, b+2.12a], 以墨西哥帽为母函数的小波变换的真正影响域是2.12a, 小波尺度a的最大值应为N/4.24(N为时间序列长度)。提出了充分利用小波函数速降性质和立方样条插值小波系数的新计算方案, 消除了该计算方案产生虚假的显著高频周期振荡问题。利用正弦函数型时间序列小波变换的解析式检验了文中提出的新计算方案的合理性。利用新计算方案对冬季Niño3.4指数进行分析, 冬季Niño3.4指数的平均整体小波功率谱表明, 存在约12年的年代际变化和准4年的年际变化, 不存在显著的准两年周期振荡。Abstract: Using the effective domain of wavelet function, the problems of cone of influence (COI) and high-frequency distortion of Mexican Hat wavelet transform in the paper of Torrence and Compo (1998) are analyzed and the solution are explored. The analytic expressions of wavelet coefficient and global wavelet power spectrum for time series of sinusoid and the relationship between period and wavelet scale are deduced. The effective domain of Mexican Hat wavelet transform is [ b-2.12a, b+2.12a] for localized time b and wavelet scale a. The authentic cone of influence for Mexican Hat wavelet transform is 2.12a. The maximum value of wavelet scale is N/4.24 (N is the length of time series). When wavelet scale approaches zero, wavelet function is singular. So the simple sum scheme replacing the integration for wavelet coefficient does not include the singularity of wavelet function, therefore the pseudo high-frequency oscillation is deduced. A new computing scheme for wavelet coefficient is designed. The major elements include improved-resolution and located time b in the new grids. The difference of effective domain of wavelet function for small different wavelet scales can be distinguished by the improved-resolution. Located time b has great contribution to its wavelet coefficient because the wavelet function is the biggest at the located point b. The pseudo significant high-frequency oscillation that is produced by the computing scheme proposed by Torrence and Compo (1998) can be eliminated by the new computing scheme of wavelet coefficient, well using sharply descending characteristics of wavelet function and cubic spline interpolation. The rationality of new computing scheme is tested by the analytic solution of wavelet transform of sinusoid time sequence. The global wavelet power spectrum from the new computing scheme indicates that the variability of winter Niño3.4 index exhibits interdecadal oscillation of about 12 years and interannual oscillation of quasi-four years, and no significant quasi-biennial oscillation. The questions of cone of influence and pseudo high-frequency oscillation have two adverse aftermaths. First is that the correct judgment of interdecadal variation can be affected and second is that pseudo quasi-biennial variation can be deduced. Therefore, the results about COI and the new computing scheme for wavelet coefficient make the conclusions induced from wavelet analysis more credible.
-
Key words:
- wavelet transform;
- cone of influence;
- computing scheme;
- Niño3.4 index
-
图 2 墨西哥帽小波函数有效定义域及新计算方案网格示意图
(a) 小波尺度a=0.47(实线) 和a=0.16(虚线) 所对应的小波函数 (浅色水平线对应到这两种小波尺度上的有效定义域) 及3个相邻时间的距平序列 (柱状), (b) 小波尺度a=0.47所对应的小波函数 (实线) 及新的高分辨率网格 (竖虚线)
Fig. 2 Schematic diagram of effective domain of Mexican wavelet function and new computing scheme
(a) wavelet functions for wavelet scale a=0.47(solid line) and a=0.16(dashed line) (light horizontal lines stand for the effective domains for these wavelet scales) and anomalous series of 3 sequential time (bar), (b) wavelet function for wavelet scale a=0.47(solid line) and new high reso lution grid (vertical dashed line)
图 3 时间序列f(t)=sin (2πt/4) (t=0, …, 40) 的平均整体小波功率谱
(点线、断线和实线分别为式 (9)、文献[7]计算方案和本文新计算方案的结果)
Fig. 3 Mean global wavelet power spectrum of time series f(t)=sin (2πt/4)(t=0, …, 40)
(dotted line, dashed line and solid line stand for the results from analytic formula (9), computing scheme in reference [7] and new computing scheme in this paper)
图 4 1951—2007年冬季Niño3.4指数距平序列及其小波变换
(a)1951— 2007年冬季Niño3.4指数距平序列, (b) 本文新计算方案计算的小波系数和影响域, (c) 文献[7]方案计算的小波系数和影响域 (图中两边的斜交叉线区域为影响域, 阴影分别表示通过α=0.1的显著性水平检验 (Monte-Carlo方法))
Fig. 4 Time evolution of winter Niño3.4 index from 1951 to 2007 and wavelet transform
(a) anomalous series of winter Niño3.4 index from 1951 to 2007, (b) wavelet coefficient computed by new computing scheme and COI in this paper, (c) wavelet coefficient computed by computing scheme and COI in reference [7] (cross-hatched regions on either end indicate the COI; shaded areas represent regions pass the test of α=0.1 level (Monte-Carlo method))
-
[1] Farge M.Wavelet transforms and their applications to turbulence.Annu Rev Fluid Mech, 1992, 24:395-457. doi: 10.1146/annurev.fl.24.010192.002143 [2] Meyers S D, Kelly B G, Brien J J O.An introduction to wavelet analysis in oceanography and meteorology:With application to the dispersion of Yanai waves.Mon Wea Rev, 1993, 121:2858-2866. doi: 10.1175/1520-0493(1993)121<2858:AITWAI>2.0.CO;2 [3] Gamage N, Blumen W. Comparative analysis of lowlevel cold fronts:Wavelet, Fourier, and empirical orthogonal function decompositions.Mon Wea Rev, 1993, 121:2867-2878. doi: 10.1175/1520-0493(1993)121<2867:CAOLLC>2.0.CO;2 [4] Weng H, Lau K M.Wavelets, period doubling, and time-frequency localization with application to organization of convection over the tropical western Pacific. J Atmos Sci, 1994, 51:2523-2541. doi: 10.1175/1520-0469(1994)051<2523:WPDATL>2.0.CO;2 [5] Gu D, Philander S G H.Secular changes of annual and interannual variability in the Tropics during the past century.J Clim, 1995, 8:864-876. doi: 10.1175/1520-0442(1995)008<0864:SCOAAI>2.0.CO;2 [6] Baliunas S, Frick P, Sokoloff D, et al. Time scales and trends in the central England temperature data (1659-1990):A wavelet analysis.Geophys Res Lett, 1997, 24(11):1351-1354. doi: 10.1029/97GL01184 [7] Torrence C, Compo G P.A practical guide to wavelet analysis.Bull Amer Meter Soc, 1998, 79(1):61-78. doi: 10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2 [8] 邓自旺, 尤卫红, 林振山.子波变换在全球气候多时间尺度变化分析中的应用.南京气象学院学报, 1997, 20(4):505-510. http://www.cnki.com.cn/Article/CJFDTOTAL-NJQX704.012.htm [9] 张存杰, 董安祥, 白虎志, 等.甘肃省河东地区伏旱的小波分析.应用气象学报, 1998, 9(3):291-297. http://qk.cams.cma.gov.cn/jams/ch/reader/view_abstract.aspx?file_no=19980342&flag=1 [10] 尤卫红, 杞明辉, 段旭.小波变换在短期气候预测模型研究中的应用.高原气象, 1999, 18(1):39-46. http://www.cnki.com.cn/Article/CJFDTOTAL-GYQX901.004.htm [11] 纪忠萍, 谷德军, 谢炯光.广州近百年气候变化的多时间尺度分析.热带气象学报, 1999, 15(1):38-47. http://www.cnki.com.cn/Article/CJFDTOTAL-RDQX901.005.htm [12] 朱益民, 孙旭光, 陈晓颖.小波分析在长江中下游旱涝气候预测中的应用.解放军理工大学学报 (自然科学版), 2003, 4 (6):90-93. http://www.cnki.com.cn/Article/CJFDTOTAL-JFJL200306021.htm [13] 李贤彬, 丁晶, 李后强.基于子波变换序列的人工神经网络组合预测.水利学报, 1999, (2):1-5. http://www.cnki.com.cn/Article/CJFDTOTAL-SLXB902.000.htm [14] 戴新刚, 汪萍, 丑纪范.准地转正压大气小波谱模式及其数值解.自然科学进展, 2004, 14(9):1012-1019. http://www.cnki.com.cn/Article/CJFDTOTAL-ZKJZ200409008.htm [15] 赵洋, 肖华勇, 李振鹏, 等.一种基于小波分析理论的灰色预测方法.西南民族大学学报 (自然科学版), 2005, 31(4):498-501. http://www.cnki.com.cn/Article/CJFDTOTAL-XNMZ200504004.htm [16] 全利红, 胡非, 程雪玲.用小波系数谱方法分析湍流湿度脉动的相干结构.大气科学, 2007, 31(1):57-63. http://www.cnki.com.cn/Article/CJFDTOTAL-DQXK200701005.htm [17] 李世雄, 刘家琦.小波变换和反演数学基础.北京:地质出版社, 1994:10-14. [18] 郑彬, 梁建茵.对流层准两年周期振荡的研究进展.热带气象学报, 2005, 21(1):79-86. http://www.cnki.com.cn/Article/CJFDTOTAL-RDQX200501009.htm