楊尚君，張 峰，石現(xiàn)峰

(西安工業(yè)大學(xué) 電子信息工程學(xué)院，西安 710021)

?

旋轉(zhuǎn)機(jī)械振動信號的小波域維納濾波去噪

楊尚君，張 峰，石現(xiàn)峰

(西安工業(yè)大學(xué) 電子信息工程學(xué)院，西安 710021)

為了去除旋轉(zhuǎn)機(jī)械振動信號采集傳輸過程中混入的噪聲干擾，文中基于維納濾波和小波閾值濾波，通過建立旋轉(zhuǎn)機(jī)械振動信號采集模型，結(jié)合振動信號濾波要求，提出了旋轉(zhuǎn)機(jī)械振動信號的小波域維納濾波算法.利用工業(yè)現(xiàn)場旋轉(zhuǎn)機(jī)械實(shí)測振動信號，對該算法進(jìn)行仿真.結(jié)果表明：該算法保持了旋轉(zhuǎn)機(jī)械振動信號的線性相位特性，濾波后信號未產(chǎn)生明顯的幅度失真；小波域維納濾波的均方誤差小于維納濾波和小波閾值濾波，去噪效果優(yōu)于維納濾波和小波閾值濾波.

旋轉(zhuǎn)機(jī)械；振動；小波域維納濾波；線性相位

旋轉(zhuǎn)機(jī)械故障檢測采用數(shù)字信號處理的方法對實(shí)際測量的振動信號進(jìn)行分析，用以參數(shù)檢測和質(zhì)量評價(jià).在采集及傳輸?shù)倪^程中，振動信號不可避免的混入噪聲干擾，對振動信號的濾波處理既要取得較好的效果，也要保證振動信號的均衡相位，以便根據(jù)這些特性,來進(jìn)行狀態(tài)檢測和故障診斷，應(yīng)用于旋轉(zhuǎn)機(jī)械振動信號的濾波處理當(dāng)中.

傳統(tǒng)的有限沖擊響應(yīng)(Finite Impulse Response，F(xiàn)IR)濾波器和無限沖擊響應(yīng)(Infinite Impulse Response，IIR)濾波器，兩者的濾波算法濾波效果和線性相位之間難以達(dá)到均衡[1]；文獻(xiàn)[2]采用維納濾波對非平穩(wěn)振動信號進(jìn)行處理，研究表明未達(dá)到預(yù)期效果；文獻(xiàn)[3]采用循環(huán)維納濾波對振動信號進(jìn)行周期性的分段處理，每段采用維納濾波方法，有效的去除了自適應(yīng)噪聲，但是循環(huán)維納濾波算法復(fù)雜度大；文獻(xiàn)[4]基于離散余弦變換(Discrete Cosine Transform，DCT)算法，保留了信號的部分離散余弦變換域的點(diǎn)數(shù)，實(shí)現(xiàn)了數(shù)據(jù)的壓縮，但是DCT濾波算法對于大數(shù)據(jù)的壓縮存在著數(shù)據(jù)丟失，濾波效果差的問題.針對振動信號濾波的敏感相位和濾波效果的問題，文中將維納濾波和小波閾值濾波相結(jié)合，通過建立旋轉(zhuǎn)機(jī)械振動信號采集模型，結(jié)合振動信號濾波要求，提出旋轉(zhuǎn)機(jī)械振動信號的小波域維納濾波算法.利用工業(yè)現(xiàn)場旋轉(zhuǎn)機(jī)械實(shí)測振動信號，對該算法進(jìn)行仿真，以期滿足濾波要求和線性相位，以適用于其他一維含噪信號的處理.

1 振動信號采集

旋轉(zhuǎn)機(jī)械振動信號的采集模型如圖1所示.圖1中ω為旋轉(zhuǎn)機(jī)械轉(zhuǎn)軸轉(zhuǎn)動的角頻率， v為轉(zhuǎn)軸轉(zhuǎn)動的線速度，理想情況下振動信號可表示為

y=A+Bcos(ωt+φ)

(1)

式中：y為理想振動信號;A為振動信號的直流分量;B為振動信號交流分量的最大值;ω為旋轉(zhuǎn)機(jī)械轉(zhuǎn)軸轉(zhuǎn)動的角頻率;φ為振動信號的初始相位;t為振動信號的時長.

圖1 傳感器測量軸振動信號的采集模型

Fig.1 Collect model of the output signal

of the transducer measuring shaft

信號處理前將其直流分量去除，隨機(jī)變量φ服從在[0,2π]區(qū)間的均勻分布，根據(jù)平穩(wěn)隨機(jī)過程的定義，有

E(y)=E[Bcos(ωt+φ)]=

(2)

Ry(t1,t2)=E[y(t1)y(t2)]=

(3)

式中：E(·)為數(shù)學(xué)期望；Ry(·,·)為自相關(guān)函數(shù)；t1，t2為不同時刻的信號時長.自相關(guān)函數(shù)僅與時間間隔(t2-t1)有關(guān)，理想狀況下振動信號為平穩(wěn)信號.由于現(xiàn)場采集過程中，振動信號不可避免會受到環(huán)境和設(shè)備自身噪聲的干擾，噪聲是非平穩(wěn)的，故含噪的振動信號是非平穩(wěn)振動信號.

根據(jù)工業(yè)現(xiàn)場的實(shí)際情況，對旋轉(zhuǎn)機(jī)械實(shí)際振動信號進(jìn)行實(shí)際采集.數(shù)據(jù)采集的相關(guān)參數(shù)如下：32倍頻采樣，采樣頻率為1 600 Hz，每通道連續(xù)采集128點(diǎn).較為理想情況下，振動信號的時域波形如圖2所示.

圖2 理想情況振動信號時域波形

Fig.2 Time domain waveform of vibration signal in ideal condition

為研究去噪性能，對信號加入非平穩(wěn)隨機(jī)噪聲，噪聲點(diǎn)數(shù)為達(dá)到與采樣點(diǎn)數(shù)匹配，故取128點(diǎn)，加噪后的振動信號如圖3所示.

圖3 加噪處理后的振動信號時域波形

Fig.3 Time domain waveform of vibration signal after adding noise

2 振動信號小波域維納濾波算法

理想狀況下，振動信號各次諧波的譜峰位置應(yīng)出現(xiàn)在50 Hz的整倍頻處，譜峰位置包含了轉(zhuǎn)軸運(yùn)行狀態(tài)的有用信息，因此譜峰的準(zhǔn)確性直接影響了后期的故障診斷.而初始相位的偏移會導(dǎo)致振動信號后期譜估計(jì)中譜峰的分裂或偏移，因此對于振動信號的濾波，單位沖擊響應(yīng)應(yīng)具有較高的線性相位[5].振動信號頻譜圖的方差和分辨力性能也直接受到噪聲的影響，因此為了獲得較好性能的方差和分辨力，應(yīng)盡可能的對振動信號的噪聲進(jìn)行去除.經(jīng)典的FIR濾波器和IIR濾波器在濾波效果和線性相位方面均難以滿足要求，故需要引入現(xiàn)代的維納濾波算法來進(jìn)行處理.根據(jù)最小均方誤差準(zhǔn)則[6],提出一種針對平穩(wěn)過程的最優(yōu)估計(jì)器.

2.1 振動信號的維納濾波算法

假定觀測信號模型為

x(n)=s(n)+w(n)

(4)

式中：s(n)為真實(shí)信號；w(n)為加性高斯白噪聲，其分布為w(n)～N(0,δ2).根據(jù)FIR濾波器準(zhǔn)則，有

(5)

式中： h(k)為FIR濾波器的單位取樣響應(yīng)；N為濾波器階次；y(n)為濾波后信號.

維納濾波算法原理如圖4所示.其中Z-1表示Z變換，e(n)表示真實(shí)信號與濾波后信號的誤差.

圖4 Wiener濾波算法原理圖

Fig.4 Wiener filtering schematic

均方誤差為

(6)

若使均方誤差最小，應(yīng)滿足維納霍夫方程

(7)

式中：φxx(m)為觀測信號x(n)的自相關(guān)函數(shù)；φsx(m)為真實(shí)信號s(n)與觀測信號x(n)的互相關(guān)函數(shù)；m為遍歷系數(shù).維納濾波算法適用于平穩(wěn)隨機(jī)信號[7]，旋轉(zhuǎn)機(jī)械振動信號在傳輸過程中不可避免的受到非平穩(wěn)噪聲的影響，因此，維納濾波算法對旋轉(zhuǎn)機(jī)械振動信號去噪具有局限性.階次為64的維納濾波實(shí)驗(yàn)結(jié)果如圖5所示.

圖5 維納64階次濾波振動信號時域波形

Fig.5 Time domain waveform of vibration signal for 64 order Wiener filtering

通過圖5可以看出，振動信號的維納濾波算法保持了信號的初始相位均衡，但濾波效果較差，達(dá)不到后期信號處理的要求.

2.2 振動信號小波閾值濾波及小波變換

根據(jù)小波閾值濾波，對含噪聲信號進(jìn)行正交小波變換.選擇合適的小波基函數(shù)和分解小波層數(shù)，對含噪信號進(jìn)行正交小波分解，得到對應(yīng)的小波分解系數(shù)，其中包含了低頻系數(shù)和高頻系數(shù).選擇合適的閾值，對分解后的系數(shù)進(jìn)行閾值處理.每一層小波系數(shù)再進(jìn)行量化處理.進(jìn)行小波反變換.將閾值處理后的小波系數(shù)進(jìn)行重構(gòu)，得到小波閾值濾波后的信號.實(shí)驗(yàn)選用軟閾值函數(shù)進(jìn)行處理[8]，數(shù)學(xué)表達(dá)式為

δ(σ)=sgn(σ)(|σ|-λ),|ω|>λ

(8)

式中：σ為小波系數(shù)；λ為閾值；δ(σ)為被估計(jì)信號小波系數(shù)的估計(jì)值.實(shí)驗(yàn)選取最小極大方差閾值，則所選的閾值為

(9)

根據(jù)選用的Coif5小波基，分解2層，小波閾值濾波的結(jié)果如圖6所示.

從圖6中可以看出，直接進(jìn)行小波閾值去噪處理的信號取得了較好的平滑特性，但濾波后信號產(chǎn)生了失真，不能作為后期的信號處理對象.

利用Haar小波作為小波基將信號從時域轉(zhuǎn)化為小波域[9].Haar小波基函數(shù)為

(10)

圖6 小波閾值濾波振動信號時域波形

Fig.6 Time domain waveform of vibration signal for wavelet threshold filtering

對于非平穩(wěn)過程，功率譜密度與頻率的冪成反比的，振動信號在經(jīng)小波變換后能夠，不同尺度間較強(qiáng)的相關(guān)性可有效去除，可以認(rèn)為非平穩(wěn)信號在經(jīng)過小波變換后起到了信號的白化作用[10]，滿足上述結(jié)果的條件須進(jìn)行正交小波變換，Haar小波作為簡單的正交函數(shù)，將振動信號從時域轉(zhuǎn)化到小波域選用Haar小波，降低了信號的非平穩(wěn)特性的同時，保留了信號的有用信息[11-12].

2.3 基于小波變換的維納濾波算法

根據(jù)維納濾波的原理，在構(gòu)建維納霍夫方程前需已知加噪信號和期望信號，利用小波閾值去噪對加噪振動信號進(jìn)行簡單的閾值去噪處理，將原始振動信號作為期望信號，兩者同時變換到小波域進(jìn)行維納濾波處理[13].小波域維納濾波算法過程如下：① 將加噪的振動信號進(jìn)行小波閾值去噪進(jìn)行預(yù)處理，得到信號，用于構(gòu)建小波域的維納霍夫方程；② 分別將原始振動信號(作為期望信號)和小波閾值去噪信號兩者分別利用Haar小波進(jìn)行小波變換，分別提取兩者的近似分量和細(xì)節(jié)分量；③ 利用兩者細(xì)節(jié)分量構(gòu)建維納霍夫方程，對小波閾值去噪信號的細(xì)節(jié)分量進(jìn)行濾波處理，利用兩者近似分量構(gòu)建維納霍夫方程，對小波閾值去噪信號的近似分量進(jìn)行濾波處理；④ 利用小波反變換函數(shù)對上述信號處理結(jié)果進(jìn)行反變換，其利用的小波基仍為Haar小波.

3 仿真及分析

小波域32階次維納濾波算法信號處理結(jié)果如圖7所示.通過圖7和圖5的對比，小波域維納濾波算法具有較好的線性相位特性，振動信號初始相位沒有發(fā)生明顯偏移.在小波域維納濾波的濾波階數(shù)小于維納濾波階數(shù)的同時，小波域維納濾波能取得更好的濾波效果.根據(jù)圖7和圖6的比較，利用小波閾值濾波算法保證了振動信號濾波后的平滑特性，該算法可用于實(shí)際的振動信號去噪環(huán)境中.由于計(jì)算量較大，該算法可做振動信號的離線分析處理.

圖7 小波域維納濾波振動信號時域波形

Fig.7 Time domain waveform of vibration signal for Wiener filtering in wavelet domain

分別計(jì)算各種算法濾波后的結(jié)果與原始振動信號做均方誤差的求解，進(jìn)行算法性能定量分析和對比，得出數(shù)據(jù)見表1.

表1 算法性能定量分析

Tab.1 Quantitative analysis of algorithm performance

濾波算法均方誤差維納濾波0.0097小波閾值濾波0.0075小波域維納濾波0.0032

從表1可看出，小波域維納濾波的均方誤差要小于維納濾波和小波閾值濾波，小波域維納濾波算法性能優(yōu)于各種單獨(dú)算法，小波域維納濾波算法濾波后的信號更接近原始的振動信號.

利用小波域維納濾波算法對振動信號進(jìn)行去噪處理，對其結(jié)果進(jìn)行頻譜分析，采取周期圖法，功率譜密度(Power Spectral Density,PSD)估計(jì)圖如圖8所示.f為頻率,DPS為功率譜密度.

圖8 濾波后振動信號的頻譜圖

Fig.8 Power spectrum of vibration signal after filtering

根據(jù)振動信號的采集模型，信號每周期采樣32點(diǎn)，采樣頻率為1 600 Hz，因此信號的固有頻率為50 Hz.從圖8中可以看出，濾波后振動信號的譜峰處于50 Hz整倍頻處，功率譜圖中的譜峰具有較好的尖銳程度，分辨率性能較好，具有較強(qiáng)辨別信號的能力，且譜峰沒有發(fā)生偏移或分裂的現(xiàn)象.因此濾波后的振動信號適用于后期的旋轉(zhuǎn)機(jī)械故障檢測.

4 結(jié) 論

1) 將振動信號從時域轉(zhuǎn)化到小波域，降低了振動信號的非平穩(wěn)特性，小波域維納濾波算法保留了維納濾波算法的線性相位特性，信號原有的初始相位未發(fā)生偏移，濾波后的信號幅值和相位未產(chǎn)生失真.

2) 小波域維納濾波去噪性能優(yōu)于維納濾波和小波閾值濾波，均方誤差小于維納濾波和小波閾值濾波.

3) 小波域維納濾波后的振動信號，功率譜圖中的譜峰處于50 Hz整倍頻處，分辨率性能好，辨別信號的能力優(yōu)于維納濾波，譜峰沒有發(fā)生偏移或分裂，適用于后期的旋轉(zhuǎn)機(jī)械故障檢測.

[1] 張峰,石現(xiàn)峰,張學(xué)智.數(shù)字信號處理原理及應(yīng)用[M].北京:電子工業(yè)出版社,2012.

ZHANG Feng,SHI Xianfeng,ZHANG Xuezhi.Principle and Application of Digital Signal Processing[M].Beijing:Electronics Industry Press,2012.

(in Chinese)

[2] 明陽.基于循環(huán)平穩(wěn)和盲源分離的滾動軸承故障特征提取方法研究[D].上海：上海交通大學(xué),2013.

MING Yang.Study on Cyclostationarity and Blind Source Separation-Based Rolling Element Bearing Fault Feature Extraction[D].Shanghai:Shanghai Jiao Tong University,2013.(in Chinese)

[3] 明陽,陳進(jìn),董廣明.基于循環(huán)維納濾波器和包絡(luò)譜的軸承故障診斷[J].振動工程學(xué)報(bào),2010,23(5):537.

MING Yang,CHEN Jin,DONG Guangming.Rolling Bearing Fault Diagnosis Based on Cyclic Wiener Filtering and Envelop Spectrum[J].Journal of Vibration Engineering,2010,23(5):537.(in Chinese)

[4] 管博,胡勁松.基于DCT的旋轉(zhuǎn)機(jī)械振動信號無相移濾波方法研究[J].汽輪機(jī)技術(shù),2007,49(4):285.

GUAN Bo,HU Jinsong.Research on the Method of the Vibration Signal of Rotating Machinery Based on DCT[J].Turbine Technology,2007,49(4):285.

(in Chinese)

[5] RINALDO R,VISETTI A.Relations Between the K -Dimensional Unconstrained and Linear Phase FIR Wiener Filters[J].Multidimensional Systems and Signal Processing,1995,6(4):347.

[6] 李楠.一種小波變換與維納濾波結(jié)合的語音抗噪研究[J].電聲技術(shù),2007,31(5):46.

LI Nan.Study on the Speech Antinoise Based on Wavelet Transform and Wiener Filtering[J].Electro Acoustic Technology,2007,31(5):46.(in Chinese)

[7] XU P,MIAO Q,TANG X.A Denoising Algorithm via Wiener Filtering in the Shearlet Domain[J].Multimedia Tools and Applications,2014,71(3):1529.

[8] 楊玲香,姚斌.基于Canny算子的小波域維納濾波去噪方法[J].兵團(tuán)教育學(xué)院學(xué)報(bào),2009,19(5):35.

YANG Lingxiang,YAO Bin.Denoisng Method via Local Wiener Filtering in Wavelet Domain Based on Canny Operator[J].Journal of Bingtuan Education Institute,2009,19(5):35.(in Chinese)

[9] ZHONG K,LEI X,LI S.Wiener Filter Based Channel Estimation for High-Speed Communication Environments[J].Wireless Personal Communications,2013,69(4):1819.

[10] 李東兵,李國平,滕國偉,等.一種基于小波和Wiener相結(jié)合圖像去噪新方法[J].電視技術(shù),2013,37(13):26.

LI Dongbing,LI Guoping,TENG Guowei,et al.New Method of De-noising Based on Wavelet and Wiener Filtering[J].Video Engineering,2013,37(13):26.(in Chinese)

[11] RAO R,MURTHY R,VENKATA R.Speech Enhancement Using Cross-correlation Compensated Multi-band Wiener Filter Combined with Harmonic Regeneration[J].Journal of Signal and Information Processing,2011(2):117.

[12] 胡耀斌,陳艾華,張春良.小波分析與維納濾波相結(jié)合的消噪方法研究[J].電力系統(tǒng)通信,2006,27(162):42.

HU Yaobin,CHEN Aihua,ZHANG Chunliang.Research of Denoising Technology about Wavelet Analysis with Wiener Filter[J].Communication of Power System,2006,27(162):42.(in Chinese)

[13] SINDELAR V,NIKOLAUS G.An Adaptation of the Wiener Filter Suitable for Analyzing Images of Isolated Single Particles[J].Journal of Structural Biology,2011,176(1):60.

(責(zé)任編輯、校對 張 超)

Wavelet Domain Denoising for Rotating Machinery Vibration Signal Based on Wiener Filtering

YANGShangjun,ZHANGFeng,SHIXianfeng

(School of Electronic Information Engineering,Xi’an Technological University,Xi’an 710021,China)

In order to remove the noise in rotating machinery vibration signal during the acquisition and transmission process, the algorithm of Wiener filtering in wavelet domain is proposed to meet the requirements of vibration signal filtering,based on Wiener filtering and wavelet threshold filtering,through the establishment of rotating machinery vibration signal acquisition model.Using the rotating machinery vibration signal which is actual measurement in the industrial field, the algorithm is simulated.The result shows: This algorithm can maintain the linear phase characteristics for rotating machinery vibration signal,the filtered signal doesn't produce the amplitude distortion; The mean square error of Wiener filtering in wavelet domain is less than that of Wiener filtering and wavelet threshold filtering.The denoising result is better than Wiener filtering and wavelet threshold filtering.

rotating machinery;vibration;Wiener filtering in wavelet domain;linear phase

10.16185/j.jxatu.edu.cn.2016.10.015

2016-01-06

國家自然科學(xué)基金項(xiàng)目(60972095；61271362)；陜西省科技廳自然基金項(xiàng)目(2012JQ8008)；陜西省教育廳科技專項(xiàng)(2012JK0545)

楊尚君(1991-)，男，西安工業(yè)大學(xué)碩士研究生.

張 峰(1979-),男,西安工業(yè)大學(xué)副教授,主要研究方向?yàn)樾盘柵c信息處理、嵌入式技術(shù)應(yīng)用、自動控制及數(shù)據(jù)庫技術(shù)運(yùn)用.E-mail:zf_zx963@163.com.

TN911.4

A

1673-9965(2016)10-0856-05