小波相关性和相干性

1、前言时域指标参数1. 均值 当观测时间趋于无穷时，信号在观测时间内取值的时间平均值就是信号的均值。均值定义为 （1）式中：是信号的观测区间。实际不可能为无穷，算出的必然包含统计误差，只能作为真值的一种估计。2. 均方值和方差 当观测时间趋于无穷时，信号在观测时间内取值平方的时间平均值就是信号的均方值，定义为： （2）如果仅对有限长的信号进行计算，则结果仅是对其均方值的估计。均方值的正平方根，为均方根值（或有效值）。 方差定义为 （3）方差反应了信号中的动态部分。方差的正平方根称为标准差。若信号的均值为零，则均方值等于方差。若信号的均值不为零时，则有下列成立 （4）3. 概率密度函数 随机信号的

2、取值落在区间内的概率可用下式表示 （5）式中：为信号取值落在区间内的总时间；为总观察时间。当时，概率密度函数定义为 （6）随机信号的取值小于或等于某一定值的概率，称为信号的概率分布函数。常用来表示。概率分布函数的定义为 （7）式中：为信号取值满足的总时间；为总的观察时间。1 相关分析1.1 相关的概念在信号分析中相关是一个非常重要的概念。所谓相关，就是指变量之间的线性联系或相互依赖关系。经典的互相关用于量化两个信号和的相关程度。两个随机信号的互相关定义为： (1)式中，。在信号平稳的假设下 若两个信号间的延迟为，则定义两个信号的互相关为 (2)式中，为信号和的观测时间。互相关函数是的函数，它完

3、整的描述了两信号之间的相关情况或取值依赖关系。式(2)说明互相关分析提供了两个信号在时域内平移时的相关度的量化方法，也提供了两信号延迟时产生的信号分量的识别方法。这种经典互相关性适用于平稳信号分析。1.2 小波相关性小波互相关类似于经典的信号互相关，有效地提供了两个信号相关性对尺度的依赖程度。设两个互相关信号和，在给定尺度和时延下，x、y的小波互相关性定义为： (3)式中，和分别为和的小波变换系数。若分离小波变换系数的实部和虚部，只讨论用实部量化给定尺度下两个信号的相关程度，则小波互相关定义为： (4)Sello和Bellazzini建议只考虑小波变换的实部，用小波交叉谱定义小波局部相关系数为

4、： (5)该小波局部相关系数定义的理论基础是经典互相关和交叉小波频谱间的关系： 。若考虑小波变换的实部和虚部提供的信息具有一定的联系。则小波互相关定义为： (6)式中和分别是式(3)定义的交叉小波相关性函数的实部和虚部。从经典互相关和小波互相关定义可以看出，小波互相关相比经典互相关引入了参数a，而正是这一参数的引入，使得将经典相关性只在时域内分析两个信号在不同延迟时的相关性，引入到在时频两域内分析信号的相似程度。也就是说小波互相关随尺度的变化而变化，是两个信号在不同延迟时的相关性，因而能够反映出信号互相关最大时，在该频率处两个信号的延迟(相差)等信息，为探测两个信号的相似程度提供更丰富的信息，

5、对分析两个信号的某一频率成分的相似程度有着重要的作用，实现了在时频两域内同时分析两个信号的相关性。1.3 小波谱为了探测信号内涵的特征信息及尺度，引入小波谱概念。小波谱首次由Hudgins和Brunet、Colloneau提出，它是基于时域和时间尺度域间的能量守恒定义的： (7)对于连续信号，其小波谱可用小波系数的模确定： （8）对于离散信号，设其二进离散小波变换为，则有 （9）说明原始信号的二范数等于小波变换各层细节信号绝对值的平方和。第j层细节的能量 称之为小波能量时谱。其相应的Fourier域的二进离散小波能量，称之为小波能量频谱。考虑到相移问题及工程实际应用中原始信号的实数性，可取实部

6、进行计算，即。由以上分析可知，小波谱是小波局部谱对平移因子的积分，是尺度（离散化表示为）的参数，即把Fourier谱的连续频谱划分为离散的频带分量，每个频带代表特定尺度下信号的频谱信息。小波谱将小波变换和Fourier谱分析结合起来，可在时域中记录信号的突变时间，又可在频域中提取信号突变的频段，通过信号在小波变换的各分解层上的小波谱探测信号内涵的特征。2. 时间序列信号小波相干性分析相干函数是计算两个随机过程频谱相关性的直接方法。如果过程是平稳的，即假设它们的频谱不随时间变化，Fourier分析能够实现频谱的精确估计，同样也能给出相干性的精确估计。但它的前提是假设信号是平稳的，而且无法辨别信号

7、的振幅和相位成分。由于实际大地测量时间序列信号受多种因素的影响，其变化过程中产生的信号是非平稳的、动态的，其内涵的特征信息的频率可能随时间而变化，表现出频率分布不均匀。因而，相干性需要作为时间函数研究，这样限制了Fourier分析的应用。为此，引入小波相干性。引入基于小波相干性的时间过程，是借鉴近来物理学中估计非平稳信号互相作用的研究。与基于Fourier的一致性相比，小波相干性能够以时间函数方式进行相干性分析，因而成为研究动态信号相互作用的有效替代方法。2.1 相干性函数对两个连续有限能量信号 和，经典的相关性是指两个信号的时间相干性。将这两个信号进行Fourier 变换，设其频率为f，则x

8、和y相干性定义为： (10)式中是两个信号和的交叉频谱密度。在信号平稳的假设下，定义为： (11)其中，E为数学期望。由Schwartz不等式，能够保证(f)值介于0和1之间。估计两个信号的相干性值要求已知和频谱及它们的交叉频谱，可用有限时间段内观测的信号求得。实际情况下，有限时间段数据的频谱和交叉频谱可由整个过程的某一段观测信号来估计。两个有限长度时间序列的交叉频谱的估计： (12)式中是时间序列在频率为时的离散Fourier系数，是的复共轭。该估计是不稳定估计，为提高估计的可靠性，又提出了平滑频谱估计。2.2 小波相干性平滑技术仅仅当各分段具有同一频谱估计时才有意义。这种情形下，要求和是平

9、稳的，即它们的频谱特征不随着时间变化而变化，也意味着(和)能分解为幅值和相位都不变的正弦波的叠加。Fourier分析适用于分析对此类信号。当x和y是非平稳时，Fourier分析和以上的相干性估计就不适用了。为解决上述非平稳信号分析中存在的问题，近年来又提出应用小波分析估计非平稳信号相干性的方法。与Fourier分析相比，小波分析能够分析具有时变频谱的信号。它能够实现信号的时频分析，即将信号的频谱特征估计为一个时间函数。在一定意义上，小波分析接近于加窗短时Fourier变换。但是，短时Fourier分析的窗的大小是不变的，而小波分析窗能够适用于信号的频率。Van Milligen et al.

10、and Santoso et al.应用Morlet小波分解信号(也可以选择其它的小波)。Morlet小波的优点在于它具有良好的时间聚集性、较高的频率分辨率、包含相位信息及其与常规信号非常相似等特点，因而用于识别两个非平稳时间序列的关联程度，进行信号的频谱估计。在频率和时刻时，Morlet小波定义如下： (13)是频率为的正弦波与以时刻为中心的高斯函数的乘积，高斯函数的的标准差与f的倒数成正比例。1小波相干性信号的小波变换是由信号与小波卷积得到的关于频率和时刻的函数： (14)Torrence 和 Webster建议通过小波频谱的光滑估计来确定小波相干性。在频率和时刻时，定义信号的光滑小波谱和

11、交叉小波谱为： (15) (16)式中*表示复共轭，是依赖于频率的一个标量。在小波相干性中，是个很重要的参数，它定义了小波相干性的时间分辨率，值越小，适应的信号频率越高，因而能够满足相干性的时间变化。与基于Fourier的相干性类似，频率和时刻的小波相干性定义为 (17)是由式(16)定义的信号的一个标量，Schwartz不等式能够保证的值介于0和1之间。当等于0时，在任意频率f和时刻t，两信号的小波相干性等于1。2小波平方相干性平方相干性用于识别两个时间序列共变过程中的频带。交叉小波能量是一般意义上的能量度量，而小波平方相干性是两个时间序列在时频空间中的协方差强度的度量。Torrence 和 Webster用平滑交叉小波谱的绝对值的平方定义小波平方相干性，并用平滑小波能量谱对其进行标准化