【摘要】傳統(tǒng)小波去噪雖然可在一定程度上去除噪聲對原始信號的干擾，但去噪效果并不理想。針對傳統(tǒng)小波去噪中存在的問題，提出一種改進(jìn)的小波去噪方法，并將改進(jìn)小波去噪與EEMD-HHT有機(jī)結(jié)合，進(jìn)而提出一種基于改進(jìn)小波去噪的EEMD—HHT信號處理新方法?；贛ATLAB軟件，分別利用EEMD-HHT方法、基于傳統(tǒng)小波去噪的EEMD-HHT信號處理方法和基于改進(jìn)小波去噪的EEMD-HHT信號處理方法對外圈故障滾動軸承進(jìn)行故障診斷試驗，試驗結(jié)果與理論計算結(jié)果對比分析表明，基于改進(jìn)小波去噪的EEMD-HHT信號處理方法最為有效。

【關(guān)鍵詞】改進(jìn)小波去噪;EEMD-HHT;改進(jìn)小波去噪與EEMD-HHT有機(jī)結(jié)合;滾動軸承故障診斷

1.引言

Huang N E等人于2007年提出基于聚合經(jīng)驗?zāi)B(tài)分解EEMD的HHT信號處理方法[1]（簡稱EEMD-HHT），該方法已在軸承故障診斷、電能質(zhì)量擾動檢測、語音信號處理和地震信號分析等不同領(lǐng)域取得了顯著成效。但在實際信號測試過程中，外界噪聲會對原始信號產(chǎn)生干擾，這必將影響所測信號的真實度，單獨(dú)使EEMD-HHT對信號進(jìn)行處理，就無法有效解決這一問題。小波變換具有去噪特性，董文智等人提出了一種基于傳統(tǒng)小波去噪的EEMD-HHT信號處理方法[2]，該方法雖然取得了一定的改進(jìn)效果，但由于采用的是基于軟閾值函數(shù)或硬閾值函數(shù)的小波去噪方法（簡稱軟閾值或硬閾值去噪），其去噪效果并不理想[3，4]。為了更好地去除外界噪聲對原始信號的干擾，張弛等人提出一種改進(jìn)的小波去噪方法，此方法在一定程度上解決了上述小波去噪方法的缺點(diǎn)[5]。將此改進(jìn)的小波去噪與EEMD-HHT有機(jī)結(jié)合，提出一種基于改進(jìn)小波去噪的EEMD-HHT信號處理新方法。分別利用EEMD-HHT方法、基于傳統(tǒng)小波去的EEMD-HHT信號處理方法和基于改進(jìn)小波去噪EEMD-HHT信號處理方法對外圈故障滾動軸承進(jìn)行故障診斷試驗，試驗結(jié)果與理論計算結(jié)果對比分析表明，基于改進(jìn)小波去噪EEMD-HHT信號處理方法最為有效。

2.EEMD-HHT的基本思想

EEMD-HHT由聚合經(jīng)驗?zāi)B(tài)分解EEMD和Hilbert變換組成，先將數(shù)據(jù)信號進(jìn)行EEMD分解，將分解得到的固有模態(tài)函數(shù)IMF分量進(jìn)行Hilbert變換，得到Hilbert譜，最后將Hilbert譜在時間軸上積分進(jìn)而得到Hilbert邊際譜[6]。對于Hilbert邊際譜，在理想的情況下，如果某個頻率對應(yīng)的幅值不為零，則表示在原信號中，此頻率存在。而在實際應(yīng)用中，如果某個頻率對應(yīng)的幅值越大，則表示此頻率存在的可能性越大，若某頻率確定存在，則其對應(yīng)的幅值必定較大[7，8]。因此，可以通過對外圈故障滾動軸承故障信號的Hilbert邊際譜進(jìn)行分析，便可判斷軸承是否發(fā)生故障。

3.改進(jìn)的小波去噪方法

在小波去噪過程中，閾值和閾值函數(shù)選取是否得當(dāng)將直接影響到小波去噪效果。傳統(tǒng)的硬閾值函數(shù)為：

（1）

軟閾值函數(shù)為：

（2）

式中，—閾值;—小波系數(shù);—處理后的小波系數(shù)。

傳統(tǒng)的軟閾值函數(shù)或硬閾值函數(shù)雖有應(yīng)用，但也存在明顯不足。具體表現(xiàn)在：硬閾值去噪，小波系數(shù)值連續(xù)性差，將使重構(gòu)的信號出現(xiàn)許多不期望的振蕩，從而失去原始信號的光滑性;軟閾值去噪，雖然去噪后的信號相對平滑，克服了小波系數(shù)不連續(xù)的缺陷，但在小波系數(shù)較大時，處理后的小波系數(shù)與處理前的小波系數(shù)總存在一定的偏差，造成了部分高頻信息的損失，可能會導(dǎo)致重構(gòu)的信號失真。因此，為了更合理的對小波系數(shù)進(jìn)行閾值量化處理，改進(jìn)的閾值函數(shù)為：

（3）

式中：

p—為平移因子

s—為收斂性因子

k—為擴(kuò)散性因子

其中m是常數(shù)，其它參數(shù)同（1）、（2）式。

當(dāng)ωj，k接近時，收斂性因子s使改進(jìn)閾值函數(shù)趨向于軟閾值函數(shù);當(dāng)ωj，k>λ時，擴(kuò)散性因子k則使改進(jìn)閾值函數(shù)趨于硬閾值函數(shù)。這樣對于不同的ωj，k，改進(jìn)閾值函數(shù)會自適應(yīng)調(diào)整為軟閾值函數(shù)或硬閾值函數(shù)，從而兼顧了傳統(tǒng)軟閾值函數(shù)和硬閾值函數(shù)的優(yōu)點(diǎn)，使閾值量化處理更為合理，小波去噪的效果也將會更佳。

也可通過調(diào)節(jié)m的大小來改善去噪效果。若去噪后的信號不夠光滑或存在尖峰時，則減小m的值;相反，當(dāng)去噪結(jié)果過于光滑而丟失了過多的細(xì)節(jié)信息時，則增大m的值。而且當(dāng)m趨向1時，改進(jìn)閾值函數(shù)趨于Garrote函數(shù);當(dāng)m趨近于0時，改進(jìn)閾值函數(shù)趨于軟閾值函數(shù)。

4.基于改進(jìn)小波去噪的EEMD-HHT

在信號測試過程中，外界噪聲信號會不可避免的對原始信號產(chǎn)生干擾，這勢必會影響所測得信號的真實度。雖然EEMD-HHT已經(jīng)具有較好的信號處理效果，但如果單獨(dú)使用EEMD-HHT對信號進(jìn)行處理，就無法有效解決噪聲信號對原始信號產(chǎn)生干擾這一問題，從而影響信號處理結(jié)果的精確度，甚至使信號處理結(jié)果失去實際意義。

基于以上考慮，由于改進(jìn)閾值函數(shù)兼顧了傳統(tǒng)軟閾值函數(shù)和硬閾值函數(shù)的優(yōu)點(diǎn)，從而使改進(jìn)小波去噪方法具有更為良好的去噪特性。為了對信號進(jìn)行更為合理的分析，現(xiàn)將改進(jìn)小波去噪方法和EEMD-HHT進(jìn)行結(jié)合，提出一種基于改進(jìn)小波去噪EEMD-HHT信號處理方法，并通過外圈故障滾動軸承故障診斷實例，來驗證了這種方法的有效性。

5.應(yīng)用實例

采用外圈故障滾動軸承做加速度信號測試試驗，采樣頻率為20KHZ，采樣長度為8192，滾動體個數(shù)z為8，軸承內(nèi)徑d為40;軸承外徑D為90;接觸角α為60度;主軸轉(zhuǎn)速為1566r/min，軸承轉(zhuǎn)動頻率?為26.1，根據(jù)外圈故障特征頻率計算公式[8]：

（4）計算得其外圈故障特征頻率為81.2HZ。傳統(tǒng)的小波去噪選用默認(rèn)閾值去噪（由MATLAB相關(guān)程序自適應(yīng)判別是使用硬閾值去噪還是使用軟閾值去噪），小波基選取sym5函數(shù)，分解層數(shù)為3，m取值為0.64，閾值計算規(guī)則同樣由MATLAB相關(guān)程序自適應(yīng)判別。

先根據(jù)所測得的信號，得到未去噪信號的時域圖1（a），再分別用改進(jìn)小波去噪和默認(rèn)閾值去噪對信號進(jìn)行去噪，得到改進(jìn)小波去噪信號的時域圖1（c）以及默認(rèn)閾值去噪信號的時域圖1（b）。為了更直觀的比較改進(jìn)小波去噪與默認(rèn)閾值去噪的效果，分別算得默認(rèn)閾值去噪信號和改進(jìn)小波去噪信號的信噪比SNT與均方根差RMSE（信號的去噪效果一般采用信噪比SNR和原始信號與去噪后的信號的均方根差RMSE來衡量，且SNR越大，RMSE越小，則去噪效果越好[9，10]），其結(jié)果如表1所示。再對去噪后的信號進(jìn)行EEMD-HHT變換，分別得到改進(jìn)小波去噪的Hilbert邊際譜圖2（c）、默認(rèn)閾值去噪的Hilbert邊際譜圖2（c），同樣得到未去噪的Hilbert邊際譜圖2（a）。

圖1 未去噪信號與去噪后信號圖

表1 改進(jìn)小波去噪與默認(rèn)閾值去噪的效果比較

衡量參數(shù) 改進(jìn)小波去噪 默認(rèn)閾值去噪

SNR

RMSE 9.18

66.21 6.43

91.12

從圖1可知，未去噪信號圖與去噪信號圖差別明顯，說明小波去噪方法確實可以明顯的減小噪聲對原信號的影響。但默認(rèn)閾值去噪信號圖和改進(jìn)小波去噪信號圖差別不明顯。不過從表1可知，改進(jìn)小波去噪的SNR大于默認(rèn)閾值去噪的SNR，而改進(jìn)小波去噪的RMSE小于默認(rèn)閾值去噪的RMSE，說明改進(jìn)小波去噪的效果比默認(rèn)閾值去噪的效果更佳。

圖2 未去噪與去噪的Hilbert邊際譜圖

比較圖2中的（a）、（b）、（c）三圖可知，在故障特征頻率81.2HZ附近，分別出現(xiàn)了一個突出的高幅值頻率86.67HZ、85.45HZ、80.57HZ，根據(jù)Hilbert邊際譜的實際意義可知，軸承確實發(fā)生了外圈故障?，F(xiàn)定義Hilbert邊際譜圖的精度值為：

（5）

計算結(jié)果表明，改進(jìn)小波去噪的Hilbert邊際譜其精度最高，為99.2%;默認(rèn)閾值去噪的Hilbert邊際譜圖其精度其次，為94.8%;未去噪的Hilbert邊際譜圖的精度最低，為93.3%。這說明基于改進(jìn)小波去噪的EEMD-HHT信號處理方法，其診斷結(jié)果更準(zhǔn)確，更有效。而且從圖2還可以看出，改進(jìn)改進(jìn)小波去噪的Hilbert邊際譜的效果更好。

6.結(jié)論

實驗結(jié)果表明，改進(jìn)小波去噪的效果比傳統(tǒng)小波去噪的效果好?；诟倪M(jìn)小波去噪的EEMD-HHT信號處理方法其精度比基于傳統(tǒng)小波去噪的EEMD-HHT信號處理方法更高。這種方法不僅可以用于軸承故障診斷，還可以用于其它領(lǐng)域的信號處理，具有普遍性。但是不管是使用改進(jìn)小波去噪還是使用傳統(tǒng)小波去噪，小波基函數(shù)的選取都是個難點(diǎn)，需要不斷的試驗才能提高其去噪效果。

參考文獻(xiàn)

[1]Huang N E，Shen Z，Long S R，et al. The empirical mode decomposi-tion and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J].Proceeding of Royal Society London，1998，454：903-995.

[2]董文智，張超.基于小波變換和EEMD分解的轉(zhuǎn)子系統(tǒng)故障診斷[J].機(jī)械科學(xué)與技術(shù)，2012，31（6）：973-976.

[3]黃玉昌，侯德文.基于改進(jìn)小波閾值函數(shù)的指紋圖像去噪[J].計算機(jī)工程與應(yīng)用，2014，50（6）：179-181.

[4]張弛，李翔，姚磊.一種改進(jìn)的小波閾值函數(shù)去噪方法[J].計算機(jī)與現(xiàn)代.2014（3）：220-222.

[5]李海東，李青.基于閾值法的小波去噪算法研究[J].計算機(jī)技術(shù)與發(fā)展，2009，19（7）：57-58.

[6]張鑫，吳亞鋒，朱帥琦.基于EEMD-HHT邊際譜的軸承故障診斷[J].科學(xué)技術(shù)與工程，2011，11（31）：7626-7629.

[7]向玲.HHT方法在轉(zhuǎn)子轉(zhuǎn)動故障診斷中的應(yīng)用[J].中國電機(jī)工程學(xué)報，2007，27（35）：85-88.

[8]張永娟，曾國英.振動環(huán)境下螺栓連接的HHT分析[J].機(jī)械設(shè)計與制造，2013（11）：161-163.

[9]王江萍.機(jī)械設(shè)備故障診斷技術(shù)及應(yīng)用[M].西安：西北工業(yè)大學(xué)出版社，2001.

[10]劉衛(wèi)東，劉尚合.小波去噪函數(shù)的改進(jìn)方法分析[J].高電壓器技術(shù)，2007，33（10）：59-63.

[11]鄭思麗，桂預(yù)風(fēng).小波閾值法在機(jī)械信號降噪研究中的應(yīng)用[J].機(jī)械設(shè)計與制造，2014（1）：223-225.

作者簡介：蔡光達(dá)（1988—），陜西理工學(xué)院機(jī)械工程學(xué)院碩士研究生在讀，研究方向：新型傳動機(jī)械的設(shè)計與研究;故障診斷。