共振稀疏分解.docVIP

共振稀疏分解：一种新的可稀疏信号的分析方法 摘要 生命和物质过程会产生大量信号，这些信号不但是不稳定的，而且是持续震荡信号和瞬态冲击信号的混合，并且这两种信号是很难线性分解的，例如声音、医疗和地理信号。因此，本文描述了一种基于信号共振的非线性信号分析方法，而这种方法不基于傅里叶变换和小波变换产生的频谱和幅值。这种方法将信号分解成一个高共振分量和一个低共振分量——高共振分量由多个同时发生的持续震荡信号组成，低共振分量由多个没有具体现状和持续时间的瞬态冲击信号组成。本文所阐述的共振稀疏分解算法使用的方法有信号稀疏表示、形态分量分析和品质因子可调小波换。 前言 频域分析法和滤波是信号处理的基础。然而，频域分析法和时频分析法并不适用于所有信号，事实上只适用于持续震荡或周期信号。那些主要由奇异点限定的分段光滑信号多数使用时域和小波变换描述、分析和处理。例如，图像扫描，眼部运动记录，潜能诱发反应，神经尖刺训练等。 然而，许多生命和物质过程产生信号不只是不稳定的，而且是持续震荡信号和瞬态冲击信号的混合，例如声音、医疗（脑电图和心电图等）和地理（海浪高度数据等）信号。这些信号既含有稳态震荡部分又含有瞬态冲击部分。脑电波包含有节奏振荡（alpha和beta波等），也包含人为测量和无节奏脑行为所产生的瞬态冲击。海浪高度数据测量的是已经流动了几百英里（100‘s）的海量的重叠高度，但是天气因素将中断这种震荡行为。当然，通过生命和物质系统测量的信号通常包含持续震荡信号和瞬态冲击信号，而这两种信号是很难线性分解的。 为了改进复杂非平稳信号的描述、分析和处理，我们阐述了一种新的基于共振的非线性信号分析方法，而这种方法不基于傅里叶变换和小波变换产生的频谱和幅值。这种方法将信号分解成一个高共振分量和一个低共振分量。其中，高共振分量由多个同时发生的持续震荡信号合成，另一方面，低共振分量由多个没有具体现状和持续时间的瞬态冲击信号合成。 这篇论文的部分内容已经出版在两个早期的会议论文中[84，85]。 图1. 单脉冲共振属性是由品质因子Q量化的，而品质因子是中心频率与频带宽度的比值。脉冲1和脉冲3在持续时间上看是一个单振荡，是低共振脉冲。脉冲2和脉冲4具有多次震荡，属于高共振分量。低品质因子小波变换（例如经典的二阶小波变换）能有效描述脉冲1和脉冲3，高品质因子小波变换能有效描述脉冲2和脉冲4。图a为时域信号，图b为频谱。 信号共振 图1用图例说明了信号共振的概念。事实上，脉冲1和脉冲3都包含了一个单振荡正弦信号。我们把这两种信号称为低共振信号是因为他们没有持续震荡。观察这两个脉冲的时域图我们发现，时域脉冲波形并不能有限体现共振程度。很明显，一个低共振脉冲既有可能是一个高频信号（脉冲1），也有可能是一个低频信号（脉冲3）。低共振脉冲不限于单边频带。因此，不能通过频率滤波的方式从一个信号中提取出低共振分量。 我们把脉冲2和脉冲4定义成高共振分量是因为他们具有持续震荡特性。两个脉冲都包含5个振荡的由哑铃函数相乘而成的正弦波（特别的，如布莱克曼窗口），如上面所说，这两个脉冲的时域波形具有相同的共振属性。同样的，一个高共振脉冲既有可能是一个高频信号（脉冲2），也有可能是一个低频信号（脉冲4），高共振分量也不能通过频率滤波的方式从一个信号中提取出来。 共振稀疏分解 正如我们所描述的，共振稀疏分解应该能够近似低分解图1中的脉冲1和脉冲2，尽管这两个信号在时域上是迭加在一起的。为了阐明共振稀疏分解算法(后面会详尽阐述)的效果，我们将这种方法应用图2的人工合成信号上。这个人工信号含有3种频率和两级共振的六个脉冲。目标是将信号的高低共振分量分离开。通过这种算法的得到的高低共振分量如图2a所示。这种算法也得到一个残余信号，从而允许随机噪声的存在。这个测试信号等同于三个信号的和：高共振分量、低共振分量和残余信号。（残余信号的幅值可以通过分解算法中的参数控制。） 图2. 共振稀疏分解和频率滤波时域波形图。(a)测试信号被稀疏表示成高低共振分量。高共振信号分量使用高品质因子RADWT变换得到。同样，低共振信号分量使用低品质因子RADWT变换得到。(b)使用LTI时间离散滤波器将测试信号分解成低、中和高频分量。(a)共振稀疏分解时域波形图(b)频率滤波时域波形图。 观察发现，线性时不变（LTI）滤波器不能实现图2a所示的信号分解，这是因为存在于高共振分量的三个频率同样也存在于地共振分量中。高低共振分量中的脉冲的不同点不是他们的频率而是他们持续振荡的程度。 当然，LTI滤波器能将测试信号分解成低、中和高频部分。使用低通、带通和高通LTI滤波器，我们能够实现测试信号的基于频率的分解，从而获得不同的频率分量，正如图2b中所示。 共振稀疏分解必然是非线性的 这里所提到的共振稀疏分级不能通过线性滤波实现