matlab 小波启发式阈值滤波,小波分析matlab工具箱简介.doc

小波分析matlab工具箱简介.doc

MATLAB的小波分析一、小波分析用于降噪的基本过程1、分解过程选定一种小波，对信号进行N层分解；2、作用阈值过程对分解得到的各层系数选择一个阈值，并对细节系数进行软阈值处理；3、重建过程降处理后的系数通过小波重建恢复原始信号；二、基本降噪模型函数一维离散小波分解命令DWTCACDDWTX,’WNAME’使用小波’WNAME’对型号X进行单层分解，求得的近似系数存放于数组CA中，细节系数存放在数组CD中；CACDDWTX,’WNAME’,’MODE’,MODE利用MODE方式进行扩展CACDDWTX,LO_D,HI_D利用指定滤波器进行小波分解WANEDECC,LWAVEDECX,N,’WNAME’使用WNAME的小波进行N层分解，C为层数，L为各层系数IDWTXIDWTCA,CD,’WNAME’利用小波WNAME把近似系数CA和CD重建为上一层近似系数XXIDWTCA,CD,’WNAME’,L重建至L层WAVERECXWAVEREC(C，L，‘WNAME‘)重建为原始信号WRCOEFXWRCOEF‘TYPE’,C,L,’WNAME’,N通过分解系数重构指定的数，TYPE为A或者DXWRCOEF‘TYPE’,C,L,’WNAME’把分解系数重建至最高层UPCOEFYUPCOEF(O,X,’WNAME’,N)用适当的滤波器作用在X上N次，求得重建系数Y，O为A表示低通滤波器，D表示高通滤波器DETCOFEDDETCOEF(C，L，N)从分解系数中提取第N层近似系数DDETCOEF(C，L，N)提取至最后一层APPCOEFAAPPCOEFC,L,’WNAME’,N用小波从分解系数中提取第N层系数WNOISESTSTDCWOISEST(C,L,S)根据传入的小波分解系数C,L对S中标识的小波层数求得其标准差，作为对噪声强度的估计；DDENCMPTHR,SORH,KEEPAPP,CRITDDENCMPIN1,IN2,X根据传入的参数IN1和IN2所指定的方式，对输入信号X求得其降噪或压缩的各级阈值。IN1DEN(用于降噪)；CMP(用于压缩)。IN2WV(使用小波变换)；WP(使用小波包变换)。THR求得的阈值；SORH‘S’软阈值；‘H’硬阈值。KEEPAPP保留的近似系数的层数。CITR为小波包变换用的熵准则WBMPENTHRWBNPENC，L,SIGMA,ALPHA根据传入分解系数用PENALTY策略确定各层阈值，并返回到数组THR中，WDCBMTHR,NKEEPWDCBM(C，L，ALPHA)用BIRGENASSART策略确定各层阈值并返回保留的系数所在层数NKEEP。THSELECTTHSTHSELET(Y，TPTR)使用指定的阈值选取方式TPTR选择用于对信号Y降噪阈值，返回到THS中，TPTRRIGSURE严格按照SURE阈值选取；TPTRSQTWOLOG对数长度阈值选择；TPTRHEUSURE(启发式SURE阈值选择)；TPTRMINIMAXI(最小极大方差阈值选择)注1、基于STEIN无偏似然估计(SURE)的软阈值估计(RIGSURE)。对于给定的阈值T，得到它的似然估计，然后将似然函数最小化，得到所需要的阈值。2、长度对数阈值(SQTWOLONG)从得到最小极大方差的阈值T乘以一个系数得到的阈值；3、启发式SURE阈值(HEURSURE)前两种方式的综合形式，因为基于SURE产生的阈值在高信号噪声比的情况下抑制噪声效果不明显，这种方法利用启发函数自动在前两种阈值中选择一个；最小极大方差阈值。4、最小极大方差阈值(MINIMAXI)是的选取的阈值产生最小的极大方差。通过统计学上聚集起的构造方法得到。以为降噪后的信号可以看成与未知回归函数的估计式相似，所以这种方法通过求得未知回归函数与原信号方差在最坏情况下的最小值来获得阈值WTHRESHYWTHRESHX,SORH,T对信号X作用阈值T作用方式为SORH‘H’,YSIGNXG|X|T硬阈值，可能会出现断点；SORH‘S’,YX1{|X|T}软阈值，不会出现断点，比较光滑。WDENXD,CXD,LXDWDENX,TPTR,SORH,SCAL,N,’WNAME’根据传入的参数，对信号X降噪，结果返回到XD中，CXD，LXD是降噪结果XD的分解系数结构。TPTR与以上选取参数和意义一致(阈值选取规则)。SORH，作用软阈值或者硬阈值。WDENCMPXC,CXC,LXC,PERF0,PERFL2WDENCMP‘GBL’,X,’WNAME’,N,THR,SORH,KEEPAPP全局阈值的降噪命令，通过小波函数WNAME指定的小波，对分解系数作用有参数THR，SORH指定的阈值，保留KEEPAPP以上的所欲近似系数。XC降噪或压缩后的信号；CXC，LXC为降噪或压缩后的小波分解系数结构。通过抑制细节函数实现降噪XNOSIDOPP赋值信号函数C,LWAVEDECX,5,’DB4’,\CAWRCOEF‘A’,C,L,’DB4’,5重建第五层近似信号INDEXL21L7从林中读取所有细节系数所在下标C1CC1INDEXCINDEX/3对细节系数进行抑制X2WAVERECC1,L,’DB4’利用抑制后的系数进行重建小波PER1NORMCA/NORMX降噪后的CA信号在原信号的能量成分PER2NORMX2/NORMXMATLAB缺省的降噪命令XNOISDOPPTHR,SORH,KEEPAPPDDENCMP‘DEN’,’WV’,X根据信号强度，给出全局阈值XDWDENCMP‘GBL’,X,’SYM4’,2,THR,SORH,KEEPAPP根据全局阈值进行降噪C,LWAVEDECX,4,’SYM4’THR,NKEEPWDCNMC,L,2得到每个层次的分层阈值；XD1,CXD,LXD,PERF0,PERFL2WDENCMP‘LVD’,C,L,WNAME,LEV,THR1,’S’根据分层阈值使用软阈值方法对信号进行降噪FUNCTIONPHASE_ALTER,PHASE_DOTWAVELET_DENOINGPHASE_DATA,PHASE_DLEV5WNAME DB4 C,LWAVEDECPHASE_D,LEV,WNAMESIGMAWNOISESTC,L,1ALPHA2THR1WBMPENC,L,SIGMA,ALPHAÀÛÓÃPENALTY¼ÔÒȶ¨½ΜÔËÃÐÖΜPHASE_DO1WDENCMP GBL ,C,L,WNAME,LEV,THR1, S ,1THR2,NKEEPWDCBMC,L,ALPHAÀÛÓÃBIRGEMASSART¼ÔÒȶ¨½ΜÔËÃÐÖΜPHASE_DO2,CXD,LXD,PERF0,PERFL2WDENCMP LVD ,C,L,WNAME,LEV,THR2, H THR,SORH,KEEPAPPDDENCMP DEN , WV ,XȱʡÃÐÖΜ½ØÐнΜÔË´¦ÀÍPHASE_DO3WDENCMP GBL ,C,L,WNAME,LEV,THR, S ,1FORII124PHASE_DATAIIPHASEIIPHASE_DII1ENDEND