时间序列的小波分æž?/p>
时间序列ï¼?/p>
Time
Series
)是地学研究中经常遇到的问题。在时间序列研究中,时域和频域是常用的两
种基本形式ã€?/p>
其中ï¼?/p>
时域分析具有时间定位能力ï¼?/p>
但无法得到关于时间序列变化的更多信息ï¼?/p>
频域分析
(如
Fourier
变换)虽具有准确的频率定位功能,但仅适合平稳时间序列分析。然而,地学中许多现象(如河å·?/p>
径流、地震波、暴雨、洪水等)随时间的变化往往受到多种因素的综合影响,大都属于非平稳序列,它们
不但具有趋势性、周期性等特征,还存在随机性、突变性以及“多时间尺度”结构,具有多层次演变规律ã€?/p>
对于这类非平稳时间序列的研究,通常需要某一频段对应的时间信息,或某一时段的频域信息。显然,æ—?/p>
域分析和频域分析对此均无能为力ã€?/p>
20
世纪
80
年代初,ç”?/p>
Morlet
提出的一种具有时
-
频多分辨功能的小波分æž?/p>
ï¼?/p>
Wavelet Analysis
)为更好
的研究时间序列问题提供了可能,它能清晰的揭示出隐藏在时间序列中的多种变化周期,充分反映系统在
不同时间尺度中的变化趋势,并能对系统未来发展趋势进行定性估计ã€?/p>
目前,小波分析理论已在信号处理、图像压缩、模式识别、数值分析和大气科学等众多的非线性科å?/p>
领域内得到了广泛的应。在时间序列研究中,小波分析主要用于时间序列的消噪和滤波,信息量系数和分
形维数的计算,突变点的监测和周期成分的识别以及多时间尺度的分析等ã€?/p>
一、小波分析基本原ç�?/p>
1.
小波函数
小波分析的基本思想是用一簇小波函数系来表示或逼近某一信号或函数。因此,小波函数是小波分æž?/p>
的关键,它是指具有震荡性、能够迅速衰减到零的一类函数,即小波函æ•?/p>
)
R
(
L
)
t
(
2
�/p>
�/p>
且满足:
�
ï‚?
�/p>
�/p>
0
dt
)
t
(
�/p>
ï¼?/p>
1
ï¼?/p>
式中ï¼?/p>
)
t
(
�/p>
为基小波函数,它可通过尺度的伸缩和时间轴上的平移构成一簇函数系ï¼?/p>
)
a
b
t
(
a
)
t
(
2
/
1
b
,
a
�/p>
�/p>
�/p>
�/p>
�/p>
其中ï¼?/p>
0
a
R,
b
a,
ï‚?/p>
�/p>
ï¼?/p>
2
ï¼?
式中ï¼?/p>
)
t
(
b
,
a
�/p>
为子小波ï¼?/p>
a
为尺度因子,反映小波的周期长度;
b
为平移因子,反应时间上的平移ã€?/p>
需要说明的是,选择合适的基小波函数是进行小波分析的前提。在实际应用研究中,应针对具体情å†?/p>
选择所需的基小波函数;同一信号或时间序列,若选择不同的基小波函数,所得的结果往往会有所差异ï¼?/p>
有时甚至差异很大。目前,主要是通过对比不同小波分析处理信号时所得的结果与理论结果的误差来判å®?/p>
基小波函数的好坏,并由此选定该类研究所需的基小波函数ã€?/p>
2.
小波变换
è‹?/p>
)
t
(
b
,
a
�/p>
是由ï¼?/p>
2
)式给出的子小波,对于给定的能量有限信号
)
R
(
L
)
t
(
f
2
�/p>
,其连续小波变换
ï¼?/p>
Continue Wavelet Transform
,简写为
CWT
)为ï¼?/p>
dt
)
a
b
t
(
f(t)
a
)
b
,
a
(
W
R
2
/
1
-
f
�/p>
�/p>
�/p>
�/p>
ï¼?/p>
3
ï¼?/p>
式中ï¼?/p>
)
b
,
a
(
W
f
为小波变换系数;
f(t)
为一个信号或平方可积函数ï¼?/p>
a
为伸缩尺度;
b
平移参数ï¼?
)
a
b
x
(
�/p>
�/p>
�/p>
)
a
b
x
(
�/p>
�/p>
的复共轭函数ã€?/p>
地学中观测到的时间序列数据大多是离散的,
设函æ•?/p>
)
t
k
(
f
�/p>
ï¼?/p>
ï¼?/p>
k=1,2,
�/p>
,N;
t
�/p>