摘 要: 為提高噪聲污染信號的檢測速度,,研究了一種融合累加平均與小波變換的方法進行信號降噪處理,并測試驗證了其有效性,。測試結(jié)果表明,,融合數(shù)字累加平均和小波變換的去噪方案在處理速度和降噪效果方面遠優(yōu)于單獨累加平均或單獨小波變換。對被測信號10 000次累加平均再小波降噪獲得的降噪效果與30 000次累加平均相同,,但耗時僅為后者的36%,,大大提高了信號檢測速度。分析比較發(fā)現(xiàn),,對帶噪信號先累加平均再小波變換方案的去噪效果優(yōu)于先小波變換再累加平均,。
關(guān)鍵詞: 信號處理;降噪,;累加平均,;小波變換
0 引言
基于光時域反射的分布式布里淵光纖傳感可以長距離感知整條光纖鏈路上被測物理參量(例如:溫度場、應(yīng)力/應(yīng)變等)的空間分布和隨時間連續(xù)變化信息,,在重要軍民大型基礎(chǔ)設(shè)施(如電網(wǎng),、油氣管道、煤礦,、隧道等)的安全監(jiān)測領(lǐng)域有廣泛的應(yīng)用前景與市場空間,,近年得到極大關(guān)注和研究[1-2]。布里淵傳感系統(tǒng)通過檢測背向散射信號的頻移和強度實現(xiàn)信息解,。通常背向布里淵散射信號包含大量噪聲,,其待測有用信號非常微弱。因此,,降低噪聲提高背向散射信號信噪比和縮短測量時間是布里淵光纖傳感系統(tǒng)研究的關(guān)鍵內(nèi)容之一,。
傳統(tǒng)的降噪方法基于單純的累加平均或者單純的小波閾值降噪。累加平均算法需要對信號進行大量的累加平均以提高信號信噪比,,因而導(dǎo)致布里淵傳感系統(tǒng)完成一次測量極為耗時,,大大降低了系統(tǒng)測量的實時性。小波閾值法能對帶噪信號進行包絡(luò)檢測,,使去噪后的信號是原始信號的近似最優(yōu)估計,,但是在噪聲類型及強度未知的情況下不同閾值的選擇對去噪效果影響極大[3]。因此,,小波閾值法通常與其他方法結(jié)合以達到更優(yōu)的降噪效果,。
本文融合累加平均與小波閾值法在信號降噪中各自優(yōu)勢,采取兩種方法相結(jié)合對噪聲信號進行降噪處理。實驗測試結(jié)果表明,,對待測信號先進行一定次數(shù)累加平均,,然后再進行小波變換降噪可以獲得更好的效果。
1 信號處理原理
1.1 數(shù)字累加平均理論
數(shù)字累加平均將一次測量的N個數(shù)據(jù)依次存儲到內(nèi)存單元中,,將下一次測量的N個數(shù)據(jù)與內(nèi)存對應(yīng)單元的數(shù)據(jù)相加再放回原內(nèi)存單元依次循環(huán)K次,,然后對各單元求其平均[4]。被測量信號中夾雜的噪聲一般為高斯白噪聲,,其期望值為零,,因此K次累加平均后改善的信噪比:
其中,SNR0為經(jīng)過累加平均處理后輸出信號的信噪比,,SNRi為未做處理的信號信噪比,。M次不同時刻采樣的累加平均可以使信噪比改善K倍。累加次數(shù)越多,,測試時間越長,,系統(tǒng)實時性越差[5-7]。
1.2 小波閾值選取原理
小波降噪法將帶噪信號經(jīng)過預(yù)處理,,利用小波變換把信號分解到各尺度中,,在每一尺度下把屬于噪聲的小波系數(shù)濾除掉,保留并增強屬于信號的小波系數(shù),,最后再經(jīng)過小波逆變換恢復(fù)信號[8-10],。對一種固定的小波函數(shù),選用的閾值不同降噪效果也不同,。針對布里淵散射信號特點,,閾值法可采取以下三種形式。
?。?)最小極大方差閾值(minimaxi threshold):屬于一種固定的閾值形式,,使得所選的閾值產(chǎn)生最小的極大方差t:
其中N為小波系數(shù)長度。
?。?)固定閾值(sqtwolog threshold):由最小極大方差的閾值t再乘以一個系數(shù)(logN)得到,。
(3)選擇啟發(fā)式閾值(heursure threshold):固定閾值和軟閾值的綜合,,根據(jù)信號信噪比大小利用啟發(fā)函數(shù)自動在固定閾值和軟閾值中選取一個作為去噪閾值,。不同信噪比條件下這兩個閾值去噪效果不同。小波閾值法降噪流程如圖1所示,。
2 實驗測試與算法設(shè)計
根據(jù)上述數(shù)字累加平均和小波閾值原理,,分別測試了基于數(shù)字累加平均、小波變換,、先累加平均再小波變換和先小波變換再累加平均等四種方案對噪聲信號的降噪性能,。其中,,以頻率為1.5 kHz、幅度為200 mV的方波與幅度為2 V的高斯白噪聲相累加模擬布里淵傳感系統(tǒng)的布里淵探測信號,,如圖2所示,。信號采集卡與實際布里淵傳感系統(tǒng)所用一致(型號為NI-PCI-5112),其最高采樣頻率為100 MHz,,精度為8位,。
信號降噪算法流程如圖3所示。將采集卡采到的數(shù)據(jù)送入計算機,,分別進行如下四種方案的測試與分析。(1)用單獨累加平均的算法對信號進行處理:在虛擬平臺上對采集的周期信號分別進行不同次數(shù)的加法運算,,之后將相對應(yīng)的加法運算結(jié)果進行平均處理,;(2)用小波閾值算法對信號進行處理:選用分解級數(shù)為5級的“sym8”小波,其閾值采取最小極大方差閾值,、固定閾值和啟發(fā)式閾值三種,;(3)先累加平均再小波閾值算法處理:不同累加平均次數(shù)后再分別采用上述三種不同的閾值對信號進行處理;(4)先小波閾值處理再累加平均運算:不同小波閾值處理后再進行不同的次數(shù)累加平均,。上述方法具體信號降噪效果將在下面分別進行分析討論,。
3 實驗結(jié)果及分析
首先,根據(jù)圖3所示的算法,,分別對上述四種方案進行驗證及分析,。圖4為累加次數(shù)分別為300次、1 000次和10 000次處理后的信號,,它們所用的測試時間分別為 0.12 s,、0.4 s和2 s。此結(jié)果表明,,隨著累加次數(shù)的增加高斯白噪聲越來越小,,但處理時間越來越長。
其次,,對被測信號分別進行基于最小極大方差閾值,、固定閾值和啟發(fā)式閾值的小波降噪,測試結(jié)果如圖5所示,。由圖可知,,最小極大方差閾值(圖5(a))和啟發(fā)式閾值(圖5(c))的選取規(guī)則較為保守,導(dǎo)致其降噪效果不理想,。經(jīng)過固定閾值處理后的信號幅度更接近原始信號(-0.2 V~0.2 V)(圖5(b)),,有效去除了隨機噪聲。由此可見,,針對本次實驗固定閾值法是最佳選擇,。從降噪處理所需時間而言,任一小波降噪法處理時間約為 0.1 Vs,僅與不到300次累加平均時間相同,。
由圖4,、圖5可知,單獨基于累加平均和小波變換法降噪存在處理時間長和降噪效果不理想等不足,。圖6為融合兩種方法的信號去噪處理結(jié)果:先進行一定次數(shù)累加平均再小波變換,。其中,對被測信號分別進行300次,、1 000次和10 000次累加平均,、再做固定閾值算法處理,測試時間分別為0.2 s,、0.5 s和2.15 s,,累加平均耗時在總耗時中占主導(dǎo)地位?;谄溆鄡煞N閾值法,,去噪效果基本相同。若僅采取累加平均法降噪達到圖6(c)所示的降噪效果,,則需約30 000次的累加平均,,其處理用時高達6 s。因此,,融合累加平均和小波變換技術(shù)的信號降噪方案大大降低信號處理時間,,提高了噪聲信號檢測速度。
最后,,將被測信號先進行小波變換降噪,、再進行不同次數(shù)的累加平均。圖7~圖9所用的閾值法分別為最小極大方差閾值,、固定閾值和啟發(fā)式閾值,。結(jié)果表明,由于被測信號的噪聲功率比較大,,經(jīng)過三種不同的閾值處理后,,固定閾值處理的降噪效果最好,最小極大方差閾值處理次之,。由于啟發(fā)式閾值法選取規(guī)則較為保守(閾值較小,,導(dǎo)致只有部分系數(shù)置零),所以降噪效果不明顯,。與圖6所示去噪效果相比較可知,,在降噪處理耗時基本相同的情況下,第四種方案(即先小波變換再累加平均)要遜于第三種方案(即先累加平均再小波變換),。
綜合以上四種方案可得出:(1)小波閾值處理用時最少,,但降噪效果差,。(2)累加平均算法的降噪效果要明顯優(yōu)于小波閾值處理,并且隨著累加次數(shù)的增加,,降噪效果會越來越好,,但運行時間會越來越長,遠遠不能滿足本研究快速處理的要求,。(3)累加平均算法與小波閾值處理相結(jié)合,,降噪效果明顯優(yōu)于單獨的數(shù)字平均濾波和小波變換方法。與單獨累加平均降噪相比較,,在相同的降噪效果基礎(chǔ)上,,該方案有效減少數(shù)字平均濾波的累加次數(shù),極大提高了系統(tǒng)噪聲信號的檢測速度,。10 000次累加平均后融合小波降噪獲得的信號與30 000次累加平均濾波效果一致,,但其耗時(2.15 s)遠小于單獨累加平均(6 s)。通過比較分析被測信號經(jīng)過小波閾值處理和累加平均算法的先后順序可知,,被測信號先經(jīng)過累加平均算法再由小波閾值處理的方法降噪效果最佳。
4 結(jié)論
為提高布里淵傳感系統(tǒng)中微弱的布里淵探測信號檢測速度,,測試了四種不同方案的噪聲信號降噪效果:累加平均,、小波變換、先累加平均再小波變換和先小波變換再累加平均,。其中,,小波變換采用三種不同的閾值選擇原理:最小極大方差閾值、固定閾值和選擇啟發(fā)式閾值,。測試結(jié)果表明,,融合數(shù)字累加平均和小波變換的信號降噪效果最佳,處理速度大大提高,。10 000次累加平均結(jié)合小波降噪方案所獲得的信號降噪效果與30 000次累加平均濾波相同,,但耗時僅為后者的36%。此外,,對帶噪信號先累加平均再小波降噪方案的去噪效果優(yōu)于先小波變換再累加平均,。此方法能夠有效應(yīng)用于信號降噪系統(tǒng)處理,可提高光纖傳感系統(tǒng)的布里淵探測信號檢測速度,。
參考文獻
[1] BAO X,,CHEN L. Recent progress in Brillouin scattering based fiber sensors[J]. Sensors,2011,,11(4):4152-4187.
[2] 趙曉東,,路元剛,胡君輝,,等.波長掃描型布里淵光時域反射儀[J].中國激光,,2012,,39(8):1-5.
[3] 宋牟平,陳翔.基于實時小波變換信號處理的相干檢測布里淵光時域反射計[J].光學(xué)學(xué)報,,2009,,29(10): 2818-2821.
[4] FARAHANI M A, WYLIE M T V,, CASTILLO-GUERRA E,, et al. Reduction in the number of averages required in BOTDA sensors using wavelet denoising techniques[J]. Journal of Lightwave Technology, 2012,,30(8):1134-1142.
[5] 李星容,,李永倩,張碩.同步疊加平均算法抑制噪聲的Labview實現(xiàn)[J].華北電力大學(xué)學(xué)報,,2009,,36(4):74-76.
[6] 周倩婷,危峻,,徐志鵬.噪聲特性對多次采集累加平均技術(shù)的影響[J].紅外與激光工程,,2010,39(5):959-962.
[7] 楊鶯,,韋育森.提取周期微弱信號的幾種軟件實現(xiàn)方法[J].測控技術(shù),,2008,27(9):22-24.
[8] 張弦,,王宏力.進化小波消噪方法及其在滾動軸承故障診斷中的應(yīng)用[J].機械工程學(xué)報,,2010,46(15):76-81.
[9] 田玉靜,,左紅偉.小波消噪閾值算法優(yōu)化[J].聲學(xué)技術(shù),, 2009,28(4):503-506.
[10] 宋牟平,,趙斌.希爾伯特變換處理的布里淵散DOFS 的研究[J].光子學(xué)報,,2005,34(9):1328-1331.