小波降噪與FFT降噪比較

郝文廣，丁常富，梁 娜

(華北電力大學 能源動力與機械工程學院，河北 保定 071003)

小波降噪與FFT降噪比較

郝文廣，丁常富，梁 娜

(華北電力大學 能源動力與機械工程學院，河北 保定 071003)

汽輪機檢測信號包含著豐富的反映其運行狀態的信息，但由于各種噪聲干擾，實測信號往往非?；靵y，以致難以從中獲得有用的信息。利用小波及小波包分解與重建以及小波變換對檢測信號消噪，并畫軸心軌跡圖與降噪前對比，結果表明：利用小波降噪對信號進行處理，效果良好。

小波分析;FFT;降噪

1 小波降噪的原理和方法

降噪是小波分析在信號處理領域的主要用途之一。由于在正交小波中，正交基的選取比傳統方法更接近實際信號本身，通過小波變換可以更容易地分離出噪聲部分。傳統降噪方法的不足之處在于使信號變換后的熵增高，無法刻畫信號的非平穩特性并且無法得到信號的相關性［1］。

振動信號的噪聲一般可以分為有色噪聲和白噪聲兩類。在降噪過程中，對于有色噪聲的消除，一般采用首先對有色噪聲白色化，然后再按照白噪聲的消噪方式進行處理。因此，消噪研究的重點是如何消除信號中的白噪聲。

一個含有噪聲的一維信號的模型可以表示成如下的形式:

式中:f(t)為真實信號;e(t)為噪聲;σ為噪聲強度;s(t)為含噪聲的信號。

認為e(t)為高斯白噪聲，則噪聲等級為1。在實際工程應用中，有用的信號通常表現為低頻信號或是一些比較平穩的信號，而噪聲信號則通常表現為高頻信號。所以消噪過程可按如下方法進行處理:首先對信號進行小波分解，根據小波分解的特點，將信號分為低頻段和高頻段，則噪聲部分通常包含在高頻段中，因此對高頻段以門限域值等形式對小波系數進行處理，這種處理符合噪聲在高頻部分均勻密集的特點。最后對信號進行重構，實現消噪的目的。對信號s(t)的消噪的目的就是要消除信號中的噪聲部分σ·e(t)，從s(t)中恢復出真實信號f(t)。

一維信號利用小波進行除噪的步驟如圖1。

圖1 小波降噪的步驟Fig．1 Steps of wavelet denoise process

2 小波降噪與FFT降噪比較

在對信號的降噪過程中，由于傅里葉分析是將信號完全在頻域中進行分析，信號在時間軸上的任何一個突變，都會影響信號的整個頻譜，不能給出信號在某個時間點的變化情況。小波分析由于能同時在時頻域中對信號進行分析，具有多分辨分析的功能，所以能在不同的分階層上有效地區分信號的突變部分和噪聲，從而實現信號的降噪。

以MATLAB的含噪信號noisdopp為例，分別以小波分析的方法和快速傅里葉變換 (FFT)的方法對其降噪，并對結果進行定量對比分析。同時為了比較兩種不同的閾值降噪的結果，進行小波分析降噪采用默認閾值和給定閾值兩種方法。小波降噪后的結果如圖2所示，FFT降噪后的結果如圖3所示。

從圖3中可以看出，在這個例子中原始信號noisdopp的初始發展階段振蕩頻率很高，但經過FFT變換后的低通濾波把這些高頻成分過濾掉了，同時很多有用的信號成分也可能丟失掉。而在圖2中，默認閾值和給定閾值方法降噪的信號都很好地保留了信號發展初期的高頻特性，且性能參數給定消噪處理后的信號較為光滑。從直觀上顯示了小波降噪較FFT降噪的優越性，下面對其降噪結果進行定量對比。對兩種方法降噪后的信號計算能量比例和標準差，這兩個數據體現了在信號降噪準則下的降噪效果，結果見表1。

表1 降噪后的信號能量百分比和標差Tab．1 After noise reduction and standard deviation of the percentage of signal energy

使用小波變換進行噪聲抑制時，降噪后的能量比例 (97.74%和96.45%)均不低于FFT消噪的結果(96.45%)，這反映了小波降噪后的信號光滑性較好;而前兩種消噪方法的標準差(28.071 4和31.054 8)則小于FFT消噪后的標準差，這體現了小波降噪后的信號與原信號保持了更高的相似度。

3 實例及結果分析

試驗采用美國本特利公司生產的ADRE for Windows 208DAIU 8通道振動數據采集分析系統。鍵相信號及振動信號均取自機組配置的3 300系統，本系統采用非接觸式電渦流傳感器測量軸徑振動，軸振測點分X方向和Y方向布置，兩方向夾角為90°。某電廠11號機組定速3 000 r/min時測得數據經降噪處理后，圖4為原信號與各層分解后波形，圖5為降噪前軸心軌跡圖，圖6為降噪后軸心軌跡圖。

該機于2009年3月10日凌晨大修后首次沖車，2 040 r/min暖機之前機組振動未見異常，剛至2 040 r/min時各瓦軸振動數據如表2所示。

圖6 降噪后軸心軌跡圖Fig．6 Orbit after denoise

表2 各軸瓦振動數據Tab．2 Bearing vibration data

隨著2 040 r/min暖機時間的延長，1X，2X方向振動開始爬升，引起變化主振頻率是基頻，但基頻相位較穩定，其中1X方向振動爬升較快，當1X方向振動爬到132 μm時，運行人員打閘停機;當轉速降至1 650 r/min時，1X方向振動最大到330 μm;轉速降到500 r/min的1X方向振動和升速至500 r/min的振動基本一樣。鑒于以上振動特征，認為造成1X，2X方向振動異常的原因是1號軸瓦處轉子和軸端汽封發生了碰摩，使高中壓轉子產生了熱彎，進而帶動2號軸瓦振動爬升。帶負荷后，4X方向振動出現爬升的原因是4號軸瓦處轉子與汽封發生了碰摩，可能碰摩部位是軸端汽封，使轉子產生了熱彎。

4 結論

通過比較小波降噪和FFT降噪對信號分析的結果及其對降噪后的軸心軌跡圖，可知小波降噪對消除噪聲對信號的影響有良好效果，進而對診斷機組故障提供依據，所以小波降噪對工程實踐中信號分析與處理起到了一定作用。

［1］胡昌華，張軍波，夏軍，等．基于Matlab的系統分析與設計——小波分析 ［M］．西安:西安電子科技大學出版社，2000．

［2］楊建國．小波分析及其工程應用［M］．北京:機械工業出版社，2005．

［3］褚福磊，彭志科，馮志鵬，等．機械故障診斷中的現代信號處理方法 ［M］．北京:科學出版社，2009．

［4］湯謹暉，歐陽美娟．小波變換在地震信號降噪中的應用 ［J］．科技廣場，2010，(5):150－152．

Tang Jinhui，Ouyang Meijuan．The application of wavelet transform for seismic signal de-noise［J］．Science Mosaic，2010，(5):150－152．

［5］賈杏托，王成儒．基于多小波變換的圖像去噪技術［J］．計算機工程與應用 2010，46(19):204－206，237．

Jia Xingtuo，Wang Chengru．Image denoising with multiwavelets ［J］．Computer Engineering and Applications．2010，46(19):204－206，237．

［6］崔彥平，傅其鳳，葛杏衛，等．機械設備故障診斷發展歷程及展望［J］．河北工業科技，2004，21(4):59－62．

Cui Yanping，Fu Qifeng，Ge Xingwei，et al．Fault diagnosis of machinery development and prospect［J］．Hebei Journalof Industrial Science ＆Technology． 2004， 21(4):59－62．

Comparison of Wavelet Denoise and FFT Denoise

Hao Wenguang，Ding Changfu，Liang Na

(School of Energy Power and Mechanical Engineering，North China Electric Power University，Baoding 071003，China)

Turbine testing signals contain rich reflects the operation condition of information，but because of the noise，the measured signal is very chaotic，so often get useful information．Use little affected wavelet packet decomposition and reconstruction and wavelet transform to detect signal de-noising，results show that by using wavelet de-noising in signal processing，the effect is good．

wavelet analysis;FFT;denoise

TK268

A

2010－10－26。

郝文廣 (1982－)，男，碩士研究生，從事發電廠系統與設備安全性、延壽與狀態檢測方面的研究，Email:344863323@qq．com。