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

1、精選優(yōu)質(zhì)文檔-傾情為你奉上1. 進(jìn)入小波分析工具包模塊，選擇1維小波包分型降噪處理模塊2. 選擇file，然后選擇載入文件，載入.mat格式的文件，載入采集的載荷時間歷程數(shù)據(jù)3. 載入數(shù)據(jù)后，會出現(xiàn)載入的信號時間歷程，如下圖所示，然后選擇分形的小波類型，已經(jīng)分形的小包等級，這里選擇db2小波，5級水平為例進(jìn)行處理4. 點(diǎn)擊analyze將會出現(xiàn)如下圖所示的小波分形處理后的結(jié)果5. 點(diǎn)擊De-noise處理模塊進(jìn)行降噪處理，將會出現(xiàn)下圖所示的操作界面， 可以根據(jù)需要選擇不同的降噪處理方式6. 選擇不同的降噪處理方式，選擇閾值選擇的方式，也可以根據(jù)需要自動調(diào)節(jié)閾值。7. 調(diào)節(jié)好降噪處理的參數(shù)后，點(diǎn)

2、擊de-noise按鈕，就會出現(xiàn)最終的小波包分析結(jié)果。最后點(diǎn)擊visual de-noise signal將會出現(xiàn)如下圖所示的降噪處理后的結(jié)果。上述就是小包分形降噪處理的簡單步驟，可以處理簡單的近似的小包處理，讀者可以根據(jù)需要進(jìn)一步擴(kuò)展，下面給出上面工具包的原始MATLAB程序，方便讀者進(jìn)一步擴(kuò)展和創(chuàng)新。clc;clear;xx=load('實(shí)驗(yàn)數(shù)據(jù).txt');x=xx(:,1)'y=xx(:,2)'x=x(end:-1:1);y=y(end:-1:1);wname ='sym6' lev=5;c,l = wavedec(y,lev,wname

3、);sigma = wnoisest(c,l,1);alpha=2;thr1 = wbmpen(c,l,sigma,alpha); thr2,nkeep = wdcbm(c,l,alpha);xd1 = wdencmp('gbl',c,l,wname,lev,thr1,'s',1); xd2,cxd,lxd,perf0,perf2= wdencmp('lvd',c,l,wname,lev,thr2,'h'); thr,sorh,keepapp = ddencmp('den','wv',x);xd3

4、= wdencmp('gbl',c,l,wname,lev,thr,'s',1);figuresubplot(4,1,1); plot(y); title('原始信號','fontsize',10);subplot(4,1,2); plot(xd1); title('使用penalty閾值降噪后的信號','fontsize',10);subplot(4,1,3); plot(xd2); title('使用Birge-Massart閾值降噪后的信號','fontsize',10);subplot(4,1,4); plot(xd2