一种改进的小波阈值信号去噪方法

张晓宁，孙丽君

(1.河南工业大学信息科学与工程学院，河南郑州 450001;2.河南工业大学电气工程学院，河南郑州 450001)

小波变换与傅里叶变换、窗口傅里叶变换相比，它是一个时间和频率的局部变换，因而能有效地从信号中提取信息，通过伸缩和平移等运算功能对函数或信号进行多尺度细化分析，解决了傅里叶变换不能解决的许多问题。

小波变换被誉为“数学显微镜”，是调和分析发展史上里程碑式的进展。小波理论被认为是对傅里叶分析的重大突破，是近年来信号处理领域的研究热点，许多学者将小波在理论上的研究成果应用到诸如图像压缩、特征提取、信号滤波和数据融合等方面。小波之所以在信号处理领域具有很大的优势，在于小波变换可以获得信号的多分辨率描述，同时，小波变换具有丰富的小波基可适应不同特征的信号。随着小波理论的发展，Mallat［1］提出了模极大值重构滤波，Xu［2］提出了空域相关滤波和Donoho［3］提出了小波域阈值滤波来消除噪声。一般地，不同性质的噪声需要采用不同的消噪方法处理。

上述算法都存在参数的选取问题，不同的参数选取对滤波的效果会有一些差异。Donoho的传统小波阈值去噪方法的实现最简单、计算量最小。但其在理论上找到的最优通用阈值，实际应用中效果并不理想。因此，文中在Donoho的传统小波阈值去噪方法的基础上，改变了小波阈值函数的部分参数得到了一种新的小波阈值函数，既避免了硬阈值函数的不连续性，又可以自适应去噪，从而有效地保存了信号的边缘信息。仿真结果表明，改进后方法有更好的去噪性能。

1 小波阈值的去噪原理

小波变换具有一种“集中”的能力。信号经小波变换后，可以认为由信号产生的小波系数包括有信号的重要信息，其幅值较大，但数目较小，而噪声对应的小波系数幅值小。通过在不同尺度上选取一合适的阈值，并将小于该阈值的小波系数置零，而保留大于阈值的小波系数，从而使信号中的噪声得到有效的抑制，最后进行小波逆变换，得到滤波后的重构信号。

小波去噪的基本思路如图1所示［4］。信号先经过预处理，然后利用小波变换讲信号分解到多尺度上，再对每一层小波系数进行阈值处理，最后对处理后的小波系数进行信号重构。

图1 小波去噪框图

设一个含噪声的一维信号模型可以表示为f(i)=s(i)+n(i)，其中s(i)为原始信号，n(i)为方差σ2的高斯白噪声，服从N(0，σ2)。

1.1 小波硬阈值去噪的步骤

(1)对信号求小波变换。

(3)利用小波变换重构，求出信号的滤波值。

1.2 小波软阈值去噪的步骤

(1)对信号求小波变换。

其中，sgn(x)为符号函数。

(3)进行小波变换重构，求出信号的滤波值。小波系数估计的软、硬阈值处理方法。如图2和图3所示。

图2 硬阈值信号

图3 软阈值信号

2 小波阈值函数的改进

硬阈值函数构造虽然简单，但其在整个小波域内是不连续的，因此重构产生的信号会产生振荡，当噪声水平较高时，这种现象尤为明显，容易出现Pseudo－Gibbs现象。这与实际应用中经常要对阈值函数进行求导运算存在矛盾，具有一定的局限性;并在信号的边缘会产生许多人为的噪声点。在实际情况下，大于阈值的小波系数中也存在噪声信号的干扰，但硬阈值函数只对小于阈值的小波系数进行处理，对大于阈值的小波系数不加处理，这与事实不相符。

软阈值函数虽然在小波域内整体连续性好，但是由于当小波系数较大时，处理过的系数与原系数之间总存在恒定的偏差，这将直接影响重构信号与真实信号的逼近程度，造成一定的高频信息损失，使信号的边缘变得模糊，给重构信号带来不可避免的误差。

在实际应用中，利用软阈值消噪信号比较光滑，但有着较大的信号失真。而利用硬阈值消噪对时变信号消噪效果是有限的。为克服软、硬阈值的缺点，文中将软阈值和硬阈值结合起来，在之前研究的基础上［5－6］增加并优化了小波阈值函数的参数，构造了一类新的阈值函数为

经验证α=0.4时，滤波的效果比较好。现在，关键在于β取何值合适。通过大量的仿真实验，当β=0.618时能获得更好的去噪效果［7］，且处理起来比较方便。

3 仿真结果与评价

(1)为检验新阈值函数的去噪效果，利用matlab软件对一段信噪比为4，并含有高斯白噪声的信号，分别用软、硬阈值方法和文中提出的改进方法进行仿真实验。方法中的 α，β分别取最优值 α=0.4，β=0.618。选用的小波函数为db3，分解层数为3，并在每分解层数上按照Sqtwolog规则进行阈值的估计。去噪结果如图4和图5所示。

图4 软硬阈值去噪结果

(2)信号去噪效果可以用信噪比(SNR)和均方误差(MSE)来描述。定义同一语音信号消噪处理后，均方误差(MSE)越小，信噪比(SNR)越大，则消噪效果越好，定义形式为

图5 本文去噪结果

其中，s(i)是原始信号;s^(i)是由小波阈值处理后的信号。

根据以上仿真和分析，3种方法对含噪声信号去噪后，信号的信噪比SNR和均方误差MSE如表1所示。

表1 信噪比SNR和均方误差MSE

从表中的数据对比可以看出，文中提出的方法比软、硬阈值方法在信噪比SNR和均方误差MSE这两个性能指标上均有明显提高。

4 结束语

小波滤波已经应用于许多理论和应用领域，尤其在信号处理中发挥着越来越重要的作用，小波去噪技术已成为信号处理研究的热点。文中介绍了小波阈值变换的原理以及软、硬阈值变换的原理和步骤，分析了软硬阈值变换的缺点，并且基于软、硬阈值方法对去噪能力的不同特点，提出了一种新的阈值函数，用该方法可以获得比软、硬阈值消噪更好的去噪效果和更高的信噪比，新的阈值函数可以避免高频信息的损失和软阈值方法中对绝对值大的小波系数总存在恒定衰减的不足，并且新的阈值可以提高消噪后信号的信噪比，减少信号的失真和震荡。

［1］MALLET S，HWANG W L.Singularity detection and processing with wavelets［J］.IEEE Trans on Information Theory，1992，38(2):617 －643.

［2］XU Yansun，JOHN B W，DENNIS M H，et al.Wavelet transform domain filters:a spatially selective noise filtration technique ［J］.IEEE Trans on Image Processing，1994，3(6):747－758.

［3］DONOHO D L.De － noising by soft－ threshold［J］.IEEE Trans on Information.Theory，1995，41(3):613 －627.

［4］潘泉，张磊，孟晋丽.小波滤波方法及应用［M］.北京:清华大学出版社，2005.

［5］崔华，宋国乡.基于小波阈值去噪方法的一种改进方案［J］.测控技术，2005(1):8 －10.

［6］张弦，李世平，孙浚清，等.一种改进的小波阈值去噪方法［J］.软件天地，2007，23(31):309 －311.

［7］谢搴，詹毅，牛聪.小波阈值去噪黄金分割法［J］.物探与化探，2006，30(3):254－257.