基于正則化核最大邊界投影維數(shù)約簡的滾動軸承故障診斷

趙孝禮, 趙榮珍, 孫業(yè)北, 何敬舉

(蘭州理工大學(xué) 機電工程學(xué)院， 蘭州 730050)

基于正則化核最大邊界投影維數(shù)約簡的滾動軸承故障診斷

趙孝禮, 趙榮珍, 孫業(yè)北, 何敬舉

(蘭州理工大學(xué) 機電工程學(xué)院， 蘭州 730050)

針對旋轉(zhuǎn)機械故障診斷中故障樣本獲取困難的現(xiàn)狀，提出一種基于正則化核最大邊界投影(Regularized Kernel Maximum Margin Projection, RKMMP)維數(shù)約簡的滾動軸承故障診斷方法。該方法首先利用RKMMP對小樣本、少標(biāo)記信息的混合故障樣本集進行訓(xùn)練降維，然后將降維后的低維敏感特征子集輸入到核極限學(xué)習(xí)機( Kernel Extreme Learning Machine, KLEM)分類器中進行故障識別。上述方法的特點是所提出的RKMMP能充分利用少量標(biāo)記樣本信息與大量無標(biāo)記樣本的故障信息，避免過學(xué)習(xí)的缺陷，同時通過添加正則化項克服小樣本問題。滾動軸承故障模擬實驗表明：該方法結(jié)合了RKMMP在維數(shù)約簡和KLEM在模式識別上的優(yōu)勢，在一定程度上能提升故障診斷的泛化能力與識別精度。該研究可為解決好故障診斷中樣本獲取困難的問題，提供理論參考依據(jù)。

故障診斷；正則化核最大邊界投影；核極限學(xué)習(xí)機分類器；維數(shù)約簡

滾動軸承作為旋轉(zhuǎn)機械中應(yīng)用廣泛且容易損壞的零部件，其運行狀態(tài)正常與否直接影響到系統(tǒng)的工作性能，因此對滾動軸承進行故障診斷有著重大的現(xiàn)實意義[1]。振動信號分析是軸承故障診斷主要手段之一，從非線性、非平穩(wěn)振動信號中提取出能綜合反映故障狀態(tài)的多域量化特征指標(biāo)的方法，已被廣泛應(yīng)用于軸承故障診斷中[2]。但多域特征的引入使得多元冗余信息滲入其中，這將導(dǎo)致特征信息之間的相關(guān)性增大、數(shù)據(jù)的分類能力被嚴重惡化的問題。因此，需用維數(shù)約簡方法對特征數(shù)據(jù)進行有效的二次特征提取，提取出低維敏感的主要特征矢量[3]。

故障特征數(shù)據(jù)的獲取是基于維數(shù)約簡的故障診斷方法的關(guān)鍵一環(huán)，故障數(shù)據(jù)常常具有“小樣本、高維數(shù)”的特點[4]。在實際問題中，無標(biāo)記的監(jiān)測數(shù)據(jù)容易獲取，而有標(biāo)記的監(jiān)測數(shù)據(jù)獲取困難。若只利用少量有標(biāo)記樣本則會導(dǎo)致學(xué)習(xí)系統(tǒng)泛化能力差；若只利用無標(biāo)記樣本，則會浪費標(biāo)記樣本中有用信息?；谟斜O(jiān)督與無監(jiān)督降維算法的不足，綜合利用少量標(biāo)記信息和大量無標(biāo)記信息，“半監(jiān)督降維算法”逐漸成為解決此類問題的關(guān)鍵[5]。

最大邊界投影(MMP)是He等[6]提出的半監(jiān)督降維算法。該算法的特點是：通過有區(qū)別的構(gòu)建不同類別標(biāo)簽的權(quán)重，最小化類內(nèi)緊湊圖的同時最大化類間分離圖；使得高維數(shù)據(jù)的低維流形嵌入后，同類數(shù)據(jù)點更加緊密，異類數(shù)據(jù)點更稀疏。但MMP本質(zhì)上是線性化的降維方法，處理非線性數(shù)據(jù)比較困難；此外，面對“小樣本問題”僅通過PCA預(yù)處理效果不明顯，且無監(jiān)督模式容易丟失大量有用信息[7]。針對這兩點不足，欲通過引入核方法，可將MMP的作用范圍引入到非線性領(lǐng)域；進一步地，通過向目標(biāo)函數(shù)添加正則化項克服小樣本問題[8]，提出正則化核最大邊界投影(RKMMP)算法，并將其應(yīng)用到“小樣本、少標(biāo)記信息”故障數(shù)據(jù)維數(shù)約簡的問題中。

故障診斷另一個核心問題是分類器的構(gòu)建。核極限學(xué)習(xí)機(KELM)是Huang等[9]根據(jù)極限學(xué)習(xí)機(ELM)[10]的隨機映射缺陷，提出的一種改進的新型網(wǎng)絡(luò)學(xué)習(xí)方法。與傳統(tǒng)分類器相比，KLEM具有學(xué)習(xí)速度快、泛化能力強、魯棒性好等優(yōu)點，在故障診斷等方面得到了廣泛應(yīng)用[11]。

鑒于上述分析，本研究擬對RKMMP與KLEM相結(jié)合的維數(shù)約簡滾動軸承故障診斷方法進行探討。欲為“小樣本、少標(biāo)記信息”的滾動軸承故障數(shù)據(jù)集的維數(shù)約簡與類別劃分，提供一種理論參考依據(jù)。

1 正則化核最大邊界投影(RKMMP)算法簡介

1.1 RKMMP算法的原理

MMP是LSDA[12]算法的半監(jiān)督擴展，半監(jiān)督學(xué)習(xí)的理論依據(jù)是根據(jù)先驗一致性假設(shè)：即同類標(biāo)簽的樣本具有相同的流形結(jié)構(gòu)，異類標(biāo)簽的樣本具有不同的流形結(jié)構(gòu)[13]。MMP既可以利用大量無標(biāo)記樣本，避免過擬合問題；又可利用少量標(biāo)記樣本刻畫數(shù)據(jù)的流形結(jié)構(gòu)。為解決MMP算法無法有效處理小樣本問題與非線性數(shù)據(jù)的缺陷，提出全新的RKMMP算法。

所提出的RKMMP算法基本思路為：首先，通過核函數(shù)將數(shù)據(jù)點映射到高維特征空間F中；然后，充分利用有標(biāo)簽與無標(biāo)簽信息，構(gòu)造有區(qū)別的權(quán)重矩陣來描述類內(nèi)與類間近鄰圖Gw、Gb，即最大化同類標(biāo)簽的數(shù)據(jù)點間的權(quán)重，最小化不同類別標(biāo)簽的數(shù)據(jù)點間的權(quán)重；通過向目標(biāo)函數(shù)中添加正則化項并優(yōu)化，使得映射后的同類樣本間的點距離更近、異類更遠。具體描述如下：

將樣本點xi的近鄰集合N(xi)分為:近鄰樣本集中與xi同類別的數(shù)據(jù)子集Nw(xi)和近鄰樣本集中與xi不同類別的數(shù)據(jù)子集Nb(xi)。定義Nw(xi)、Nb(xi)分別為：

(1)

Nb(xi)=N(xi)-Nw(xi)

(2)

式中:c(xi)表示xi的類別標(biāo)簽，k表示xi的近鄰個數(shù)。

(3)

(4)

式中，γ為衡量權(quán)值數(shù)，參考文獻[6]取值50。為使高維特征空間投影到低維特征空間后，類內(nèi)近鄰圖數(shù)據(jù)點更加緊湊，同時類間圖數(shù)據(jù)點更加離散，建立兩個優(yōu)化目標(biāo)函數(shù),即：

(5)

(6)

式中:Dw、Db分別為類內(nèi)、類間對角矩陣，Lb=Db-Wb。對式(5)施加一個約束條件aTXDwXTa=1，進一步地將式(5)簡化為：

(7)

綜合兩目標(biāo)函數(shù)式(6)、(7)，得到MMP的優(yōu)化目標(biāo)函數(shù)，即：

(8)

將式(8)轉(zhuǎn)換成求解廣義特征值問題見式(9)，即

X(βLb+(1-β)Ww)XTa=λXDwXTa

(9)

當(dāng)訓(xùn)練樣本數(shù)目小于樣本的特征維數(shù)時,MMP算法會出現(xiàn)病態(tài)奇異矩陣問題，也就是“小樣本問題”。這容易導(dǎo)致MMP無法計算出最佳投影矩陣或計算結(jié)果不穩(wěn)定。為使奇異矩陣變成正定矩陣，通過添加一個正則化常數(shù)δ與單位矩陣I乘積所得常數(shù)矩陣δI作為正則因子，即將式(9)轉(zhuǎn)變?yōu)镽MMP的廣義特征值問題：

X(βLb+(1-β)Ww)XTa=λ(XDwXT+δI)a

(10)

RKMMP則是通過引入式(11)所示的RBF核函數(shù)，將RMMP的應(yīng)用擴展到非線性領(lǐng)域。

Ki,j=<φ(xi),φ(xj)>=exp(-‖xi-xj‖/t)

(11)

式中:φ表示非線性函數(shù);t表示核參數(shù)。在Hilbert空間F中式(10)變形為

φ(X)[βLb+(1-βWw](X)Taφ=
λ[φ(X)Dwφ(X)T+δI]aφ

(12)

式中，aφ是由F空間所有向量線性組合而成的解向量，則存在向量aφ=[a1,a2,..,an]滿足：

(13)

結(jié)合式(12)、(13)，對式(12)兩邊做內(nèi)積，可得KRMMP優(yōu)化目標(biāo)式(14)：

K[βLb+(1-β)Ww]Ka=λ(KDwK+δI)a

(14)

1.2 RKMMP算法的實現(xiàn)流程

RKMMP算法具體步驟如下：

輸入：高維空間數(shù)據(jù)集X={(x1,c1), (x2,c2),..(xl,cl),xl+1,...xl+u},xi∈Rd,

輸出：低維特征向量Y，投影映射矩陣A。

(1)將原始數(shù)據(jù)集X進行數(shù)據(jù)歸一化，得到Xo。

(2)將Xo通過核函數(shù)K(·)映射到高維特征空間F中，得到φ(Xo)。

(3)利用有標(biāo)記樣本和未標(biāo)記樣本，根據(jù)式(5)、(6)通過k近鄰法對數(shù)據(jù)集構(gòu)造類內(nèi)、類間近鄰圖。

(4)根據(jù)權(quán)值矩陣，求解式(14)的廣義特征值和特征向量。得前m個最小非零特征值對應(yīng)的特征向量，構(gòu)成投影向量aφ={a1，a2,…，am}，并將φ(Xo)通過a進行降維投影，得到映射后的低維特征集Y=aTφ(Xo)。

2 核極限學(xué)習(xí)機(KELM)算法簡介

2.1 KELM分類器

KELM作為ELM算法的擴展，是由Huang等根據(jù)SVM與ELM的優(yōu)化角度，將核函數(shù)引入到極限學(xué)習(xí)機中，克服了隨機初始值導(dǎo)致的不穩(wěn)定問題。

(15)

式中：β表示輸出權(quán)重矩陣；h(xp)表示輸入xp經(jīng)隱含層映射后的輸出矩陣；若ELM的隱含層輸出函數(shù)G(w,x,b)已知，則計算h(xp)形式如下

h(xp)=[G(w1·xp+b1),…,G(wL·xp+bL)]

(16)

式中:w為輸入權(quán)值，b為偏值；w、b在訓(xùn)練中隨機賦值。

由Karush-Kuhn-Tucker (KKT)定理可知，ELM的訓(xùn)練學(xué)習(xí)等效為求解如式(17)的對偶優(yōu)化問題：

(17)

式中：C為懲罰系數(shù)，用以權(quán)衡結(jié)構(gòu)風(fēng)險和經(jīng)驗風(fēng)險之間的比例；τp表示理論輸出tp與實際輸出f(xp)的誤差；αp為拉格朗日算子，求解式(17)可得：

(18)

式中：H為隱含層輸出矩陣，T=[t1,t2,...,tN]表示輸入樣本集的目標(biāo)向量。HHT在求解廣義逆過程中，因數(shù)據(jù)存在復(fù)共線性導(dǎo)致HHT非奇異，容易引起不穩(wěn)定現(xiàn)象。故將核函數(shù)思想引入ELM中，根據(jù)Mercer條件定義核矩陣：

Ψ=HHT∶Ψi,j=h(xi)·(xj)=K(xi,xj)

(19)

通過核矩陣K代替HHT，由式(18)、(19)代入式(15)中得KLEM輸出模型：

f(xp)=[K(xp,x1),…,K(xp,x1)]ζ

(20)

(21)

式中，ζ即為KELM網(wǎng)絡(luò)的輸出權(quán)值。

由以上過程可知，通過向KELM添加偏值常量1/C，可增加系統(tǒng)的穩(wěn)定性與泛化能力；同時KLEM無需設(shè)置眾多參數(shù)，只需選擇合適的核函數(shù)。相對于ELM，KELM解決分類與回歸問題效果更好，具有更強大的非線性逼近能力。

2.2 KELM分類器故障辨識流程

KLEM算法故障辨識流程圖見圖1。它的具體實現(xiàn)步驟如下：

圖1 KELM分類器的故障辨識流程圖Fig.1 The flow chart of fault identification of the KELM classifier

(1)將輸入樣本分為訓(xùn)練樣本集與測試樣本集，確定核參數(shù)及懲罰因子。

(2)求解出輸出權(quán)值β，寄存KELM參數(shù)，得訓(xùn)練學(xué)習(xí)后的KELM分類器。

(3)將測試樣本與訓(xùn)練樣本輸入KLEM分類器中進行故障診斷。

3 設(shè)計的RKMMP維數(shù)約簡的故障診斷方法

3.1 混合故障數(shù)據(jù)集、RKMMP與KELM相結(jié)合的故障診斷方法的探討

為解決故障診斷過程中，訓(xùn)練樣本與有標(biāo)記信息樣本較少導(dǎo)致故障診斷困難的現(xiàn)狀。擬以軸承故障數(shù)據(jù)集為實驗對象，提出RKMMP維數(shù)約簡與KELM分類器相結(jié)合的故障診斷方法。該方法所提出的RKMMP半監(jiān)督降維方法，一方面能綜合利用少量有標(biāo)記故障樣本與大量無標(biāo)記故障樣本，提取出最佳分類的低維特征子集，另一方面正則化與核方法的引入能有效克服小樣本問題，增強其處理非線性問題的能力。故障診斷的實質(zhì)是模式識別，方法中的KELM分類器能快速、穩(wěn)定、精確地實現(xiàn)故障類別的量化表征。

所設(shè)計的故障診斷方法主要思路為：首先，以滾動軸承的振動信號為對象，從時域、頻域、時頻域等多方面構(gòu)建能全面反映故障信息的高維、小樣本、少標(biāo)記信息的混合故障數(shù)據(jù)集，然后將數(shù)據(jù)集輸入到RKMMP中降維處理，得到低維敏感的特征子集；最后，將低維特征子集輸入到KELM分類器中進行訓(xùn)練學(xué)習(xí)與故障辨識，得出故障診斷結(jié)果。

3.2 故障診斷方法實現(xiàn)的流程與步驟

設(shè)計的基于RKMMP維數(shù)約簡的故障診斷方法流程，見圖2。整個故障診斷方法由如下幾個步驟實現(xiàn)：

圖2 故障診斷方法流程圖Fig.2 The flow chart of fault diagnosis method

步驟1 對機械振動信號進行采集，從時域、頻域和時頻域方面做特征提取，得到原始數(shù)據(jù)特征集X。

步驟2 根據(jù)故障數(shù)據(jù)的具體特性確定核參數(shù)t，并將特征集經(jīng)高斯核函數(shù)K(·)映射到高維特征空間F，在此空間進行RMMP降維，得到低維特征集Y。

步驟3 將低維特征集Y輸入KELM進行模式識別，得到測試樣本的故障類型。

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

4.1 實驗數(shù)據(jù)的來源

實驗中所采用的實驗數(shù)據(jù)來源于美國凱斯西儲大學(xué)(CWRU)電氣工程實驗室的滾動軸承實驗數(shù)據(jù)[14]。實驗軸承的型號為6205-RS JEM SKF深溝球軸承。軸承的具體參數(shù)為：內(nèi)徑25 mm，外徑52 mm，厚度15 mm，滾動體直徑8.18 mm。實驗所選取的數(shù)據(jù)是在如下條件下測得：電機負載為2 hp，在采樣頻率為48 kHZ，轉(zhuǎn)速為1 750 r/min情況下，通過驅(qū)動端軸承座上的加速度傳感器拾取振動信號。實驗中軸承的損傷是用電火花加工模擬，選取故障尺寸為直徑0.53 mm、深度為0.28 mm(故障等級為嚴重)的三種典型故障類型：外圈故障、內(nèi)圈故障、滾動體故障。以2 048個采樣點為一組樣本，分別截取軸承外圈、軸承內(nèi)圈、滾動體故障以及正常四種狀態(tài)下的振動信號各50組，以其中每種狀態(tài)的20組作為訓(xùn)練樣本(其中前5組為有標(biāo)記故障樣本，后15組為無標(biāo)記樣本)，余下的30組作為測試樣本。為建立高維、小樣本混合故障特征集，需從時域、頻域、時頻域多個角度構(gòu)建如表1的特征集合，共得到108個特征參數(shù)(得到4×20×108維的小樣本故障數(shù)據(jù)特征集)。其中，因軸承信號包含四種故障類型，故EMD最少能分解出4個固有模態(tài)函數(shù)(IMF)和一個余項信號，為充分把握原數(shù)據(jù)的特征信息，提取前5層IMF分量及復(fù)雜度作為能量特征，即35～44號特征；同時采用db4小波包函數(shù)對振動信號進行6層正交小波包分解，分解出的各頻帶的能量特征，即45～108號特征。

表1 特征參數(shù)

4.2 應(yīng)用情況與分析

從實驗數(shù)據(jù)中構(gòu)建高維混合故障特征集，輸入到RKMMP中進行維數(shù)約簡，再通過KELM分類器進行故障識別。根據(jù)文獻[6]，設(shè)置參數(shù)近鄰值k=5、調(diào)節(jié)因子設(shè)置為0.5、正則化因子設(shè)置為1；為充分涵蓋原始數(shù)據(jù)集的本質(zhì)信息，目標(biāo)維數(shù)統(tǒng)一降維到3維(目標(biāo)維數(shù)=故障類別數(shù)-1)；為驗證所提方法的有效性，分別采用PCA、KPCA、KLSDA、MMP、KMMP算法對原始故障數(shù)據(jù)集進行維數(shù)約簡(其中，KLSDA、MMP、KMMP經(jīng)過PCA預(yù)處理)，作對比試驗。為方便后續(xù)工作，分別記六種算法{PCA, KPCA, KLSDA, MMP, KMMP, RKMMP}={R1, R2, R3, R4, R5, R6}。對核函數(shù)參數(shù)選取，采用五折交叉法選取最優(yōu)高斯核參數(shù)，得到KPCA、KLSDA、KMMP、RKMMP的核參數(shù)分別為t=88、75、20、60。在降維過程中，半監(jiān)督算法只利用有標(biāo)記故障樣本，無監(jiān)督算法利用所有樣本，但不利用故障標(biāo)簽信息；有監(jiān)督算法只利用有標(biāo)記樣本進行降維學(xué)習(xí)。為評估方法的有效性，實驗從以下三個角度進行驗證：

4.2.1可視化降維與故障辨識結(jié)果對比

將故障數(shù)據(jù)集輸入到降維算法中，得六種方法的降維后測試樣本三維嵌入圖，見圖3。

(a)PCA (b)KPCA (c)KLSDA

(d)MMP (e)KMMP (f)RKMMP圖3 測試樣本基于不同方法的降維結(jié)果Fig.3 The test sample results based on different methods of dimensionality reduction

從圖3可以看出：RKMMP降維后的聚類與分類效果最好，其降維后的各故障類內(nèi)距離小、類間間距明顯；PCA的降維后測試樣本聚類與分類效果最差；其他四種算法聚類與分類效果介于二者之間。除RKMMP外，其他五種降維方法的故障間都存在著一定程度的交叉混疊現(xiàn)象，其中外圈與滾子故障降維效果六種算法都不理想；正常狀態(tài)的各類算法降維結(jié)果都較為成功。

同時，將降維后的低維特征子集輸入ELKM分類器中進行故障模式識別。KLEM核函數(shù)取RBF型，交叉優(yōu)化參數(shù)得C=100，t=40。訓(xùn)練得到六種算法降維后各類故障的識別準確率，見圖4所示(其中，外圈、軸內(nèi)圈、滾動體以及正常四類狀態(tài)分別記為F1、F2、F3、F4，而F5為四類故障的總體平均識別率)。圖4中，六種降維方法的總體平均識別率逐漸上升，平均識別率RKMMP最高，PCA最低；針對單個故障，滾子故障的識別率六類算法都未能準確劃分出；正常狀態(tài)識別率最好，六種算法都能準確識別出來。

圖4 六類算法的故障診斷結(jié)果Fig.4 Fault diagnosis result of six algorithms

由圖3、圖4得如下結(jié)論：

1) PCA的降維效果及識別正確率最差，這由于PCA是基于二階統(tǒng)計信息，只能保持數(shù)據(jù)的最大方差方向投影，忽略了原始特征空間的非線性特征結(jié)構(gòu)，故效果較差；KPCA降維效果及識別率一般，這是因為KPCA雖然具有較強的非線性提取能力，但其保持局部信息的能力較低；PCA、KPCA無監(jiān)督方式導(dǎo)致其無法在類別信息的指導(dǎo)下進行維約簡，易丟失大量利于分類的信息，整體識別率較低。

2) KLSDA、MMP的降維效果及識別正確率一般。這是由于KLSDA雖然能挖掘出局部流形結(jié)構(gòu)信息，但過度依賴于少量有標(biāo)記的故障樣本，導(dǎo)致過學(xué)習(xí)問題，無法取得理想的故障診斷結(jié)果；MMP雖然能夠綜合利用有類別標(biāo)簽與無類別標(biāo)簽信息，但MMP本質(zhì)是線性降維方法，故降維與辨識結(jié)果未能取得最優(yōu)。

3) RKMMP降維與故障識別結(jié)果比KMMP等其它五種方法都好的原因是，KMMP無法克服小樣本問題，PCA處理后容易丟失大量有用信息，導(dǎo)致無法取得最優(yōu)。RKMMP既能利用無標(biāo)記故障樣本中的故障信息，又能克服KLSDA等有監(jiān)督算法出現(xiàn)的過學(xué)習(xí)問題，核方法及正則化的改進提高了其應(yīng)用范圍及泛化能力，故取得最優(yōu)的故障診斷結(jié)果，圖3的降維情況與圖4的結(jié)果證明了基于RKMMP維數(shù)約簡的故障診斷方法的有效性。

4.2.2 RKMMP解決小樣本與少標(biāo)記信息的性能驗證

為測試RKMMP算法對小樣本問題處理能力，通過改變各類故障的訓(xùn)練樣本與測試樣本的比例，設(shè)置每組10/40、15/35、20/30、25/25、30/20、40/10(其中訓(xùn)練樣本小于25時為小樣本情況，有標(biāo)記故障樣本為前5個不變)。并得六類降維算法低維測試樣本輸入KELM分類器所得的平均識別率，見圖5。

圖5 不同訓(xùn)練樣本數(shù)對應(yīng)的平均識別準確率Fig.5 The average recognition accuracy of different training sample

從圖5中可以看出：①總體上，除KLSDA只利用有標(biāo)記故障樣進行訓(xùn)練不變外，其他五種降維方法的識別準確率都隨著訓(xùn)練樣本數(shù)的增加而增加，這是因為訓(xùn)練樣本越多，就越能表征數(shù)據(jù)的整體情況，故障識別率得到一定提升;②PCA、KPCA在少訓(xùn)練樣本情況下，準確率下降明顯；MMP、KMMP雖然在一定程度能夠克服小樣本問題，但效果不明顯，這是因為PCA預(yù)處理作為無監(jiān)督方式，容易丟失大量分類信息，不能有效處理“小樣本問題”，導(dǎo)致故障診斷效果下降;③其中，RKMMP的穩(wěn)定性能最好，即使在小樣本情況下也能取得較高的識別率，這是因為正則化的改進對小樣本問題處理效果明顯，從而說明RKMMP能夠有效處理小樣本問題，推廣適應(yīng)能力較強。

為進一步地驗證方法中RKMMP綜合利用有標(biāo)記與無標(biāo)記樣本的能力，即測試算法中標(biāo)記樣本與無標(biāo)記樣本比例對降維效果的影響。設(shè)置如下實驗選取訓(xùn)練樣本與測試樣本比例為30/20不變(去除小樣本因素影響)，分別依次增加各類故障有標(biāo)記樣本數(shù)目，設(shè)置各狀態(tài)有標(biāo)記訓(xùn)練樣本分別為2、5、10、15、20、25、30個。經(jīng)RKMMP等算法訓(xùn)練降維后，并得低維測試樣本輸入KELM分類器所得的平均識別率，見圖6。

從圖6中可以看出：

(1)隨著訓(xùn)練類別標(biāo)簽樣本數(shù)的增加，KLSDA、MMP、KMMP、KRMMP的分類識別率隨之增加，這是因為隨著故障標(biāo)記的增加，先驗信息越多，過擬合現(xiàn)象隨之減少，故障識別率增加。而PCA、KPCA無監(jiān)督算法無法利用有標(biāo)記故障樣本，故無影響。

(2)其中，有監(jiān)督學(xué)習(xí)算法KLSDA受故障標(biāo)記影響較大，識別率變化比較明顯，說明有監(jiān)督算法嚴重依賴故障標(biāo)記信息；而半監(jiān)督算法對少量類別標(biāo)記的訓(xùn)練樣本也有著良好的故障分類性能，說明半監(jiān)督算法能夠同時利用樣本的有標(biāo)記與無標(biāo)記故障信息。

(3)RKMMP的穩(wěn)定性能最好，識別準確率一直高于其他算法；尤其是在標(biāo)記樣本極其少量的情況下，RKMMP、KMMP、MMP等半監(jiān)督算法就能獲得較高的故障診斷精度，進而說明半監(jiān)督學(xué)習(xí)算法的在處理少量有標(biāo)記故障樣本時的優(yōu)越性。

圖6 有故障標(biāo)記數(shù)目對故障診斷的影響Fig.6 The influence of fault labeled number for fault diagnosis

4.2.3 方法中KELM分類器的性能驗證

為進一步驗證分類器KELM優(yōu)異的魯棒性與穩(wěn)定性，訓(xùn)練與測試樣本之比為20/30等參數(shù)不變的情況下。參考文獻[15]向原始混合故障數(shù)據(jù)集中添加系數(shù)a=0.1、0.2、0.4、0.8的隨機干擾噪聲，將經(jīng)RKMMP降維后的結(jié)果輸入到KELM、BP、SVM、ELM分類器進行的故障辨識結(jié)果作對比，得到故障識別平均準確率，如表2所示。

由表2中可以看出：添加隨機噪聲后，ELM的故障診斷結(jié)果波動較大，而KELM變化?。徽f明KELM的識別率的穩(wěn)定性要優(yōu)于ELM，原因是ELM初始值的隨機性導(dǎo)致其訓(xùn)練不穩(wěn)定且泛化能力差，KELM在核映射條件下穩(wěn)定性能更好；由此證明了經(jīng)RKMMP降維后的KELM分類器具有較好的魯棒性能，穩(wěn)定性好、泛化能力強的特點。

表2 不同分類器的抗干擾能力對比

5 結(jié) 論

為解決故障診斷過程中故障樣本獲取困難的現(xiàn)狀，提出基于改進的半監(jiān)督降維算法RKMMP維數(shù)約簡的滾動軸承故障診斷方法。首先，新提出的RKMMP算法具有良好的維數(shù)約簡能力，能充分利用有標(biāo)簽與無標(biāo)簽信息，克服小樣本問題，提取出低維敏感、具有鑒別能力的特征子集；然后，將低維特征子集輸入到KELM分類器中，KELM分類器能更精確、穩(wěn)定地實現(xiàn)故障類別的劃分，適應(yīng)多類樣本問題。實驗結(jié)果表明，所提方法能夠?qū)ι儆?xùn)練樣本、少標(biāo)記信息的軸承故障數(shù)據(jù)集進行有效地維數(shù)約簡與故障辨識，為滾動軸承故障診斷提供了一種新的解決思路。該研究后續(xù)工作可以從以下幾個方面進行：

(1)只對軸承的嚴重故障進行診斷，尚未涉及對早期、輕度等故障類型進行診斷。

(2)KELM分類器的核參數(shù)選取只采用交叉驗證獲得，需要進一步的優(yōu)化。

[1] 張文義，于德介，陳向民. 基于信號共振稀疏分解與能量算子解調(diào)的軸承故障診斷方法[J]. 中國電機工程學(xué)報，2013, 33(20):111-118.

ZHANG Wenyi, YU Dejie, CHEN Xiangmin. Fault diagnosis of rolling bearings based on resonance-based sparse signal decomposition and energy operator demodulating[J]. Proceedings of the CSEE, 2013, 33( 20): 111-118.

[2] 黃宏臣，韓振南，張倩倩，等. 基于拉普拉斯特征映射的滾動軸承故障識別[J]. 振動與沖擊，2015,34(5): 128-134.

HUANG Hongchen, HAN Zhennan, ZHANG Qianqian，et al. Method of fault diagnosis for rolling bearings based on Laplacian eigenmap[J]. Journal of Vibration and Shock, 2015,34(5): 128-134.

[3] 李鋒，湯寶平，陳法法. 基于線性局部切空間排列維數(shù)化簡的故障診斷[J]. 振動與沖擊，2012, 31(13): 36-40.

LI Feng, TANG Baoping, CHEN Fafa. Fault diagnosis model based on dimension reduction using linear local tangent space alignment[J]. Journal of Vibration and Shock, 2012, 31(13): 36-40.

[4] 韋佳，楊創(chuàng)新，馬千里，等. 基于局部重構(gòu)與全局保持的半監(jiān)督判別分析方法[J]. 華南理工大學(xué)學(xué)報(自然科學(xué)版)，2010, 38(7):45-49.

WEI Jia, YANG Chuangxin, MA Qianli, et al. Semi-supervised discriminant analysis method based on local reconstruction and global preserving[J]. Journal of South China University of Technology(Nature Science Edition), 2010, 38(7):44-49.

[5] ZHU X. Semi-supervised learning literature survey[J]. Computer Science, 2008, 37(1):63-77.

[6] HE X, CAI D, HAN J. Learning a maximum margin subspace for image retrieval[J]. IEEE Transactions on Knowledge & Data Engineering, 2008, 20(2):189-201.

[7] JIANG L, SHI T, XUAN J. An intelligent fault diagnosis method of rolling bearings based on regularized kernel Marginal Fisher analysis[C]// Journal of Physics: Conference Series, 2012:597-604.

[8] PIMA I, ALADJEM M. Regularized discriminant analysis for face recognition[J]. Pattern Recognition, 2004, 37(9):1945-1948.

[9] HUANG G B, ZHOU H, DING X, et al. Extreme learning machine for regression and multiclass classification[J]. IEEE Transactions on Systems Man & Cybernetics, Part B: Cybernetics A Publication of the IEEE Systems Man & Cybernetics Society, 2012, 42(42):513-529.

[10] HUANG G B, ZHU Q Y, SIEW C K. Extreme learning machine: Theory and applications[J]. Neurocomputing, 2006, 70(1/2/3):489-501.

[11] 陳紹煒，柳光峰，冶帥. 基于核極限學(xué)習(xí)機的模擬電路故障診斷研究[J]. 西北工業(yè)大學(xué)學(xué)報，2015, 33(2):290-294.

CHEN Shaowei, LIU Guangfeng, YE Shuai. A method of fault diagnosis for analog circuit based on KELM[J]. Journal of Northwestern Polytechnical University, 2015,33(2):290-294.

[12] CAI D, HE X, ZHOU K, et al. Locality sensitive discriminant analysis[C]. IJCAI, 2007:708-713.

[13] 楊望燦，張培林，吳定海，等. 基于改進半監(jiān)督局部保持投影算法的故障診斷[J]. 中南大學(xué)學(xué)報(自然科學(xué)版)，2015, 47(6):2059-2064.

YANG Wangcan, ZHANG Peilin, WU Dinghai, et al. Fault diagnosis based on improved semi-supervised locality preserving projections[J]. Journal of Central South University(Science and Technology), 2015, 47(6): 2059-2064.

[14] LOPARO K A. Bearings vibration data set, Case Western Reserve University [DB/OL]. http://csegroups.case.edu/bearingdata-center/home.

[15] CHEN F, TANG B, CHEN R. A novel fault diagnosis model for gearbox based on wavelet support vector machine with immune genetic algorithm[J]. Measurement, 2013, 46(1):220-232.

Fault diagnosis of rolling bearings based on the dimension reduction using the regularized kernel maximum margin projection

ZHAO Xiaoli, ZHAO Rongzhen, SUN Yebei, HE Jingju

(School of Mechanical and Electronical Engineering, Lanzhou University of Technology, Lanzhou 730050, China)

In view of the difficulty of fault samples’ acquisition in the fault diagnosis of rotating machinery, a novel fault diagnosis method for rolling bearings based on the dimension reduction using the regularized kernel maximum margin projection (RKMMP) was proposed. In the method, the RKMMP was used to reduce the dimension of the mixed fault data set of small samples and less labeled information. Then the low-dimensional sensitive feature subset was input into a kernel extreme learning machine (KLEM) classifier for training and fault identification. The characteristics of the method is that the proposed RKMMP can make full use of the labeled information of small samples and the fault information of numerous unlabeled samples, and avoid the problem of over-fitting. At the same time, a regularization term was added to overcome the disadvantage of small samples. The simulations on rolling bearing fault diagnosis show that the method combines the advantages of RKMMP in dimension reduction and of KLEM in pattern recognition. To a certain extent, it can improve the generalization ability of fault diagnosis and the accuracy of recognition. The study solves the problem of samples’ acquisition and provides a theoretical support to the rolling bearing fault diagnosis.

dynamic response ; multiaxial random excitation; vibration fatigue; modal shape

國家自然科學(xué)基金(51675253)；高等學(xué)校博士學(xué)科點專項科研基金(20136201110004)

2016-04-20 修改稿收到日期： 2016-06-12

趙孝禮 男，碩士生，1991年生

趙榮珍 女，博士，教授，博士生導(dǎo)師，1960年生

TP18; TH165

A

10.13465/j.cnki.jvs.2017.14.016