第一作者鄭近德男，博士生，1986年生

通信作者程軍圣男，博士，教授，博士生導(dǎo)師，1968年生

廣義經(jīng)驗(yàn)?zāi)B(tài)分解性能分析與應(yīng)用

鄭近德,程軍圣,曾鳴,羅頌榮

(湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國(guó)家重點(diǎn)實(shí)驗(yàn)室，長(zhǎng)沙410082)

摘要：針對(duì)經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition，EMD)的均值曲線采用三次樣條擬合而容易引起包絡(luò)過(guò)沖和不足等缺陷，相關(guān)學(xué)者提出了許多改進(jìn)均值曲線的變種EMD方法，取得了一定的效果。廣義經(jīng)驗(yàn)?zāi)B(tài)分解(Generalized EMD，GEMD)方法綜合了多種改進(jìn)EMD方法，通過(guò)定義不同的均值曲線對(duì)信號(hào)進(jìn)行逐階篩分，從得到的每一階分量中選取最優(yōu)作為最終的廣義內(nèi)稟模態(tài)函數(shù)(Generalized Intrinsic Mode Function，GIMF)，由于每一階的GIMF分量都是最優(yōu)的，因此相較于EMD等單一均值曲線篩分方法，GEMD分解結(jié)果也是最優(yōu)的。論文對(duì)GIMF分量準(zhǔn)則進(jìn)行了改進(jìn)以及對(duì)GEMD性能進(jìn)行了分析，并將GEMD應(yīng)用于仿真和實(shí)測(cè)信號(hào)的分析，結(jié)論表明GEMD分解是完備的和正交的，有比EMD更強(qiáng)的分解能力，而且適合機(jī)械振動(dòng)信號(hào)的處理和故障診斷。

關(guān)鍵詞：經(jīng)驗(yàn)?zāi)B(tài)分解；廣義經(jīng)驗(yàn)?zāi)B(tài)分解；局部特征尺度分解；分解能力；故障診斷

收稿日期：2013-09-12修改稿收到日期：2014-01-09

中圖分類號(hào):TH165+.3; TN911文獻(xiàn)標(biāo)志碼：A

基金項(xiàng)目：973項(xiàng)目

基金項(xiàng)目：北京理工大學(xué)基礎(chǔ)研究基金(20120142009)

Performance analysis and application of generalized empirical mode decomposition

ZHENGJin-de,CHENGJun-sheng,ZENGMing,LUOSong-rong(State key Laboratory of Advanced Design and Manufacture for Vehicle Body, Hunan University, Changsha 410082, China)

Abstract:Aiming at that the mean curve defined in empirical mode decomposition (EMD) fitted with cubic spline may cause envelope overshoot and undershoot, many improved EMDs for improving the mean curve are proposed and some good effects are achieved. Generalized empirical mode decomposition (GEMD) integrates several improved EMD methods and selects the best component from components obtained by sifting different mean curves in each rank as the final generalized IMF (GIMF). Since the GIMF is the best in each rank, the corresponding results of GEMD are also the best. Here, GEMD was introduced firstly and then an improved criterion of GIMF was developed. Furthermore, GEMD was employed to analyze simulated and actual mechanical vibration signals. The results showed that GEMD is complete and orthogonal, it has a better capacity of decomposition than that of EMD, it is suitable for mechanical fault diagnosis as well.

Key words:EMD; GEMD; local characteristic-scale decomposition; capacity of decomposition; fault diagnosis

經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition，EMD)是由Huang等[1-3]提出的一種自適應(yīng)的數(shù)據(jù)處理方法。EMD方法通過(guò)定義瞬時(shí)頻率具有物理意義的內(nèi)稟模態(tài)函數(shù)(Intrinsic Mode Function，IMF)，同時(shí)定義一種基于包絡(luò)的均值曲線，通過(guò)從原始信號(hào)中不斷分離出包絡(luò)均值曲線，即執(zhí)行篩分過(guò)程，將原始信號(hào)分解為若干個(gè)IMF分量和一個(gè)剩余項(xiàng)之和。

由于原EMD方法中均值曲線定義為上、下包絡(luò)線的均值，其中上、下包絡(luò)線分別由三次樣條擬合極大值和極小值而成，因此，不可避免地會(huì)產(chǎn)生欠包絡(luò)和過(guò)包絡(luò)以及擬合誤差等問(wèn)題。Chen等[4]提出了用B樣條代替三次樣條，以及Peel等[5-6]等提出用有理樣條取代三次樣條，這兩種改進(jìn)三次樣條插值的方法雖然克服了三次樣條帶來(lái)的問(wèn)題，但同時(shí)又引入了其它參數(shù)，且由于極值點(diǎn)是固定的，對(duì)大部分信號(hào)上述兩種分解方法的效果不是很明顯。蓋強(qiáng)和鄒巖崑等[7-8]提出的局域波分解方法提出了兩種新的均值曲線的定義方法，新的均值曲線避開(kāi)了上下包絡(luò)線的定義，直接由數(shù)據(jù)的極值點(diǎn)定義均值曲線，且只采用一次三次樣條擬合，因此節(jié)約了計(jì)算時(shí)間和提高了計(jì)算速度。FREI等提出了一種新的時(shí)頻能分析方法——本征時(shí)間尺度分解(Intrinsic Time-scale Decomposition， ITD)，ITD基于極值的特征尺度和數(shù)據(jù)本身線性變換，提高了計(jì)算速度，避免了擬合誤差以及端點(diǎn)效應(yīng)向數(shù)據(jù)內(nèi)部的擴(kuò)散影響，然而ITD分解由于采用線性變換，會(huì)導(dǎo)致得到的分量出現(xiàn)波形失真。最近，筆者[9]提出了另一種時(shí)頻分析方法——局部特征尺度分解方法(Local Characteristic-scale Decomposition， LCD)，LCD克服了ITD的不足，給出了一種新的均值曲線定義，然而LCD方法本質(zhì)上也是EMD和ITD的變形方法。

上述改進(jìn)的EMD方法本質(zhì)上是對(duì)EMD均值曲線的改進(jìn)，可能對(duì)某些信號(hào)而言，上述改進(jìn)的方法在某一方面有一定的提高，而對(duì)另外一些信號(hào)原EMD方法仍具有其獨(dú)特的優(yōu)勢(shì)；即使對(duì)同一個(gè)信號(hào)不同頻段的分量，可能改進(jìn)的方法提取出的第一個(gè)分量較接近實(shí)際值，而對(duì)于第二個(gè)分量則可能是EMD提取的分量更接近實(shí)際值。當(dāng)用上述方法分析實(shí)際信號(hào)時(shí)，無(wú)法判定哪種信號(hào)處理方法是最合適的和最優(yōu)的。

基于上述考慮，綜合各種改進(jìn)EMD方法的特點(diǎn)，筆者[10]提出了一種新的時(shí)頻分析方法——廣義經(jīng)驗(yàn)?zāi)B(tài)分解(Generalized Empirical Mode Decomposition, GEMD)。GEMD基于不同的均值曲線而定義，從高頻到低頻進(jìn)行篩分，從每階得到的分量中選擇最優(yōu)分量作為最終的廣義IMF(Generated IMF，GIMF)。再原始信號(hào)中減去該GIMF分量，并將剩余信號(hào)視為原始信號(hào)，不斷重復(fù)上述篩分過(guò)程，直到剩余信號(hào)是一單調(diào)或常函數(shù)，或者能量與原始信號(hào)相比可以忽略的函數(shù)。

本文詳細(xì)研究GEMD方法中GIMF分量的判據(jù)準(zhǔn)則，在原帶寬準(zhǔn)則的基礎(chǔ)上提出了一種瞬時(shí)頻率波動(dòng)能量差準(zhǔn)則；同時(shí)對(duì)GEMD分解的正交性和完備性進(jìn)行了驗(yàn)證，結(jié)果表明GEMD分解是正交的和完備的。最后，以兩個(gè)正弦信號(hào)疊加的混合信號(hào)模型為例，對(duì)GEMD方法的頻率分辨能力進(jìn)行了研究，結(jié)果表明GEMD可分解信號(hào)范圍要大于EMD，能分解出EMD不能分解的兩頻率較為接近的信號(hào)的疊加信號(hào)，是一種具有良好性能的信號(hào)分解方法。

1廣義經(jīng)驗(yàn)?zāi)B(tài)分解

1.1均值曲線的定義

GEMD方法的關(guān)鍵是定義信號(hào)合理的均值曲線。然而難點(diǎn)在于，到目前為止，大部分關(guān)于均值曲線的定義都是經(jīng)驗(yàn)性的。綜合多種參考文獻(xiàn)，論文采用以下六種具有代表性的均值曲線的定義方法：

(1)包絡(luò)均值曲線(Envelop Mean，EM)，定義詳見(jiàn)EMD方法[1]。

(2)固有時(shí)間尺度均值(Intrinsic Time-scale Mean，ITM)定義參見(jiàn)文獻(xiàn)[11]。

(3)局部特征尺度均值(Local Characteristic-scale Mean，LCM)，定義參見(jiàn)文獻(xiàn)[9]。

(4)極值點(diǎn)中值均值(Extremum- midpoint mean，EMM)。

假設(shè)待分析信號(hào)為Yt(t≥0)，所有極值點(diǎn)為(τk,Yk)(k=1,2,…,M)，相鄰兩個(gè)極大值(或極小值)的中點(diǎn)可表示為

(1)

定義

(2)

再用三次樣條函數(shù)擬合和插值所有的(tk+1,Bk+1)即得到信號(hào)Yt的均值曲線。

(5)極值域均值(Extremum Field Mean，EFM)，定義參見(jiàn)文獻(xiàn)[7]。

(6)局部積分均值(Local Integral Mean，LIM)，定義參見(jiàn)文獻(xiàn)[12]。

(3)

其中:w(t) 是權(quán)重因子。 特別地，二者可定義為：

(4)

上述六種均值曲線的定義方法從不同的角度定義信號(hào)的均值曲線，大部分是經(jīng)驗(yàn)性的，它們的物理模型還有待于進(jìn)一步驗(yàn)證。

1.2GEMD

設(shè)待分析信號(hào)為x(t)，廣義經(jīng)驗(yàn)?zāi)B(tài)分解(GEMD)的主要步驟如下[10]。

(1)u0(t)=x(t)，j=1；

(2)分別用上述六種方法提取信號(hào)uj-1(t)均值曲線：m1(t),m2(t), …,m6(t)；

(3)將上述均值曲線mi(t)從uj-1(t)中分離出，得到剩余信號(hào)ri(t)，即

ri(t)=uj-1(t)-mi(t)(i=1, 2, … 6)

(5)

uj(t)=uj-1(t)-Ij(t)

(6)

(6)j=j+1，重復(fù)上述步驟(2)～(5)，直到剩余分量un(t)為一常函數(shù)，或單調(diào)函數(shù)或能量與u0(t)相比可以忽略的函數(shù)，停止上述過(guò)程。

原始信號(hào)x(t)可重寫(xiě)為

(7)

GEMD方法綜合了EMD，ITD和LCD等這類分解方法的優(yōu)勢(shì)，繼承了EMD的篩分過(guò)程，通過(guò)逐層選優(yōu)，選出GIMF分量，再對(duì)剩余信號(hào)重復(fù)篩分。與EMD等單一的篩分過(guò)程相比，GEMD得到的分量在每一層的分量中是最優(yōu)的，因此，理論上GEMD分解的結(jié)果也是最優(yōu)的。

1.3GIMF判據(jù)

GEMD的原文獻(xiàn)采用頻率帶寬最小來(lái)選擇GIMF分量。頻率帶寬反映瞬時(shí)頻率的調(diào)制和波動(dòng)，頻率帶寬越小，表明瞬時(shí)頻率受調(diào)制越小，波動(dòng)越小。

當(dāng)瞬時(shí)頻率為常函數(shù)時(shí)，帶寬準(zhǔn)則能夠很好反映瞬時(shí)頻率的聚集性。若φ(t)為非常函數(shù)時(shí)，帶寬準(zhǔn)則失去了預(yù)期的結(jié)果。此時(shí)，采取如下方案來(lái)選擇GIMF：

(4)第i個(gè)pre-GIMF作為最終的GIMF。

原GIMF判據(jù)中通過(guò)帶寬準(zhǔn)則選擇GIMF，而改進(jìn)后的判據(jù)彌補(bǔ)了帶寬準(zhǔn)則的不足，從不同的IMF分量中選擇瞬時(shí)頻率更精確，帶寬更小的分量作為最優(yōu)，相鄰的GIMF分量帶寬無(wú)重疊或者重疊更小，因此分解正交性更好。

此外，由于GEMD是從每層不同的IMF分量中選擇最優(yōu)，如果EMD能將待分解信號(hào)分解出，則GEMD也一定能分解出；如果EMD不能分解出，而其它方法定義的均值曲線，如LCM或LIM可以分解出，則說(shuō)明GEMD也能分解出；只有所有的均值定義方法不能分解出時(shí)，GEMD才不能分解出，因此，GEMD的分解能力要比EMD強(qiáng)。

2仿真分析

考慮仿真信號(hào)

si(t)=x0(t)+xi(t), i=1,2。 x0(t)=

[1+0.2cos(2π2t)]cos(2πf0t+2πf0/10t2),

x1(t)=cos(2πa1f0t),

x2(t)=cos(2πa2f0t)，t∈[0, 1]

(8)

si(t)由調(diào)幅調(diào)頻信號(hào)x0(t)和余弦信號(hào)疊加而成，其中:f0=30,a1=0.5,a2=0.7。為比較，分別采用EMD和GEMD對(duì)s1(t)和s2(t)進(jìn)行分解，分解結(jié)果分別如圖1～2和圖4～5所示。

C i表示EMD分解得到的第i個(gè)IMF分量， res表示剩余項(xiàng)，下同。 圖1　混合信號(hào)s 1(t)的EMD分解結(jié)果 Fig.1 The EMD results of mixed signal s 1(t)

G i表示GEMD分解得到的 第i個(gè)GIMF分量，res表示剩余項(xiàng)，下同。 圖2　混合信號(hào)s 1(t)的GEMD分解結(jié)果 Fig.2 The GEMD results of mixed signal s 1(t)

從上到下依次為：第一個(gè)和第二個(gè)分量與真實(shí)分量的絕對(duì)誤差(實(shí)線表示EMD分解結(jié)果的絕對(duì)誤差，虛線表示GEMD分解結(jié)果的絕對(duì)誤差)、EMD分解的重構(gòu)誤差、GEMD分解的重構(gòu)誤差。 圖3　混合信號(hào)s 1(t)的EMD和 GEMD分解分量的絕對(duì)誤差和重構(gòu)誤差 Fig.3 The absolute error and reconstruction error of mixed signal s 1(t) generated by EMD and GEMD

圖4　混合信號(hào)s 2(t)的EMD分解結(jié)果 Fig.4 The EMD results of mixed signal s 2(t)

圖5　混合信號(hào)s 2(t)的GEMD分解結(jié)果 Fig.5. The GEMD results of mixed signal s 2(t)

從上到下依次為：第一個(gè)和第二個(gè)分量與對(duì)應(yīng)真實(shí)分量的絕對(duì)誤差，EMD分解的重構(gòu)誤差，GEMD分解的重構(gòu)誤差。 圖6　混合信號(hào)s 2(t)的EMD 和GEMD分解分量的絕對(duì)誤差和重構(gòu)誤差 Fig.6 The absolute error and reconstruction error of mixed signal s 2(t) generated by EMD and GEMD

分析圖1至圖6可得到以下結(jié)論。首先，由圖1和圖2看出，由于兩個(gè)分量成分的瞬時(shí)頻率(30+6tHz，15 Hz)相差較大，EMD和GEMD都能將二者完全分開(kāi)，且與真實(shí)分量絕對(duì)誤差較小(圖3第1, 2行)，但GEMD分解效果更好，GIMF分量與真實(shí)分量的絕對(duì)誤差幅值在大部分時(shí)間上要小于IMF分量與真實(shí)分量的絕對(duì)誤差幅值。其次，由圖4和圖5可以看出，由于兩個(gè)分量成分的瞬時(shí)頻率(30+6tHz，21 Hz)較為接近，EMD已不能將二者完全分開(kāi)，分解結(jié)果與真實(shí)分量相差較大，而GEMD則能夠?qū)⒍咄耆珠_(kāi)，與真實(shí)分量的絕對(duì)誤差較小(圖6第1, 2行)。再次，由圖3和圖6也可以看出，對(duì)于兩個(gè)混合信號(hào)，EMD和GEMD分解都是完備的，絕對(duì)重構(gòu)誤差數(shù)量級(jí)在10-15。

3GEMD分解能力分析

文獻(xiàn)[17]以簡(jiǎn)化的兩個(gè)正弦信號(hào)的疊加模型，對(duì)EMD方法的分解能力進(jìn)行了詳細(xì)的研究。為了對(duì)比分析，仍考慮兩個(gè)正弦信號(hào)疊加模型：

x(t)=cos2πt+acos(2πft+φ)

(9)

其中:φ表示兩信號(hào)的初相位差，a表示兩信號(hào)幅值比，f(0

定義評(píng)價(jià)指標(biāo)

δ(a,f)=

(10)

圖7　兩個(gè)諧波疊加模型的EMD的分解能力圖 Fig.7 The capacity of EMD decomposition of two harmonic superposition model

EMD和GEMD的分解能力分別如圖7和圖8所示。需要說(shuō)明的是，圖7中EMD可分解的分界頻率(設(shè)為fc)與文獻(xiàn)[17]有所差異，這主要是因?yàn)椴蓸宇l率fs以及fs與最大頻率fmax的關(guān)系對(duì)分解效果的影響[19-20]。由圖7和圖8可以發(fā)現(xiàn)，在相同的采樣頻率和最大頻率，相同的最大迭代次數(shù)(8次)等條件下，GEMD的分解能力范圍要明顯大于EMD的分解能力范圍。在幅值比a≥1時(shí)，曲線af2=1可視為不可分解部分的下邊界，af=1可視為可分解部分的上邊界。但在二者曲線之間的部分，GEMD的可分解部分仍大于EMD的可分解部分。此外，在幅值比0.1a1時(shí)，GEMD能分別出頻率更接近的信號(hào)。更詳細(xì)地，①當(dāng)0.1≤a<0.2時(shí)，GEMD可分解邊界頻率fc≈0.6；EMD可分解邊界頻率fc≈0.4；②當(dāng)0.2≤a<0.6時(shí)，GEMD可分解邊界頻率fc≈0.7；而EMD可分解邊界頻率fc≈0.5；③當(dāng)0.6≤a≤1時(shí)，GEMD的可分解邊界頻率fc≈0.77. 而EMD可分解邊界頻率fc≈0.55。

圖8　兩個(gè)諧波疊加模型的GEMD的分解能力圖 Fig.8 The capacity of GEMD decomposition of two harmonic superposition model

4應(yīng)用分析

為了說(shuō)明GEMD的有效性和相對(duì)于EMD的優(yōu)越性，論文將其應(yīng)用于具有故障的機(jī)械振動(dòng)信號(hào)分析。實(shí)驗(yàn)數(shù)據(jù)采用美國(guó)凱撒西儲(chǔ)大學(xué)公開(kāi)的軸承數(shù)據(jù)[21]，采用電火花技術(shù)對(duì)6205-2RS JEM SKF型深溝球軸承設(shè)置大小為0.017 78 cm，深度0.027 94 cm的單點(diǎn)故障；在轉(zhuǎn)速為1 750 r/min,負(fù)載為0 HP，采樣頻率為12 kHz的條件下，采集振動(dòng)加速度信號(hào)。不失一般性，以外圈故障軸承的振動(dòng)加速度信號(hào)為例，經(jīng)計(jì)算，轉(zhuǎn)頻：fr≈29.2 Hz，外圈故障特征頻率：fo≈104.6 Hz。振動(dòng)信號(hào)時(shí)域波形如圖9所示。

圖9　具有外圈故障的軸承振動(dòng)加速度信號(hào) Fig.9 Vibrationacceleration signal of rolling bearing with out race fault

采用GEMD對(duì)其進(jìn)行分解，結(jié)果如圖10所示。為了比較和識(shí)別故障頻率，圖11畫(huà)出了前五個(gè)GIMF的包絡(luò)譜。由圖11可以看出，前四個(gè)高頻分量被故障特征頻率調(diào)制，第五個(gè)分量被故障特征頻率和轉(zhuǎn)頻調(diào)制，解調(diào)后的包絡(luò)譜很明顯地反映了外圈故障特征，與實(shí)際相符。

其中:G i表示第i個(gè)GIMF，圖中只畫(huà)出了前五個(gè)GIMF及剩余項(xiàng) 圖10　圖9所示信號(hào)的GEMD分解結(jié)果 Fig.10 The results of signal shown in fig.9 decomposed by GEMD; G istands for the ith GIMF and only the first five GIMFs and its residue are shown.

圖11　從上到下依次為前五個(gè)GIMF的包絡(luò)譜 Fig.11 Envelope spectrum of the first five GIMFs

其中:c i表示第i個(gè)IMF，圖中只畫(huà)出了前五個(gè)IMF及剩余項(xiàng) 圖12　圖9所示信號(hào)的EMD分解結(jié)果 Fig.12 The results of signal shown in fig.9 decomposed by EMD; c i stands for the ith IMF and only the first five IMFs and its residue are shown

圖13　從上到下依次為前五個(gè)IMF的包絡(luò)譜 Fig.13 Envelope spectrum of the first five IMFs

采用EMD對(duì)圖9所示的信號(hào)進(jìn)行分解，結(jié)果如圖12所示，為了與GEMD結(jié)果比較，其前五個(gè)分量的包絡(luò)譜如圖13所示。由圖中可以看出，前三個(gè)IMF的包絡(luò)譜與前三個(gè)GIMF的包絡(luò)譜非常接近，從中也可以識(shí)別故障，而第四和第五個(gè)IMF的包絡(luò)譜包含的信息則不明顯，故障特征頻率和轉(zhuǎn)頻被低頻干擾較嚴(yán)重。因此，比較而言，GEMD分解有一定的優(yōu)越性。

5結(jié)論

對(duì)最新提出的廣義經(jīng)驗(yàn)?zāi)B(tài)分解(GEMD)方法的完備性和正交性進(jìn)行了數(shù)值驗(yàn)證，依據(jù)GEMD按照頻率從高頻到低頻的分解特征，提出了一種基于瞬時(shí)頻率能量波動(dòng)差的GIMF分量最優(yōu)評(píng)價(jià)準(zhǔn)則，并對(duì)其分解能力進(jìn)行了詳細(xì)研究。得出如下結(jié)論：

(1)GEMD與EMD方法篩分流程相似，分解是完備的和正交的，但GEMD有比EMD更好的正交性。

(2)與EMD相比，GEMD分解能力更強(qiáng)，可分解信號(hào)的范圍更廣。

(3)將GEMD應(yīng)用于機(jī)械故障振動(dòng)信號(hào)的分析，結(jié)果表明GEMD不僅適合于機(jī)械振動(dòng)信號(hào)的分析，而且有比EMD更好的故障診斷效果。

綜上，GEMD是一種有效的信號(hào)分解方法。當(dāng)然，GEMD也有其不足之處，如GIMF分量的評(píng)價(jià)準(zhǔn)則直接決定了分解的效果，而最優(yōu)GIMF選取準(zhǔn)則還有待于進(jìn)一步研究。相信隨著研究的深入和理論的完善，GEMD將會(huì)是一種優(yōu)秀的非線性和非平穩(wěn)信號(hào)處理方法。

參考文獻(xiàn)

[1]Huang N E, Zheng S, Long S R,et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proc. Roy. Soc. London, 1998, 454: 903-995.

[2]Wu Z, Huang N E. A study of the characteristics of white noise using the empirical mode decomposition method[J].Proc. R. Soc. London A, 2004, 460: 1597-1611.

[3]Hou T Y, Yan M P, Wu Z. A variant of the EMD method for multi-scale data, Adv. Adapt. Data Anal.1,2009:483-516.

(下轉(zhuǎn)第155頁(yè)) N E, Zheng S, Long S R,et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proc. Roy. Soc. London, 1998, 454: 903-995.

[2]Wu Z, Huang N E. A study of the characteristics of white noise using the empirical mode decomposition method[J].Proc. R. Soc. London A, 2004, 460: 1597-1611.

[3]Hou T Y, Yan M P, Wu Z. A variant of the EMD method for multi-scale data, Adv. Adapt. Data Anal.1,2009:483-516.

[4]Chen Qiu-hui, Huang Norden, Riemenschneider Sherman, et al. A B-spline approach for empirical mode decompositions[J].Advances in Computational Mathematics, 2006, 24: 171-195.

[5]Peel M C, McMahon T A, Pegram G G S. Assessing the performance of rational spline-based empirical mode decomposition using a global annual precipitation dataset[J]. Proc. R. Soc. A, 2009, 465: 1919-1937.

[6]Pegram G G S, Peel M C,McMahon T A. Empirical mode decomposition using rational splines: an application to rainfall time series[J]. Proc. R. Soc. A., 2008, 464:1483-1501.

[7]蓋強(qiáng). 局域波時(shí)頻分析方法的理論研究與應(yīng)用[D]. 大連:大連理工大學(xué), 2001.

[8]鄒巖崑. 局域波分析的理論方法研究及應(yīng)用[D]. 大連:大連理工大學(xué), 2004.

[9]程軍圣，鄭近德，楊宇. 一種新的非平穩(wěn)信號(hào)分析方法—局部特征尺度分解法[J].振動(dòng)工程學(xué)報(bào)，2012, 25(2): 215-220.

CHENG Jun-sheng, ZHENG Jin-de, YANG Yu.A new nonstationary signal analysis approach-the local characteristic-scale decomposition method[J]. Journal of Vibration Engineering, 2012, 25(2):215-220.

[10]Zheng Jin-de, Cheng Jun-sheng, Yang Yu. Generalized empirical mode decomposition and its applications to rolling element bearing fault diagnosis[J]. Mechanical Systems and Signal Processing, 2013, 40(1):136-153.

[11]Frei M G,Osorio I.Intrinsic time-scale decomposition:Time-frequency-energy analysis and real-time filtering of non-stationary signals[J].Proc. Royal Soc. A, 2007, 463:321-342.

[12]Hong Hong, Wang Xin-long, Tao Zhi-yong. Local integral mean-based sifting for empirical mode decomposition[J]. IEEE signal processing letters, 2009, 16(10): 841-844.

[13]Cohen L. Time frequency analysis: theory and applications [M]. Prentice Hall Press, 1995.

[14]Xuan B, Xie Q, Peng S.EMD sifting basedon band width [J]. IEEE Signal Process. Lett, 2007,14(8):537-540.

[15]鄭近德，程軍圣，楊宇. 一種新的估計(jì)瞬時(shí)頻率的方法—經(jīng)驗(yàn)包絡(luò)法[J].振動(dòng)與沖擊，2012, 31 (17): 86-90.

ZHENG Jin-de, YANG Yu, CHENG Jun-sheng,A new instantaneous frequency estimation approach-empirical envelope method[J].Journal of Vibration and Shock, 2012, 31 (17): 86-90.

[16]Huang N E, Wu Z H, Long S R, et al. On the frequency [J]. Advances in Adaptive Data Analysis, 2009, 1(2): 177-229.

[17]Rilling G,Flandrin P. One or two frequencies the empirical mode decomposition answers[J]. IEEE Trans. Signal Process., 2008, 56(1): 85-95.

[18]楊宇，曾鳴，程軍圣.局部特征尺度分解方法及其分解能力研究[J].振動(dòng)工程學(xué)報(bào)，2012, 25(5)：602-609.

YANG Yu, ZENG Ming, CHENG Jun-sheng.Research on local characteristic-scale decomposition and its capacities[J].Journal of Vibration Engineering, 2012, 25(5)：602-609.

[19]Rilling G, Flandrin P. On the influence of sampling on the empirical mode decomposition[J]. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), 2006:444-447.

[20]Rilling G, Flandrin P. Sampling effects on the empirical mode decomposition [J]. Advances in Adaptive Data Analysis,2009, 1(1): 43-59.

[21]Xu H B, Chen G H, An intelligent fault identification method of rolling bearings based on LSSVM optimized by improved PSO[J]. Mechanical Systems and Signal Processing,2013,35:167-175. Qiu-hui, Huang Norden, Riemenschneider Sherman, et al. A B-spline approach for empirical mode decompositions[J].Advances in Computational Mathematics, 2006, 24: 171-195.

[5]Peel M C, McMahon T A, Pegram G G S. Assessing the performance of rational spline-based empirical mode decomposition using a global annual precipitation dataset[J]. Proc. R. Soc. A, 2009, 465: 1919-1937.

[6]Pegram G G S, Peel M C,McMahon T A. Empirical mode decomposition using rational splines: an application to rainfall time series[J]. Proc. R. Soc. A., 2008, 464:1483-1501.

[7]蓋強(qiáng). 局域波時(shí)頻分析方法的理論研究與應(yīng)用[D]. 大連:大連理工大學(xué), 2001.

[8]鄒巖崑. 局域波分析的理論方法研究及應(yīng)用[D]. 大連:大連理工大學(xué), 2004.

[9]程軍圣，鄭近德，楊宇. 一種新的非平穩(wěn)信號(hào)分析方法—局部特征尺度分解法[J].振動(dòng)工程學(xué)報(bào)，2012, 25(2): 215-220.

CHENG Jun-sheng, ZHENG Jin-de, YANG Yu.A new nonstationary signal analysis approach-the local characteristic-scale decomposition method[J]. Journal of Vibration Engineering, 2012, 25(2):215-220.

[10]Zheng Jin-de, Cheng Jun-sheng, Yang Yu. Generalized empirical mode decomposition and its applications to rolling element bearing fault diagnosis[J]. Mechanical Systems and Signal Processing, 2013, 40(1):136-153.

[11]Frei M G,Osorio I.Intrinsic time-scale decomposition:Time-frequency-energy analysis and real-time filtering of non-stationary signals[J].Proc. Royal Soc. A, 2007, 463:321-342.

[12]Hong Hong, Wang Xin-long, Tao Zhi-yong. Local integral mean-based sifting for empirical mode decomposition[J]. IEEE signal processing letters, 2009, 16(10): 841-844.

[13]Cohen L. Time frequency analysis: theory and applications [M]. Prentice Hall Press, 1995.

[14]Xuan B, Xie Q, Peng S.EMD sifting basedon band width [J]. IEEE Signal Process. Lett, 2007,14(8):537-540.

[15]鄭近德，程軍圣，楊宇. 一種新的估計(jì)瞬時(shí)頻率的方法—經(jīng)驗(yàn)包絡(luò)法[J].振動(dòng)與沖擊，2012, 31 (17): 86-90.

ZHENG Jin-de, YANG Yu, CHENG Jun-sheng,A new instantaneous frequency estimation approach-empirical envelope method[J].Journal of Vibration and Shock, 2012, 31 (17): 86-90.

[16]Huang N E, Wu Z H, Long S R, et al. On the frequency [J]. Advances in Adaptive Data Analysis, 2009, 1(2): 177-229.

[17]Rilling G,Flandrin P. One or two frequencies the empirical mode decomposition answers[J]. IEEE Trans. Signal Process., 2008, 56(1): 85-95.

[18]楊宇，曾鳴，程軍圣.局部特征尺度分解方法及其分解能力研究[J].振動(dòng)工程學(xué)報(bào)，2012, 25(5)：602-609.

YANG Yu, ZENG Ming, CHENG Jun-sheng.Research on local characteristic-scale decomposition and its capacities[J].Journal of Vibration Engineering, 2012, 25(5)：602-609.

[19]Rilling G, Flandrin P. On the influence of sampling on the empirical mode decomposition[J]. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), 2006:444-447.

[20]Rilling G, Flandrin P. Sampling effects on the empirical mode decomposition [J]. Advances in Adaptive Data Analysis,2009, 1(1): 43-59.

[21]Xu H B, Chen G H, An intelligent fault identification method of rolling bearings based on LSSVM optimized by improved PSO[J]. Mechanical Systems and Signal Processing,2013,35:167-175.