小波分析整理.pdfVIP

青海湖地区近50 年降水量周期变化分析 杨沈斌，张弥，吕开龙 1．问题 青海湖位于青海省东北部的青海湖盆地内，既是中国最大的内陆湖泊，也是中国最大 的咸水湖。近年来有多篇关于青海湖面积受气候变化影响的研究报告，认为温度升高，降水 减少和蒸发量大是造成面积下降的几个重要原因。为此，本实验拟采用连续小波分析方法对 青海湖地区年降水量周期变化进行分析，探讨降水量变化与青海湖面积变化的关系。 2. 资料 考虑到刚察气象观测站正好处于青海湖边上，所得数据对反映青海湖周边降水变化具有 较好的代表性。因此，以青海湖西北位置的刚察气象台站1961-2007 年年降水量资料为例， 对该站年降水资料进行小波分析，获取其周期变化特征。同时，获取了青海湖1961-2007 年 的面积变化资料。从资料发现，青海湖面积近50 年呈现下降趋势，平均下降3 km2/a 。图1 和图2 分别显示了1961-2007 年刚察气象站年降水量和青海湖面积的变化，以及对应的趋势 线。 图1 1961-2007 年刚察气象台站年降水量及趋势线 图2 1961-2007 年青海湖面积及趋势线 3. 小波分析方法 小波变换具有多分辨率分析的特点，并且在时频两域都具有表征信号局部特征的能力。 小波变换通过将时间系列分解到时间频率域内，从而得出时间系列的显著的波动模式，即周 期变化动态，以及周期变化动态的时间格局（Torrence and Compo, 1998 ）。小波（Wavelet ）， 即小区域的波，是一种特殊的、长度有限，平均值为零的波形。它有两个特点：一是“小”， 二是具有正负交替的“波动性”，即直流分量为零。小波分析是时间（空间）频率的局部化 分析，它通过伸缩平移运算对信号（函数）逐步进行多尺度细化，能自动适应时频信号分析 的要求，可聚焦到信号的任意细节。小波分析将信号分解成一系列小波函数的叠加，而这些 小波函数都是由一个母小波（mother wavelet ）函数经过平移与尺度伸缩得来的。用这种不 规则的小波函数可以逼近那些非稳态信号中尖锐变化的部分，也可以去逼近离散不连续具有 局部特性的信号，从而更为真实的反映原信号在某一时间尺度上的变化。小波分析这种局部 分析的特性使其成为对非稳态、不连续时间序列进行量化的一个有效工具（Stoy et al., 2005）。 小波是一个具有零均值且可以在频率域与时间域内进行局部化的数学函数（Grinsted et al., 2004 ）。一个小波被称为母小波（mother wavelet ），母小波可沿着时间指数经过平移与尺 度伸缩得到一系列子小波。子小波可以通过尺度（s，频率的反函数）函数和时间（n ）位置 或平移来描述。利用一系列子小波，一个信号可以在不同的时间尺度上进行计算并显示出详 细的特征尺度。拉伸更大的小波窗口，使其宽度更大便可以分析时间系列中波动较大的部分 并捕捉大尺度（低频）事件的特征。相反，压缩较小的窗口将包含小尺度（高频）的事件信 息。当信号被子小波相乘，被s 与n 唯一的表达，我们可以计算出信号在时间频率域一个具 体位置的系数。如果信号在时间n 上的谱成分可以与小波s 比较，那么计算的小波系数具有 相对较大的值。在其它n 与s 的组合（如其它的子小波）上都进行这样的计算，那么将会产 生一系列系数（小波变换）来表达信号在时间频率域内的分解。通过这样的变化便可得到时 间系列的波动模式（周期变化模式）以及这些模式随时间的变化（Furon et al., 2008; Jevrejeva et al., 2003 ）。 小波变换可以分为连续小波变换（the Continuous Wavelet 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 尺度化以及转