李兵，徐榕，賈春寧，郭清晨

(1.軍械工程學(xué)院 四系，河北 石家莊050003;2.總裝備部駐上海地區(qū)軍事代表室，上海201109)

0 引言

振動(dòng)信號(hào)分析是對(duì)發(fā)動(dòng)機(jī)進(jìn)行結(jié)構(gòu)、故障分析和狀態(tài)監(jiān)測(cè)最為廣泛、也是最行之有效的方法之一。其中，振動(dòng)信號(hào)處理與特征參數(shù)提取是實(shí)現(xiàn)發(fā)動(dòng)機(jī)故障診斷的關(guān)鍵與核心［1］。

小波分析具有多尺度特性和“數(shù)學(xué)顯微”特性，因此近年來(lái)小波分析被越來(lái)越廣泛地應(yīng)用于發(fā)動(dòng)機(jī)故障診斷中并取得了很好的效果［2-7］。然而傳統(tǒng)的小波分析是基于頻率的線性分解，對(duì)于具有特定非線性形態(tài)特征的信號(hào)，應(yīng)用傳統(tǒng)的線性小波分析方法有時(shí)得不到很好的效果。

形態(tài)小波變換是以數(shù)學(xué)形態(tài)學(xué)為基礎(chǔ)的一種小波變換，是小波理論非線性擴(kuò)展研究的一個(gè)方向，它是由Goutsias 和Heijmans 在2000年時(shí)首先提出的［8-9］，它成功地將大多數(shù)線性小波和非線性小波統(tǒng)一起來(lái)，并形成了多分辨分析的統(tǒng)一框架。作為一種非線性小波，形態(tài)小波兼顧數(shù)學(xué)形態(tài)學(xué)的形態(tài)特性與小波的多分辨率特性，具有良好的細(xì)節(jié)保留和抗噪聲性能。目前，形態(tài)小波已經(jīng)在圖像壓縮、圖像融合以及電力系統(tǒng)故障檢測(cè)等方面得到了廣泛的應(yīng)用［10-13］，在機(jī)械設(shè)備振動(dòng)信號(hào)中的分析也開始嶄露頭角［14-16］。

盡管形態(tài)小波具有很多優(yōu)良的性質(zhì)，但也存在一個(gè)明顯的限制，那就是濾波器結(jié)構(gòu)是固定的，即濾波算子對(duì)信號(hào)全局來(lái)講是固定不變的;而信號(hào)的突變處往往要求具有不同于平緩處的濾波特性，需要能夠根據(jù)信號(hào)局部特性的不同而對(duì)濾波算子進(jìn)行相應(yīng)的改進(jìn)［17-20］。

因此，本文在形態(tài)小波和自適應(yīng)小波理論的基礎(chǔ)上，提出了一種自適應(yīng)形態(tài)梯度提升小波分解方案(AMGLW)，在有效保留反映發(fā)動(dòng)機(jī)故障特征的信號(hào)分量的前提下可以同時(shí)抑制噪聲;在此基礎(chǔ)上，提出采用改進(jìn)非負(fù)矩陣分解技術(shù)(INMF)對(duì)信號(hào)進(jìn)行特征提取，計(jì)算用于發(fā)動(dòng)機(jī)故障診斷的特征參數(shù)集。采用實(shí)測(cè)的發(fā)動(dòng)機(jī)在5 種狀態(tài)下的振動(dòng)信號(hào)對(duì)本文提出的信號(hào)處理與特征參數(shù)計(jì)算方法進(jìn)行了驗(yàn)證，并與傳統(tǒng)的方法進(jìn)行了對(duì)比。

1 自適應(yīng)形態(tài)提升小波

1.1 形態(tài)小波的概念

將線性小波中的線性濾波器用非線性形態(tài)濾波器代替，就可以構(gòu)成一類非線性形態(tài)小波，主要分為對(duì)偶小波和非對(duì)偶小波。對(duì)偶小波包含2 個(gè)分析算子和一個(gè)合成算子，分析算子一個(gè)作用于信號(hào)，另一個(gè)作用于細(xì)節(jié)信號(hào)。非對(duì)偶小波是對(duì)偶小波的一個(gè)特殊情況，它包含2 個(gè)合成算子，信號(hào)合成算子與細(xì)節(jié)信號(hào)合成算子，線性正交小波屬于非對(duì)偶小波，其信號(hào)分析算子與細(xì)節(jié)分析算子分別對(duì)應(yīng)于低通和高通濾波器;有關(guān)形態(tài)小波的具體概念和詳細(xì)闡述可參見文獻(xiàn)［9］。

1.2 自適應(yīng)形態(tài)提升小波

Piella 等［17-20］提出了一種根據(jù)信號(hào)局部梯度的自適應(yīng)更新方案，其框架如圖1所示。

圖1中D 為判決函數(shù)，λd為更新算子，π 為預(yù)測(cè)算子;φ0為信號(hào)合成算子.⊕d、Θd分別為廣義加法和減法算子。在標(biāo)準(zhǔn)提升方法的更新過程中，更新算子和加、減法算子固定不變。

雖然Piella 等提出的自適應(yīng)提升小波具有非冗余性和完美重構(gòu)等優(yōu)點(diǎn)，但對(duì)濾波器的限制十分嚴(yán)重，而且其本質(zhì)仍然是線性小波變換，很難采用更加靈活的方式構(gòu)造非線性小波，比如形態(tài)小波。

圖1 自適應(yīng)更新提升方案Fig.1 Adaptive update lifting scheme

本文提出了一種自適應(yīng)形態(tài)梯度提升小波(AMGLW)方案，其結(jié)構(gòu)圖與圖1一樣，與Piella 等提出的方法不同，本文采用形態(tài)梯度提升算子和平滑濾波器構(gòu)造自適應(yīng)形態(tài)提升小波變換。

在判斷信號(hào)奇異點(diǎn)上采用與Piella 等相同的方法，即梯度法，信號(hào)的局部梯度值定義為

判別函數(shù)定義為

式中T 為閾值。

當(dāng)d(n)=1 時(shí)，采用形態(tài)梯度更新算子，即

提升后的信號(hào)和細(xì)節(jié)系數(shù)分別為

當(dāng)d(n)=0 時(shí)，采用平滑算子，令提升后的信號(hào)和細(xì)節(jié)系數(shù)分別為

由上述各式可以看出，當(dāng)信號(hào)局部梯度值較大時(shí)即信號(hào)出現(xiàn)沖擊特征時(shí)，采用梯度提升算子，保留信號(hào)的沖擊特征;當(dāng)信號(hào)較為平滑時(shí)，采用均值濾波，以有效地濾除噪聲。

2 改進(jìn)非負(fù)矩陣分解

雖然AMGLW 能夠有效地提取出原始信號(hào)中包含的故障信息，但處理后的空間維數(shù)仍然非常高，不能直接用來(lái)進(jìn)行分類。為節(jié)省存儲(chǔ)空間，降低計(jì)算復(fù)雜度，避免“維數(shù)災(zāi)難”，必須對(duì)信號(hào)進(jìn)行進(jìn)一步的特征提取。

非負(fù)矩陣分解(NMF)算法是1999年由Lee 和Seung［21］在《Nature》上提出的一種新的特征提取方法。非負(fù)矩陣分解的心理學(xué)和生理學(xué)構(gòu)造依據(jù)是對(duì)整體的感知由對(duì)組成整體的部分的感知構(gòu)成的(純加性的)，這也符合直觀的理解:整體是由部分組成的，因此它在某種意義上抓住了智能數(shù)據(jù)描述的本質(zhì)。此外，這種非負(fù)性的限制導(dǎo)致了相應(yīng)描述在一定程度上的稀疏性，稀疏性的表述已被證明是介于完全分布式的描述和單一活躍分量的描述之間的一種有效數(shù)據(jù)描述形式。非負(fù)矩陣分解算法實(shí)現(xiàn)簡(jiǎn)便，而且具有可解釋性和明確的物理意義，已經(jīng)引起許多科學(xué)家和研究人員的廣泛重視。

2.1 非負(fù)矩陣分解主要思想

非負(fù)矩陣分解的主要思想為:已知非負(fù)矩陣V，尋找適當(dāng)?shù)姆秦?fù)矩陣因子W 和H，使得

即給定數(shù)據(jù)向量集合Vn×m，其中n 為數(shù)據(jù)樣本的維數(shù)，m 為集合中數(shù)據(jù)樣本的個(gè)數(shù)，這個(gè)矩陣可以近似的分解為矩陣Wn×r和矩陣Hr×m的乘積。一般情況下，r 的選擇要滿足(n+m)r ＜nm，從而W 和H的維數(shù)都會(huì)小于原始矩陣V，由此可得實(shí)現(xiàn)對(duì)原始數(shù)據(jù)維數(shù)壓縮。

2.2 非負(fù)矩陣分解的算法實(shí)現(xiàn)

2001年Lee 和Seung 提出了簡(jiǎn)化的計(jì)算非負(fù)矩陣分解的2 種算法［22］。一種算法以最小化剩余的Frobenius-Norm 為目標(biāo)函數(shù);另一種算法則是最小化修正的Kullback-Laebler 散度為目標(biāo)函數(shù)。

本文采用最小化修正的Kullback-Laebler 散度為目標(biāo)函數(shù)，給定非負(fù)矩陣Vn×m，尋找矩陣Wn×r和矩陣Hr×m，使得矩陣V 和矩陣WH 的Kullback-Laebler 散度最小，因此非負(fù)矩陣分解算法可轉(zhuǎn)化為如下的帶約束的優(yōu)化問題:

Lee 和Seung 等提出了采用乘性迭代規(guī)則對(duì)上述優(yōu)化問題進(jìn)行求解，為了消除尺度對(duì)于基矩陣的影響，對(duì)基矩陣W 限制其L-1 范數(shù)為1.

對(duì)應(yīng)于(10)式的乘性迭代規(guī)則為迭代公式如下所示:

由(11)式～(12)式可以看出，在上述算法的每一步迭代過程中，W 和H 的新值可以通過其當(dāng)前值與一些因子的乘積來(lái)獲得，(13)式保證了基矩陣的范數(shù)為1.

2.3 非負(fù)矩陣分解的改進(jìn)

原始的非負(fù)矩陣分解技術(shù)是一種非監(jiān)督的學(xué)習(xí)算法，因此由非負(fù)矩陣分解技術(shù)直接得出的特征參數(shù)并不一定取得很好的分類效果，這已經(jīng)在某些研究中得到了驗(yàn)證。文獻(xiàn)［23 -24］提出了一種鑒別非負(fù)矩陣分解方法，該方法將鑒別分析直接引入到NMF 的優(yōu)化過程中，取得了比原始NMF 更好的分類精度，但該方法存在著算法無(wú)法收斂的問題;文獻(xiàn)［25］將NMF 矩陣分解后得到的編碼系數(shù)進(jìn)一步進(jìn)行Fisher 線性鑒別分析，同樣取得了較好的分類效果。文獻(xiàn)［26］提出了一種更為簡(jiǎn)單的改進(jìn)方法，該方法在NMF 分解后采用Fisher 判別準(zhǔn)則對(duì)分解后的基向量進(jìn)行選擇，在提高分類精度的同時(shí)進(jìn)一步壓縮了特征維數(shù)。本文采用這種方法對(duì)NMF 進(jìn)行改進(jìn)。

該方法的主要思想為，NMF 算法收斂后得到基向量矩陣Wn×r和系數(shù)矩陣Hr×m(特征矩陣)，計(jì)算系數(shù)矩陣中每個(gè)特征的Fisher 判別度:

式中:C 為樣本類別數(shù);mik和mjk分別為第i 類和第j類樣本第k 個(gè)特征的均值;σik和σjk分別為第i 類和第j 類樣本第k 個(gè)特征的方差。Fisher 判別度越大，證明該特征的區(qū)分性能越高。

因此，根據(jù)上式計(jì)算的F 值對(duì)所有r 個(gè)特征進(jìn)行排序，選擇前d 個(gè)特征所對(duì)應(yīng)的基向量組成改進(jìn)后的基向量矩陣W'n×r，由此可以去除原基向量矩陣中與分類信息無(wú)關(guān)的基向量，不但進(jìn)一步壓縮了特征向量的維數(shù)，而且可以提高特征子集的分類精度。

3 發(fā)動(dòng)機(jī)故障診斷

3.1 發(fā)動(dòng)機(jī)試驗(yàn)

實(shí)驗(yàn)在某型發(fā)動(dòng)機(jī)上進(jìn)行，加速度傳感器粘貼在第一缸缸蓋處，實(shí)驗(yàn)在空間較大的車間內(nèi)進(jìn)行，計(jì)算機(jī)、電荷放大器等裝置在車間里的控制間內(nèi)放置。實(shí)驗(yàn)時(shí)發(fā)動(dòng)機(jī)轉(zhuǎn)速為1 200 r/min，采樣頻率為40 kHz，采樣點(diǎn)數(shù)為發(fā)動(dòng)機(jī)一個(gè)工作周期。為結(jié)合發(fā)動(dòng)機(jī)故障診斷研究的需要，實(shí)驗(yàn)中除采集正常運(yùn)行工況的數(shù)據(jù)外，還根據(jù)實(shí)際的實(shí)驗(yàn)條件，采用人為手動(dòng)設(shè)置故障的方式，對(duì)部分發(fā)動(dòng)機(jī)典型的故障模式進(jìn)行模擬。實(shí)驗(yàn)中設(shè)置發(fā)動(dòng)機(jī)在5 種狀態(tài)下工作，分別為正常、1 缸失火、2 缸失火、進(jìn)氣門間隙過大、排氣門間隙過大。實(shí)驗(yàn)中在發(fā)動(dòng)機(jī)每種運(yùn)行狀態(tài)下采集20 個(gè)樣本，共100 個(gè)樣本。

圖2 發(fā)動(dòng)機(jī)在正常狀態(tài)下的振動(dòng)信號(hào)Fig.2 Vibration signal from engine with normal state

圖2為發(fā)動(dòng)機(jī)在正常狀態(tài)下缸蓋振動(dòng)信號(hào)的時(shí)域波形，圖中給出了信號(hào)中主要沖擊分量所應(yīng)的發(fā)動(dòng)機(jī)工作過程。由于傳感器安裝在發(fā)動(dòng)機(jī)1 缸體上，因此信號(hào)中的主要分量為1 缸燃爆信號(hào)、1 缸進(jìn)、排氣門落座信號(hào)和2 缸燃爆信號(hào)。

3.2 基于AMGLW 的發(fā)動(dòng)機(jī)信號(hào)處理

本文分別采用傳統(tǒng)的線性小波變換(采用db5小波)、極大形態(tài)提升小波、Piella 提出的自適應(yīng)提升小波(ALW)和本文提出的自適應(yīng)形態(tài)梯度提升小波(AMGLW)對(duì)發(fā)動(dòng)機(jī)正常信號(hào)進(jìn)行分解，分解層數(shù)為3.在分解中，設(shè)定判別函數(shù)的閾值為信號(hào)梯度最大值的1/2.

圖3 圖2中信號(hào)4 種小波分解的第3 層近似系數(shù)Fig.3 Third level approximation coefficient of signal in Fig.2 deposed by four type wavelet schemes

圖3給出了采用前述的4 種小波分解方法對(duì)信號(hào)進(jìn)行3 層分解后的第3 層近似信號(hào)?？梢园l(fā)現(xiàn)，本文提出的自適應(yīng)形態(tài)梯度提升小波具有最好的沖擊信號(hào)提取效果。

圖4為發(fā)動(dòng)機(jī)在5 種狀態(tài)下的振動(dòng)信號(hào)。圖5為采用AMGLW 對(duì)發(fā)動(dòng)機(jī)5 種狀態(tài)下振動(dòng)信號(hào)進(jìn)行3 層分解后近似系數(shù)的波形圖，從中可以看出，AMGLW 在大幅壓縮信號(hào)長(zhǎng)度的前提下，不僅有效地保留了反映發(fā)動(dòng)機(jī)工作狀態(tài)地沖擊特征信號(hào)，而且大大降低了噪聲的干擾，為準(zhǔn)確有效地判斷發(fā)動(dòng)機(jī)工作狀態(tài)提供了堅(jiān)實(shí)的基礎(chǔ)。

3.3 特征提取

采用AMGLW 對(duì)發(fā)動(dòng)機(jī)故障信號(hào)進(jìn)行3 層分解后，可以得到得到一個(gè)512 ×100 的樣本集。

每類狀態(tài)下選擇10 個(gè)樣本作為訓(xùn)練樣本，這樣可以得到訓(xùn)練矩陣X512×50，表1給出了發(fā)動(dòng)機(jī)故障信號(hào)經(jīng)AMGLW 處理后的部分訓(xùn)練樣本向量集，表中每行對(duì)應(yīng)發(fā)動(dòng)機(jī)的一種工作狀態(tài)。

圖4 發(fā)動(dòng)機(jī)在5 種狀態(tài)下的振動(dòng)信號(hào)Fig.4 Vibration signals from engine with five states

在對(duì)發(fā)動(dòng)機(jī)訓(xùn)練樣本集X 采用NMF 技術(shù)進(jìn)行分解之前，首先采用PCA 技術(shù)確定基向量秩r，經(jīng)PCA 分析后發(fā)現(xiàn)，前40 個(gè)特征值的累積貢獻(xiàn)率超過了99%，因此在對(duì)發(fā)動(dòng)機(jī)故障信號(hào)進(jìn)行非負(fù)矩陣分解時(shí)選擇基向量秩為r=40.

在訓(xùn)練過程開始之前，采用PCA 對(duì)非負(fù)矩陣分解的基向量和系數(shù)向量進(jìn)行初始化。采用PCA 對(duì)X 進(jìn)行分解后，取前40 個(gè)特征值對(duì)應(yīng)的特征向量的絕對(duì)值作為初始化的基矩陣W，取X 在W 上的投影系數(shù)的絕對(duì)值作為初始化的系數(shù)矩陣H.非負(fù)矩陣分解訓(xùn)練迭代次數(shù)設(shè)置為200.

訓(xùn)練結(jié)束后可得到基向量矩陣W512×40和編碼矩陣H40×50，表2給出了由非負(fù)矩陣分解得到的40 個(gè)基向量中部分基向量，表3給出了由非負(fù)矩陣分解得到的部分編碼系數(shù)，即特征參數(shù)，表3中的25 個(gè)數(shù)據(jù)與表1中25 個(gè)樣本為一一對(duì)應(yīng)關(guān)系。

圖5 采用AMGLW 對(duì)圖4信號(hào)進(jìn)行3 層分解后近似系數(shù)Fig.5 Third level approximation coefficients of signals in Fig.4 deposed by AMGLW

圖6給出了由(14)式計(jì)算出的由非負(fù)矩陣分解后得到的40 個(gè)特征值的Fisher 判別度，可以看出部分特征值的Fisher 判別度明顯高于其他特征值，因此，本文選擇前20 個(gè)Fisher 判別度較大的特征值所對(duì)應(yīng)的基向量作為改進(jìn)后的基向量矩陣W'.

對(duì)于訓(xùn)練樣本來(lái)說，可用編碼矩陣對(duì)其進(jìn)行描述，對(duì)于剩余的測(cè)試樣本，可采用下式得到編碼矩陣:

由此可得到描述發(fā)動(dòng)機(jī)故障信號(hào)的特征參數(shù)集。本文將基于AMGLW 和NMF 技術(shù)得到的特征參數(shù)集記為FAMGLW-NMF，將基于改進(jìn)NMF 技術(shù)所得到的特征參數(shù)集記為FAMGLW-INMF.

表1 非負(fù)矩陣分解部分訓(xùn)練樣本集示意圖(每行代表一種狀態(tài))Tab.1 Part of the training samples for NMF (Each row represents one engine state)

表2 非負(fù)矩陣分解得到的基向量(部分)示意圖Tab.2 Part of the basic vectors acquired by NMF

表3 非負(fù)矩陣分解得到的編碼向量示意圖(每行代表一種狀態(tài))Tab.3 Coding vectors acquired by NMF (Each row represents one engine state)

圖6 非負(fù)矩陣分解40 個(gè)特征值的Fisher 判別度Fig.6 Fisher discriminate degree of the forty features acquired by NMF

3.4 發(fā)動(dòng)機(jī)故障分類

為證明提出的特征提取方法的有效性和優(yōu)越性，本文還采用廣泛應(yīng)用的統(tǒng)計(jì)特征來(lái)提取信號(hào)的特征，本文采用了4 個(gè)有量綱統(tǒng)計(jì)特征參數(shù)(均值、標(biāo)準(zhǔn)差、均方根和峰-峰值)和6 個(gè)無(wú)量綱特征參數(shù)(偏度、峭度、峰值指標(biāo)、脈沖指標(biāo)、波形指標(biāo)和裕度指標(biāo))作為統(tǒng)計(jì)特征參量，分別計(jì)算原始信號(hào)、包絡(luò)信號(hào)和小波3 層分解后各層系數(shù)的統(tǒng)計(jì)特征參數(shù)作為一個(gè)特征參數(shù)集，其維數(shù)為60，記該特征參數(shù)集為FSTAT.此外，為驗(yàn)證AMGLW 對(duì)發(fā)動(dòng)機(jī)信號(hào)處理的效果，本文還計(jì)算了AMGLW 3 層分解后各層系數(shù)的統(tǒng)計(jì)特征參數(shù)，其維數(shù)為40，記為FAMGLW-STAT.

采用前述各種特征參數(shù)提取方法，可得到描述發(fā)動(dòng)機(jī)故障信號(hào)的4 個(gè)特征參數(shù)集，即FSTAT、FAMGLW-STAT、FAMGLW-NMF和FAMGLW-INMF.為驗(yàn)證結(jié)果的通用性，采用3 種常見的分類器，即最近鄰分類器(KNNC)、人工神經(jīng)網(wǎng)絡(luò)(BP-NN)和支持向量機(jī)分類器(SVM)對(duì)發(fā)動(dòng)機(jī)的5 種狀態(tài)進(jìn)行分類。

在分類實(shí)驗(yàn)中，從發(fā)動(dòng)機(jī)每類狀態(tài)下隨機(jī)選擇10 個(gè)樣本作為訓(xùn)練樣本，剩余的10 個(gè)樣本作為測(cè)試樣本，為保證結(jié)果的有效性，將此隨機(jī)選擇過程重復(fù)50 次，并將50 次運(yùn)算結(jié)果的平均值作為評(píng)價(jià)標(biāo)準(zhǔn)。表4給出了前述4 個(gè)特征子集的分類精度。

表4 各特征子集的發(fā)動(dòng)機(jī)故障診斷精度Tab.4 Engine fault diagnosis accuracy of different feature subsets

由表4可以看出，在發(fā)動(dòng)機(jī)故障診斷中，本文所提出的基于AMGLW 和改進(jìn)非負(fù)矩陣分解的特征子集FAMGLW-INMF取得了最好的分類精度，而且特征子集的維數(shù)也比傳統(tǒng)的統(tǒng)計(jì)特征參數(shù)FSTAT、FAMGLW-STAT和FAMGLW-NMF低。這一方面是由于AMGLW 能夠有效地提取振動(dòng)信號(hào)中能夠反映發(fā)動(dòng)機(jī)工作狀態(tài)的有用分量，另一方面是因?yàn)楦倪M(jìn)的非負(fù)矩陣分解不僅可以去除原基向量矩陣中與分類信息無(wú)關(guān)的基向量，還進(jìn)一步降低了特征向量的維數(shù)，因此可以以較少的特征參數(shù)獲得較高的故障診斷精度。

4 結(jié)論

1)提出了一種自適應(yīng)形態(tài)梯度提升小波(AMGLW)方法。針對(duì)傳統(tǒng)小波變換濾波器結(jié)構(gòu)固定的缺點(diǎn)，采用信號(hào)的局部梯度作為判斷信號(hào)奇異性的度量指標(biāo)，在信號(hào)突變處采用提出的形態(tài)梯度提升算子以保留信號(hào)的沖擊特征，在信號(hào)緩變處采用平滑算子以抑制噪聲。

在發(fā)動(dòng)機(jī)故障信號(hào)中的應(yīng)用表明，AMGLW 能在強(qiáng)噪聲和強(qiáng)干擾條件下用較少的小波系數(shù)有效地提取信號(hào)中的沖擊特征，因此可有效地提取振動(dòng)信號(hào)中能夠反映發(fā)動(dòng)機(jī)工作狀態(tài)的有用分量。

2)在采用AMGLW 對(duì)發(fā)動(dòng)機(jī)振動(dòng)信號(hào)處理的基礎(chǔ)上，采用改進(jìn)非負(fù)矩陣分解技術(shù)對(duì)信號(hào)進(jìn)行特征參數(shù)提取，實(shí)際的故障診斷結(jié)果表明，與傳統(tǒng)的信號(hào)處理與特征提取技術(shù)相比，本文提出的基于AMGLW 與改進(jìn)非負(fù)矩陣分解的特征子集不僅維數(shù)較低，而且具有更高的分類精度。

本文提出的信號(hào)處理與特征提取方法為準(zhǔn)確判斷發(fā)動(dòng)機(jī)工作狀態(tài)提供了一種新的行之有效的技術(shù)手段，并可推廣至其他機(jī)械設(shè)備比如齒輪、軸承等設(shè)備的故障診斷中去。

References)

［1］ Wu J D，Chen J C.Continuous wavelet transform technique for fault signal diagnosis of internal combustion engines［J］.NDT and E International，2006，39(4):304 -311.

［2］ 王仲生，何紅，陳錢.小波分析在發(fā)動(dòng)機(jī)早期故障識(shí)別中的應(yīng)用研究［J］.西北工業(yè)大學(xué)學(xué)報(bào)，2006，24(1):68 -71.WANG Zhong-sheng，HE Hong，CHEN Qian.Exploring effective early identification of aero-engine rotor faults ［J］.Journal of Northwestern Polytechnical University，2006，24(1):68 - 71.(in Chinese)

［3］ 曹建軍，張培林，張英堂，等.基于提升小波包變換的發(fā)動(dòng)機(jī)缸蓋振動(dòng)信號(hào)特征提?。跩］.振動(dòng)與沖擊，2008，27(2):34-38.CAO Jian-jun，ZHANG Pei-lin，ZHANG Ying-tang，et al.Feature extraction of an engine cylinder head vibration signal based on lifting wavelet package transformation ［J］.Journal of Vibration and Shock，2008，27(2):34 -38.(in Chinese)

［4］ 吳定海，張培林，任國(guó)全，等.基于雙樹復(fù)小波包的發(fā)動(dòng)機(jī)振動(dòng)信號(hào)特征提取研究［J］.振動(dòng)與沖擊，2010，29(4):160 -165.WU Ding-hai，ZHANG Pei-lin，REN Guo-quan，et al.Feature extraction of an engine vibration signal based on dual-tree wavelet package transformation ［J］.Journal of Vibration and Shock，2010，29(4):160 -165.(in Chinese)

［5］ Wu J D，Liu C H.Investigation of engine fault diagnosis using discrete wavelet transform and neural network［J］.Expert Systems with Applications，2008，35(3):1200 -1213.

［6］ Wu J D，Liu C H.An expert system for fault diagnosis in internal combustion engines using wavelet packet transform and neural network［J］.Expert Systems with Applications，2009，36(3 Part 1):4278 -4286.

［7］ Zhou R，Bao W，Li N，et al.Mechanical equipment fault diagnosis based on redundant second generation wavelet packet transform［J］.Digital Signal Processing:a Review Journal，2010，20(1):276 -288.

［8］ John Goutsias，Henk J A M Heijmans.Nonlinear multiresolution signal decomposition schemes-Part I:morphological pyramids［J］.IEEE Transactions on Image Processing，2000，9(11):1862 -1876.

［9］ Henk J A M Heijmans，John Goutsias.Nonlinear multiresolution signal decomposition schemes-Part II:morphological wavelets［J］.IEEE Transactions on Image Processing，2000，9(11):1897 -1913.

［10］ Flouzat G，Amram O，Laporterie F，et al.Multiresolution analysis and reconstruction by a morphological pyramid in the remote sensing of terrestrial surfaces［J］.Signal Processing，2001，81(10):2171 -2185.

［11］ Lin P L，Huang P Y.Fusion methods based on dynamic-segmented morphological wavelet or cut and paste for multifocus images［J］.Signal Processing，2008，88(6):1511 -1527.

［12］ Zhang J F，Smith J S，Wu Q H.Morphological undecimated wavelet decomposition for fault location on power transmission lines［J］.IEEE Transactions on Circuits and Systems I:Regular Papers，2006，53(6):1395 -1402.

［13］ Ji T Y，Lu Z，Wu Q H.Detection of power disturbances using morphological gradient wavelet［J］.Signal Processing，2008，88(2):255 -267.

［14］ 章立軍，陽(yáng)建宏，徐金梧，等.形態(tài)非抽樣小波及其在沖擊信號(hào)特征提取中的應(yīng)用［J］.振動(dòng)與沖擊，2007，26(10):56-60.ZHANG Li-jun，YANG Jian-hong，XU Jin-wu，et al.Morphological undecimated wavelet and its application to feature extraction of impulsive signal［J］.Journal of Vibration and Shock，2007，26(10):56 -60.(in Chinese)

［15］ Hao R，Chu F.Morphological undecimated wavelet decomposition for fault diagnostics of rolling element bearings［J］.Journal of Sound and Vibration，2009，320(4 -5):1164 -1177.

［16］ 沈路，周曉軍，劉莉，等.形態(tài)小波降噪方法在齒輪故障特征提取中的應(yīng)用［J］.農(nóng)業(yè)機(jī)械學(xué)報(bào)，2010，41(4):217 -221.SHEN Lu，ZHOU Xiao-jun，LIU Li，et al.Application of morphological wavelet de-noising in extracting gear fault feature［J］.Transactions of the Chinese Society for Agricultural Machinery，2010，41(4):217 -221.(in Chinese)

［17］ Piella G，Heijmans H J A M.Adaptive lifting schemes with perfect reconstruction［J］.IEEE Transactions on Signal Processing，2002，50(7):1620 -1630.

［18］ Piella G，Pesquet-Popescu B，Heijmans H.Adaptive update lifting with a decision rule based on derivative filters［J］.IEEE Signal Processing Letters，2002，9(10):329 -332.

［19］ Heijmans H J A M，Pesquet-Popescu B，Piella G.Building nonredundant adaptive wavelets by update lifting［J］.Applied and Computational Harmonic Analysis，2005，18(3):252 -281.

［20］ Piella G，Pesquet-Popescu B，Heijmans H J A M.Gradient-driven update lifting for adaptive wavelets［J］.Signal Processing:Image Communication，2005，20(9 -10):813 -831.

［21］ Lee D D，Seung H S.Learning the parts of objects by non-negative matrix factorization［J］.Nature，1999，401(6755):788 -791.

［22］ Lee D D，Seung H S.Algorithms for non-negative matrix factorization［J］.Advances in Neural Information Processing Systems，2001，13:556 -562.

［23］ Buciu I，Pitas I.A new sparse image representation algorithm applied to facial expression recognition［C］∥Machine Learning for Signal Processing XIV - Proceedings of the 2004 IEEE Signal Processing Society Workshop.Sao Luis:IEEE，2004:539 -548.

［24］ Buciu I，Pitas I.NMF，LNMF，and DNMF modeling of neural receptive fields involved in human facial expression perception［J］.Journal of Visual Communication and Image Representation，2006，17(5):958 -969.

［25］ Guimet F，Boqué R，F(xiàn)erré J.Application of non-negative matrix factorization combined with fisher's linear discriminant analysis for classification of olive oil excitation-emission fluorescence spectra［J］.Chemometrics and Intelligent Laboratory Systems，2006，81(1):94 -106.

［26］ Cho Y C，Choi S.Nonnegative features of spectro-temporal sounds for classification［J］.Pattern Recognition Letters，2005，26(9):1327 -1336.