基于小波分析的信号去噪.doc

基于小波分析的信号去噪 一、实验目的 1、掌握小波分析的原理； 2、利用小波分析进行信号去噪，并编写Matlab程序。 二、实验内容 1、使用不同小波函数对信号去噪，比较消噪效果； 2、采取不同分解层数对信号去噪，比较消噪效果； 3、阈值设定方法对信号去噪的影响； 三、实验原理 小波分析方法是一种窗口大小（即窗口面积）固定但其形状可改变，时间窗和频率窗都可改变的时频局部化分析方法。即在低频部分具有较高的频率分辨率和较低的时间分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率，所以被誉为数学显微镜。正是这种特性，使小波变换具有对信号的自适应性。原则上讲，传统上使用傅里叶分析的地方，都可以用小波分析取代。小波分析优于傅里叶变换的地方是，它在时域和频域同时具有良好的局部化性质。 小波函数的定义：设为平方可积函数，即，若其傅里叶变换（是的傅里叶变换）满足 称为一个基本小波或母小波（Mother Wavelet），并称上式为小波函数的允许条件。与标准的傅立叶变换相比，小波分析中用到的小波函数不具有唯一性，对于一个时频分析问题，如何选者最佳的小波基函数是一个重要的问题。常用的小波函数有Haar小波、dbN小波、Morl小波、Mexh小波、Meyer小波等，不同的小波函数对应不同的尺度函数和性能。从下图中可以看出小波变换与傅立叶变换在时频窗口特性上有很大的不同，更显示了上述小波变换的特点。 图-1 小波变换的时频分析窗 小波变换的多分辨率分析实际上就是对一个频带信号进行低频分解，对每一步分解出来的低频部分在分解，使频率分辨率越来越高，其目的是构造一个理想的正交小波基。小波包分析实际上就是对与多分辨率分析没有分解的高频信号也进行逐层分解，进一步提高时频分辨率。小波分析地这些原理与特点与测控领域中的滤波原理非常相似，常常被用于信号噪声的消除。 运用小波分析进行一维信号消噪和识别信号中的发展趋势是小波分析的一个重要应用之一。在实际的工程应用中，所分析的信号可能包含噪声和一些不应有的趋势项，对这种信号，首先需要作信号的预处理，将信号的噪声和趋势项去处，提取有用信号。 由于小波变换具有对信号的自适应性，可以识别和处理信号在任意局部所包含的频率成分，能识别信号中的奇异信号、非平稳随机信号。 下面以小波分析在信号消噪方面的应用，介绍小波分析在信号的应用处理的数据为本教研室正在进行的捷联惯导当中加速度传感器的输出。 对于小波消噪的算法主要分为三步： (1)用小波变换或小波包变换将含噪数据变换到小波域； (2)确定一个阈值标准并用作用到小波域； (3)重构数据，得到消噪后的信号；小波函数种类繁多，不同的小波函数适用不同的情况，下面以常用的五种小波函数Daubechies小波、Coiflets小波、Symlets小波、正交小波、双正交小波对本实验数据处理，选择最合适的小波函数。 clear all; fid = fopen(动态测试一周1.txt); A = fscanf(fid,%f,[1 inf]); meg_sensor=A;meg_sensor=meg_sensor(1:40000);lengtha=size(meg_sensor); %%%%%%%%%%%%%%%%%%%%%%%% 小波分解 [c_db,l_db]=wavedec(meg_sensor,6,db5); [c_coif,l_coif]=wavedec(meg_sensor,6,coif5); [c_sym,l_sym]=wavedec(meg_sensor,6,sym5); [c_bior,l_bior]=wavedec(meg_sensor,6,bior5.5); [c_rbio,l_rbio]=wavedec(meg_sensor,6,rbio5.5); %%%%%%%%%%%%%%%%%%%%%%% 强制消噪 meg_sensor_db=wrcoef(a,c_db,l_db,db5,6); meg_sensor_coif=wrcoef(a,c_coif,l_coif,coif5,6); meg_sensor_sym=wrcoef(a,c_sym,l_sym,sym5,6); meg_sensor_bior=wrcoef(a,c_bior,l_bior,bior5.5,6); meg_sensor_rbio=wrcoef(a,c_bior,l_bior,rbio5.5,6); %%%%%%%%%%%%%%%%%%%%%% 结果显示 %%%%%%%%%%%%%%%%%%%%%%% 原始信号和强制滤波之后的信号比较 figure(1);k=1:1:lengtha(2); subplot