小波方差制作步骤.docVIP

第六章 时间序列的小波分析 时间序列（Time Series）是地学研究中经常遇到的问题。在时间序列研究中，时域和频域是常用的两种基本形式。其中，时域分析具有时间定位能力，但无法得到关于时间序列变化的更多信息；频域分析（如Fourier变换）虽具有准确的频率定位功能，但仅适合平稳时间序列分析。然而，地学中许多现象（如河川径流、地震波、暴雨、洪水等）随时间的变化往往受到多种因素的综合影响，大都属于非平稳序列，它们不但具有趋势性、周期性多时间尺度结构，具有多层次演变规律20世纪80年Morlet提出的具有时-频多分辨功能的它清的揭示出隐藏在时间序列中的种变化周期，反映系统在不同时间尺度中的变化趋势，并能对系统未来发展趋势进行定性估计。小波分析理论已在信号处理、图像压缩、模式识别、数值分析和大气科学等众多的非线性科学领域内且满足： （1） 式中，为基小波函数，它可通过尺度的伸缩和时间轴上的平移构成一簇函数系： 其中， （2） 式中，为子小波；a为尺度因子，反映小波的周期长度；b为平移因子，反应时间上的平移。 需要说明的是，选择合适的基小波函数是进行小波分析的前提。在实际应用研究中，应针对具体情况选择所需的基小波函数；同一信号或时间序列，若选择不同的基小波函数，所得的结果往往会有所差异，有时甚至差异很大。目前，主要是通过对比不同小波分析处理信号时所得的结果与理论结果的误差来判定基小波函数的好坏，并由此选定该类研究所需的基小波函数。 2. 小波变换 若是由（2）式给出的子小波，对于给定的能量有限信号，其连续小波变换（Continue Wavelet Transform，简写为CWT）为： （3） 式中，为小波变换系数；f(t)为一个信号或平方可积函数a为伸缩尺度b平移参数为的复共轭，（k=1,2,…,N; 为取样间隔），则式（3）的离散小波变换形式为： （4） 由式（3）或（4）可知小波分析的基本原理，即通过增加或减小伸缩尺度a来得到信号的低频或高频信息，然后分析信号的概貌或细节，实现对信号不同时间尺度和空间局部特征的分析。 实际研究中，最主要的就是要由小波变换方程得到小波系数，然后通过这些系数来分析时间序列的时频变化特征。 3. 小波方差 将小波系数的平方值在b域上积分，就可得到小波方差，即 () 小波方差随尺度a的变化过程，称为小波方差图。由式()可知，它反映波动的能量随尺度a的分布。因此，小波方差图确定中种尺度扰动的相对强度存在的主要时间尺度，即主周期。 ⑴计算小波系数；⑵绘制小波系数图（实部、模和模方）、小波方差图和主周期变化趋势图，并分别说明各图在分析径流多时间尺度变化特征中的作用。 表1 某流域某水文观测站1966-2004年实测径流数据（×108m3） 年份 径流量 年份 径流量 年份 径流量 年份 径流量 年份 径流量 1966 1.438 1974 2.235 1982 0.774 1990 1.806 1998 1.709 1967 1.151 1975 4.374 1983 0.367 1991 0.449 1999 0.000 1968 0.536 1976 4.219 1984 0.562 1992 0.120 2000 0.000 1969 1.470 1977 2.590 1985 3.040 1993 0.627 2001 2.104 1970 3.476 1978 3.350 1986 0.304 1994 1.658 2002 0.009 1971 4.068 1979 2.540 1987 0.728 1995 1.025 2003 3.177 1972 2.147 1980 0.807 1988 0.492 1996 0.955 2004 0.921 1973 3.931 1981 0.573 1989 0.007 1997 1.341 分析 1. 选择合适的基小波函数是前提 在运用小波分析理论解决实际问题时，选择合适的基小波函数是前提。只有选择了适合具体问题的基小波函数，才能得到较为理想的结果。目前，可选用的小波函数很多，如Mexican hat小波、Haar小波、Morlet小波和Meyer小波等。在本例中，我们选用Morlet连续复小波变换来分析径流时间序列的多时间尺度特征。原因如下： 1.1 径流演变过程中包含“多时间尺度”变化特征且这种变化是连续的，所以应采用连续小波变换来进行此项分析。 1.2实小波变换只能给出振幅和正负