小波阈值去噪算法设计及其应用

1、-. z.北方民族大学学士学位论文 论文题目: 小波阈值去噪算法的设计及其应用 院(部)名 称: 数学与信息科学学院 学 生 姓 名: 黄 慧 东 专 业: 信息与计算科学 学 号: 20100433 指导教师: 黄 永 东 论文提交时间: 2013年5月14日 论文辩论时间:学位授予时间:北方民族大学教务处制-. z小波阈值去噪算法的设计及其应用摘要本文主要阐述了小波阈值去噪算法的设计及其应用.对小波进展了初步的介绍，小波分析是分析未经过任何处理的信号所含有的不同的性质，进而用于图像处理、小波滤波、数据隐藏等.比方声音信号频率的上下，发声时间的长短、振幅、旋律等各个方面.从平稳的波形之中发现

2、突变的尖峰.小波分析是依照各种小波基函数对分解原始信号的一种分析方法. 介绍了小波滤波并列举了几种常用的小波滤波算法.时至今日，小波滤波成为了一种新的滤波思路，其功能除了去噪、降噪以外，还兼有平滑、锐化和保存信号特征的功能.第三章则较为详细介绍了小波阈值去噪算法并进展了算法设计，最后还给出了小波阈值去噪算法的应用实例.小波阈值去噪就是将经过小波分解后的信号通过选取适当的阈值过滤掉带噪信号，再用小波逆变换进展小波重构.关键字：小波分析，小波变换，小波滤波，小波阈值去噪.designofwaveletthresholddenoisingalgorithmanditsapplicationabstr

3、actthisarticlefocusesonthewaveletthresholdingalgorithmdesignanditsapplication.thefirstchapterintroducesthewaveletconductedpreliminary, waveletanalysis isananalysisofvariouschangesinthecharacteristicsoftheoriginalsignal , andfurtherusedindatapression, noiseremoval , featureselection. fore*amplesingin

4、gsignal: thetrebleorbass, soundduration , undulatingmelodyandsoon. waveletanalysisistheuseofavarietyof waveletfunction on rawsignal deposition.thesecondchapterintroducesthewaveletfilteringandlistsseveralmonlyusedwaveletfilteringalgorithms. today, waveletfilteringhasbeeanewfilterideas, inadditiontoit

5、sfunctionnoising , noisereduction , italsobinessmooth, sharpenandretainthefunctionofthesignalcharacteristics .thethirdchapterisamoredetaileddescriptionofthewaveletthresholdingalgorithmandalgorithmdesign , andfinallygivesthewaveletthresholdingalgorithme*amples . waveletthresholdingisbasedontheeffecti

6、vesignalandnoisehavedifferentpropertiesatdifferentdepositionscale , constructedusingmathematicaltoolsappropriatethreshold , andthetargetsignalwaveletcoefficientsthresholdingkeywords: waveletanalysis, wavelettransform, waveletfiltering, waveletthresholding . TOC o 1-3 h z u 目录 HYPERLINK l _Toc9289 第一

7、章 小波初步介绍1 HYPERLINK l _Toc22826 1.1小波分析的开展1 HYPERLINK l _Toc18766 1.2小波分析的应用1 HYPERLINK l _Toc31857 第二章 小波滤波算法5 HYPERLINK l _Toc30698 2.1 小波滤波5 HYPERLINK l _Toc5681 2.2 常用的小波滤波算法6 HYPERLINK l _Toc1666 第三章 小波阈值去噪算法8 HYPERLINK l _Toc18633 3.1小波阈值去噪算法8 HYPERLINK l _Toc28480 3.2小波阈值去噪的具体步骤:8 HYPERLINK l

8、 _Toc8228 3.3 阈值函数9 HYPERLINK l _Toc7185 3.4小波阈值去噪算法的算法设计9 HYPERLINK l _Toc18995 3.5小波阈值去噪算法实例： PAGEREF _Toc1899512 HYPERLINK l _Toc28913 总结 PAGEREF _Toc28913 18-. z第一章 小波初步介绍1.1小波分析的开展什么是小波？小波就是具有衰减性和波动性，并且它波动时振幅有正有负的波.最早提出小波分析的是法国人j.morlet，他擅长于信号处理.在1986年，数学家y.meyer在一个意外的情况下得到了一个较为标准意义上的的小波基，在那之后，

9、他和s.mallat一起创立用多尺度分析来构造小波基的方法，从那以后，小波分析迈上了不断开展的道路.在这段时期，当属小波十讲真正使得小波分析得到了广泛的推广,到了近年来，小波域阈值滤波方法被donoho等人所建立起来，取得了非常巨大的成功.小波分析与fourier变换和窗口fourier变换两者作比拟，小波得益于其时间和频率的局域变换的特点，从而能有效并且准确的从信号之中提取讯息，通过伸缩和平移等功能，对函数或信号进展多尺度细化分析multiscaleanalysis，许多复杂问题不能通过傅里叶变换解决，但小波分析却可以办到，也是由于这点，小波变换被高度称赞，有的人们还称其为数学显微镜，成为了

10、调和分析开展历史上一座耀眼的里程碑.小波分析是一个正在蓬勃开展的新领域，而且经由无数学者数年来的不断探寻钻研，小波分析根底的数学形式系统已然创立好了，理论根底得到了进一步.如今它在信号分析、声音处理、数据隐藏、图像水印、地震勘察、水文分析等领域中的表现非凡，取得了成功，得到了广泛认可. 小波分析 是分析未经过任何处理的信号所含有的不同的性质，进而用于图像处理、小波滤波、数据隐藏等.比方声音信号频率的上下，发声时间的长短、振幅、旋律等各个方面.从平稳的波形发现突变的尖峰.小波分析是依照各种小波基函数对分解原始信号的一种分析方法. 1.2小波分析的应用在现实的应用场景中，人们所接收到的信号往往夹杂

11、着各种各样不同的噪声，这些噪声对信号形成了一定干扰，信号的有效性不高，有时甚至能使信号失效，因此信号去噪是信号处理过程中的重中之重.然而传统的根据信号的频谱分析进展带通或低通滤波方法虽然比拟简单易操作，却对有效频带的噪声束手无策，也无法实现高分辨率.小波分析的出现正好提供了一种新的分析方法.它具有时频分析，多分辨率分析等长处，这些长处与人类的视觉系统的特征一样，因此小波分析的应用领域非常广泛.图像处理图像是二维的信号，其小波变换就相当于二次一维信号的小波变换：第一次的一维信号的小波变换就相当于图像的行变换.第二次的一维信号的小波变换就相当于图像的列变换.小波变换用于图像压缩有良好的效果，已经形

12、成图像压缩的标准，如jpeg2000.图像压缩小波变换用于图像压缩：采用小波进展压缩.在经过小波变换之后，统计特性变得更好，行和列之间的相关性得以消失.有损压缩：依照视觉成像机制，把各个不同分辨率的小波系数进展比特匹配.接着转换到一维作熵编码，算术编码和霍夫曼编码等. 无损压缩：采取整数小波变换，舍入误差并不存在.缺乏的是不能进展比特分配.小波变换用于图像压缩图 1小波变换用于图像压缩数据隐藏左边为原始图像1024*768 右边为信息隐藏后的图像1024*768此为隐藏的5320*280交通图像，可完全恢复数字水印数字水印是一种新的加密方法，有效弥补了传统方法的缺乏，这些年来人们将其应用了到方

13、方面面.数字水印技术就是通过特定的算法将所需要隐藏的水印信息直接嵌入到多媒体容之中，但原容的价值和使用并不会因此而受到影响，并且人的知觉系统发觉不到水印的存在，只有用指定的的检测器或阅读器才能够为人们所看到. 银行提款码嵌入指纹，网上进展身份认证:左边为指纹原始图像，右为嵌入水印的图像小波滤波算法2.1 小波滤波小波去噪是指在小波域上处理小波系数以抑制或去除噪声.在20世纪90年代初，小波去噪以及降噪开场出现在了文献中，随着小波去噪方法的多样化和研究的不断深入，近年来更是开展迅速，时至今日，小波滤波成为了一种新的滤波思路，其功能除了去噪、降噪以外，还兼有平滑、锐化和保存信号特征的功能.小波滤波

14、方法有很多，首先是mallat等人提出了基于小波变换模极大值原理的滤波方法，后由*u提出了基于信号尺度间相关性的空域相关滤波算法，随后donoho等人另辟蹊径，又提出了小波阈值滤波方法.2.1.1小波滤波原理及模型： 首先，我们来了解下什么事小波变换平方可积的函数*t的小波变换定义为：如果满足容许性条件则w*a,b= 其中.称为*(t)的小波变换.其中为小波基函数，a为尺度因子，b为平移因子，其作用是确定对*(t)分析的时间位置，也即时间中心.假设参数a和b都是连续变量，则称其为连续小波变换，此时连续小波变换存在逆变换，形如：利用此公式可以恢复原信号. 选择适当的母小波，可以使得小波基函数及其

15、傅丽叶变换都能够有着较为不错的局部性，从这里可以得出来的结论是小波分析是时-频分析7.2.1.2根本原理：小波变换具有以下性质：时频局部化.小波基函数在时、频领域都具有有限或近似有限的定义域，所以在其伸缩平移后依然在时、频域是局部性的；多分辨率.原始信息分辨率最高，逐层分解后分辨率降低，这样的分解方法可以很好的表达信号的非平稳特征；选基灵活.小波变换能够灵活选择变换基，可以根据不同应用场合选择恰当的小波函数；窗口固定但形状可变.尺度因子的伸缩变化能引起小波奇函数的窗口变化，时间窗口随a的增大而增大，频域窗口则相反； 以上特点就是小波去噪方法之所以能够成功的关键所在.可以看出，小波变换的时频局部

16、化特性和多分辨率性使得小波滤波方法能在去噪的很好地保存信号的非平稳特征.2.2 常用的小波滤波算法迄今为止，常用的小波滤波方法有三类：基于小波变换模极大值原理的滤波方法2基于信号尺度空间相关性的空域相关滤波算法；3根据幅值较大的小波系数由重要的信号产生这一个假设的阈值滤波方法.2.2.1空域相关滤波最先提出利用尺度空间相关性进展信号滤波思想的是witkin，即从尺度上由粗到细搜索经过子带分解后的信号边缘，再得到去噪后的信号，随后，*u在1994年进一步提出了空域相关滤波法.经过小波变换后得到的信号的小波系数即使在不同的尺度间里也具有很强的相关性，噪声的小波系数却是弱相关的或者是不相关的.将相邻

17、尺度上的系数做乘法运算得到空域相关的系数，假设*点经过归一化后的相关系数的幅值大于其小波系数的幅值，我们则认为该点为信号的突变点，将该点提取出来.对提取出的小波系数进展逆变换后，就能得到滤波信号.*u 在提出空域相关滤波方法的时候，并没有给出来噪声方差的估计以及终止迭代过程的阈值.pan 在计算过各分解尺度上的噪声方差之后，提出了噪声能量阈值的理论计算方法，还给出终止迭代过程的准则，首次提出了理论相对完整的空域相关滤波算法.并对此算法进展了改良，得到了具有令人满意滤波效果的自适应性空域相关滤波算法.这里值得注意的是，相关系数的定义将会直接影响滤波的结果7.2.2.2模极大值重构滤波小波系数的模

18、极大值表达了小波的奇异性，这种方法则是利用了这个性质来重构信号，即通过处理各个尺度上小波系数的模极大值来改变小波的奇异性，而描述信号和噪声在多尺度空间上不同特性的lipschitz指数就可以度量函数的奇异性.mallat认为噪声的 lipschitz指数小于零，根据这点可以去除掉噪声，留下模极大值点处的值，再用交替投影法alternative projection,ap来重构信号.基于模极大值的重构滤波算法已经有了很好的理论根底，滤波性能也相对稳定，对噪声并无特殊的要求，只要信号与噪声的 lipschitz 指数不一样，就可以利用此方法来滤波.其缺点就是：由模极大值重构信号的算法不但复杂，而且

19、计算的难度大，在实际的应用中的滤波效果并不佳.2.2.3小波域阈值滤波小波域阈值滤波算法最先由 donoho 等人提出，其主要思想就是把小波分解后各层小波系数的模值与一特定的阈值进展比拟，小于阈值的视为由噪声产生，设为零，小波系数大于阈值时将保存或增强，然后用信号阈值函数处理小波系数，最后，再进展小波逆变换来重构信号，从而到达消除噪声的目的.该方法较为重要的是对阈值和阈值函数的选择，这两个要素决定了去噪后的信号的准确度.第三章 小波阈值去噪算法3.1小波阈值去噪算法小波阈值去噪即是依照有效信号和噪声在不同的分解尺度上有着不一样的性质 ，使用数学的工具构造出适宜的阈值，并且对目标信号的小波系数进

20、展相应的阈值处理.其过程，就是利用小波分解将信号按尺度进展多层分解，然后根据以往的经历和学到的知识设置一个阈值，小于的小波系数认为主要由噪声引起，大于的小波系数认为是由有效信号引起，接着消除每个尺度之中由噪声所引发的小波系数，根据情况保存或是加强由有效的信号引发的小波系数，最后，再进展小波逆变换将经过阈值处理后的小波系数重构得出去噪后的信号5.大致流程图如下：含噪信号小波分解阈值去噪重构信号小波逆变换3.2小波阈值去噪的具体步骤:利用小波分解处理带噪信号.确定n层的小波分解，然后将信号分解为n层.阈值处理.采取恰当的阈值和阈值函数对每一层的高频系数进展处理，将噪声信号过滤出来.小波的重构.依照

21、小波分解的第n层低频的系数和每一层的高频的系数，来进展信号重构. 在以上描述的3个步骤之中，最为关键重要的便是如何应用恰当的方法选取阈值和阈值函数，信号中噪声去除的好不好，全靠这两个点.阈值的选取方式大致分为两种，全局阈值和局部适应阈值.全局阈值适用于每一层的小波系数，局部阈值是为各层小波系数单独设置的，因而选取更加灵活.1 全局统一阈值dj统一阈值：.其中的是噪声的标准方差,n为信号的长度.该阈值是在高斯模型下针对多维独立的正态变量联合分布所得出的5;2基于零均值的正态分布的置信区间阈值:=3，其中为噪声的标准方差.该阈值是假设0均值的正态分布变量必定会处在区间3，3之，因此我们通常都会认为

22、绝对值大于3的系数是由信号产生的，否则就是由噪声所产生的5.3.3 阈值函数 主要有这么两种阈值函数选取方式，硬阈值法和软阈值法硬阈值法：其中为选取的阈值下同.即将信号的绝对值与阈值比照，信号的绝对值小于等于阈值则设其为0，假设信号的绝对值大于所选阈值，则信号将保持不变.软阈值法：其中这里将信号的绝对值与设置的阈值作比拟，信号的绝对值小于等于阈值则设其为0，假设信号的绝对值大于阈值，则将其设置为信号绝对值与阈值的差.比照这两种阈值函数的选取.硬阈值方法能很好的反映信号的非平稳特征，但是它是不连续的，会导致重构之后的信号有所波动，所以去噪后的信号残留的噪声会很明显；软阈值法具有较好的连续性，所以

23、去噪结果更平滑，然则它的导数不具有连续性，其预计小波系数与初始小波系数值相比有一个固定不变的差值，而且对大于阈值的系数进展定值压缩与噪声随着小波系数增大而减少的事实不符合，所以软阈法处理后的小波系数较好处理，但可能会使图像边缘模糊出现失真.3.4小波阈值去噪算法的算法设计 小波阈值算法以其简单,运算量小的特点,得以大量的运用在信号去噪方面.美中缺乏之处有二:分解层数越大，白噪声信号的模极大值会跟着下降,这就会导致在用全局阈值处理不同的分解层时会在低频系数中错误的去掉太多的有效信息，在高频系数反而使大量的噪声残留.在进展阈值处理时,硬阈值处理能较好的保存有效信号的局部特征,但由于硬阈值法的不具有

24、连续性,去除噪声后的结果不能稳定；软阈值处理的信号尽管在去除噪声后变得更为光滑,但估计的小波系数与原小波系数之间存在固定的偏差.本文采用的改良方法如下:3.4.1各层采用不同的阈值1其中k为一个常数，j为分解层数，a为lipschitz指数，为第j层的小波系数.lipschitz指数能衡量信号的奇异性,即是说a越小,信号突起越多.对于一般信号a0,即有用信号对应小波系数随分解尺度j的增大而变大.而白噪声的lipschitz指数则为负值,即噪声对应的小波系数随分解尺度j的增大而减小.而白噪声的lip指数满足：2有1式和2式可知：其中为噪声对应的第j层小波系数,由以上式子不难看出噪声说对应的第j+

25、1层小波系数的最大值比第j层小波系数的最大值要小,所以,此文在经过阈值处理时每层系数都会使用不一样的阈值,选取恰当的阈值计算方法就可以得到第一层阈值,往下的每层阈值为上一层阈值的倍，即.3.4.2采用新的阈值函数由于软阈值和硬阈值的缺乏之处明显,在这里选用一种新的阈值函数.该阈值函数要尽可能的结合软、硬阈值的优点，即是说，既要保证新阈值函数的连续性,又要尽可能防止软阈值函数中的固定偏差.拟阈值函数如下:其中，式中n是一个为正的实数，是经过处理前的初始小波系数,是经过处理后的小波系数.该阈值函数具有连续性,并且当的时时候，函数在高阶上可导,这时增大，t的值就会随之减小,这样，处理后的小波系数也不

26、会改变.当n的取值非常大时,新阈值函数类似于软阈值函数;当n接近于0的时候,新阈值函数则类似于硬阈值函数.这也即是说新的函数能够通过调节正实数n的大小从而得以在软阈值函数和硬阈值函数之间切换,比两种经典的阈值函数灵活了许多.图1 硬阈值函数、软阈值函数及改良法制函数的曲线 图2 block信号及去噪结果图1是硬阈值函数、软阈值函数和新阈值函数曲线,阈值=1，n=0.5.从图中能够看出硬阈值函数的不连续性和软阈值函数存在的固定偏差.而新的阈值函数在阈值点处连续,而且在这里的渐近线是硬阈值曲线.3.4.3实验仿真 在这里使用matlab软件首先把高斯白噪声加到block信号中去,接着再使用dbl小

27、波将初始信号利用小波分解分为5层.通过最小极大方差法得到阈值,分别使用一般的小波阈值去噪方法和本文里设计的小波阈值去噪方法处理小波系数.仿真的结果如图2所示. 从图中可以看出,硬阈值法不能完全去除噪声信号,处理后的信号包含了一些锯齿,不光滑.软阈值已经失去了很多有效的信号,使信号在奇异点变化缓慢.而本文设计的小波阈值去噪方法恰恰防止了软、硬阈值两者的缺乏之处,既成功去除掉了更多的噪声,又还很好的保存了更多的有效信号.3.5小波阈值去噪算法实例：心电信号去噪3.5.1心电信号的噪声特点心电信号中混杂着的噪声是不平稳的，在其采集过程中会受到来自多类型噪声的干扰，其中以下三种噪声影响较大:1当来自电

28、源的磁场作用于心电图机时，会与周围人体之间形成一个环形电路;2人体肌肉紧也会产生电路而造成影响;3人体呼吸系统也会产生噪声影响. 这些额外的磁场与电路所产生的噪声会严重地干扰心电信号，心电信号会因为这些干扰而变得畸形，还会影响之后对信号的处理,特别是计算时机因此而得出错误的判断，导致不良的后果.可以看出，准确有效的去除心电信号中的噪声是非常有必要的.(心电信号)3.5.2小波分析的去噪原理 在现实的应用中，信号一般波动较大，不具有线性还有不少的奇异点.一个原始的一维信号用公式可表示为：其中，f(t)为真实的信号，s(t)为含噪的信号，e(t)为噪声，为噪声的标准偏差.有效信号通常处于低频地带，

29、而高频信号一般是噪声信号. 在使用小波分析分解原始信号之后，在高频小波系数中包含了大量的噪声信号，有用信号往往集中在数量较少的低频小波系数；而噪声对应的小波系数则相反. 综合上面的表达，我们能够利用阈值去噪法来处理小波系数，然后应用小波逆变换进展信号重构，这样就能够消除信号中的噪声.3.5.3小波分解的示意图：小波分解系数示意图3.5.4一维信号利用小波进展除噪的步骤预处理小波变换多尺度分解各尺度小波系数除噪小波逆变换重构信号除噪后的信号含噪信号1.利用小波变换来去噪的流程示意图：在小波变换域上进展阈值处理：多层小波分解阀值操作多层小波重构3.5.5阈值函数和阈值的选取 阈值函数可以采用硬阈值

30、法或者是软阈值法，阈值则可选择固定阈值、stein无偏似然估计阈值、启发式阈值和极大极小阈值中的一个.3.5.6小波函数的选择小波变换并不只是由正弦函数所唯一确定的，小波基分为了很多种，不同的小波对应不同的信号去噪.对于特定的信号来说，假设小波的选择不当，去噪的结果很有可能相差很远，甚至还有可能丧失一些有用的信息.面对如此之多小波，必须经过大量的仿真研究结果来筛选出能出色处理心电信号的小波.依照不少文献记载，对心电去噪最适宜的函数是b样条函数： b样条函数是一种非紧支撑正交的对称小波，光滑性比拟高，频率特性不错，分频能力也强，还有频带相干较小的特性.在信号处理之中小波的作用就相当于带通滤波器，

31、并且它的对称性和反对称性分别等价为线性相位和广义线性相位.我们知道的是，当一个带通滤波器是非线性相位或广义线性相位的时候，它将会使通过它的信号发生生畸变.本文从理论的和实践的两个角度来看，紧支集的小波则是最有力的一个.b样条函数是一种根本的样条函数，并且它的支撑区域最小据上所述，b样条小波确实是一个很适宜的选择.3.5.7小波去噪效果的评价式中的yi表示标准的原始信号, *i 则表示经处理后的估计的信号.其中，信噪比mse越大越好, 多尺度熵snr则越小越好.3.5.8小波去噪程序 输出评价结果及去噪后的心电信号去噪效果的评价 根据选取的阈值去噪及重构求各层的阈值 提取各尺度小波系数 对被选的心电信号进展小波分解 选取其中一个导联的数据 调用心电数据开场matlab程序：对100.dat from mit-bih-db的单导联数据进展去噪处理c l=wavedec(e,3,db5);% 从c中提取尺度3下的近似小波系数ca3=appcoef(c,l,db5,3);程序运行结果tr =0.0708 0.3636 1.0552snr =121.6977mse =0.00243.5.9 仿真分析 通过仿真结果来看，对应小波函数和阈值的不同选择，会导致去噪的效果也大相径庭.选择db5小波，默认阈值对心电去噪处理效果较好.总结小波分析开展至今已经