分形降噪處理的小波分析MATLAB工具箱教程

1、1. 進(jìn)入小波分析工具包模塊，選擇1維小波包分型降噪處理模塊wavelet toolbox main menufile window helpiwavelet 1-daveletpaccemcontinuous wavelet 1-dcomplex continuous wavelet 1-dspecialized tools 1dswt oenoising vodensity estimation vdregression estimation vowavelet coefficients selection bdfractional brownian generation 1-dtwodl

2、menslo 門 alwavelet 2-dwavelet packet 2-dspecialized tools 2-0swt denoising 2-0wavelet coefficients selection 2-0image fusionmultiple 1dmultisignal analysis 1-dmultivariate denoisingmultiscale pnnc. comp. analysiswavelet designnew wavelet for cwtdisplayiwavelet displaywavelet packet displayextensions

3、ignal extensionimage extensionclose2. 選擇file，然后選擇載入文件，載入.mat格式的文件，載入采集的 載荷時間歷程數(shù)據(jù)3. 載入數(shù)據(jù)后，會出現(xiàn)載入的信號時間歷程，如下圖所示，然后選擇分形的小波類型，己經(jīng)分形的小包等級，這里選擇db2小波，5 級水平為例進(jìn)行處理4. 點擊analyze將會出現(xiàn)如下圖所示的小波分形處理后的結(jié)果wavelet 1-d2 o-25cj o-24202-442024decomposition at level 5 s:a5*d5*d4*d3*<22*dl産,iki -l.iiti11 1v、垂重盡垂量1*-240 iii1

4、11詹耋1iiikiiirii2003004005006007008009001000匚v*xv*x- 1(7 y-c m.100l x yx-y«history<<rview axesdatav 0024)wavelet 肋 2 level 5statisticscompresshistogra. 一de-noiseanalyzefull decompositiat level5 disoiavmodeshow synthestzclose5. 點擊de-noise處理模塊進(jìn)行降噪處理，將會出現(xiàn)下圖所示的操作界面，可以根據(jù)需要選擇不同的降噪處理方式6. 選擇不同的降噪處

5、理方式，選擇閾值選擇的方式，也可以根據(jù)需 要自動調(diào)節(jié)閾值。select thresholding method fixed form thresholdo softhardselect noise structure unsealed white noiselev int select threshi nt. dependent threshold settingsde-noiseresiduals i7. 調(diào)節(jié)好降噪處理的參數(shù)后，點擊de-noise按鈕，就會出現(xiàn)最終的 小波包分析結(jié)果。最后點it visual de-noise signal將會出現(xiàn)如下圖所示的降噪處理后的結(jié)果。上述就是小包

6、分形降噪處理的簡單步驟，可以處理簡單的近似的小包處理，讀者可以根據(jù)需要進(jìn)一步擴展，下面給出上面工具包的原始matlab程序，方便讀者進(jìn)一步擴展和創(chuàng)新。clc;clear;xx=load(，實驗數(shù)據(jù).txt)； x=xx(:, 1廣； y=xx(:，2) ； x=x(cnd：-l: 1); y=y (end：-l:1); wname = sym6 ； lev=5;c，1二 wavedec(y, lev, wname); sigma = wnoisest (c, 1, 1);alpha=2;thrl = wbmpen(c, 1, sigma, alpha):thr2, nkeep = wdcbm(

7、c，1，alpha);xdl 二 wdencmpc gbl，c，1，wname, lev, thrl, s，1);xd2, cxd, lxd, perfo, perf2= wdencmpc lvd，c, 1，wname, lev, thr2, h)； thr, sorh, keepapp = ddencmpc den， wv，x);xd3 二 wdencmpc gbl，c，1, wname, lev, thr, s，1);plot(y); titlec 原始信號， fontsize，10);plot (xdl); title ( 使用penalty闊值降噪后的信號， font size，10); plot(xd2): title( 使用birge-massart閾值降噪后的信號， fontsizc，10); plot(xd