




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、小波分析张晓楠 始于2015-6-25目录小波分析11. 重要学习网址:12. 小波分析的要点12.1 目的12.2 方法1(1)连续小波变换CWT2(2)小波功率谱3(3)小波功率谱边缘效应及影响锥3(4)小波功率谱的显著性检验4(5)尺度选择42.3 具体步骤43. Matlab: Crosswavelet and Wavelet Coherence54. 文献举例111. 重要学习网址:小波分析:/research/wavelets/交叉小波:/wavelet-coherence/wtc-exa
2、mplehttp:/noc.ac.uk/using-science/crosswavelet-wavelet-coherence/home.php/home.php?mod=space&uid=1148346&do=blog&id=882008/home.php/home.php?mod=space&uid=1148346&do=blog&id=8820082. 小波分析的要点2.1 目的小波分析是一个强有力的统计工具,最早使用在信号处理与分析领域中,通过对声音、图像、地震等信号进行降噪、重建、提取
3、,从而确定不同信号的震动周期出现在哪个时间或频域上。现在广泛的应用于很多领域。在地学中,各种气象因子、水文过程、以及生态系统与大气之间的物质交换过程都可以看作是随时间有周期性变化的信号,因此小波分析方法同样适用于地学领域,从而对各种地学过程复杂的时间格局进行分析。如,温度的日变化周期、年变化周期出现在哪些事件段上,在近100年中,厄尔尼诺-拉尼娜现象的变化周期及其出现的时间段,等等。2.2 方法小波变换具有多分辨率分析的特点,并且在时频两域都具有表征信号局部特征的能力。小波变换通过将时间系列分解到时间频率域内,从而得出时间系列的显著的波动模式,即周期变化动态,以及周期变化动态的时间格局(Tor
4、rence and Compo, 1998)。小波(Wavelet),即小区域的波,是一种特殊的、长度有限,平均值为零的波形。它有两个特点:一是“小”,二是具有正负交替的“波动性”,即直流分量为零。小波分析是时间(空间)频率的局部化分析,它通过伸缩平移运算对信号(函数)逐步进行多尺度细化,能自动适应时频信号分析的要求,可聚焦到信号的任意细节。小波分析将信号分解成一系列小波函数的叠加,而这些小波函数都是由一个母小波(mother wavelet)函数经过平移与尺度伸缩得来的。用这种不规则的小波函数可以逼近那些非稳态信号中尖锐变化的部分,也可以去逼近离散不连续具有局部特性的信号,从而更为真实的反映
5、原信号在某一时间尺度上的变化。小波分析这种局部分析的特性使其成为对非稳态、不连续时间序列进行量化的一个有效工具(Stoy et al., 2005)。小波是一个具有零均值且可以在频率域与时间域内进行局部化的数学函数(Grinsted et al., 2004)。一个小波被称为母小波(mother wavelet),母小波可沿着时间指数经过平移与尺度伸缩得到一系列子小波。子小波可以通过尺度(s,频率的反函数)函数和时间(n)位置或平移来描述。利用一系列子小波,一个信号可以在不同的时间尺度上进行计算并显示出详细的特征尺度。拉伸更大的小波窗口,使其宽度更大便可以分析时间系列中波动较大的部分并捕捉大尺
6、度(低频)事件的特征。相反,压缩较小的窗口将包含小尺度(高频)的事件信息。当信号被子小波相乘,被s与n唯一的表达,我们可以计算出信号在时间频率域一个具体位置的系数。如果信号在时间n上的谱成分可以与小波s比较,那么计算的小波系数具有相对较大的值。在其它n与s的组合(如其它的子小波)上都进行这样的计算,那么将会产生一系列系数(小波变化)来表达信号在时间频率域内的分解。通过这样的变化便可得到时间系列的波动模式(周期变化模式)以及这些模式随时间的变化(Furon et al., 2008; Jevrejeva et al., 2003)。小波变化可以分为连续小波变化(the Continuous Wa
7、velet Transform, CWT)与离散小波变换(Discrete Wavelet Transform, DWT)。离散小波变化DWT是数据的紧凑表示,长用于降噪与数据压缩。连续小波变化CWT更适合于信号特征的提取(Grinsted et al., 2004)。CWT作为时间系列间歇式波动特征提取的工具被广泛的应用的地球物理学研究中(Grinsted et al., 2004; Furon et al., 2008)。(1)连续小波变换CWT可以将具有等时间步长t的离散时间系列xn (n=1, N)的连续小波变换定义为小波函数0尺度化以及转换下的xn的卷积: (1)式中*表示共轭复数,
8、N是时间系列的总数据个数,(t/s)1/2是一个用于小波函数标准化的因子从而使得小波函数在每个小波尺度s上具有单位能量。通过转换小波尺度s并沿着时间指数n进行局部化,最终可得到一幅展示时间系列在某一尺度上波动特征及其随时间变化的图谱,即小波功率谱(Torrence and Compo, 1998; Torrence and Webster, 1999; Grinsted et al., 2004)。对一个时间系列进行小波转换时,母小波的选择显得尤为重要,Farge(1992)曾经讨论过母小波选择时需要考虑的因素,例如正交与非正交、负值与实值、母小波的宽度与图形等等。正交小波函数一般用于离散小波
9、变换,非正交小波函数即可用于离散小波变换也可用于连续小波变换(Torrence and Compo, 1998)。通常在对时间系列进行分析时,希望能够得到平滑连续的小波振幅,因此非正交小波函数较为合适。此外,要得到时间系列振幅和相位两方面的信息,就要选择复值小波,因为复值小波具有虚部,可以对相位进行很好的表达(Torrence and Compo, 1998)。Morlet小波不但具有非正交性而且还是由Gaussian调节的指数复值小波。 (2)式中t为时间,0是无量纲频率。当0=6,小波尺度s与傅里叶周期(period)基本相等(, = 1.03s)(Torrence and Webster
10、, 1999),所以尺度项与周期项可以相互替代。由此可见,Morlet小波在时间与频率的局部化之间有着很好的平衡(Grinsted et al., 2004)。此外,Morlet小波中还包含着更多的振动信息,小波功率可以将正、负峰值包含在一个宽峰之中(Torrence and Compo, 1998)。(2)小波功率谱为使计算更为快捷,公式5-1的卷积在傅里叶域内执行(Torrence and Compo, 1998; Grinsted et al., 2004)。定义为小波功率谱(wavelet power spectrum),该功率谱表达了时间系列在给定小波尺度和时间域内的波动量级(Laf
11、renire and Sharp, 2003)。由于我们采用的Morlet母小波为复值小波,因此也为复数,其复值部分可以解释为局部相位(Torrence and Compo, 1998)。将小波功率谱在某一周期上进行时间平均,我们可以得到小波全谱(global wavelet spectrum), (3)小波全谱能够表明时间系列真实功率谱的无偏、一致估计(Torrence and Compo, 1998)。由于小波全谱可以显示出背景谱量度,所以局部小波谱的峰值可以得到验证。因为该特性,通过小波全谱图中可以清晰的辨别时间系列的周期波动特征及其强度。(3)小波功率谱边缘效应及影响锥由于小波变换假设
12、数据是循环的,所以当我们处理有限长度的时间系列时,在小波功率谱中会出现边缘效应,即在功率谱的起始及末端部分出现误差。由于该原因,需要我们在时间系列的末尾补零从而使得分析的时间系列的总长度N大于2m而小于2m+1。但是,当我们采取这样的措施时会在小波功率图谱边缘引起端点不连续以及谱振幅下降的现象。在这种情况下,需要明确一个概念,即影响锥(Cone of Influence, COI),影响锥COI表示小波谱区域以及相应的边缘效应。在COI的边缘小波谱值会下降e-2(Torrence and Compo, 1998; Grinsted et al., 2004; Furon et al., 200
13、8)。(4)小波功率谱的显著性检验小波功率谱的统计显著性可以对照一个原假设进行评价,该原假设为假设信号由一个给定背景功率谱(Pk)的稳定过程产生,通常背景功率谱为白噪声或红噪声(Torrence and Compo, 1998; Lafrenire and Sharp, 2003)。由于许多地球物理时间系列具有红噪声特征(即方差随着尺度的增加或频率的下降而增加),所以常采用红噪声作为背景谱对小波谱进行检验。红噪声过程可以很好的由一阶自回归过程(AR1)来模拟(Torrence and Compo, 1998; Grinsted et al., 2004)。一个由lag-1自相关处理的AR1的傅
14、里叶功率谱可以定义为: (4)式中k为傅里叶频率指数。通常在研究中,每个尺度上用COI以外的值以5%的显著水平进行估计。(5)尺度选择在进行小波变换时,还需要选择一系列尺度s。本研究使用非正交小波变换,我们可以使用任意一组的尺度来构建较完整的图像。一系列尺度可以用2的分数幂来表达:,j = 0, 1, J (5) (6)式中,s0为可分辨的最小尺度,J为确定的最大尺度。s0应该被选择恰当以便使相等的傅里叶周期近似于2t。 一个足够小的j的选择依赖于小波方程谱空间的宽度。2.3 具体步骤(1) 数据预处理数据时间系列必须是连续等时间步长。进行标准化处理(2) 母小波选择可选择Mexican ha
15、t 小波,或Morlet小波。通常在对时间系列进行分析时,希望能够得到平滑连续的小波振幅,因此非正交小波函数较为合适。此外,要得到时间系列振幅和相位两方面的信息,就要选择复值小波,因为复值小波具有虚部,可以对相位进行很好的表达(Torrence and Compo, 1998)。Morlet小波不但具有非正交性而且还是由Gaussian调节的指数复值小波。(3) 尺度选择如时间序列为47年的年降水数据,时间系列长度N=47,为了减小功率谱的边缘效应,在进行交互小波变换时选择26个数据。时间步长dt=1,即一年一个数据。j可选择0.125。(4) 显著性检验由于许多地球物理时间系列具有红噪声特征
16、(即方差随着尺度的增加或频率的下降而增加),所以常采用红噪声作为背景谱对小波谱进行检验(在程序中lag1=0.72)。在计算中,每个尺度上用COI以外的值以5%的显著水平进行估计。参考文献(1) Farge M. Wavelet transforms and their applications to turbulence. Annual Review of Fluid Mechanics. 1992. 24: 395-457.(2) Furon A, Wagner Riddle C, Smith C R, et al. Wavelet analysis of wintertime and s
17、pring thaw CO2 and N2O fluxes from agricultural fields. Agricultural and Forest Meteorology. 2008. 148, 1305-1317.(3) Grinsted A, Jevrejeva S, Moore J. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Processes in Geophysics. 2004. 11: 561-566.(4
18、) Jevrejeva S, Moore J C, Grinsted A. Influence of the Arctic oscillation and El Nino -southern oscillation (ENSO) on ice conditions in the Baltic Sea: the wavelet approach. Journal of Geophysical Research. 2003. 108 (D21), 4677.(5) Lafrenire M, Sharp M. Wavelet analysis of inter-annual variability
19、in the runoff regimes of glacial and nival stream catchments, Bow Lake, Alberta. Hydrological Process. 2003. 17, 1093-1118.(6) Stoy P C, Katual G G, Siqueira M B S, et al. Variablity in net ecosystem exchange from hourly to inter-annual time scale at adjacent pine and hardwood forests: a wavelet ana
20、lysis. Tree Physiology. 2005. 25: 887-902.(7) Torrence C, Compo G P. A practical guide to wavelet analysis. Bulletin of the American Meteorological Society. 1998. 79(1): 61-78.(8) Torrence C, Webster P J. Interdecadal Changes in the ENSO-Monsoon system. 1999. Journal of climate. 1999. 12: 2679-2690.
21、3. Matlab: Crosswavelet and Wavelet CoherenceCross wavelet analysis and wavelet coherence are powerful methods for testing proposed linkages between two time series.Check the histograms of the time series to ensure that they are not too far from normally distributed. Consider transforming the time s
22、eries, if the pdfs of the time series are far from Gaussian. When choosing a transformation, it is preferable to choose an analytic transformation such as taking the logarithm if the data is log-normally distributed. In other cases a simple percentile transformation might be useful. An advantage of
23、using that particular transformation is that it does not have any outliers.Consider what the expectations are for the outcome of the analysis given the proposed linking mechanism. We caution against blindly applying these methods to randomly chosen data sets. Like other statistical tests some datase
24、t sets will display highly statistically significant links simply by chance.When a wavelet has been chosen the CWTs of both time series are calculated. We suggest a scale resolution of 10 scales peroctaveand use of the Morlet wavelet unless there are good grounds to do otherwise. For geophysical tim
25、e series an AR1 red noise assumption is often suitable and can be used to calculate the significance level of the wavelet power. Remember to take special care not to misinterpret results inside the COI.From the two CWTs the XWT is calculated. The XWT exposes regions with high common power and furthe
26、r reveals information about the phase relationship. If the two series are physically related we would expect a consistent or slowly varying phase lag that can be tested against mechanistic models of the physical process. The circular mean of the phase angles can be used to quantify the phase relatio
27、nship.Also, from two CWTs the WTC can be calculated which can be thought of as the local correlation between the time series in time frequency space. Where XWT unveils high common power, WTC finds locally phase locked behavior. The more desirable features of the WTC come at the price of being slight
28、ly less localized in time frequency space. The significance level of the WTC has to be determined using Monte Carlo methods.Wavelet coherence toolbox example举例:The time series we will be analyzing are the winter Arctic Oscillation index (AO) and the maximum sea ice extent in the Baltic (BMI).首先,通过直方
29、图检视两组变量的分布情况AO和BMI,图1显示AO的分布近似于正态分布,而BMI存在着明显的双峰分布(highly bi-modal)情况。这种情况下,我们考虑将BMI数据转换为一组百分位数来代表。转换后的序列相对于气候呈现更好地线性相关性(The transformed series probably reacts more linearly to climate.)。图1百分位数转换函数为boxpdf(),转换之后的分布情况如图2所示。图2随之而来是连续小波变换(Continuous wavelettransform,CWT),其目的在于时间序列在时域空间的表达(The CWT expan
30、ds the time series into time frequency space.)。图3再之后是交叉小波变换(Cross wavelet transform,XWT)。The XWT finds regions in time frequency space where the time series show high common power.图4再下来是小波相干(Wavelet coherence,WTC)。The WTC finds regions in time frequency space where the two time series co-vary共变(but
31、does not necessarily have high power).图5Waveletcoherencebetween the winter Arctic Oscillation index and the sea ice extent in the Baltic. Arrows show thephaserelationship(相位关系)between the two series. Here the arrows point left which means that BMI is in anti-phase with AO. Light shading shows region influenced by edge effects.AcknowledgementsSpecial thank to Dr. Aslak Grinsted, I learned more from his website(/Home).附上
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025届广东省惠州市第三中学物理高二下期末考试试题含解析
- 冰心的课件教学课件
- 宣传宪法的课件
- 二零二五年玻璃密封件生产与销售合同
- 2025年物流运输合同范本790字:物流与家具行业配送协议
- 二零二五年度叉车操作员实习实训及就业合同
- 二零二五年度玻璃钢化粪池智能化改造升级合同
- 二零二五年度补充协议书:供应链管理补充协议
- 2025版购物中心装修工程转让合同协议
- 二零二五年度智能电子产品代工合作协议
- 最新-伤口愈合新进展和美容缝合课件
- 中压交联电缆电缆正、负和零序计算
- 优衣库商业模式分析
- 调度系统介绍课件
- tpo41阅读听力部分参考答案
- 华联学院日语能力考试N5试题二及参考答案
- 大小中型矿山的划分,矿山生产建设规模分类一览表
- 黑布林The Clever Woman 聪明的妇人公开课课件
- 采购年中工作总结汇报PPT(24P)
- 施耐德ATV31变频器说明书
- 房屋建筑构造(地基与基础)课件
评论
0/150
提交评论