基于小波变换微弱生命信号提取的研究
一般含噪的一维信号的模型可表示为：
s(k)=f(k)+εe(k)，k=0，1，…n-1 (3)
式中，s(k)为含噪信号，f(k)为有用信号，e(k})为噪声信号。
利用小波检测微弱生命信号的实质是提取强噪声背景下的生命信号，这个过程即去噪，在小波去噪的方法中比较常用的是阈值去噪法。
小波阈值去噪可分为3部分：
(1)信号的小波分解选择一个小波函数对信号进行分解计算。
(2)小波分解高频系数的阂值量化 对各分解尺度下的高频系数选择一个阈值进行阈值量化处理。
(3)小波重构 根据小波分解的最底层低频系数和各层高频系数进行小波重构。
最关键的是阈值的选择以及阈值的量化，该步骤完成的好坏决定信号消噪的质量。在阈值去噪中，阈值函数体现了对超过和低于阈值的小波系数模的不同处理策略以及不同估计方法。设ω是原始小波系数，η(ω)表示阈值化后的小波系数，T是阈值，I(x)为示性函数。
常见阈值函数有：(a)硬阈值函数(图3a)，η(ω=ωI(|ω|>T)；(b)软阈值函数(图3b)，η(ω)=(ω-sgn(ω)T)I(|ω|>T)。本文引用地址：http://www.eepw.com.cn/article/195600.htm
小波阈值去噪方法除阈值函数的选取外，另一个关键因素是阈值估计。如果阈值太小，去噪后的信号仍然有噪声存在；阈值太大，重要的信号特征又被过滤掉，引起偏差。从直观上看，对于给定的小波系数，噪声越大，阈值就越大。常见的阈值估计方法有Visushrik阈值、SUREShrink阈值、GCV阈值等。其中，SUREShrink阈值估计方法是在SURE准则下得到的阈值，该准则是在均方差准则的无偏估计，专门针对软阈值函数得出的结论，且阈值趋近于理想值。设原始信号小波系数估计通过软阈值函数萎缩得到，即
阈值的选择可通过下面的风险函数定义：
由于小波变换的正交性，风险函数可以写成：
可以证明，当V服从Guass分布时，有下面的等式成立
式中，P(|Yi|>t)服从二项分布，其概率可用|Yi|>t出现的频率近似，可得到风险函数的表达式如下：
式中，I是示性函数，^表示两数取小。
则最佳闽值选择可以通过最小化风险函数得到，即
，对于最佳阈值的选择可以在一个有限的范同内，即t*∈{Y1，Y2，…，YN}。在实际应用中，SUREShrink阈值去噪法能获得较为满意的去噪效果，这是一种误差较低的阈值去噪方法。
小波分解生命信号和噪声所产生的高频系数叠加构成信号的高频系数向量，一般的阈值选取方法会将生命信号的高频部分当作噪声信号滤除掉，而采用基于Stein无偏似然估计SUREShrink阈值选取规则，可以保留生命信号的高频部分，这对于有极少一部分高频信息在噪声范围内的信号很有用。
4 小波去噪的MATLAB仿真
一般检测到的微弱生命信号的背景强噪声主要是工频干扰信号，因此采用正弦信号模拟人体心跳信号频率为0．7 Hz、幅度是1，模拟的工频干扰信号频率为50 Hz、幅度是心跳信号的10倍，和Matlab提供的噪声noissin信号叠加，可近似组成强噪声背景下的生命信号，采用db3小波进行信号分解，并对信号进行SUREShrink阈值估计，并采用heursure函数实现。
叠加信号去噪仿真图如图4所示，叠加信号经过小波阈值去噪法去噪后，可得到较好的生命信号，小波分解和重构的细节，如图5和图6所示。根据Mallat算法的基本思想，高频信号和低频信号分别可以从图中反映出来，其中a1和d1分别反映模拟生命信号的正弦信号，和强噪声干扰的工频信号，这就说明对微弱生命信号的提取小波可以取得很好的效果，由于这里所使用的是模拟的生命信号，在实际应用时还应进行改进。
5 结束语
生命信号由于本身的特点，传统的傅里叶变换对其消噪和提取显得无能为力，因为傅里叶变换对信号的分析只是在频域中进行，不能反映信号某一点的变化情况，而小波变换可以对信号在时频两域进行分析，很适合探测信号的瞬时状态，对微弱生命信号可以进行有效去噪和提取。通过仿真表明，小波变换很适合微弱生命信号的检测，可以在这一领域发挥重要作用。
评论