鄭近德 蘇繆涎 潘海洋

摘要： 自適應(yīng)噪聲輔助集成經(jīng)驗(yàn)?zāi)B(tài)分解（Complete Ensemble Empirical Mode Decomposition with Adaptive Noise， CEEMDAN）解決了集成經(jīng)驗(yàn)?zāi)B(tài)分解在集成平均過(guò)程中的分解不完備問(wèn)題，但噪聲殘留和虛假分量問(wèn)題仍然存在。針對(duì)CEEMDAN的不足，提出了自適應(yīng)噪聲加權(quán)優(yōu)選經(jīng)驗(yàn)?zāi)B(tài)分解（Weighted Mean?optimized Empirical Mode Decomposition with Adaptive Noise， WMEMDAN）。該方法用改進(jìn)的均值曲線構(gòu)造方式提取內(nèi)稟模態(tài)函數(shù)（IMF），以正交性最小為依據(jù)，從不同權(quán)重的迭代篩分結(jié)果中選取出最優(yōu)IMF，改善了CEEMDAN的分解能力，同時(shí)通過(guò)對(duì)不同權(quán)重下的分解結(jié)果進(jìn)行篩選，確保每一階的IMF分量都是整體最優(yōu)，減少虛假分量和殘留噪聲。仿真和實(shí)驗(yàn)信號(hào)分析結(jié)果表明，WMEMDAN在減少虛假分量和提高分解精度等方面具有優(yōu)勢(shì)。將所提方法應(yīng)用于滾動(dòng)軸承和齒輪的故障診斷，分析結(jié)果表明了方法的有效性和優(yōu)越性。

關(guān)鍵詞： 故障診斷; 軸承; 齒輪; CEEMDAN; 均值優(yōu)選

引 言

Huang等[1]提出的經(jīng)驗(yàn)?zāi)B(tài)分解（Empirical mode decomposition， EMD）是一種基于數(shù)據(jù)驅(qū)動(dòng)的自適應(yīng)分解方法，可以將信號(hào)從高頻到低頻自適應(yīng)地分解為一系列內(nèi)稟模態(tài)函數(shù)（Intrinsic mode function， IMF）和趨勢(shì)項(xiàng)之和[2?4]。EMD被廣泛應(yīng)用于諸多領(lǐng)域，但也存在著模態(tài)混疊[5?7]、端點(diǎn)效應(yīng)[8?9]、包絡(luò)過(guò)沖和不足[10?11]等缺陷。為了克服模態(tài)混疊，WU等[12]利用高斯白噪聲輔助分解，提出了集成經(jīng)驗(yàn)?zāi)B(tài)分解（Ensemble empirical mode decomposition， EEMD）。EEMD通過(guò)添加高斯白噪聲能夠抑制模態(tài)混疊，但其也帶來(lái)了噪聲殘留的問(wèn)題[13?15]。Yeh等[16]針對(duì)EEMD殘余噪聲問(wèn)題提出了補(bǔ)充總體平均經(jīng)驗(yàn)?zāi)B(tài)分解（Complementary Ensemble Empirical Mode Decomposition， CEEMD）。CEEMD通過(guò)添加符號(hào)為正負(fù)的高斯白噪聲后，對(duì)分解結(jié)果集成平均，可以達(dá)到減少噪聲殘余的目的。

然而，不論是EEMD還是CEEMD，在信號(hào)加入不同高斯白噪聲后，分解結(jié)果的數(shù)目都很難保持一致，致使在集成平均過(guò)程中不能保證分解的完備性。鑒于此，自適應(yīng)噪聲完整集成經(jīng)驗(yàn)?zāi)B(tài)分解[17?21]被提出，其解決了分解結(jié)果數(shù)量不一致的問(wèn)題，并具有分解完備性。但是，CEEMDAN仍然存在一些問(wèn)題有待解決：（1）端點(diǎn)效應(yīng)和包絡(luò)過(guò)沖與不足問(wèn)題;（2）易得到虛假分量;（3）殘余噪聲影響分解結(jié)果的精度。

CEEMDAN對(duì)極值點(diǎn)進(jìn)行三次樣條曲線插值擬合來(lái)構(gòu)造均值曲線，然后將均值曲線從信號(hào)中迭代篩分，且不斷更新，而每次篩分采用的均值曲線都是基于上一次迭代的剩余信號(hào)[22]，從而自適應(yīng)地分解信號(hào)。鑒于均值曲線的構(gòu)造和迭代篩分在CEEMDAN算法中起著決定性作用，考慮對(duì)這兩個(gè)核心部分加以改進(jìn)，以解決CEEMDAN存在的問(wèn)題。針對(duì)插值擬合導(dǎo)致的端點(diǎn)效應(yīng)和包絡(luò)過(guò)沖與不足等問(wèn)題，考慮對(duì)信號(hào)進(jìn)行延拓，同時(shí)對(duì)相鄰三個(gè)極值點(diǎn)加權(quán)平均，均勻極值點(diǎn)的分布，來(lái)抑制包絡(luò)的過(guò)沖和不足問(wèn)題。在迭代篩分過(guò)程中，考慮在均值曲線還未迭代更新之前，將均值曲線從當(dāng)前剩余信號(hào)中盡可能完全篩分出[22]，在篩分過(guò)程中引入不同的權(quán)重，以正交性作為衡量篩分是否徹底的依據(jù)，從不同權(quán)重下的篩分結(jié)果中選取出最優(yōu)結(jié)果，最大程度地減少篩分結(jié)果中的噪聲殘留。對(duì)每次篩分過(guò)程都加以優(yōu)選，在經(jīng)過(guò)多次迭代更新后，可以確保IMF分量達(dá)到整體最優(yōu)，以減少虛假分量的存在，提高分解的精度。

1 自適應(yīng)噪聲加權(quán)優(yōu)選經(jīng)驗(yàn)?zāi)B(tài)分解

1.1 CEEMDAN方法

CEEMDAN在每一階段添加特定的白噪聲，對(duì)每一階段的分解結(jié)果進(jìn)行集成平均，解決了EEMD集成平均過(guò)程中分解結(jié)果數(shù)量不一致的問(wèn)題，在不增加計(jì)算量的前提下，降低了信號(hào)的重構(gòu)誤差，具有更好的分解完備性。在CEEMDAN方法中，設(shè)原始信號(hào)為x（t），ωm （t）為第m次向信號(hào)添加的白噪聲序列，定義算子Ep（·）為EMD的第p個(gè)IMF分量，具體方法描述如下：

1.2 WMEMDAN方法

WMEMDAN的關(guān)鍵在于均值曲線的構(gòu)造和優(yōu)選。在均值曲線構(gòu)造方面，不同于CEEMDAN直接對(duì)極值點(diǎn)進(jìn)行包絡(luò)擬合，WMEMDAN先對(duì)相鄰極值點(diǎn)加權(quán)平均，再對(duì)加權(quán)平均后的極值點(diǎn)進(jìn)行包絡(luò)擬合。在均值曲線優(yōu)選方面，為了盡可能從信號(hào)中分離出均值曲線，在不同權(quán)重下，以最小正交性為依據(jù)選取整體最優(yōu)的分解結(jié)果。正交性越小意味著分解結(jié)果的分解正交性越好，每一階分量之間的混淆程度越低，越能準(zhǔn)確地反映真實(shí)的物理過(guò)程。此外，WMEMDAN還對(duì)原始信號(hào)進(jìn)行端點(diǎn)延拓，在分解出結(jié)果后去除延拓部分，避免端點(diǎn)效應(yīng)的影響。該方法具體步驟如下：

WMEMDAN引入了不同權(quán)重β，將不同權(quán)重的均值曲線從信號(hào)中分離，并計(jì)算分解結(jié)果的正交性指標(biāo)。正交性指標(biāo)越小，說(shuō)明分解結(jié)果中每一階IMF分量彼此之間的混淆程度越低，越能反映信號(hào)中的故障沖擊特征，有利于提取故障特征，因此從不同的權(quán)重分解結(jié)果中優(yōu)選出正交性最小的分解結(jié)果作為WMEMDAN方法的最終IMF分量。引入權(quán)重的目的是為了從不同權(quán)重下獲得整體最優(yōu)的IMF分量，如果權(quán)重β過(guò)大，均值曲線被過(guò)度分離，導(dǎo)致信號(hào)X（t）中包含的有用信息被剔除，IMF分量呈現(xiàn)失真情況。如果權(quán)重β個(gè)數(shù)太少的話，難以確保所得IMF分量的整體最優(yōu)性。綜合考慮之下，令權(quán)重β在0到2之間以0.1為步長(zhǎng)取值。在此范圍內(nèi)，權(quán)重β最大值為2，不容易出現(xiàn)過(guò)度分離的情況。以一定步長(zhǎng)取值，權(quán)重值連續(xù)變化，不會(huì)因?yàn)闄?quán)重值的缺少而導(dǎo)致分解結(jié)果的整體最優(yōu)性無(wú)法保證，當(dāng)然步長(zhǎng)越小，權(quán)重個(gè)數(shù)越多，理論上能夠獲得更好的整體最優(yōu)性，但增加的計(jì)算量和提高的整體最優(yōu)性的實(shí)際性價(jià)比不高，因此，論文以0.1為步長(zhǎng)進(jìn)行取值。

在WMEMDAN方法中，改進(jìn)的均值曲線構(gòu)造方式抑制了過(guò)沖和不足，利用端點(diǎn)延拓減少端點(diǎn)效應(yīng)的影響，優(yōu)選均值曲線能使結(jié)果達(dá)到整體最優(yōu)的狀態(tài)，因此，理論上能夠提高分解能力，解決CEEMDAN在噪聲殘留和虛假分量方面的缺陷。

1.3 加權(quán)因子α

式（5）中引入了加權(quán)因子α，本節(jié)將分析α對(duì)分解結(jié)果的影響。在0.5?4的范圍內(nèi)以0.05為步長(zhǎng)取值作為α的樣本，采用下式的疊加信號(hào)為例進(jìn)行分析。

信號(hào)x（t）由幅值為2的隨機(jī)噪聲x1（t）、調(diào)幅信號(hào)x2（t）以及余弦信號(hào)x3（t）疊加而成。計(jì)算不同α下分解結(jié)果的正交性指標(biāo)和前三個(gè)分量的誤差能量比，如圖1和2所示。

不同加權(quán)因子α下的正交性指標(biāo)大致呈現(xiàn)先下降再上升的趨勢(shì)，大約在2范圍內(nèi)，正交性指標(biāo)數(shù)值較小，分解正交性更好。對(duì)比不同分量的誤差能量比，當(dāng)加權(quán)因子α<2.5時(shí)，各個(gè)分量的誤差能量比基本趨于平穩(wěn)，波動(dòng)不大，而加權(quán)因子α>3時(shí)，各個(gè)分量的誤差能量比有了明顯地增長(zhǎng)趨勢(shì)，和真實(shí)分量的吻合程度開(kāi)始下降。綜合正交性指標(biāo)和誤差能量比，可以認(rèn)為加權(quán)因子α不適宜過(guò)小或者過(guò)大，在2附近是一個(gè)比較合適的選取范圍，既能夠?qū)崿F(xiàn)平緩極值點(diǎn)波動(dòng)，其正交性指標(biāo)和誤差能量比都能夠獲得較好的數(shù)值，確保分解結(jié)果的分解正交性和分解精度。

經(jīng)過(guò)多次試驗(yàn)分析都驗(yàn)證了該范圍是一個(gè)合適的選取區(qū)域，所以論文在2附近選取α，同時(shí)為了獲得最好的分解效果，從中以正交性最小依據(jù)進(jìn)行選取。

2 仿真分析

2.1 疊加仿真信號(hào)分析

為了驗(yàn)證WMEMDAN的優(yōu)越性，不失一般性，仍考慮式（10）所示的疊加信號(hào)，其時(shí)域波形圖如圖3所示。用WMEMDAN對(duì)疊加信號(hào)進(jìn)行分解，添加的噪聲幅值為0.25，集成次數(shù)為50，分解個(gè)數(shù)設(shè)置為5，加權(quán)因子取2，其結(jié)果如圖4所示。為了便于對(duì)比，圖5中給出相同參數(shù)條件下CEEMDAN的分解結(jié)果。從圖4和5可知，WMEMDAN和CEEMDAN都能夠分解疊加信號(hào)，WMEMDAN前三個(gè)IMF分量分別對(duì)應(yīng)真實(shí)分量x1（t），x2（t），x3（t），而CEEMDAN中對(duì)應(yīng)真實(shí)分量的是IMF1，IMF3和IMF5，中間存在著x1（t）和x2（t）的虛假分量IMF2與IMF4，兩者的幅值很大。對(duì)于疊加信號(hào)而言，相比于CEEMDAN，所提出的方法有效地減少了虛假分量。再考慮兩種方法的分解精度，表1中給出了IMF分量的正交性指標(biāo)（IO）和相關(guān)系數(shù)（Cor）。在相同參數(shù)條件下，WMEMDAN具有更小的IO，這意味著分解結(jié)果的模態(tài)混疊程度更低，分解效果更優(yōu)。WMEMDAN前三個(gè)分量的相關(guān)系數(shù)Cor分別為0.8576，0.9511和0.9835，都大于CEEMDAN的相關(guān)系數(shù)0.8101，0.8617和0.9833，這意味著WMEMDAN的分解結(jié)果和真實(shí)分量相關(guān)性大，更與真實(shí)分量吻合。兩種方法分解結(jié)果和分解精度的對(duì)比，驗(yàn)證了WMEMDAN在減少虛假分量的同時(shí)，能夠獲得更優(yōu)的分解精度。

2.2 滾動(dòng)軸承故障仿真信號(hào)

再考慮構(gòu)造滾動(dòng)軸承內(nèi)圈故障模型來(lái)驗(yàn)證WMEMDAN的優(yōu)勢(shì)。滾動(dòng)軸承內(nèi)圈故障信號(hào)S（t）可以看做是實(shí)際故障特征信號(hào)Sfault（t）和噪聲信號(hào)Snoise（t）的疊加，即

在沖擊信號(hào)Sfault（t）中，δ為阻尼系數(shù)，fn為固有頻率，ff為故障特征頻率，F(xiàn)s為采樣頻率，frf為轉(zhuǎn)頻，Q為載荷，G為滾珠質(zhì)量系數(shù)，N為采樣點(diǎn)數(shù)，u為采樣頻率Fs和故障特征頻率ff的整數(shù)比。

設(shè)置采樣頻率Fs=12 kHz，固有頻率fn=1500 Hz，采樣點(diǎn)數(shù)N=4096，載荷Q=5 kN，滾珠質(zhì)量系數(shù)G=2，轉(zhuǎn)頻frf和內(nèi)圈故障特征頻率ff分別為5 Hz和27 Hz，如圖6所示。分別采用CEEMDAN和WMEMDAN兩種方法分解仿真信號(hào)，添加的噪聲幅值為0.25 dB，集成次數(shù)為200次，設(shè)置分解結(jié)果為10，加權(quán)因子為2，結(jié)果如圖7所示。在前4個(gè)IMF分量中，WMEMDAN峭度最大的分量是IMF3，而CEEMDAN是IMF4。分別對(duì)兩個(gè)分量進(jìn)行幅值譜和功率譜分析，如圖8所示。CEEDMAN中幅值最大的譜線為202 Hz的干擾信號(hào)，而故障特征頻率淹沒(méi)在大量的干擾頻率中。WMEMDAN中故障特征頻率所在譜線最為明顯，干擾頻率被很好地抑制，故障特征頻率的識(shí)別效果更好。

通過(guò)上述分析，對(duì)于疊加信號(hào)和軸承故障仿真信號(hào)而言，WMEMDAN能夠有效地分解信號(hào)，同時(shí)減少虛假分量、提高分解精度，該方法的優(yōu)勢(shì)得到初步的驗(yàn)證。

3 實(shí)測(cè)信號(hào)分析

為了進(jìn)一步驗(yàn)證所提出方法的優(yōu)越性，將WMEMDAN應(yīng)用于軸承故障診斷和齒輪故障診斷中。所用軸承故障信號(hào)由圖9所示的ID?25/30型軸承全壽命試驗(yàn)臺(tái)采集，測(cè)試軸承為6206?2RSJ/C3深溝球軸承。齒輪故障診斷信號(hào)由圖10所示裝置采集，采用模數(shù)為2.5，齒數(shù)為37的標(biāo)準(zhǔn)齒輪作為測(cè)試齒輪。

3.1 軸承故障診斷

3.1.1 軸承內(nèi)圈故障

軸承故障直徑為0.4 mm，主軸轉(zhuǎn)速為1500 r/min，轉(zhuǎn)頻?ri=25 Hz，負(fù)載7 kN，采樣頻率為8192 Hz，內(nèi)圈故障特征頻率為?i=135 Hz。內(nèi)圈故障信號(hào)的時(shí)域圖如圖11所示。先采用WMEMDAN方法對(duì)信號(hào)進(jìn)行分解，添加的噪聲幅值為0.25，集成次數(shù)為50，分解個(gè)數(shù)設(shè)置為6，加權(quán)因子取1.98，其分解結(jié)果如圖12（a）所示，如無(wú)特殊說(shuō)明，分解結(jié)果只取前4個(gè)分量（下同），對(duì)分量做包絡(luò)譜分析，如圖12（b）所示。圖13中是相同參數(shù)條件下CEEMDAN的分解結(jié)果和包絡(luò)譜分析。

從圖12（b）和13（b）中可以看出，WMEMDAN和CEEMDAN前兩個(gè)分量中都存在故障特征頻率和轉(zhuǎn)頻，能夠有效地識(shí)別，但是IMF3和IMF4分量卻有效地說(shuō)明WMEMDAN方法的優(yōu)越性。WMEMDAN的IMF3中轉(zhuǎn)頻和故障特征頻率所在譜線非常明顯，干擾信號(hào)被很好地抑制，能夠準(zhǔn)確地識(shí)別，而CEEMDAN的IMF3中故障特征頻率幅值小，無(wú)法識(shí)別，且存在39.3，59，83.6 Hz等明顯的干擾頻率。WMEMDAN和CEEMDAN的IMF4分量中都存在轉(zhuǎn)頻，但是前者低頻段的無(wú)關(guān)頻率少，轉(zhuǎn)頻譜線明顯，后者出現(xiàn)了14.7和34.4 Hz的兩個(gè)干擾頻率。上述結(jié)果表明，WMEMDAN能夠有效地分解內(nèi)圈故障信號(hào)，同時(shí)抑制干擾信號(hào)的效果優(yōu)于CEEMDAN方法。為了量化說(shuō)明WMEMDAN方法在分解精度方面的優(yōu)越性，用下式所示的SNRd來(lái)衡量信噪比的大小。

IMF3和IMF4的SNRd如表2所示。不論是IMF3還是IMF4中，WMEMDAN的SNRd都小于CEEMDAN的SNRd，這意味著干擾信號(hào)所占的比重少，抑制效果更佳。

3.1.2 軸承外圈故障

低速軸承故障信號(hào)因其振動(dòng)頻率低，診斷故障的難度更大，以低速軸承外圈故障信號(hào)為例加以分析。軸承故障直徑為0.3 mm，主軸轉(zhuǎn)速為150 r/min，轉(zhuǎn)頻?ro=2.5 Hz，負(fù)載為5 kN，采樣頻率為8192 Hz，故障特征頻率為?o=8.92 Hz。低速外圈故障信號(hào)的時(shí)域圖如圖14所示。在與內(nèi)圈故障設(shè)置相同參數(shù)的條件下，仍采用WMEMDAN和CEEMDAN方法對(duì)低速故障信號(hào)加以分解，其結(jié)果和相應(yīng)的包絡(luò)譜分析如圖15和16所示。在圖16（b）中只存在兩倍轉(zhuǎn)頻，故障特征頻率無(wú)法辨識(shí)，這意味著CEEMDAN對(duì)于低速故障信號(hào)的分解能力不夠，而WMEMDAN的IMF4中出現(xiàn)了故障特征頻率和兩倍轉(zhuǎn)頻，和干擾頻率的區(qū)分度高，能夠準(zhǔn)確診斷出外圈故障。仍用式（15）的SNRd來(lái)量化地說(shuō)明WMEMDAN在低速外圈軸承故障診斷中的優(yōu)勢(shì)。CEEMDAN的包絡(luò)譜中并未出現(xiàn)明顯的故障特征頻率，只是出現(xiàn)了兩倍轉(zhuǎn)頻，WMEMDAN中在IMF4分量中出現(xiàn)了故障特征頻率和兩倍的轉(zhuǎn)頻，為此，分別計(jì)算IMF4中故障頻率和轉(zhuǎn)頻的SNRd，如表3所示。不論是故障特征頻率還是轉(zhuǎn)頻，WMEMDAN的SNRd都小于CEEMDAN，因此，所提出方法的分解效果更好。

3.2 齒輪故障分析

為了驗(yàn)證WMEMDAN方法的有效性，以齒輪斷齒故障為例進(jìn)行分析。齒輪故障分析中主軸轉(zhuǎn)速為720 r/min，采樣頻率為2048 Hz，齒輪故障特征頻率fg= 12 Hz。故障信號(hào)如圖17所示。添加的噪聲幅值為0.3，集成次數(shù)為100，分解個(gè)數(shù)為7，加權(quán)因子取2，兩種方法的分解結(jié)果和包絡(luò)譜分析如圖18?19所示。對(duì)比圖18（b）和19（b）可知，CEEMDAN的故障特征頻率隱藏于干擾頻率中，準(zhǔn)確地提取故障特征頻率比較困難。雖然IMF2中出現(xiàn)了2倍故障特征頻率，但是該倍頻和干擾頻率的區(qū)分度不高。在圖18（b）的IMF2中，出現(xiàn)了故障特征頻率，其譜線幅值明顯高于其余頻率，有利于判斷故障類型。如表4中所示，前三個(gè)IMF分量中，WMEMDAN的SNRd略低于CEEMDAN，尤其是IMF2分量中，兩者SNRd的差值最大，與包絡(luò)譜分析結(jié)果相吻合。雖然IMF4中，WMEMDAN的優(yōu)勢(shì)不夠明顯，但并不影響WMEMDAN方法識(shí)別齒輪故障特征頻率，因此，可以認(rèn)為所提出的方法在有效地識(shí)別齒輪故障的同時(shí)，抑制了干擾信號(hào)，能夠獲得更好的分解精度。

綜上所述，軸承內(nèi)圈故障、低速軸承外圈故障以及齒輪故障實(shí)測(cè)分析都驗(yàn)證了，相比于CEEMDAN，所提出的WMEMDAN方法在抑制虛假分量和干擾信號(hào)方面具有優(yōu)勢(shì)。

4 結(jié) ?論

（1） WMEMDAN在改善CEEMDAN分解能力的同時(shí)，減少分解過(guò)程中產(chǎn)生的虛假分量，抑制殘留噪聲對(duì)分解結(jié)果的影響，分解精度更高。

（2） 將所提出的方法應(yīng)用于軸承內(nèi)圈故障、低速軸承外圈故障和齒輪故障診斷中，相比CEEMDAN，具有更好的故障診斷能力。不論是方法理論分析還是實(shí)驗(yàn)數(shù)據(jù)分析，都驗(yàn)證了所提出方法的優(yōu)越性，有效地改進(jìn)了CEEMDAN的缺陷。

參考文獻(xiàn)：

[1] 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， A Mathematical Physical & Engineering Sciences，1998， 454（1971）： 903-995.

[2] Penga Z K， Tseb P W， Chua F L. An improved Hilbert-Huang transform and its application in vibration signal analysis[J]. Journal of Sound and Vibration， 2005， 286： 187?205.

[3] Yu Dejie， Cheng Jnusheng， Yang Yu. Application of EMD method and Hilbert spectrum to the fault diagnosis of roller bearings[J]. Mechanical Systems and Signal Processing， 2005， 19（2）：259-270.

[4] Yang Yu， Yu Dejie， Cheng Junsheng. A roller bearing fault diagnosis method based on EMD energy entropy and ANN[J]. Journal of Sound & Vibration， 2006， 294（1）：269-277.

[5] Deering R， Kaiser J F. ?The use of a masking signal to improve empirical mode decomposition[C]. Proceedings， IEEE Int. Conf. Acoust. Speech， Signal Process. （ICASSP305）， Philadelphia， USA， 2005， 4： iv/485-iv/488.

[6] 何 劉， 林建輝， 丁建明， 等. 調(diào)幅-調(diào)頻信號(hào)的經(jīng)驗(yàn)?zāi)B(tài)分解包絡(luò)技術(shù)和模態(tài)混疊[J]. 機(jī)械工程學(xué)報(bào)，2017，53（2）：1-10.

He Liu， Lin Jianhui， Ding Jianming， et al. Empirical mode decomposition envelope technique and mode mixing problem in amplitude modulation-frequency modulation signals[J]. Journal of Mechanical Engineering， 2017，53（2）：1-10.

[7] 曹 瑩， 段玉波， 劉繼承. Hilbert-Huang變換中的模態(tài)混疊問(wèn)題[J]. 振動(dòng).測(cè)試與診斷，2016，36（3）：518-523+605-606.

Cao Ying， Duan Yubo， Liu Jicheng. Research and application of mode-mixing in Hilbert-Huang transform[J]. Journal of Vibration， Measurement and Diagnosis， 2016，36（3）： 518-523+605-606.

[8] 王學(xué)敏， 黃方林. EMD端點(diǎn)效應(yīng)抑制的一種實(shí)用方法[J]. 振動(dòng).測(cè)試與診斷，2012，32（3）：493-497.

Wang Xuemin， Huang Fanglin. Practical method to restrain the end effect of EMD[J]. Journal of Vibration， Measurement and Diagnosis， 2012，32（3）：493-497.

[9] 徐力彬， 宋余慶， 劉 毅. 基于鏡像延拓和窗函數(shù)的端點(diǎn)效應(yīng)抑制方法[J]. 計(jì)算機(jī)工程，2015，41（4）：112-116.

Xu Libin， Song Qingyu， Liu Yi. Restraining method for end effect based on mirror extension and window function[J]. Computer Engineering， 2015，41（4）：112-116.

[10] Zheng Jinde， Pan Haiyang， Liu Tao， et al. Extreme-point weighted mode decomposition[J]. Signal Processing， 2018， 142： 366-374.

[11] Jia Linshan， Zhang Qing， Zheng Xiang， et al. The empirical optimal envelop and its application to local mean decomposition[J]. Digital Signal Processing， 2019，87： 166-177.

[12] WU Zhaohua， Huang N E. Ensemble empirical mode decomposition： A noise-assisted data analysis method [J]. Advances in Adaptive Data Analysis， 2009， 1（1）： 1-41.

[13] 卓仁雄， 肖金鳳. 基于改進(jìn)的集合經(jīng)驗(yàn)?zāi)B(tài)分解的電動(dòng)機(jī)滾動(dòng)軸承故障診斷研究[J]. 機(jī)械制造與自動(dòng)化，2019，48（1）：36-39.

Zhuo Renxiong， Xiao Jinfeng. Research on fault diagnosis method of motor bearing based on improved EEMD and SVM[J]. Machinery Manufacturing and Automation， 2019，48（1）：36-39.

[14] 程軍圣， 王 ?健， 桂 ?林. 改進(jìn)的EEMD方法及其在滾動(dòng)軸承故障診斷中的應(yīng)用[J]. 振動(dòng)與沖擊，2018，37（16）：51-56.

Cheng Junshegn， Wang Jian， Gui Lin. An improved EEMD method and its application in rolling bearing fault diagnosis[J]. Journal of Vibration and Shock， 2018，37（16）：51-56.

[15] 孔德同， 劉慶超， 雷亞國(guó)， 等. 一種改進(jìn)的EEMD方法及其應(yīng)用研究 [J]. 振動(dòng)工程學(xué)報(bào)，2015，28（6）：1015-1021.

Kong Detong， Liu Qingchao， Lei Yaguo， et al. The improved EEMD method and its application[J]. Journal of Vibration Engineering， 2015，28（6）：1015-1021.

[16] Yeh J R， Shieh J S， Huang N E. Complementary ensemble empirical mode decomposition： A novel noise enhanced data analysis method[J]. Advances in Adaptive Data Analysis， 2014， 2（2）：1000042.

[17] Torres Maria E， Colominas M A， Gastón Schlotthauer， et al. A complete ensemble empirical mode decomposition with adaptive noise [C]. 2011 IEEE International Conference on Acoustics， Speech and Signal Processing （ICASSP）， Prague， Czech Republic， 2011： 4144?4147.

[18] Colominas M A， Gastón Schlotthauer， Torres María E. Improved complete ensemble EMD： A suitable tool for biomedical signal processing[J]. Biomedical Signal Processing and Control， 2014， 14： 19-29.

[19] Wang Lei， Liu Zhiwen， Miao Qiang， et al. Complete ensemble local mean decomposition with adaptive noise and its application to fault diagnosis for rolling bearings [J]. Mechanical Systems and Signal Processing， 2018， 106： 24-39.

[20] Wang Jun， Luo Yuyan， Tang Lingyu， et al. A new weighted CEEMDAN-based prediction model： An experimental investigation of decomposition and non-decomposition approaches[J]. Knowledge-Based Systems， 2018，160： 188-199.

[21] Hassan A R， Bhuiyan M I H. Computer-aided sleep staging using complete ensemble empirical mode decomposition with adaptive noise and bootstrap aggregating [J]. Biomedical Signal Processing and Control， 2016， 24： 1-10.

[22] 鄭近德， 潘海洋， 程軍圣. 均值優(yōu)化經(jīng)驗(yàn)?zāi)B(tài)分解及其在轉(zhuǎn)子故障診斷中的應(yīng)用[J]. 機(jī)械工程學(xué)報(bào)，2018，54（23）：93-101.

Zheng Jinde， Pan Haiyang， Cheng Junsheng. Mean-optimized empirical mode decomposition and its application in rotor fault diagnosis[J]. Journal of Mechanical Engineering， 2018，54（23）：93-101.

[23] 鄭近德， 程軍圣， 楊宇. 基于改進(jìn)的ITD和模糊熵的滾動(dòng)軸承故障診斷方法[J]. 中國(guó)機(jī)械工程， 2012， 23（19）： 2372-2377.

Zheng Jinde， Cheng Junsheng， Yang Yu. A rolling bearing fault diagnosis method based on improved ITD and fuzzy entroy[J]. China Mechanical Engineering， 2012， 23（19）： 2372-2377.

[24] Lu Siliang， He Qingbo， Wang Jun. A review of stochastic resonance in rotating machine fault detection[J].Mechanical Systems and Signal Processing， 2019， 116：230-260.