第六章小波分析方法在濾波和消噪方面的應(yīng)用

第六章小波分析方法在濾波和消噪方面的應(yīng)用16.1小波分析在常規(guī)濾波方面的應(yīng)用低通濾波要求：要求保留原信號(hào)中特定的低頻范圍信號(hào)，去掉高頻分量。方法：使用小波包算法，保留需要低頻部分，高頻部分置零。高通濾波要求：保留原信號(hào)中的高頻范圍信號(hào)，去掉低頻分量。方法：使用小波包算法，保留需要高頻部分，低頻部分置零。帶通濾波要求：保留原信號(hào)中的某個(gè)特定頻帶，去掉其他頻率成分。方法：使用小波包算法，保留需要頻率部分，其他的部分置零。6.1小波分析在常規(guī)濾波方面的應(yīng)用低通濾波26.2小波分析在消噪方面的應(yīng)用1白噪聲的特點(diǎn)白噪聲是隨機(jī)性的。不同的和不相關(guān)。不同時(shí)刻的采樣值不相關(guān)；零均值且能量無(wú)限；時(shí)域表現(xiàn)是均勻密集的；包含有全部頻譜。白噪聲的小波分析特點(diǎn)1）實(shí)際工程中有用信號(hào)通常為低頻信號(hào)或一些比較平穩(wěn)的信號(hào)，噪聲通常為高頻信號(hào)。通過(guò)小波分解，噪聲表現(xiàn)在高頻部分，且衰減嚴(yán)重。各尺度上不相關(guān)；6.2小波分析在消噪方面的應(yīng)用1白噪聲的特點(diǎn)32）消噪方法可分為強(qiáng)制消噪和門限消噪強(qiáng)制消噪直接將小波分解的高頻系數(shù)置零，然后進(jìn)行信號(hào)重構(gòu)。門限消噪根據(jù)經(jīng)驗(yàn)和某種依據(jù)設(shè)定門限值（閾值），對(duì)高頻部分系數(shù)用門限值處理，大于門限的保留，低于門限的置零。門限消噪又可分為硬閾值和軟閾值消噪，前者設(shè)定固定閾值，后者根據(jù)估計(jì)計(jì)算自動(dòng)獲取。2）消噪方法可分為強(qiáng)制消噪和門限消噪4第六章小波分析方法在濾波和消噪方面的應(yīng)用詳解ppt課件5一維信號(hào)的消噪步驟一維信號(hào)的小波分解。選擇一個(gè)小波并確定一個(gè)小波分解的層次，然后對(duì)信號(hào)進(jìn)行N層小波分解；小波分解的高頻系數(shù)的閾值量化。對(duì)第一到第N層的每一個(gè)高頻系數(shù)，選擇一個(gè)閾值進(jìn)行軟閾值量化處理；一維小波的重構(gòu)。根據(jù)小波分解的第N層低頻系數(shù)和經(jīng)過(guò)量化的各層高頻系數(shù)，進(jìn)行一維信號(hào)的小波重構(gòu)。一維信號(hào)的消噪步驟一維信號(hào)的小波分解。選擇一個(gè)小波并確定一個(gè)6信號(hào)模型：f(i)為真實(shí)信號(hào)，e(i)為噪聲。一般的工程應(yīng)用中，有用信號(hào)常為低頻信號(hào)，或比較平穩(wěn)的信號(hào)，噪聲常為高頻。信號(hào)模型：f(i)為真實(shí)信號(hào)，e(i)為噪聲。7用于信號(hào)消噪的主要小波指令函數(shù)名功能cwt一維連續(xù)小波變換dwt單尺度一維離散小波變換dwtper單尺度一維離散小波變換（周期性）wavedec多尺度一維小波分解一維小波分解函數(shù)用于信號(hào)消噪的主要小波指令函數(shù)名功能cwt一維連續(xù)小波變換d8一維小波重構(gòu)函數(shù)函數(shù)名功能idwt單尺度一維離散小波逆變換idwtper單尺度一維離散小波重構(gòu)（周期性）waverec多尺度一維小波重構(gòu)upwlev單尺度一維小波分解的重構(gòu)wrcoef對(duì)一維小波系數(shù)進(jìn)行單支重構(gòu)upcoef一維系數(shù)的直接小波重構(gòu)一維小波重構(gòu)函數(shù)函數(shù)名功能idwt單尺度一維離散小波逆變換i9一維小波消噪和壓縮函數(shù)函數(shù)名功能thselect信號(hào)消噪的閾值選擇wthresh進(jìn)行軟閾值或硬閾值處理wthcoef一維信號(hào)的小波系數(shù)閾值處理wden用小波進(jìn)行一維信號(hào)的自動(dòng)消噪ddencmp獲取在消噪或壓縮過(guò)程中的默認(rèn)值閾值（軟或硬）、熵標(biāo)準(zhǔn)wdencmp用小波進(jìn)行信號(hào)的消噪和壓縮一維小波消噪和壓縮函數(shù)函數(shù)名功能thselect信號(hào)消噪的閾10thselect

THR=THSELECT(X,TPTR)returnsthresholdX-adaptedvalueusingselectionruledefinedbystringTPTR.TPTR選項(xiàng)閾值選擇規(guī)則rigrsure采用史坦（stein)的無(wú)偏似然估計(jì)（UnbiasedRiskEstimate)原理（SURE)進(jìn)行自適應(yīng)閾值選擇sqtwolog固定的閾值形式，等于sqrt(2*log(length(s)))heursure啟發(fā)式閾值選擇minimaxi用極大極小原理選擇的閾值thselect

THR=THSELECT(X,TPTR11舉例init=2055415866;rand('seed',init);x=randn(1,1000);plot(x);thr1=thselect(x,'rigrsure');thr2=thselect(x,'sqtwolog');thr3=thselect(x,'heursure');thr4=thselect(x,'minimaxi');執(zhí)行程序后輸出：[thr1thr2,thr3,thr4]＝2.73163.71693.71692.2163舉例init=2055415866;12wthreshY=wthresh(X,SORH,T)returnssoft(ifSORH='s')orhard(ifSORH='h')T-thresholdingoftheinputvectorormatrixX.Tisthethresholdvalue.

Y=wthresh(X,‘s’,T)returnsY=SIGN(X).(|X|-T)+,把信號(hào)的絕對(duì)值與閾值比較，小于和等于閾值的點(diǎn)變?yōu)?，大于閾值的點(diǎn)變?yōu)橹蹬c閾值之差。

Y=wthresh(X,'h',T)returnsY=X.1-(|X|>T),與閾值比較，大于閾值等于原值，小于等于閾值置零，hardthresholdingiscruder.wthreshY=wthresh(X,SORH,T)13舉例y=linspace(-1,1,100);figure(1);subplot(311);plot(y)title('原始信號(hào)')grid;thr=0.4;ythard=wthresh(y,'h',thr);subplot(312);plot(ythard);title('硬閾值信號(hào)');grid;ytsoft=wthresh(y,'s',thr);subplot(313);plot(ytsoft);title('軟閾值信號(hào)');grid;舉例y=linspace(-1,1,100);14第六章小波分析方法在濾波和消噪方面的應(yīng)用詳解ppt課件15wthcoef進(jìn)行一維信號(hào)小波系數(shù)閾值處理。nc=wthcoef(‘d’,c,l,n,p)，返回小波分解結(jié)構(gòu)nc，n為尺度向量，p是一個(gè)包含把較小系數(shù)置零的百分比信息向量，與n同長(zhǎng)度。nc=wthcoef(‘d’,c,l,n)將指定尺度的高頻系數(shù)全部置零；nc=wthcoef(‘a(chǎn)’,c,l)將低頻系數(shù)全部置零。nc=wthcoef(‘t’,c,l,n,t,sorh)返回經(jīng)閾值處理后的分解向量。n為指定尺度向量，t為對(duì)應(yīng)的閾值向量，sorh用來(lái)指定硬或軟閾值。wthcoef進(jìn)行一維信號(hào)小波系數(shù)閾值處理。16loadnoissin;s=noissin(1:1000);subplot(421)plot(s);title('原始信號(hào)')[c,l]=wavedec(s,3,'db3');n=[1,2,3];%設(shè)置尺度向量p=[98,99,97];%設(shè)置閾值向量nc1=wthcoef('d',c,l,n,p);%對(duì)高頻系數(shù)進(jìn)行閾值處理nc2=wthcoef('d',c,l,n);%對(duì)n指定尺度的高頻系數(shù)全部置零nc3=wthcoef('a',c,l);%低頻系數(shù)置零%nc4=wthcoef('t',c,l,n);ss1=waverec(nc1,l,'db3');ss2=waverec(nc2,l,'db3');ss3=waverec(nc3,l,'db3');subplot(422)plot(ss1);title('消噪后信號(hào)1')subplot(425)plot(ss2);title('消噪后信號(hào)2');subplot(426)plot(ss3);title('消噪后信號(hào)3')ss4=wrcoef('a',c,l,'db3',3);subplot(529)plot(ss4);title('消噪后信號(hào)4')loadnoissin;subplot(422)17

18wden用小波進(jìn)行一維信號(hào)的自動(dòng)消噪[XD,CXD,LXD]=WDEN(X,TPTR,SORH,SCAL,N,'wname')[XD,CXD,LXD]=WDEN(C,L,TPTR,SORH,SCAL,N,'wname')TPTR閾值選擇規(guī)則；見(jiàn)前面。Scal定義所乘的閾值是否需要重新調(diào)整：Scal＝‘one’時(shí)不用調(diào)整；Scal＝‘sln’時(shí)，根據(jù)第一層的系數(shù)進(jìn)行一次噪聲層的估計(jì)來(lái)調(diào)整閾值；Scal＝‘mln’時(shí)，在不同層估計(jì)噪聲，以此調(diào)整閾值。XD消噪后的信號(hào)。CXD,LXD消噪后的小波分解結(jié)構(gòu)。wden用小波進(jìn)行一維信號(hào)的自動(dòng)消噪19snr=3;%設(shè)置信噪比init=2055615866;%設(shè)置隨機(jī)數(shù)的初始值[xref,x]=wnoise(3,11,snr,init);lev=5;xd=wden(x,'heursure','s','one',lev,'sym8');figure(1)subplot(321);plot(xref);axis([1,2048,-10,10]);title('原始信號(hào)')subplot(322);plot(x);axis([1,2048,-10,10]);title('含噪信號(hào)，信噪比3');subplot(323);plot(xd);axis([1,2048,-10,10]);title('用heursure閾值去噪后的信號(hào)');xd2=wden(x,'rigrsure','s','sln',lev,'sym8');subplot(324);plot(xd2);axis([1,2048,-10,10]);title('用rigrsure閾值去噪后的信號(hào)');xd3=wden(x,'sqtwolog','s','sln',lev,'sym8');subplot(325);plot(xd2);axis([1,2048,-10,10]);title('用sqtwolog閾值去噪后的信號(hào)');xd4=wden(x,'minimaxi','s','sln',lev,'sym8');subplot(326);plot(xd4);axis([1,2048,-10,10]);title('用minimaxi閾值去噪后的信號(hào)');snr=3;%設(shè)置信噪比xd2=wden(x,'rigrsu20第六章小波分析方法在濾波和消噪方面的應(yīng)用詳解ppt課件21小波分析在平穩(wěn)信號(hào)消噪中的應(yīng)用平穩(wěn)信號(hào)往往表現(xiàn)為低頻信號(hào)。1）周期信號(hào)中混有白噪聲一般使用自相關(guān)過(guò)程消噪即可2）時(shí)頻受限信號(hào)進(jìn)行小波包分解，剔除頻限之外部分，再進(jìn)行門限消噪。小波分析在平穩(wěn)信號(hào)消噪中的應(yīng)用平穩(wěn)信號(hào)往往表現(xiàn)為低頻信號(hào)。22小波分析在非平穩(wěn)信號(hào)消噪中的應(yīng)用工程實(shí)際應(yīng)用中存在大量的非平穩(wěn)信號(hào)，一方面需要消除噪聲，同時(shí)需要保留突變成分，而突變成分可能就在高頻中。一般而言，噪聲尤其是白噪聲在高頻部分具有均勻的表現(xiàn)，而突變部分只在某些分量中幅度較大。在多個(gè)尺度下同時(shí)進(jìn)行觀察，以減少誤判和漏判。消噪時(shí)可以只保留各尺度與突變點(diǎn)對(duì)應(yīng)的模極大值，其他用零代替；也可以保留突變點(diǎn)及其附近的模極大值，效果更好。小波分析在非平穩(wěn)信號(hào)消噪中的應(yīng)用工程實(shí)際應(yīng)用中存在大量的非平23正交小波變換對(duì)應(yīng)的小波變換模極大值點(diǎn)在不同尺度上存在移位現(xiàn)象，是由于不對(duì)稱造成的。需要使用對(duì)稱小波進(jìn)行分析。如果將相鄰兩個(gè)尺度層上的小波變換對(duì)應(yīng)相乘，能進(jìn)一步增強(qiáng)突變點(diǎn)的表現(xiàn)，并能更好的消除噪聲。正交小波變換對(duì)應(yīng)的小波變換模極大值點(diǎn)在不同尺度上存在移位現(xiàn)象24第六章小波分析方法在濾波和消噪方面的應(yīng)用詳解ppt課件25第六章小波分析方法在濾波和消噪方面的應(yīng)用詳解ppt課件26第六章小波分析方法在濾波和消噪方面的應(yīng)用詳解ppt課件27第六章小波分析方法在