向 丹，葛 爽

(1.廣東技術(shù)師范學(xué)院 自動化學(xué)院，廣州 510635;2.華南理工大學(xué) 機械與汽車工程學(xué)院，廣州 510640)

利用振動信號對機械狀態(tài)進行監(jiān)測是機械故障診斷技術(shù)中的一種有效手段，而故障特征的提取是故障診斷的關(guān)鍵，關(guān)系到故障診斷的準(zhǔn)確性和早期預(yù)報的可能性。因此，人們從各種角度提取振動信號的特征參數(shù)，包括時域、頻域和時頻域等，取得了一定的效果［1－2］。但由于機械故障信號的非線性、非平穩(wěn)特性，使得傳統(tǒng)的基于線性系統(tǒng)的特征提取方法不能有效提取故障的非線性特征。因此，人們開始探索非線性動力學(xué)參數(shù)的特征提取方法。

近年來，人們將熵的概念，包括樣本熵、能量熵、奇異熵和近似熵等，作為特征參數(shù)提取的方法引入到故障診斷領(lǐng)域。YAN等［3］將近似熵概念引入到軸承狀態(tài)監(jiān)測上，并取得了很好的效果。Pincus［4］提出了近似熵的改進算法樣本熵，并廣泛用于生理信號處理上。趙志宏等［5］提出了一種基于經(jīng)驗?zāi)B(tài)分解與樣本熵的軸承故障診斷方法。流形學(xué)習(xí)作為一種非線性的數(shù)據(jù)

維數(shù)約簡方法，能有效地挖掘高維非線性數(shù)據(jù)內(nèi)在幾何分布特征，具有很好的非線性復(fù)雜信息處理能力，已被用于機械故障診斷的提取領(lǐng)域。粟茂林等［6］運用局部切空間排列算法對最優(yōu)小波系統(tǒng)矩陣進行優(yōu)化，有效地提取了引發(fā)故障的沖擊成分。葛爽等［7］提出一種基于擴散映射與支持向量機的能量耗損信號分析，對于能量耗損信號分析非常有效。

雖然故障診斷理論和方法的研究取得了很大的進步，但對于不確定、非線性的復(fù)雜故障診斷問題，采用傳統(tǒng)的統(tǒng)計和信號處理的特征提取方法的效果十分有限?；诖耍疚奶岢隽嘶谛〔ò鼧颖眷睾土餍螌W(xué)習(xí)的故障特征提取模型，并通過試驗驗證該模型的有效性。

1 小波包變換與樣本熵

1.1 小波包變換

根據(jù)多分辨率分析理論，定義小波包變換兩個函數(shù)［8］。

其中φ(t)與φ(t)分別為尺度函數(shù)與小波函數(shù)。

在正交的情況下，函數(shù)wm(t)(m=0，1，2)分別為

式中，j為尺度參數(shù)，wj，m，n為小波包原子。與小波變換不同，小波包變換將信號的高頻部分與低頻部分分別進行分解。小波包分解的結(jié)構(gòu)圖如圖1所示。

圖1 3層小波包分解結(jié)構(gòu)圖Fig.1 3 layer wavelet packet decomposition structure diagram

圖中S(K)表示原始信號，a表示低頻，d表示高頻，末尾數(shù)表示小波包分解的層數(shù)，即尺度數(shù)。

1.2 樣本熵

設(shè)N個數(shù)據(jù)組成的時間序列為x(n)，樣本熵的計算方法如下［4］:

(1)將時間序列按序號組成一個維數(shù)為m的向量;

(2)定義向量 Xm(i)與 Xm(j)之間的距離d［Xm(i)，Xm(j)］是兩個對應(yīng)元素中最大差值的絕對值，即

(3)統(tǒng)計Xm(i)與Xm(j)之間距離小于等于r的j(1≤j≤N －m，j≠i)的數(shù)目，記作 Bi。對于 1≤i≤N －m，定義

(4)增加維數(shù)到m+1，計算Xm+1(i)與Xm+1(j)小于等于r的距離的個數(shù)，記為Ai

定義:

(5)樣本熵定義為

從樣本熵定義可以看出，所求樣本熵與m，r的取值有關(guān)，因此，確定m，r兩個參數(shù)的選取對于樣本熵的計算很重要。根據(jù)文獻(xiàn)［9］的研究成果，當(dāng)m=1或2，r=0.1－0.25 std(std為原始數(shù)據(jù)的標(biāo)準(zhǔn)差)時所得到的樣本熵的統(tǒng)計特征比較合理，因此本研究取m=2，r=0.15 std。

2 流形學(xué)習(xí)與局部切空間排列

2.1 流形學(xué)習(xí)

流形學(xué)習(xí)法是一種基于微分幾何與拓?fù)鋵W(xué)的非線性高維數(shù)據(jù)處理方法，采用該方法能將高位數(shù)據(jù)在保持其主特征的基礎(chǔ)上降低數(shù)據(jù)維數(shù)。假定給定數(shù)據(jù)集X={Xi，i=1，…，n}，Xi∈RD為獨立同分布的，M 為嵌入在D維歐氏空間中的d維流形，定義嵌入映射f∶M?Rd→RD，其中D?d，流形學(xué)習(xí)的任務(wù)就是在給定的觀測數(shù)據(jù)X的前提下，獲取其低維流形表達(dá)式Y(jié)={Yi，i=1，…，n}∈RD，構(gòu)造出從高維空間到低維空間的非線性映射 f(·)［10］。

2.2 局部切空間排列

局部切空間排列(LTSA)算法通過逼近每一樣本點的切空間來構(gòu)建低維流形的局部幾何結(jié)構(gòu)，觀測數(shù)據(jù)點在局部切空間的投影獲取局部低維坐標(biāo)，然后通過仿射變換得到全局坐標(biāo)。給定一個樣本點集{x1，x2，…，xn}，xi∈Rm，LTSA 的具體算法如下［11－13］:

第一步，尋找每個樣本點的領(lǐng)域，不妨設(shè)xi={xi1，xi2，…，xik］為樣本點xi包括自身在內(nèi)的最近的K個鄰域點所構(gòu)成的矩陣。

第二步，局部線性投影。LTSA計算一個d維的仿射子空間來逼近Xi中的點，即:

其中 Θ =［θ1，…，θk］且 Q 的列數(shù)為 d，記=xilk為鄰域矩陣 xi的中心點，Qi∑vi為中心化鄰域矩陣Xi－=的奇異值分解，即 Qi，Vi分別為對應(yīng)于最大的d個奇異值的左右奇異向量所構(gòu)成的矩陣。這樣很容易求出式(12)的最優(yōu)解為:

從而可以得到局部坐標(biāo)系統(tǒng)為:

第三步，將所有這些有交疊的局部坐標(biāo)Θi=［θ(i)1，…)］排列起來以得到一個全局坐標(biāo)系統(tǒng)T=［τ1，…，τN］。

LTSA能較好的恢復(fù)出流形等距的低維空間子集，而且LTSA并不要求這個低維空間子集是凸的，也就是說，LTSA方法對于帶有空洞的流形，也具有很好的適用性。

3 基于小波包樣本熵與流形學(xué)習(xí)的特征提取模型

由于利用原始信號的樣本熵只能在一個尺度域進行分析，無法準(zhǔn)確區(qū)分機械運行狀況，利用小波包分析的多分辨率特性，將小波包分解與樣本熵結(jié)合起來在多個尺度對原始信號進行特征的提取，然后利用流形學(xué)習(xí)的非線性復(fù)雜信息處理能力進一步提取高維數(shù)據(jù)的低維特征，該模型不僅增大了故障診斷率，而且降低了特征數(shù)據(jù)的復(fù)雜度，從而增強了故障模式別的分類性能。該特征提取模型如圖2所示，具體步驟如下:

(1)信號分解:利用DB2小波對信號進行3層小波包多分辨率分解。

(2)信號重構(gòu):對經(jīng)過分解得到的8個頻帶的序列進行重構(gòu)，得到8個小波包重構(gòu)信號，每個重構(gòu)信號分別包含了原始信號各個不同頻段的信息。

(3)計算小波包樣本熵特征向量:由式(9)分別計算8個小波包樣本熵，構(gòu)成一個8維的特征向量:

(4)特征值的歸一化處理。

(5)流形學(xué)習(xí)進一步特征提取:利用局部切空間排列對初提取的小波包樣本熵特征向量進行低維流形特征的提取。

4 試驗與結(jié)果分析

4.1 試驗數(shù)據(jù)采集

本文的試驗研究是在軸承振動試驗臺進行完成的，試驗對象為JZQ－250型變速箱的滾動軸承，通過人為制造一些滾動軸承微小故障(故障直徑都約為0.2 mm)來模擬真實故障進行試驗和數(shù)據(jù)的采集。圖3為滾動軸承故障診斷試驗臺。

圖2 故障特征提取模型Fig.2 Fault feature extraction model

圖3 軸承振動試驗臺Fig.3 Bearing vibration test bench

在齒輪箱滾動軸承的試驗中，電機的轉(zhuǎn)速為1 750 r/min，采用了正常滾動軸承和具有內(nèi)圈故障、外圈故障和滾動體故障的滾動軸承，分別進行振動信號的采集。數(shù)據(jù)采集卡選用的是研華PCI－1711數(shù)據(jù)采集卡，5通道振動數(shù)據(jù)同步采集，采樣頻率為12 kHz。

4.2 特征提取分析

小波包特征熵特征提取方法基本思路是在某一分解尺度上將信號在不同的頻帶內(nèi)分解，分別提取不同頻帶內(nèi)的分解系數(shù)作為相應(yīng)的特征向量。圖4為試驗軸承內(nèi)圈故障信號的第3層小波包分解結(jié)果。每個頻帶的頻率范圍為750 Hz。

由圖4可知，出現(xiàn)故障的齒輪箱在旋轉(zhuǎn)過程中會產(chǎn)生周期性沖擊的激振力，同時激發(fā)不同的頻率成分，其振動信號為非平穩(wěn)的復(fù)雜調(diào)制信號，細(xì)節(jié)信號中基本都出現(xiàn)了明顯的有規(guī)律的沖擊信號。對小波包分解重構(gòu)的8個細(xì)節(jié)信號分別求出m=2時的樣本值，圖5是滾動軸承四種不同狀態(tài)信號的小波包樣本熵值的分布圖。

由圖5可以看出，滾動軸承不同狀態(tài)下小波包樣本熵值是不同的，不同頻段樣本熵值的大小反映了該頻段信號分布的確定性。表1為不同頻段小波包樣本熵的平均值，可以看出，正常滾動軸承的小波包熵值大于其他3種工作狀態(tài)的熵值，這是因為在正常狀態(tài)下，振動信號的樣本分布相對平均和不確定。當(dāng)出現(xiàn)內(nèi)圈或外圈或滾動體故障時，在相應(yīng)的頻帶內(nèi)會出現(xiàn)相應(yīng)的共振頻率，此時，在這個頻帶內(nèi)的樣本值將會增大，樣本分布的不確定性減少，從而使得小波包熵值減少。由于滾動體故障樣本集中更為嚴(yán)重，所以其小波包樣本值最小。

圖4 內(nèi)圈故障信號及其小波包分析Fig.4 Fault signal of inner race and its wavelet packet analysis

圖5 各頻帶小波包樣本熵值分布圖Fig.5 Wavelet packet sample entropy distribution map of each band

表1 各頻帶小波包樣本熵平均值Tab.1 Average value of each band wavelet packet sample entropy

因此，可以通過小波包對振動信號進行不同頻段的分解和重構(gòu)，得到包含不同的頻段信息的樣本熵，將之作為故障特征向量(此為特征提取第一階段)，再用流形學(xué)習(xí)進一步提取低維的故障特征(此為特征提取第二階段)。由于不同的流形學(xué)習(xí)法對不同數(shù)據(jù)的特征提取效果不同，作為比較，本文將用一種典型的非線性流形學(xué)習(xí)(LTSA)和一種典型的線性流形學(xué)習(xí)(PCA)對提取的小波包樣本熵進行降維。針對滾動軸承四種不同工作狀態(tài)信號，每個狀態(tài)信號取1 024個點作為一個樣本，采用db2小波包對信號進行3層次分解，分別提取第3層各節(jié)點的樣本熵作為特征向量，得到8維的列向量，再將8維的特征向量用流形學(xué)習(xí)法降到3維。圖6為LTSA和PCA降維的效果圖，鄰域參數(shù)K=7，目標(biāo)維數(shù)為3。

圖6 該模型的特征提取效果圖Fig.6 Features extraction effect map of this model

由圖6可以看出，對于提取的小波包樣本熵特征數(shù)據(jù)，LTSA比PCA的特征提取效果更好，從圖6(a)可以看出，滾動軸承故障的四種狀態(tài)特征已經(jīng)基本完全分開，證明了該模型具有很好的特征提取效果。為了量化分類效果，將兩個階段提取的特征用SVM進行分類。四種狀態(tài)信號先各提取10個樣本作為訓(xùn)練樣本，然后再各提取10個樣本作為測試樣本。將訓(xùn)練樣本輸入到SVM進行訓(xùn)練，SVM算法分類設(shè)置為多分類，核函數(shù)選擇徑向基核函數(shù)，表2給出了特征提取兩個階段基于SVM的故障狀態(tài)識別結(jié)果。

表2 支持向量機識別結(jié)果Tab.2 Support vector machine recognition results

第一階段特征SVM訓(xùn)練所花時間為6.3 s，第二階段中用LTSA提取的特征SVM訓(xùn)練所花時間為2.02 s，而用PCA提取的特征SVM訓(xùn)練所花時間為3.3 s。在第一階段，由于小波包分解頻帶中含有大量信息在特征空間中產(chǎn)生各種重疊，使得后面分類器(SVM)訓(xùn)練時間加長，并且降低了特征的分辨率，所以其故障識別率不是很好。而第二階段提取的特征的故障識別率達(dá)到100%，訓(xùn)練時間比第一階段短，對于相似故障，文獻(xiàn)［5］提出的特征提取方法其故障識別率只有95.33%，本文提出的特征提取方法識別效果明顯好于文獻(xiàn)［5］提出的方法。從而證明了該特征提取模型具有很好的故障模式分類性能。

5 結(jié)論

基于小波包、樣本熵和流形學(xué)習(xí)理論，提出了一種新的故障特征提取模型，首先選擇容錯性好的樣本熵作為特征參數(shù)，通過對振動信號進行基于db2小波的小波包分解，提取信號的小波包樣本熵作為特征向量，然后利用流形學(xué)習(xí)進一步提取低維特征。試驗結(jié)果表明，基于小波包樣本熵和流形學(xué)習(xí)的故障特征提取模型不僅提高了特征分辨率，而且簡化了特征數(shù)據(jù)，增強了故障模式識別的分類性。

［1］于德介，程軍圣，成瓊.基于復(fù)小波變換相位譜的齒輪故障診斷［J］.振動、測試與診斷，2004，24(4):275－277.YU De-jie，CHENG Jun-sheng，CHENG qiong.Gear fault diagnosis of complex wavelet transform based on phase spectrum［J］. JournalofVibration， Measurement ＆Diagnosis，2009，29(2):201 －204.

［2］鞠萍華，秦樹人，秦毅，等.多分辨EMD方法與頻域平均在齒輪早期故障診斷中的研究［J］.振動與沖擊，2009，28(5):97－101.JU Ping-hua，QIN shu-ren，QIN Yi，et al.Study of the average in the early fault diagnosis of gear multiresolution and frequency domain EMD method［J］.Journal of Vibration and Shock，2009，28(5):97 －101.

［3］Yan R Q，Gao R X.Approximate entropy as a diagnostic tool for machine health monitoring［J］.Mechancial Systems and Signal Processing，2007，21(2):824 －839.

［4］Pincus S M.Approximate entropy(ApEn)as a complexity measure［J］.Chaos，1995，5(1):110 －117.

［5］趙志宏，楊紹普.一種基于樣本熵的軸承故障診斷方法［J］.振動與沖擊，2012，31(6):136－140.ZHAO Zhi-hong，YANG Shao-pu.A bearing fault diagnosis method based on sample entropy［J］.Journal of Vibration and Shock，2012，31(6):136 －140.

［6］粟茂林，梁霖，王孫安，等.基于連續(xù)小波系數(shù)非線性流形學(xué)習(xí)的沖擊特征提取方法［J］.振動與沖擊，2012，31(1):106－111.LI Mao-Lin，LIANG Lin，WANG Sun-an，et al.Mechanical impact feature extraction method based on nonlinear manifold learning of continuous wavelet coefficients［J］.Journal of Vibration and Shock，2012，31(1):106 －111.

［7］葛爽，謝小鵬，黃博，等.基于能量耗損的故障診斷信號分析方法［J］.潤滑與密封，2012，37(8):7－10.GE Shuang，XIE Xiao-peng，HUANG Bo，et al.New algorithm for manifold-learning-based energy dissipation signal analysis［J］.Lubrication Engineering，2012，37(8):7 －10.

［8］Nikolaou N G，Antoniadis I A.Rolling element bearing fault diagnosis using wavelet packets［J］.NDT ＆ E International，2002，35(3):197－205.

［9］Pincus S M.Assessing serial irregularity and its implications for health［J］.Annals of the New York Academy of Sciences，2002，954:245－267.

［10］Silva V D，Tenenbaum J B.Global versus local methods in monlinear dimensionality reduction［J］.Advances in Neural Information Processing Systems，2002:705 －712.

［11］張熠卓，徐光華，梁霖，等.利用增量式非線性流形學(xué)習(xí)的狀態(tài)監(jiān)測方法［J］.西安交通大學(xué)學(xué)報.2011，45(1):64－68.ZHANG Yi-zhuo，XU Guang-hua，LIANG Lin，etal.Condition monitoring method for mechanical equipments based on incremental nonlinear manifold learning［J］.Journal of Xi＇an Jiao Tong University.2011，45(1):64 －68.

［12］ZHANG Zhen-yue，ZHA Hong-yuan.Principal manifolds and nonlineardimension reduction via localtangentspace alignment［J］.SIAM Journal on Scientific Computing，2005，26(1)313－338.

［13］梁霖，徐光華，栗茂林，等.沖擊故障特征提取的非線性流形學(xué)習(xí)方法［J］.西安交通大學(xué)學(xué)報，2009，43(11):95－99.LIANG Lin，XU Guang-hua，LI Mao-lin，et al.Nonlinear manifold learning method ofmechanlcalimpulse faults extraction［J］.Journal of Xi＇an Jiao Tong University，2009，43(11):95－99.