基于MATLAB小波去噪方法的研究

《基于MATLAB小波去噪方法的研究》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、!!第!2!0!卷!!第!!7期!!!!!!!!!!!!!!!计!!!算!!!机!!!仿!!!!真!!!!!!!!!!!!!!!!!2!00!3!年!7!月!!!!文章编号：1006-934（82003）07-0119-04基于MATLAB小波去噪方法的研究杜浩藩，丛爽（中国科学技术大学自动化系，安徽合肥230027）摘要：在实际的计算机控制系统中，采样信号不可避免的受到各种噪声和干扰的污染，使得由辨识采样信号得到的系统模型存在偏差而妨碍了系统控制精度的提高。Donoho和Johnstone提出的小波收缩去噪算法对去除叠加性高斯白噪声非常有效

2、。我们对此在MATLAB环境下做了详尽的仿真研究，并利用滤波器、傅立叶变换/反变换和小波收缩三种方法对实际电机随动系统的辨识采样数据进行了去噪对比实验，总结了应用小波去噪的一些实际经验。关键词：小波变化；滤波；去噪中图分类号：TP301.6文献标识码：A1引言其中di为含噪信号，fi为“纯净”采样信号，zi为独立同分布iid在计算机控制中，需要从控制系统中获得采样信号（观的高斯白噪声N（0，1）（用符号可表示为：zi~N（0，1）），!为测数据）用于系统的模型辨识或实时控制，而这些采样信号噪声水平，信号长度为N。将不可避免的受到噪声的污染。如

3、何从这些受噪声干扰的为了从含噪信号di中还原出真实信号fi，可以利用信号信号中估计得到“纯净”的信号是建立系统高精度模型和实和噪声在小波变换下的不同的特性，通过对小波分解系数进现高性能控制的关键。行处理来达到信号和噪声分离的目的。在实际工程应用中，滤波器去噪是实际应用最广泛的一种方法，但时常在滤有用信号通常表现为低频信号或是一些比较平稳的信号，而除噪声的同时导致了有用信号的失真，它是从纯频域的角度噪声信号则通常表现为高频信号，所以我们可以先对含噪信来分析应该消除哪些频率范围内的噪声。1995年Donoho和号进行小波分解（如进行三层分解）：［

4、1］、［2］提出了小波收缩去噪的技术，他们研究的是在JohnstoneS=OA1+OD1=OA2+OD2+OD1叠加性高斯白噪声环境下检测出真实信号的情况，利用正交小波变换和高斯随机变量的性质对信号的小波分解系数做=OA3+OD3+OD2+OD1（1）阈值量化（threshokking），无失真的还原出真实信号。文章［3］其中OAi为分解的近似部分，ODi为分解的细节部分，i=1，中对Donoho-Johnstone的去噪方法做了总结推广，研究了在2，3，则噪声部分通常包含在OD1，OD2，OD3中，用门限阈值对高斯白噪声情况下选择不同的小波

5、变换形式（如非正交小波小波系数进行处理，重构信号即可达到去噪的目的［4］。变换）的去噪效果，并公式化了实际数据中的几种更复杂的总结去噪过程，可以分成以下三个步骤：噪声模型。1）对观测数据作小波分解变化［3］：本文对Donoho-Johnstone的收缩去噪方法在MATLAB环W0d=W0f+!·W0z（2）境下做了较为详尽的研究，总结了应用小波收缩去噪的一些实际经验。文章组织如下：第一节是小波去噪原理；第二节其中d表示观测数据向量f1，f2，⋯fN，f是真实信号向量f1，是小波去噪过程中两种门限阈值处理方式和几种阈值选取f2，⋯，fN，z是高

6、斯随机向量z1，z2，⋯，zN，其中用到了小波分方法；第三节是MATLAB中的小波去噪对比实验研究；最后解变换是线性变换的性质（小波分解变换本质上是一种积分是小结。变换）。2）对小波系数W0d作门限阈值处理（根据具体情况可以2小波去噪原理使用软阈值处理或硬阈值处理，而且可以选择不同的阈值形叠加性高斯白噪声是最常见的噪声模型，受到叠加性高式，这将在后面作详细讨论），比如选取最著名的阈值形式：斯白噪声“污染”的观测信号可以表示为：tN=!#2kogN（3）di=fi+!·zi，i=1，⋯，N（"）门限阈值处理可以表示为"tW0d，可以证明当N趋于

7、N无穷大时使用阈值形式（3）对小波系数作软阈值处理可以几基金项目：中科院优秀青年学者奖基金项目乎完全去除观测数据中的噪声［5］。收稿日期：2002-06-173）对处理过的小波系数作逆变换W0-1重构信号：—119—f!=W0-l!tW0d（4）值，计算公式为：N即可得到受污染采样信号去噪后的信号。0，N#32t={N（7）0.3936+0.l829!Iog23阈值的选取和阈值量化3）rigrsure：采用史坦（Stein）的无偏似然估计原理进行阈Donoho-Johnstone的小波收缩去噪方法的关键步骤是如值选择，首先得到一个给定阈值t的

8、风险估计，选择风险最小何选择阈值和如何进行门限阈值处理，在这我们将作较为详的阈值t!作为最终选择。细的讨论。4）heursure：选择启发式阈值，它是sgtwoIo