強噪聲背景下頻率加權(quán)能量算子和變分模態(tài)分解在軸承故障提取中的應(yīng)用

徐元博， 蔡宗琰， 胡永彪， 丁 凱

(長安大學(xué)道路施工技術(shù)與裝備教育部重點實驗室， 陜西 西安 710064)

引 言

對于大部分機械零部件來說，無論是軸承、傳動軸或者其他零部件，當(dāng)它們出現(xiàn)故障時，都會出現(xiàn)異常振動。機械設(shè)備正常工作時表現(xiàn)出的振動量在時域和頻域都具有一定的數(shù)據(jù)特征并且具有一定的規(guī)律性，因此一旦檢測出實際振動數(shù)據(jù)的特征量不符合正常工作情況下的規(guī)律性時，即可判定設(shè)備的運行有異常情況或故障隱患[1]，因此利用振動信號診斷法來判斷故障特征已經(jīng)得到了廣泛應(yīng)用。機械故障產(chǎn)生的振動信號一般具有很強的非線性與非平穩(wěn)性特征，同時伴隨著噪聲。經(jīng)驗?zāi)B(tài)分解(Empirical Mode Decomposition，EMD)作為一種非線性和非平穩(wěn)信號的處理方法[2]，得到了很多學(xué)者的關(guān)注，并在機械故障診斷領(lǐng)域[3-5]，得到了較好的應(yīng)用，但由于該算法的篩分特性，使其很大程度上要依賴于極值點的找尋和停止法則，同時對背景噪聲有著較強的敏感性，所以易出現(xiàn)模態(tài)混疊現(xiàn)象，使該算法魯棒性較差。為了克服這一缺點，N E Huang等在其基礎(chǔ)上提出了集成經(jīng)驗?zāi)B(tài)分解[6](Ensemble Empirical Mode Decomposition，EEMD)，該算法是在整個時頻空間中加入均勻分布的白噪聲。利用這一特點，EEMD在克服模態(tài)混疊現(xiàn)象上取得良好效果，同樣在機械故障領(lǐng)域得到了廣泛應(yīng)用[7-8]。EEMD作為噪聲輔助型算法，在分解后的信號中必然會留下殘余白噪聲，理論上經(jīng)過無限次集成平均后，可以完全消除殘余白噪聲，但在實際應(yīng)用中這是不可能實現(xiàn)的，并且在加入白噪聲數(shù)量大的情況下，平均集成次數(shù)也隨之增加，這就造成數(shù)據(jù)處理要花費大量時間。基于此，本文采用一種新的自適應(yīng)信號模態(tài)分解方法，變分模態(tài)分解[9](Variational Mode Decomposition, VMD)，該方法由Dragomiretski-y等于2014年提出，其最大的特點就是摒棄了EMD和EEMD算法中的遞歸篩分思想，而是將對模態(tài)的估計轉(zhuǎn)變?yōu)樽兎謫栴}的求解, 在頻域不斷更新各模態(tài)及其中心頻率, 最后各模態(tài)經(jīng)傅里葉逆變換到時域[10]。同時VMD算法融合了維納濾波組，不僅可以有效地處理背景噪聲、增強該算法的魯棒性，而且縮短了程序的運行時間。文獻[11-13]已經(jīng)將該方法應(yīng)用到旋轉(zhuǎn)機械的故障診斷當(dāng)中，取得了良好的效果。因此本文將VMD和一種改進的能量算子方法—頻率加權(quán)能量算子結(jié)合(Frequency-Weighted Energy Operator, FWEO)應(yīng)用于振動機械的故障診斷中。FWEO作為能量算子中的一種方法，不同于希爾伯特變換(Hilbert Transform, HT)，它不屬于積分變換，而是基于求導(dǎo)運算，因此計算復(fù)雜性低，最主要的是，相對于傳統(tǒng)的Teager能量算子(Teager Energy Operator, TEO)，在面對噪聲與其他頻率成分干擾時，它具有更強的魯棒性[14]。最后本文將該方法應(yīng)用在仿真信號和真實軸承故障信號上，并與其他方法進行對比，驗證了該方法具有更加優(yōu)秀的性能，也為其他科研領(lǐng)域提供了新的參考。

1 VMD基本原理

VMD算法的目標(biāo)就是將一個實數(shù)信號f(t)分解為多個互不相關(guān)的稀疏性子信號，也就是本征模函數(shù)，假設(shè)每個模態(tài)uk都具有中心頻率和有限帶寬，約束條件為各模態(tài)之和等于輸入信號f(t)，且每個模態(tài)的估計帶寬之和最小，步驟如下[9]：

(1) 通過對每個模態(tài)函數(shù)uk進行希爾伯特變換(Hilbert Transform)得到解析信號的單邊譜，如下

(1)

式中δ(t)表示狄利雷克分布函數(shù)；*表示卷積。

(2) 每個經(jīng)過希爾伯特變換后的模態(tài)函數(shù)同ejωkt進行混合，使每個模態(tài)函數(shù)的頻譜移頻到基帶上，式子如下

(2)

式中ωk代表中心頻率。

(3) 通過梯度的二范數(shù)的平方對其模態(tài)的帶寬進行估計，最終的公式如下

(3)

(4)

式中 {uk}={u1,u2,…,uk}代表K個模態(tài)函數(shù);{ωk}={ω1,ω2,…,ωk}代表每個模態(tài)函數(shù)的中心頻率。

為了使約束問題轉(zhuǎn)化為非約束問題，引入二次懲罰因子和拉格朗日乘子，二次懲罰因子用來保證信號重建后的保真度，朗格朗日乘子用來保證其約束的嚴格性，其擴展的拉格朗日表達式如下

(5)

使用Parseval/Plancherel傅里葉等距變換, 將式(5)中的{uk}和{ωk}的二次優(yōu)化問題轉(zhuǎn)變到頻域, 分別得到各模態(tài)的頻域上的更新，具體表達式如下:

(6)

(7)

2 頻率加權(quán)能量算子基本原理

傳統(tǒng)的TEO能量算子一般稱之為Teager能量算子[15]，方法理念易于理解，假設(shè)有一組連續(xù)信號

x(t)=acos(ωt)

(8)

式中a代表幅值；ω代表瞬時頻率。

則一個能量算子定義如下

(9)

它的離散形式可以表示如下

ψ(x(n))=x2(n)-x(n-1)x(n+1)

(10)

將式(8)和其導(dǎo)數(shù)形式帶入式(9)得到以下式子：

(11)

將式(11)中兩個式子進行簡單運算可以得到瞬時頻率和幅值：

(12)

對于一般的解調(diào)信號，通常使用解析信號的形式來對其進行定義，如下式

X(t)=x(t)+jH[x(t)]=Aejφ(t)

(13)

式中H[x(t)]代表信號x(t)的HT；A(t)代表信號的瞬時幅值；φ(t)代表信號的瞬時相位。而信號的能量則是由信號的幅值的平方?jīng)Q定[16-17]，其表達式為

(14)

(15)

令x(t)=Acos(ωt+φ)，則：

(16)

(17)

將式(15)和(16)帶入式(14)，得出

Γ[x(t)]=A2ω2

(18)

可以看到新的能量算子也具有相類似的性質(zhì)。

對于一個離散信號x(n)，則它的能量算子為

h2(n+1)+h2(n-1)]+

h(n+1)+h(n-1)]

(19)

3 實驗與討論

為了體現(xiàn)VMD-FWEO在強背景噪聲干擾下仍能較好地提取出故障頻率，本部分采用一組軸承故障模擬信號來進行驗證，并與較多文獻中提出的EMD-HT和EEMD-TEO等[18-22]故障提取方法進行對比，從而體現(xiàn)其優(yōu)越性。

實驗信號采用一組軸承故障模擬仿真信號，如下式

(20)

為了體現(xiàn)FWEO在強噪聲和有其他頻率干擾下的優(yōu)越性能，模擬環(huán)境采用強噪聲和其他頻率干擾下的復(fù)雜背景環(huán)境，混入噪聲為信噪比SNR=-10 dB的高斯白噪聲，4個干擾頻率分別為5，27，43和810 Hz。圖1分別為沖擊信號和混合后的復(fù)雜信號時域圖。

現(xiàn)分別用HT，TEO和FWEO三種方法做出各自的包絡(luò)譜和能量譜，如圖2所示。從圖2(a)中可以看到原始信號的頻譜圖中只含有干擾頻率而無故障頻率100Hz。通過圖2(b)，(c)和(d)對比發(fā)現(xiàn)，雖然HT可以勉強提取出故障特征頻率，但其幅值非常小，也未出現(xiàn)故障頻率的倍頻，且在其包絡(luò)譜中仍可以看到干擾頻率，因此它的提取效果可以說是失敗的；而TEO的能量譜，在噪聲和其他頻率成分干擾下，其提取故障特征頻率的能力受到很大限制，可以看出，在能量譜中幾乎都是噪聲頻率而無故障特征頻率；顯然FWEO的效果要強于前兩種方法，從其能量譜中可以看出，不僅譜線頂點幅值非常高，而且沒有干擾頻率的譜線，基本都為故障頻率及其倍頻，因此說明在強噪聲背景下FWEO的提取故障特征頻率的能力更加突出。值得注意的是此次實驗是在沒有前處理的情況下進行。

圖1 沖擊信號和混合信號時域圖Fig.1 Time domain of impact signal and composite signal

圖2 HT, TEO和FWEO性能對比Fig.2 Comparison of the performances among HT, TEO and FWEO

圖3為經(jīng)過前處理后，EMD-HT，EEMD-TEO與VMD-FWEO的對比結(jié)果，從圖3(a)和(b)可以看出，在經(jīng)過EMD，EEMD的前處理，過濾掉大量噪聲和干擾頻率后，HT和TEO可以將故障特征頻率提取出來，從圖3(c)可以看出，VMD-FWEO的效果還是要優(yōu)于前兩種方法。同時對比圖3(c)和圖2(d)后發(fā)現(xiàn)，在無需前處理的情況下， FWEO雖然可以提取出故障特征頻率，但仍然可以看到有較多的噪聲頻率存在于FWEO能量譜中。而經(jīng)過VMD處理后，可以看出FWEO能量譜更加的清晰干凈，噪聲頻率干擾明顯減少。因此，可以看出VMD-FWEO在強噪聲和多頻率干擾的背景環(huán)境下對軸承故障特征頻率有著較好的提取能力。

圖3 EMD-HT, EEMD-TEO和VMD-FWEO的性能對比Fig.3 Comparison of the performances among EMD-HT, EEMD-TEO and VMD-FWEO

4 振動篩軸承故障應(yīng)用

4.1 試驗?zāi)康呐c安排

為了進一步驗證該方法在實際工況下的優(yōu)越性和實用性，將其應(yīng)用在振動機械中常用的振動篩的軸承故障中。該實驗的目的有二：一、通過實驗驗證該種方法不僅可以在旋轉(zhuǎn)機械上使用，而且可以在振動機械上應(yīng)用；二、振動機械不同于旋轉(zhuǎn)機械，在相同的工況下，它的背景噪聲和其他振動干擾要明顯強于旋轉(zhuǎn)機械，對算法的魯棒性要求更高。因此可以在實際工況下驗證該種方法相較于其他傳統(tǒng)方法的優(yōu)越性。

實驗采用的振動篩如圖4所示，分別對滾動軸承進行內(nèi)環(huán)點蝕、外環(huán)點蝕，點蝕尺寸參數(shù)如表1所示。軸承參數(shù)如表2所示。軸承的內(nèi)外圈故障特征頻率可以通過下式得到：

(21)

(22)

式中fo，fi分別代表外圈和內(nèi)圈故障特征頻率；d代表滾子直徑；D代表滾道的節(jié)徑；Z代表滾動體數(shù)量；α代表接觸角；f代表軸的轉(zhuǎn)動頻率。根據(jù)上式計算出的故障頻率和相關(guān)頻率如表3所示。

圖4 振動篩Fig.4 Vibrating screen

型號位置 寬度深度長度1308內(nèi)環(huán)1.00.2軸承寬度/21308外環(huán)0.70.2軸承寬度/2

表2 軸承參數(shù)

表3 本次實驗相關(guān)特征頻率

針對振動篩的工作特點，本次實驗專門準(zhǔn)備了4塊大偏心塊和4塊小偏心塊，配以單電機和雙電機運轉(zhuǎn)，當(dāng)雙電機運轉(zhuǎn)時兩電機轉(zhuǎn)速相同。將它們進行交叉組合來獲得不同的振動量級，隨著振動量級的增大，背景噪聲也隨之增強，從而采集不同振動量級下的故障數(shù)據(jù)，實驗安排如表4所示。

表4 實驗安排

實驗安排中，同樣采用EMD-HT，EEMD-TEO和VMD-FWEO三種方法進行對比。小偏心單電機振動量級最小，大偏心雙電機振動量級最大。為了更好體現(xiàn)該方法的優(yōu)越性，本次試驗采用外圈振動量級最大(實驗6)工況下和內(nèi)圈振動量級最大(實驗3)工況下進行對比實驗。

同時，為了更好地體現(xiàn)VMD-FWEO及其FWEO單獨使用時的優(yōu)越性，本次實驗還采用了美國西儲大學(xué)軸承實驗驗證中心提供的公共數(shù)據(jù)作為對比，該數(shù)據(jù)由旋轉(zhuǎn)機械設(shè)備中采集，因此噪聲較小。旋轉(zhuǎn)實驗臺如圖5所示。該設(shè)備的軸承參數(shù)如表5所示。相關(guān)頻率參數(shù)如表6所示。

表5 軸承參數(shù)

表6 相關(guān)特征頻率

圖5 旋轉(zhuǎn)實驗臺Fig.5 Experiment rig

4.2 外圈實驗

振動篩與模擬實驗臺的外圈故障信號如圖6和7所示，從圖中可以看出從兩種設(shè)備中提取出的外圈故障頻率都已經(jīng)被噪聲與其他干擾淹沒，只有各自的基頻頻率15.17和29.95 Hz。

現(xiàn)用EMD，EEMD和VMD分別對兩組信號進行分解，分解后的模態(tài)分量，利用K-L散度[23]從中選取與原信號相似度最大的模態(tài)信號，兩種工況下的K-L散度如表7和8所示。

表7 K-L實驗臺數(shù)據(jù)散度值

表8 K-L振動篩數(shù)據(jù)散度值

圖6 振動篩外圈原始信號Fig.6 Original signal of outer ring of vibrating screen

圖7 實驗臺外圈原始信號Fig.7 Original signal of outer ring of experimental rig

從中分別合適的模態(tài)分量做HT包絡(luò)譜，TEO能量譜和FWEO能量譜，其對比如圖8所示，左邊為實驗臺的故障頻譜圖，右邊為振動篩故障頻譜圖。

圖8 實驗臺和振動篩外圈故障特征對比Fig.8 Comparison of faulty feature of outer ring between the vibrating screen and the experimental rig

故障信號經(jīng)過EMD，EEMD和VMD分解處理后，去除了大量的高頻噪聲和系統(tǒng)中的干擾頻率，因此從圖中可以看到，任何一種方法，都可將兩種工況下的外圈故障頻率提取出來。但從細節(jié)上依然可以看出，不管是在模擬實驗臺上還是振動篩上，通過VMD-FWEO得出的能量譜中的故障特征譜線幅值明顯高于其他兩種方法得出的譜線幅值，說明提取故障特征能力較強。

同時，文獻[24-25]提出可將瞬時頻率作為評價算法魯棒性是否優(yōu)劣的一種標(biāo)準(zhǔn)，瞬時頻率譜線波動越大，說明魯棒性越差，易受噪聲干擾，反之亦然。因此可以采用瞬時頻率來判斷三種方法的優(yōu)劣。EMD，EEMD和VMD分解后的效果，瞬時頻率對比如圖9所示。通過對比可以明顯發(fā)現(xiàn)兩個特點：第一，經(jīng)過VMD分解后，不管從模擬實驗臺上還是振動篩上提取的軸承故障信號，其模態(tài)分量的故障頻率譜線的波動都明顯小于其他兩種方法；其次，從圖中可以清楚看到，當(dāng)使用模擬實驗臺上提取的故障信號時，EMD和EEMD分解的模態(tài)分量的頻率譜線波動較小，而采用振動篩上提取的故障信號時，其波動明顯增大，而VMD分解的模態(tài)分量的故障頻率譜線在兩種工況下，除個別時點外，波動并無大幅度的變化。因此，從這點可以說明基于EMD和EEMD提取軸承故障的方法易受噪聲影響，而基于VMD的方法有更好的魯棒性，抗噪的能力較強。

圖9 瞬時頻率對比Fig.9 Comparison of instantaneous frequency

4.3 內(nèi)圈實驗

從外圈實驗可以看出，在EMD，EEMD和VMD的前處理幫助下，HT，TEO和FWEO三種方法都能將故障特征頻率提取出來。因此，為了比較三種方法在真實數(shù)據(jù)下單獨使用時，提取故障特征頻率能力的優(yōu)劣，內(nèi)圈故障實驗在不進行前處理的工況下進行。

現(xiàn)用Hilbert，TEO和FWEO的三種方法對兩種工況下的軸承內(nèi)圈故障信號進行處理，如圖10所示。從圖10中可以看出，三種方法在沒有前處理的情況下，均能將旋轉(zhuǎn)模擬實驗臺的內(nèi)圈故障特征頻率及其倍頻提取出來。而對于在強噪聲背景下提取出的振動篩軸承故障信號，僅有FWEO能提取出故障特征頻率及其倍頻，而其他兩種方法都是失敗的。因此通過這個對比實驗，可以發(fā)現(xiàn)，傳統(tǒng)的HT和TEO能量算子易受噪聲和其他干擾頻率的影響，在復(fù)雜的背景環(huán)境下不適合單獨使用，必須與有效的前處理方法結(jié)合使用才能提取出軸承故障特征頻率。而對于FWEO能量算子，在兩種工況下都可將故障特征頻率提取出來，說明該方法具有更好的魯棒性，不易受到外部環(huán)境的影響。同時，在單獨使用FWEO時，可以提高運算效率，大大節(jié)省程序運行時間。

5 結(jié) 論

通過仿真信號和真實的振動篩軸承故障實驗，有幾點發(fā)現(xiàn)：

(1) 實驗結(jié)果對比說明，傳統(tǒng)的HT和TEO能量算子可以在背景噪聲和振動干擾較小的情況下單獨使用而不需要前處理，但噪聲一旦加大，其提取故障特征的能力明顯下降甚至是失敗的，因此需要同一些有效的前處理方法聯(lián)立使用。但對于FWEO能量算子來說，當(dāng)噪聲增大時，它的故障特征提取能力依然較強，再不需要前處理的情況下，依然可以較好的將故障特征頻率提取出來。

(2) EMD-HT，EEMD-TEO和VMD-FWEO三種方法通過實驗看出，在兩種工況下，都可有效提取出故障特征頻率，但通過其瞬時頻率的波動程度，可以看出，VMD-FWEO對于噪聲有著更強的抗干擾性能。因此將其應(yīng)用到軸承故障領(lǐng)域或其他機械故障診斷方向是可行的。

圖10 內(nèi)圈故障特征對比Fig.10 Comparison of fault characteristic of inner ring

(3) 因本次實驗通過改變不平衡量大小獲得振動干擾，所以未模擬出實際工況下的復(fù)雜的沖擊振動；另外將其應(yīng)用于同時存在內(nèi)外圈故障的復(fù)合故障情況時，都將是我們的下一步研究工作。

[1] 普亞松，郭德偉，張文斌．故障診斷技術(shù)在煤礦機械設(shè)備中的應(yīng)用[J]．工礦自動化，2015，41(4)：36—39.

Pu Yasong, Guo Dewei, Zhang Wenbin. Application of fault diagnosis technologies in coal mine machinery[J]. Industry and Mine Automation, 2015, 41(4): 36—39.

[2] Huang N E, Shen Z, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings of the Royal Society of London, 1998, 454(A): 903-995.

[3] 吳虎生，呂建新，吳廬山，等．基于EMD 和 SVM的柴油機氣閥機構(gòu)故障診斷[J]．中國機械工程，2010，21(22)：2710—2714.

Wu Husheng, Lü Jianxin, Wu Lushan, et al.Fault diagnosis for diesel valve train based on SVM and EMD[J]. China Mechanical Engineering, 2010, 21(22): 2710—2714.

[4] 冷軍發(fā)，荊雙喜，禹建功，等．EMD與能量算子解調(diào)在提升機齒輪箱故障診斷中的應(yīng)用[J]．煤炭學(xué)報，2013，38(2)：531—535．

Leng Junfa, Jin Shuangxi, Yu Jiangong, et al.Application of fault diagnosis method based on EMD and energy operator demodulation to hoist gearbox[J]. Journal of China Coal Society, 2013, 38(2): 531—535.

[5] 石瑞敏，楊兆建．基于改進 EMD 的多繩摩擦提升機載荷信息特征提取[J]．煤炭學(xué)報，2014，39(4)：782—788．

Shi Ruimin,Yang Zhaojian. Feature extractionfor hoisting load of multiple rope friction hoist based on improved EMD[J]. Journal of China Coal Society, 2014, 39(4): 782—788.

[6] Zhaohua Wu, NORDEN E Huang. Ensemble empirical mode decomposition: a noise-assisted data analysis method[J]. Advances in Adaptive Data Analysis, 2011, 1(01): 1—41.

[7] 張淑清，黃文靜，胡永濤，等．基于總體平均經(jīng)驗?zāi)Ｊ椒纸饨旗睾突旌螾SO-BP算法的軸承故障診斷方法[J]．中國機械工程，2016，27(22)：3048—3054．

Zhang Shuqing, Huang Wenjing, Hu Yongtao, et al. Bearing fault diagnosis method based on EEMD approximate entropy and hybrid PSO-BP algorithm[J]. China Mechanical Engineering, 2016, 27(22): 3048—3054.

[8] ?vokelj M, Zupan S, Prebil I. Multivariate and multiscale monitoring of large-size low-speed bearings using ensemble empirical mode decomposition method combined with principal component analysis[J]. Mechanical Systems & Signal Processing, 2010, 24(4):1049—1067.

[9] Dragomiretskiy K, Zosso D. Variational mode decomposition[J]. IEEE Transactions on Signal Processing, 2013, 62(3): 531—544.

[10] 謝 平，楊芳梅，李欣欣，等．基于變分模態(tài)分解-傳遞熵的腦肌電信號耦合分析[J]．物理學(xué)報，2016，65(11)：(118701-1)—(118701-9)．

Xie Ping,Yang Fangmei, Li Xinxin,et al. Functional coupling analyses of electroencephalogram and electromyogram based on variational mode decomposition transfer entropy[J].Acta Phys. Sin., 2016, 65(11): (118701-1)—(118701-9).

[11] 馬增強，李亞超，劉 政，等.基于變分模態(tài)分解和Teager能量算子的滾動軸承故障特征提取[J].振動與沖擊，2016，35(13)：134—140．

Ma Zengqiang, Li Yachao, Liu Zheng. Rolling bearings’ fault feature extraction based on variational mode decomposition and Teager energy operator[J]. Journal of Vibration and Shock. 2016, 35(13): 134—140.

[12] 劉長良，武英杰，甄成剛．基于變分模態(tài)分解和模糊C均值聚類的滾動軸承故障診斷[J].中國電機工程學(xué)報，2015，35(13)：3358—3366．

Liu Changliang, Wu Yingjie, Zhen Chenggang. Rolling bearing fault diagnosis based on variational mode decomposition and fuzzy C-means clustering[J]. Proceedings of the CSEE, 2015, 35(13): 3358—3366.

[13] 唐貴基，王曉龍．變分模態(tài)分解方法及其在滾動軸承早起故障診斷中的應(yīng)用[J].振動工程學(xué)報，2016，29(4)：639—648．

Tang Guiji, Wang Xiaolong. Variational mode decomposition method and its applicationon incipient-fault diagnosis of rolling bearing[J].Journal of Vibration Engineering, 2016, 29(4): 639—648.

[14] O′Toole JM, Temko A, Stevenson N. Assessing instantaneous energy in the EEG: A non-negative, frequency-weighted energy operator[C]. 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBC, 2014:3288—3291.

[15] Kaiser J F. On a simple algorithm to calculate the 'energy' of a signal[C]. International Conference on Acoustics, Speech, and Signal Processing， IEEE, 2002:381—384.

[16] Loutridis S J. Instantaneous energy density as a feature for gear fault detection[J]. Mechanical Systems and Signal Processing, 2006, 20(5): 1239—1253.

[17] Moustafa W, Cousinard O, Bolaers F, et al. Instantaneous energy and circular statistics applied on ultrasonic signals for very low speed bearing fault diagnosis[C].International Congress on the Engineering of Industrial Risks, Ciri，2014.

[18] 高 強，杜小山，范 虹，等．滾動軸承故障EMD診斷方法研究[J]．振動工程學(xué)報，2007，20(1)：15—18．

Gao Qiang, Du Xiaoshan, Fan Hong, et al. An empirical mode decomposition based method for rolling bearing fault diagnosis[J]. Journal of Vibration Engineering, 2007, 20(1): 15—18.

[19] Tang Y N, Lu C, Hu J M. Rolling bearing fault diagnosis and health assessment based on EMD-TEO[J]. Advanced Materials Research, 2014, 998-999: 470—475.

[20] 王天金，馮志鵬，郝如江，等．基于Teager能量算子的滾動軸承故障診斷研究[J]. 振動與沖擊，2012，31(2)：1—5.

Wang Tianjin, Feng Zhipeng, Hao Rujiang, et al. Fault diagnosis of rolling element bearings based on Teager energy operator[J]. Journal of Vibration and Shock, 2012, 31(2): 1—5.

[21] Feng Z, Zuo M J, Hao R, et al. Ensemble empirical mode decomposition based teager energy spectrum for bearing fault diagnosis[J]. Journal of Vibration and Acoustic-Transaction of the ASME, 2013, 135(3): 031013.

[22] 佘 博，田福慶，湯 健，等．基于Teager能量算子和EEMD的滾動軸承故障診斷方法[J/OL]． 北京工業(yè)大學(xué)學(xué)報，http://kns.cnki.net/kcms/detail/11.2286.T.20170414.0853.030.html．

She Bo, Tian Fuqing, Tang Jian, et al. Fault diagnosis of rolling bearing based on Teager energy operator and EEMD[J]. Journal of Beijing University of Technology, http://kns.cnki.net/kcms/detail/11.2286.T.20170414.0853.030.html．

[23] 韓中和，李文華．基于K-L散度的EMD虛假分量識別方法[J]．中國電機工程學(xué)報，2012，32(11)：112—117．

Han Zhonghe, Li Wenhua. A false component identification method of EMD based on Kullback-leibler Divergence[J].Proceeding of the CSEE, 2012, 32(11): 112—117.

[24] Norden E Huang, Zhaohua Wu, Steven R Long, et al. On instantaneous frequency[J]. Advances in Adaptive Data Analysis, 2011, 1(02):177—229.

[25] Daubechies I, Lu J, Wu H T. Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool[J]. Applied and Computational Harmonic Analysis, 2011, 30(2): 243—261.