小波相关性和相干性.docx

小波相关性和相干性前言 时域指标参数 1. 均值 当观测时间T趋于无穷时，信号在观测时间T内取值的时间平均值就是信号x(t)的均值。均值定义为 mx=lim1TT®¥òx(t)dt0T式中：T是信号的观测区间。实际T不可能为无穷，算出的mx必然包含统计误差，只能作为真值的一种估计。 2. 均方值和方差 当观测时间T趋于无穷时，信号在观测时间T内取值平方的时间平均值就是信号x(t)的均方值，定义为： f2x=lim1TT®¥òT0x(t)dt 2如果仅对有限长的信号进行计算，则结果仅是对其均方值的估计。均方值的正平方根，为均方根值xmax。 方差定义为 s2x=lim1TT®¥2()xt-mdtxò0T方差反应了信号x(t)中的动态部分。方差的正平方根sx称为标准差。若信号x(t)的均值为零，则均方值等于方差。若信号x(t)的均值不为零时，则有下列成立 s3. 概率密度函数 随机信号x(t)的取值落在区间内的概率可用下式表示 Pprb=x<x(t)£x+Dx=limDTT2x=fx-mx 22T®¥式中：DT为信号x(t)取值落在区间(x,x+Dx内的总时间；T为总观察时间。 当Dx®0时，概率密度函数定义为 p(x)=lim1éDTù limêúT®¥DxëTûDx®¥随机信号x(t)的取值小于或等于某一定值d的概率，称为信号的概率分布函数。常用P(x)来表示。概率分布函数的定义为 P(x)=Pprb=x(t)£d=limDTdTT®¥式中：DTd为信号x(t)取值满足x(t)£d的总时间；T为总的观察时间。 1 相关分析 1.1 相关的概念 在信号分析中相关是一个非常重要的概念。所谓相关，就是指变量之间的线性联系或相互依赖关系。经典的互相关用于量化两个信号x(t)和y(t)的相关程度。两个随机信号的互相关RXY定义为： RXY=CXYCXXCYY (1) 式中，CXY=Ex(t)y(t)，CXX=Ex(t)2，CYY=Ey(t)2。 在信号平稳的假设下 CXY=Ex(t)y(t)=1limT®¥Tòx(t)y(t)dt若两个信号间的延迟为t，则定义两个信号的互相关为 RXY(t)=CXY(t)CXXCYY (2) 式中，CXY(t)=1limT®¥Tòx(t)y(t-t)dt，T为信号x(t)和y(t)的观测时间。互相关函数RXY(t)是t的函数，它完整的描述了两信号之间的相关情况或取值依赖关系。 式(2)说明互相关分析提供了两个信号在时域内平移时的相关度的量化方法，也提供了两信号延迟时产生的信号分量的识别方法。这种经典互相关性适用于平稳信号分析。 1.2 小波相关性 小波互相关类似于经典的信号互相关，有效地提供了两个信号相关性对尺度的依赖程度。设两个互相关信号x(t)和y(t)，在给定尺度a和时延u下，x、y的小波互相关性定义为： WC(a,u)=EWXX(a,t)WYY(a,t+u) (3) XY式中，WXX(a,t)和WYY(a,t+u)分别为x(t)和y(t)的小波变换系数。 若分离小波变换系数的实部RWXX(a,t)和虚部IWYY(a,t+u)，只讨论用实部量化给定尺度a下两个信号的相关程度，则小波互相关定义为： WRXY(a,u)=RWCRWCXXXY(a,u)YY(a,0)RWC(a,0) (4) Sello和Bellazzini建议只考虑小波变换的实部，用小波交叉谱WXY(a,u)=WXX(a,u)WYY(a,u)定义小波局部相关系数为： WLCC(a,u)=RWXY(a,u)WXX(a,u)WYY(a,u) (5) 该小波局部相关系数定义的理论基础是经典互相关和交叉小波频谱间的关系： ò+¥-¥x(t)y(t)dt=1Cyòò0+¥+¥-¥RWXY(a,t)dadt。 若考虑小波变换的实部和虚部提供的信息具有一定的联系。则小波互相关定义为： WRXY(a,t)=RWC(a,t)XYXX2+IWC(a,t)XY2WC(a,0)WCYY(a,0) (6) 式中RWCXY和IWCXY分别是式(3)定义的交叉小波相关性函数的实部和虚部。 从经典互相关和小波互相关定义可以看出，小波互相关相比经典互相关引入了参数a，而正是这一参数的引入，使得将经典相关性只在时域内分析两个信号在不同延迟时的相关性，引入到在时频两域内分析信号的相似程度。也就是说小波互相关随尺度的变化而变化，是两个信号在不同延迟时的相关性，因而能够反映出信号互相关最大时，在该频率处两个信号的延迟(相差)等信息，为探测两个信号的相似程度提供更丰富的信息，对分析两个信号的某一频率成分的相似程度有着重要的作用，实现了在时频两域内同时分析两个信号的相关性。 1.3 小波谱 为了探测信号内涵的特征信息及尺度，引入小波谱概念。小波谱首次由Hudgins和Brunet、Colloneau提出，它是基于时域和时间尺度域间的能量守恒定义的： òx-¥+¥2(t)dt=1KyòòCX(a,t)0-¥+¥+¥2dadta (7) 2对于连续信号x(t)，其小波谱Wx(a,t)可用小波系数的模确定： WX(a,t)=CX(a,t)C*X(a,t)=CX(a,t) 2对于离散信号x(t)，设其二进离散小波变换为Cx(j)，则有 说明原始信号的二范数等于小波变换各层细节信号绝对值的平方和。第j层细节的能量 jÎzx2£2åCX(j)2Wx(j)=Cx(j) 称之为小波能量时谱。其相应的Fourier域的二进离散小波能量2Wx(w)=åW(j)exp(-x称之为小波能量频谱。考虑到相移问题及工程实际应用中原始jwN)，信号的实数性，可取实部进行计算，即Wx(w)=åW(j)cos(wtxN)。 由以上分析可知，小波谱是小波局部谱对平移因子的积分，是尺度a的参数，即把Fourier谱的连续频谱划分为离散的频带分量，每个频带代表特定尺度a下信号的频谱信息。小波谱将小波变换和Fourier谱分析结合起来，可在时域中记录信号的突变时间，又可在频域中提取信号突变的频段，通过信号在小波变换的各分解层上的小波谱探测信号内涵的特征。 2. 时间序列信号小波相干性分析 相干函数是计算两个随机过程频谱相关性的直接方法。如果过程是平稳的，即假设它们的频谱不随时间变化，Fourier分析能够实现频谱的精确估计，同样也能给出相干性的精确估计。但它的前提是假设信号是平稳的，而且无法辨别信号的振幅和相位成分。由于实际大地测量时间序列信号受多种因素的影响，其变化过程中产生的信号是非平稳的、动态的，其内涵的特征信息的频率可能随时间而变化，表现出频率分布不均匀。因而，相干性需要作为时间函数研究，这样限制了Fourier分析的应用。为此，引入小波相干性。 引入基于小波相干性的时间过程，是借鉴近来物理学中估计非平稳信号互相作用的研究。与基于Fourier的一致性相比，小波相干性能够以时间函数方式进行相干性分析，因而成为研究动态信号相互作用的有效替代方法。 2.1 相干性函数 对两个连续有限能量信号 x(t)和y(t)，经典的相关性是指两个信号的时间相干性。将这两个信号进行Fourier 变换，设其频率为f，则x和y相干性定义为： r(f)=SSXXXY(f)(f)SYY(f) (10) 式中SXY(f)是两个信号x(t)和y(t)的交叉频谱密度。在信号平稳的假设下，SXY(f)定义为： SXY(f)=òR-¥+¥XY(t)exp(-i2pft)dt (11) 其中，RXY(t)=E(x(t)y(t-t)，E为数学期望。由Schwartz不等式，能够保证(f)值介于0和1之间。 估计两个信号的相干性值要求已知x和y频谱及它们的交叉频谱，可用有限时间段内观测的信号求得。实际情况下，有限时间段数据的频谱和交叉频谱可由整个过程的某一段观测信号来估计。两个有限长度时间序列的交叉频谱的估计： *=x(f)×y(f) (12) Sxy%*(f)是y式中x%(f)是时间序列x(t)在频率为f时的离散Fourier系数，y%(f)的复共轭。该估计是不稳定估计，为提高估计的可靠性，又提出了平滑频谱估计。 2.2 小波相干性 平滑技术仅仅当各分段具有同一频谱估计时才有意义。这种情形下，要求x和y是平稳的，即它们的频谱特征不随着时间变化而变化，也意味着x(和y)能分解为幅值和相位都不变的正弦波的叠加。Fourier分析适用于分析对此类信号。当x和y是非平稳时，Fourier分析和以上的相干性估计就不适用了。 为解决上述非平稳信号分析中存在的问题，近年来又提出应用小波分析估计非平稳信号相干性的方法。与Fourier分析相比，小波分析能够分析具有时变频谱的信号。它能够实现信号的时频分析，即将信号的频谱特征估计为一个时间函数。在一定意义上，小波分析接近于加窗短时Fourier变换。但是，短时Fourier分析的窗的大小是不变的，而小波分析窗能够适用于信号的频率。 Van Milligen et al. and Santoso et al.应用Morlet小波分解信号(也可以选择其它的小波)。Morlet小波的优点在于它具有良好的时间聚集性、较高的频率分辨率、包含相位信息及其与常规信号非常相似等特点，因而用于识别两个非平稳时间序列的关联程度，进行信号的频谱估计。在频率f和时刻t时，Morlet小波定义如下： yyt,ft,f(u)=f×exp(i2pf(u-t)×exp(-(u-t)2s2) (13) (u)是频率为f的正弦波与以时刻为中心的高斯函数的乘积，高斯函数的的标准差与f的倒数成正比例。 1小波相干性 信号x(u)的小波变换是由信号与小波卷积得到的关于频率f和时刻的函数： +¥W信号的光滑小波谱SWXX(t,f)=ò-¥x(u)×y*t,f(u)du (14) Torrence 和 Webster建议通过小波频谱的光滑估计来确定小波相干性。在频率f和时刻t时，定义(t,f)和交叉小波谱SWXXXY(t,f)为： (t,f)dt (15) SWXX(t,f)=òt+d/2t-d/2W*XX(t,f)WXXSWXY(t,f)=òt+d/2t-d/2W*XX(t,f)WYY(t,f)dt (16) 式中*表示复共轭，d是依赖于频率的一个标量。在小波相干性中，d是个很重要的参数，它定义了小波相干性的时间分辨率，d值越小，适应的信号频率越高，因而能够满足相干性的时间变化。与基于Fourier的相干性类似，频率f和时刻t的小波相干性WC(t,f)定义为 WC(t,f)=SW|SWXXXY(t,f)(t,f)|SWYY(t,f)| (17) SWXY(t,f)是由式(16)定义的信号x和y的一个标量，Schwartz不等式能够保证WC(t,f)的值介于0和1之间。当d等于0时，在任意频率f和时刻t，两信号的小波相干性WC(t,f)等于1。 2小波平方相干性 平方相干性用于识别两个时间序列共变过程中的频带。交叉小波能量是一般意义上的能量度量，而小波平方相干性是两个时间序列在时频空间中的协方差强度的度量。Torrence 和 Webster用平滑交叉小波谱的绝对值的平方定义小波平方相干性Rn(t,f)，并用平滑小波能量谱对其进行标准化 22R2W(t,f)=-1XY2(t,f)sWY(t,f)2sWX(t,f)-12 (18) <*>表示在时间和尺度上均进行平滑。因子S-1用于将其转换为能量密度，0Rn(t,f)1。 3. 小波包能量 采用离散小波包变换正交基分解信号的优点是，信号总能量可以分割成各种时频单元。每个单元能量贡献率的数学表达是一个小波函数，其伸缩系数如下： Ef=ò¥-¥f2(t)dt2j2j=ååòm=1n=1¥-¥fj(t)fj(t)dt mn当小波基正交时，公式就变为 2j2fijEf=åEi=1=åòi=1¥-¥fij(t)2dt fj(t)是离散小波变换系数或信号公式中的sj(k)Yj,k(t)或wj(k)jij,k(t)，Eif是信号能量。