有色噪声在物理小波变换下功率谱分布的研究.pdfVIP

有色噪声在物理小波变换下 功率谱分布研究 郑庆庆高静怀 (西安交通大学电子与信息工程学院波动与信息研究所．西安710049) 1．引言 实际问题中的噪声多为有色噪声(如红噪声等)。在用Fourier变换、小波变换等进行信 号分析时，掌握有色噪声在该变换下的分布规律，是辨识有效信号和噪声的基础。这方面的 研究是目前很活跃的研究领域…之一。本文研究有色噪声在匹配地震子波的物理小波(以下 简称物理小波)下的分布规律，推导出利用物理小波分析信号时必须的参数，最后把物理小 波用于有色噪声中信号的辨识。 2．原理 2．1红噪声的傅立叶谱与小波谱 一个简单的模型就是单变量滞后为1的自回归过程【AR(1)】[11： 工^=锨^．1-I-Z^ 其中口为设定的延迟因子为1(1ag-1)的自相关系数，Xo=0，z。为高斯白噪声，它 的平均功率谱(归一化后)为： 最=雨万瓦1--孤0[2瓦厕 ‘2) 其中k=0…N12是频率域样点序号。通过选择—个恰当的lag一1自相关系数，我们可以用 (2)式来模拟红噪声谱。注意到a=0时便是白噪声谱。我们利用自适应滤波器来解决该 参数的估计问题(详细讨论省略)。 小波变换定义为： ．8 矗 一 lI ^X ∞ 3 ％L：， 、一、 七 ．y LV七 M∑脚 式中，YcI为信号的傅立叶变换，k=O…N一1，s为尺度因子，w：G)为信号的小波变换。 文献【l】以Morlet小波、Paul小波等为基本小波，研究结果表明：红噪声的小波功率谱服从z； 分布，在给定的时间f，平均小波功率谱满足(2)式。 2．2与物理小波有关的几个物理量表达式 163 子波的幅值，f为能量衰减因子，卢为能量延迟时间，‰为视频率。该小波的傅立叶变换 为小波能量下降到它的最大值的e-2时的时间点。可以证明(这里省略)，对于物理小波， 影响锥因子长度为：△r=(√丽+卢)s。 为了将以尺度为变量的小波谱转换为时一频谱，必须找到给定的小波函数中尺度与频率 (周期的倒数)的对应关系。设允为傅立叶周期，它与尺度的关系为(推导过程略)： 4z 允= ， 对物理小波归一化， 令 门g。(f)12=1得： exp(邮2)，相应的归一化的物 go(f)=(等)1Hexp(--z(t一卢)2)exp(itOot)，go(0)=(等)¨4 1 理小波的傅氏变换为：莒o(co)=(≤【．)¨4 Z巧I ，·r“ 对誊(s∞)归一化得：营(sc0)=(竿)¨2营o(sco)。 讲 3有色噪声在物理小波下功率谱分布 与Morlet小波情况类似，以物理小波为基本小波的小波变换也是线性变换，我们猜测 在这种小波变换下，有色噪声的小波功率谱应服从z；分布，局部小波谱(指定时间处的谱) 可能满足(2)式。为了验证上述猜想，我们进行MonteCarlo仿真。图l是仿真结果，从 图可见，理论曲线与仿真试验结果是一致的。 图2和图3是用不同的基本小波对同一含噪信号进行小波变换后得到得小波功率谱，结 果表明，用物理小波作为基本小波，有效信号能量在时一频域的分布更加集中。图4中有效 信号为两个，其中含有有色噪声(其它参数及说明见图注)。图6为能量分布空间中切出的 四个局部谱，从这些局部谱可见，只有在有效信号处，局部谱的能量超过95％的置信线。 利用此特征，我们可辨识有效信号和噪声。 由以上