37小波分析支噪与压缩.docVIP

３７小波分析支噪与压缩 知识背景 小波分析中最重要的两个运用是信号去噪和信号压缩，使用小波分解可以将原信号分解为一系列的近似分量和细节分量，信号的噪声集中表现在细节分量。使用一定的与阈值处理细节分量，再经过小波重构后可以得到平滑的信号。信号的小波压缩与小波去噪在本质上是相同的，即经过小波重构后可以得到平滑的信号。小波与小波去噪的主要区别在于阈值的选择算法不同，可以使用全局阈值，也可以使用自适应值。 主要内容 本练习主要讲述小波分析中信号去噪和信号压缩的基本思想，以及进行相应处理的相关函数，主要通过函数的用法实例和MATLAB自带的实例来揭示小波分析中信号的处理方法。由于小波分析中所用到的函数众多，本练习只是给出常用和比较重要的分析函数和工具包。对于具体的需要，读者可以通过MATLAB的帮助文件来获得详细资料。 练习过程 （一）小波去噪和压缩阈值选取方案可以通过函数ddencmp,其用法为： 〔THR,SORH,KEEPAPP,CRIT〕=ddencmp(dorc,worwp,x) 可以有三种使用形式: １　[THR,SORH,KEEPAPP,CRIT]=DDENCMP(IN1,IN2,X) [THR,SORH,KEEPAPP]=DDINCMP(IN1,’wv’,X) [THR,SORH,KEEPAPP,CRIT]=DDINCMP(IN1,’wp’,X) 根据信号x和参数IN1和IN2自动生成小波去噪和压缩的阈值选取方案。其中IN1可选值为“den”为信号去噪，”cmp”为信号压缩。IN2可选值为“wv”使用小波分解，取值为“cmp”时使用小波包分解。而THR为函数选择的阈练习３７　小波去噪与压缩 值，SORH为ddencmp函数选择的阈值使用方式，取值为“s”时为软阈值，取值为“h”时为硬阈值。而KEEPAPP决定了是否对近似分量进行阈值处理，可选１或者０。CRIT为使用小波包分解时选取的熵函数类型。 例： init=21564148656;randn(‘seed’,init); x=randn(1,800); [thr,sorh,keepapp]=ddencmp(‘den’,’wv’,x) [thr,sorh,keepapp]=ddencmp(‘cmp’,’wv’,x) [thr,sorh,keepapp]=ddencmp(‘den’,’wp’,x) [thr,sorh,keepapp]=ddencmp(‘cmp’,’wp’,x) 在不同的信号处理方式下所得结果如下： thr= 3.6851752973086 sorh= s keepapp= 1 thr= 0. 67980588564722 sorh= h keepapp= 1 thr= 4.23105962519793 sorh= h keepapp= 1 thr= 0. 67980588564722 sorh= h keepapp= 1 （二）当选取用于小波去噪的阈值时，可以使用函数thselect,其用法为： thr=thsrlect(x,tptr) 其中thr为返回用于小波去噪的阈值，x为有噪信号，参数tptr则指定了阈值选取算法，可选值分别为：rigrsure, heursure, sptwolog, minimaxi.举例如下： init=21564148656;randn(‘seed’,init); x=randn(1,800); thr=thselect(x,,’rigrsure’) thr=thselect(x.,’sqtwolog’) thr=thselect(x,’minimaxi’) 在参数取值不同的情况下所得结果分别为： thr= 1.99908916578119 thr= 3.65639487136385 thr= 3.65639487136385 thr= 2.15746129710980 (三)一维信号的小波去噪处理的函数主要为wden,函数开式为： [xd,cxd,lxd]=wden(in1,in2,in3,in4,in5,in6,in7) 主要运用方式为：[XD,CXD,LXD]=WDEN(X,TPTR,SORH,SCAL,N,’wname’) [XD,CXD,LXD]=WDEN(C,LTPTR,SORH,SCAL,N,’wname’) 函数处理后返回去噪的信号xd,及其小波分解结构〔cxd,lxd〕,X为有噪信号。参数tptr指定阈值选取算法，其可选值为：rigrsure, heursure, sptwoolog, minimaxi. sorh为函数选择的阈值使用方式，当取值为‘s’时，为软阈值；取值为‘h’时，为硬阈值。Scal可以取的参数值为：one, sln, mln;