邱愛中,崔曉,張海泉,邱大為

(1.鄭州師范學院,鄭州 450044；2.中國科學技術(shù)大學,合肥 230000)

機械故障診斷中的傳感器往往處于惡劣的環(huán)境中，直接獲取的原始信號必然含有大量噪聲，必須通過降噪才能獲取有用的信號信息特征。弱故障信號的提取和恢復一直是機械故障診斷的難點和要點，例如，動力裝備的轉(zhuǎn)子、齒輪、軸承等部件在故障初期，不明顯的故障信號往往比較微弱，有時完全淹沒在強噪聲中，而且信號中往往同時混雜了白噪聲和脈沖噪聲。

為提高降噪效果，文獻[1-2]指出常規(guī)的離散小波變換會產(chǎn)生偽Gibbs效應(yīng)，改進措施是采用atrous算法進行非降采樣的小波變換(Undecimated Discrete Wavelet Transform，UDWT)。然而該變換會使計算量猛增，且冗余信息不易消除，增加后續(xù)數(shù)據(jù)處理的難度。文獻[3-5]指出復數(shù)小波在數(shù)據(jù)處理上相對冗余且計算量較少，不足之處是要構(gòu)造完全重構(gòu)的逆濾波器在實踐中非常困難。而雙樹復小波變換(Dual-Tree Complex Wavelet Transform，DT-CWT)[6-7]完全重構(gòu)的濾波器較易實現(xiàn)，且具有近似平移不變性、能完全重構(gòu)、相對較小的計算量等優(yōu)點。文獻[8-9]認為中值濾波能較好抑制脈沖噪聲，且具有良好的邊緣保持特性，可作為非線性濾波的典型代表。

因此，研究了復數(shù)小波的構(gòu)造，簡化了濾波器的設(shè)計，在此基礎(chǔ)上提出一種基于雙樹復小波和中值濾波相結(jié)合的消噪方法。

1 基于DT-CWT和中值濾波的降噪法

1.1 雙樹復小波構(gòu)造的數(shù)學模型

雙樹復小波的實現(xiàn)是通過一對濾波器組同時處理輸入數(shù)據(jù)[4-6]，其包含2個平行的小波分支(小波樹)，2棵樹中的實數(shù)濾波器h0(n)和h1(n)，g0(n)和g1(n)分別表示1個共軛正交濾波器對。

H(0)=1。

(1)

(3)

設(shè)φ為正交尺度函數(shù)，h={h0,h1,…,hN}是對應(yīng)的雙尺度方程濾波器，則構(gòu)造這樣的正交小波時濾波器應(yīng)滿足

(4)

(5)

其中，當n=0時，δ=1；n≠0時，δ=0。

1)小波函數(shù)ψ(x)需具有m階消失矩特性。

2)低通濾波器函數(shù)系數(shù)要滿足

(6)

3)H(ω)在ω=π處存在m重零點。

滿足以上條件可得

(7)

得到濾波器后，令

(8)

在0～2π范圍內(nèi)求出|H0(ω)|的上界值，如果其不大于2m-1，則H就為所求濾波器函數(shù)。聯(lián)立(4)～(8)式就可以構(gòu)造符合要求的各種互為Hilbert變換的正交小波。

1.2 中值濾波

中值濾波的原理是設(shè)計一種滑動窗口處理抽樣數(shù)據(jù)，用滑動窗口中心位置的信號抽樣值取代當前窗口內(nèi)所有抽樣的對應(yīng)中位值。例如，假設(shè)數(shù)據(jù)處理的窗口抽樣長度n=2k+1(k為正整數(shù))，那么第i時刻輸入信號序列在窗口內(nèi)的抽樣樣本就為(xi-k，xi-k+1，…，xi，…，xi+k-1，xi + k)，則此時中值濾波器的輸出為

yi=med(xi-k，xi-k+1，…，xi，…，xi+k-1，xi + k)，

(9)

式中：med()表示取此時處理數(shù)據(jù)窗口內(nèi)所有抽樣數(shù)的中位數(shù)。由此可知，窗口寬度若取為n=2k+1，如果信號中混雜有脈沖信號噪聲，信號抽樣序列中脈沖寬度小于或等于k的脈沖噪聲將被去除；而當脈沖噪聲寬度大于或等于k+1抽樣間隔時，脈沖噪聲未得到抑制，仍被保留。所以中值濾波對脈沖信號寬度滿足上述條件的脈沖噪聲有較好的去除抑制作用，但對隨機白噪聲的抑制能力卻微乎其微[7]。

1.3 基于DT-CWT和中值濾波的降噪方法

在離散正交復小波變換構(gòu)造的基礎(chǔ)上，與中值濾波數(shù)據(jù)處理原理有機結(jié)合，提出了該降噪思想，充分考慮了DT-CWT的完全重構(gòu)降噪特性，并可同時抑制白噪聲和脈沖噪聲，其具體步驟為：

1)將被檢測信號用DT-CWT進行變換，得到2棵樹的正交離散復小波每層的小波系數(shù)。

2)由于噪聲主要集中在小波分解細節(jié)系數(shù)中，故需對每層的細節(jié)系數(shù)首先進行中值濾波處理，再實施噪聲估計和閾值選擇，進行閾值降噪濾波，而讓尺度系數(shù)保持不變。

3)利用小波變換軟閾值估計方法，獲取降噪閾值，對細節(jié)系數(shù)實施軟閾值降噪處理。

4)將2棵樹的系數(shù)按照上述方法進行閾值修改后，對系數(shù)再進行逆DT-CWT，重構(gòu)得到降噪后的被分析信號。

2 仿真分析

2.1 非平穩(wěn)、非線性信號

為了驗證降噪效果，進行MATLAB仿真檢驗，利用該軟件產(chǎn)生具有代表性的3種非平穩(wěn)信號Heavysine，Bloeks，Bumps，并在其中加入白噪聲和有色噪聲，然后用上述算法及一些常規(guī)方法對染噪信號進行降噪處理，并對比降噪效果。

利用降噪處理后重構(gòu)信號的信噪比Rsn和均方差誤差Rmse對這些降噪算法進行對比評判。評判的標準是降噪后信號的信噪比越高，均方差越小，則降噪效果越好；反之，說明降噪效果相對較差。評判的公式為

(10)

(11)

不同方法對應(yīng)的信噪比和均方差見表1。由表可知：基于DT-CWT和中值濾波的降噪方法能夠獲得比其他降噪法更高的信噪比，顯示出較好的對白噪聲降噪效果。

表1 不同方法消噪處理后的結(jié)果

某非平穩(wěn)、非線性染噪信號及不同降噪方法處理后的結(jié)果如圖1所示。從圖中可以看出：基于DT-CWT和中值濾波的降噪方法降噪效果相對較好，可以更好地保留強高斯白噪聲背景下的信號特征的細節(jié)信息，更好地再現(xiàn)原始無噪信號，其他方法降噪后仍含有較大噪聲，原始信號細節(jié)信息有所丟失。

圖1 不同方法對染噪信號的降噪結(jié)果

2.2 沖擊信號

混雜白噪聲和脈沖噪聲的染噪信號及不同方法對應(yīng)的降噪結(jié)果如圖2所示。從圖中可以看出：其他方法只能抑制白噪聲，對去除脈沖噪聲無能為力；而基于DT-CWT和中值濾波的降噪方法可以更好地去除強噪聲，還原原始的無噪信號，保留強脈沖噪聲背景下的弱信號特征信息，實現(xiàn)更好的降噪效果。

圖2 不同方法去除脈沖噪聲的效果比較

3 水輪機軸承故障監(jiān)測的試驗實例

以某水電站水輪發(fā)電機軸承為試驗對象，長期不間斷監(jiān)測其工作情況。水輪機型號為HL220-LJ250，采用29492型推力調(diào)心滾子軸承,內(nèi)徑460 mm，外徑800 mm。在水導軸承附近加裝渦流傳感器測量軸承振動信號，采樣頻率為512 Hz，采樣長度為1 024。不斷記錄得到一列信號數(shù)據(jù)，發(fā)現(xiàn)一水輪機軸承的振動峰值及均方根值等觀察監(jiān)測指標從2014年6月5日后逐漸增大，測得故障振動信號如圖3a所示，由于測量環(huán)境惡劣，測量受強烈水流沖擊的影響，故障信息基本無法識別。

圖3 水輪機軸承沖擊故障信號的提取

為了確定故障源并驗證本研究效果，采用基于DT-CWT和中值濾波的降噪方法對此振動信號進行降噪處理，結(jié)果如圖3b所示，原始振動信號中摻雜的由水流沖擊產(chǎn)生的強噪聲被有效地剔除。與該軸承正常振動信號(圖3c)對比可以看出，該水輪機出現(xiàn)了周期性的振動沖擊。

將軸承取出并拆套，發(fā)現(xiàn)其內(nèi)圈產(chǎn)生了局部損傷，葉輪轉(zhuǎn)動時滾動體轉(zhuǎn)動到損傷部位時產(chǎn)生明顯的沖擊力，導致該機械故障信號的產(chǎn)生。故障信號的沖擊頻率為9.51 Hz，與圖中監(jiān)測的頻率正好吻合。更換軸承后重新安裝，發(fā)電機恢復正常運轉(zhuǎn)，信號恢復正常。

4 結(jié)束語

仿真試驗和應(yīng)用實例分析表明，基于DT-CWT和中值濾波的降噪方法可獲得比常規(guī)小波降噪以及CWT降噪更高的信噪比，不僅能抑制高斯白噪聲，去除偶然的短暫脈沖噪聲，還能有效地提取有用的周期性沖擊信號，但對單個偶然短暫的有用脈沖信號的提取卻無能為力。

機械故障通常表現(xiàn)為周期性的沖擊特征，基于DT-CWT和中值濾波的降噪方法可以消噪后保留和提取軸承、齒輪等故障的周期性沖擊信號特征信息，根據(jù)其周期性沖擊所對應(yīng)的特征頻率和信息對其進行簡便故障識別，比常用的頻譜分析方法更為簡便，為進一步的頻譜分析提供基礎(chǔ)。