第一作者王波男，博士生，1982年生

通信作者劉樹林男，教授，博士生導(dǎo)師，1963年生

基于量子遺傳算法優(yōu)化RVM的滾動(dòng)軸承智能故障診斷

王波1,2,劉樹林1,蔣超1,張宏利1

(1.上海大學(xué)機(jī)電工程與自動(dòng)化學(xué)院,上海200072； 2.滁州學(xué)院機(jī)械與汽車工程學(xué)院,安徽滁州239000)

摘要：提出了基于量子遺傳算法(QGA)優(yōu)化相關(guān)向量機(jī)(RVM)核函數(shù)參數(shù)的方法，通過仿真比較了量子遺傳算法與其它方法在核函數(shù)參數(shù)優(yōu)化方面的性能，結(jié)果表明基于量子遺傳算法優(yōu)化出的算法性能優(yōu)于其它方法的優(yōu)化性能。將基于量子遺傳算法優(yōu)化的相關(guān)向量機(jī)(QGA-RVM)應(yīng)用于滾動(dòng)軸承的故障診斷；采用總體平均經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)將滾動(dòng)軸承故障信號(hào)自適應(yīng)地分解成多個(gè)內(nèi)稟模態(tài)函數(shù)(IMF)，將IMF能量作為故障特征輸入到QGA-RVM進(jìn)行最終的故障診斷。結(jié)果表明，該方法能夠快速準(zhǔn)確地診斷出滾動(dòng)軸承故障，驗(yàn)證了該方法的有效性和穩(wěn)定性；此外，通過與支持向量機(jī)(SVM)的對(duì)比分析，顯示了RVM在智能故障診斷應(yīng)用中的優(yōu)越性。

關(guān)鍵詞：量子遺傳算法;故障診斷;相關(guān)向量機(jī);EEMD

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目(51175316);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(20103108110006)資助項(xiàng)目;滁州學(xué)院規(guī)劃研究項(xiàng)目(2014GH20)資助

收稿日期：2014-05-05修改稿收到日期：2014-10-11

中圖分類號(hào):TH212；TH213.3

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2015.17.034

Abstract:A novel method to optimize relevance vector machine (RVM)’s kernel function parameters based on the quantum genetic algorithm (QGA) was proposed. It was compared with other optimization algorithms with simulations. The results showed that the optimization method based on QGA is superior to other optimization methods. The model of RVM optimized with QGA (QGA-RVM) was applied in fault diagnosis of rolling bearings. Fault signals were decomposed adaptively into some intrinsic mode functions (IMFs) with the ensemble empirical mode decomposition (EEMD). The IMF energy as fault features was inputted into QGA-RVM for final fault diagnosis. Experimental results showed that the proposed method can diagnose rolling bearings’ faults rapidly and accurately, its validity and stability are verified; moreover, the superiority of RVM in intelligent fault diagnosis is revealed through the comparative analysis between QGA-RVM and SVM.

Rolling bearings’ intelligent fault diagnosis based on RVM optimized with Quantum genetic algorithm

WANGBo1,2,LIUShu-lin1,JIANGChao1,ZHANGHong-li1(1. School of Mechatronic Engineering and Automation,Shanghai University, Shanghai 200072, China;2. School of Mechanical and Automotive Engineering, Chuzhou Institute, Chuzhou 239000, China)

Key words:quantum genetic algorithm; fault diagnosis; relevance vector machine; ensemble empirical mode decomposition(EEMD)

在工程實(shí)際中，獲取機(jī)械設(shè)備故障樣本是一個(gè)難點(diǎn)。由于缺乏豐富的故障樣本，無法使用大量故障樣本訓(xùn)練診斷模型，在機(jī)械故障智能診斷領(lǐng)域涉及的通常是如何處理小樣本數(shù)據(jù)問題。支持向量機(jī)(Support Vector Machine, SVM) 具有處理小樣本問題的優(yōu)越性、良好的泛化能力等優(yōu)點(diǎn)，作為智能識(shí)別器在機(jī)械故障診斷領(lǐng)域得到了廣泛的關(guān)注和研究，同時(shí)也取得了良好的故障診斷應(yīng)用效果。然而，SVM也存在著一些缺點(diǎn)：①雖然所需的支持向量遠(yuǎn)小于訓(xùn)練樣本數(shù)，但會(huì)隨著訓(xùn)練樣本集的增大線性增加，造成過擬合和計(jì)算時(shí)間的增加；②核函數(shù)必須滿足Mercer 條件，即核函數(shù)必須是正定連續(xù)對(duì)稱函數(shù)；③無法獲得概率式的預(yù)測，無法獲取預(yù)測中的不確定性，同時(shí)也不符合實(shí)際應(yīng)用需求；④需要估計(jì)折中系數(shù)C(回歸時(shí)需要估計(jì)不敏感參數(shù)ε)，通常采用交叉驗(yàn)證的方法確定參數(shù)的值，這樣會(huì)增加計(jì)算時(shí)間，計(jì)算復(fù)雜度也會(huì)隨著訓(xùn)練樣本數(shù)的增加而增加。

Tipping[1-2]提出了一種基于貝葉斯理論的統(tǒng)計(jì)學(xué)習(xí)算法: 相關(guān)向量機(jī)(Relevance Vector Machine, RVM)，其相對(duì)于SVM有以下一些優(yōu)點(diǎn)：①可以輸出預(yù)測值的后驗(yàn)概率分布；②比SVM具有更加稀疏的決策函數(shù)，即需要相關(guān)向量比支持向量更少便可獲得與SVM類似的分類準(zhǔn)確率;③核函數(shù)不需要滿足Mercer 條件，可以構(gòu)造更多的核函數(shù)。近年來，在機(jī)械故障診斷領(lǐng)域，也有一些RVM的應(yīng)用與研究成果涌現(xiàn)[3-6]。然而，還未得到相關(guān)研究者的足夠重視。理論上RVM對(duì)于處理小樣本問題更具有優(yōu)勢，更加適合處理工程實(shí)際問題。同時(shí)，RVM相對(duì)簡單的決策函數(shù)所需較少的計(jì)算時(shí)間更加適合在線診斷。

基于SVM故障診斷模型參數(shù)(折中系數(shù)C與核函數(shù)參數(shù))的設(shè)定，很大程度影響著最終的故障診斷準(zhǔn)確率。RVM中沒有涉及折中參數(shù)C，僅需選擇合適的核參數(shù)值，緩解了參數(shù)對(duì)其診斷結(jié)果的影響。然而，RVM的分類結(jié)果對(duì)核參數(shù)依然十分敏感。傳統(tǒng)的方法常采用網(wǎng)格搜索和交叉驗(yàn)證方法相結(jié)合以選取最佳的核參數(shù)，這種枚舉方法消耗了大量的時(shí)間和資源。

為此，提出了一種基于RVM的故障診斷模型，利用量子遺傳算法(QGA)優(yōu)化RVM，自適應(yīng)地選擇最優(yōu)的核函數(shù)參數(shù)。同時(shí)，分別采用遺傳算法(GA)和粒子群算法(PSO)優(yōu)化RVM的核參數(shù)，比較分析了不同算法之間的優(yōu)劣，顯示了QGA在核參數(shù)優(yōu)化方面的優(yōu)勢。提取滾動(dòng)軸承振動(dòng)信號(hào)的EEMD能量作為故障特征，通過優(yōu)化后的RVM診斷模型進(jìn)行故障識(shí)別。結(jié)果表明了所提出方法的可行性和有效性。

1相關(guān)向量機(jī)分類理論模型

(1)

其中：ω=(w1,w2,…,wN)T為權(quán)重的系數(shù)向量，Kn(x)為選用的核函數(shù)。RVM希望預(yù)測出輸入樣本xi屬于某個(gè)類別的后驗(yàn)概率，因此可通過Sigmoid函數(shù)

(2)

將y(x)映射到(0,1)區(qū)間進(jìn)行分類[1]。

假設(shè)每個(gè)樣本獨(dú)立同分布，并且每次的預(yù)測為獨(dú)立事件，假設(shè)p(t|x)采用Bernoulli分布，可得預(yù)測結(jié)果t的后驗(yàn)概率的似然函數(shù)為

P(t|ω)=

(3)

根據(jù)貝葉斯理論的推理過程，為了解決SVM由于過多支持向量造成的過擬合現(xiàn)象，需給定權(quán)重ω某個(gè)先驗(yàn)概率分布，因此合理的假設(shè)權(quán)重ω符合Gaussian正態(tài)分布，則

(4)

其中：α=(α0,α1,α2,…,αN)為引入控制權(quán)重wi的向量，稱為超參數(shù)；N(·)為正態(tài)分布函數(shù)。

RVM的核心特點(diǎn)就是引入超參數(shù)α對(duì)決策函數(shù)的權(quán)重進(jìn)行控制，當(dāng)模型訓(xùn)練結(jié)束時(shí)，如果超參數(shù)很大時(shí)，對(duì)應(yīng)權(quán)重的值為零，極少數(shù)非零權(quán)重對(duì)應(yīng)的基函數(shù)構(gòu)成了最終的決策函數(shù)，實(shí)現(xiàn)了模型的稀疏化[2]。非零權(quán)重對(duì)應(yīng)的訓(xùn)練樣本的向量稱為相關(guān)向量。

對(duì)于給定的α值，如果已知訓(xùn)練樣本的概率分布已知，根據(jù)貝葉斯定理可推導(dǎo)出ω的后驗(yàn)概率分布

(5)

其中：p(t|α)為證據(jù)函數(shù)，p(t|ω)為似然函數(shù)。

由于權(quán)重P(t|ω)不符合高斯分布，沒有解析解，因此采用Laplace逼近方法來近似。又由于

(6)

因此，為了求解在已知α條件下的ω的最可能的權(quán)重ωMP，即等價(jià)于求

(7)

(8)

其中：yi=σ{y(xi;ω)}，A=diag(α0,α1,α2,…,αN)。

為了求解ωMP，需要對(duì)其進(jìn)行不斷的迭代求以取最大值，可采用迭代加權(quán)最小二乘方法來求解。在求解過程中，需要首先求得對(duì)數(shù)后驗(yàn)分布的梯度向量g及Hessian矩陣H，由式(8)微分可得：

g=ωlog{p(t|ω)p(ω|α)}=ΦT(t-y)-Aω(9)

H=ωωlog{p(t|ω)p(ω|α)}=-(ΦTBΦ+A) (10)

其中：B=diag(β1,β2,…,βN)，βn=σ{y(xn)}[1-σ{y(xn)}],Φ稱之為設(shè)計(jì)矩陣，且Φnm=K(xn,xm-1)。

利用Laplace逼近方法中Gaussian正態(tài)分布逼近后驗(yàn)概率P(t|ω)，是對(duì)后驗(yàn)概率分布的眾數(shù)位置處函數(shù)的二次逼近，此時(shí)的后驗(yàn)協(xié)方差

Σ′=(-H|ωMP)-1=(ΦTBΦ+A)-1

(11)

對(duì)于K(K>2)類多分類問題，似然函數(shù)(3)可以表達(dá)為

(12)

此時(shí)的輸出值為多值輸出yk(xn;ωk)，每一個(gè)輸出都有各自獨(dú)立的權(quán)重向量ωk和超參數(shù)向量αk。如果直接求解多分類問題時(shí)，K類的Hessian矩陣H及Σ′都會(huì)非常大，導(dǎo)致計(jì)算量也很大，因此可將多分類問題分解成多個(gè)二分類問題來解決。

2基于量子遺傳算法的核函數(shù)參數(shù)優(yōu)化

2.1過程描述

核函數(shù)的選擇對(duì)于最終診斷率有著很大的影響，然而至今也未形成完整的理論解決該問題。在相關(guān)的研究中，高斯徑向基核函數(shù)的表現(xiàn)通常較好。因此，選用高斯徑向基核函數(shù)作為RVM的基函數(shù)。

群智能優(yōu)化算法被引入到核參數(shù)的選擇中，利用群智能算法卓越的尋優(yōu)能力，自適應(yīng)的選取最優(yōu)的核函數(shù)參數(shù)，如GA[7]和PSO[8]。然而GA存在著收斂過慢、易陷入局部最優(yōu)等缺陷。PSO也存在著易陷入局部最優(yōu)、局部尋優(yōu)能力較差等缺點(diǎn)。QGA是在傳統(tǒng)GA中引入量子計(jì)算的概率和機(jī)制形成的新型算法，充分利用量子計(jì)算具有天然的并行性，極大地加快了處理信息的速度。與傳統(tǒng)GA不同的是，QGA中的染色體用量子比特取代GA中的二進(jìn)制串，采用量子比特編碼，使得一個(gè)染色體可以同時(shí)表達(dá)多個(gè)量子態(tài)的疊加，即僅通過一個(gè)小數(shù)量種群的量子個(gè)體來代替?zhèn)鹘y(tǒng)數(shù)量較大的個(gè)體。此外，量子門操作這一獨(dú)特的處理步驟，使得QGA擁有良好的全局尋優(yōu)能力。QGA比GA擁有更好的多樣性特征以及更佳的收斂性。因此，采用GA的改進(jìn)算法QGA自適應(yīng)地選取最優(yōu)核參數(shù)。圖1為QGA優(yōu)化RVM核參數(shù)的流程圖。

圖1　基于QGA優(yōu)化RVM核參數(shù)的流程圖 Fig.1 The flow chart of the parameter optimization of RVM based on QGA

具體的核參數(shù)優(yōu)化步驟如下：

步驟1QGA將核參數(shù)看做一個(gè)染色體。隨機(jī)產(chǎn)生N條以量子比特為編碼的染色體對(duì)核參數(shù)進(jìn)行編碼，構(gòu)成初始種群

其中一個(gè)量子比特表示為|ψ〉=α|0〉+β|1〉，α、β為幅常數(shù)，m為一條量子染色體的長度。

步驟2對(duì)初始的種群P(0)中的每個(gè)個(gè)體進(jìn)行一次二進(jìn)制的測量，將種群P(0)轉(zhuǎn)化為長度m的二進(jìn)制串構(gòu)成的種群Q(0)。

步驟3對(duì)種群Q(0)進(jìn)行適應(yīng)度評(píng)估。對(duì)隨機(jī)產(chǎn)生的待選核參數(shù)值的二進(jìn)制編碼轉(zhuǎn)化成十進(jìn)制，即形成初始的核參數(shù)數(shù)組，將其每一個(gè)核參數(shù)值分別代入核函數(shù)中，利用RVM對(duì)待測數(shù)據(jù)集進(jìn)行分類，對(duì)適應(yīng)度(此時(shí)為分類準(zhǔn)確率)進(jìn)行評(píng)估，保留當(dāng)前最佳適應(yīng)度(最佳分類準(zhǔn)確率)下的核參數(shù)值和對(duì)應(yīng)的適應(yīng)度，即當(dāng)前最優(yōu)個(gè)體。

步驟4通過量子門進(jìn)化操作來更新P(0)獲得新的種群P(1)。選擇量子旋轉(zhuǎn)門使算法朝著最優(yōu)解的方向搜索，此時(shí)

步驟5令t=t+1,并返回步驟2。以此不斷迭代更新，最終獲取最優(yōu)的RVM核函數(shù)參數(shù)。

2.2關(guān)鍵參數(shù)的設(shè)置

在利用QGA自適應(yīng)地獲取最優(yōu)RVM核函數(shù)參數(shù)之前，需要對(duì)優(yōu)化算法中的一些關(guān)鍵的參數(shù)進(jìn)行初始化設(shè)置，主要包括以下參數(shù)的設(shè)置：

初始化種群。P(0)中各個(gè)(α,β)均初始化為(1/sqrt(2),1/sqrt(2))，表示在初始進(jìn)化時(shí)一個(gè)染色體所有狀態(tài)均以相同的概率進(jìn)行疊加。

適應(yīng)度函數(shù)Fitness。采用QGA優(yōu)化RVM的目的是為了獲得最優(yōu)的分類準(zhǔn)確率。因此，選取分類的準(zhǔn)確率作為QGA優(yōu)化核函數(shù)參數(shù)中的適應(yīng)度函數(shù)，即

算法的結(jié)束條件。當(dāng)最大的迭代數(shù)超過初始的設(shè)置代數(shù)或當(dāng)連續(xù)10代的最佳適應(yīng)度(分類準(zhǔn)確率)之差的絕對(duì)值小于0.001，此時(shí)核參數(shù)的優(yōu)化程序終止。

2.3仿真分析

為了驗(yàn)證QGA在RVM核參數(shù)優(yōu)化中的有效性，采用高斯徑向基函數(shù)作為核函數(shù)，選取UCI標(biāo)準(zhǔn)數(shù)據(jù)庫中的Wine、Iris以及Cancer三個(gè)標(biāo)準(zhǔn)數(shù)據(jù)集作為實(shí)驗(yàn)數(shù)據(jù)，分別采用QGA、GA和PSO選取RVM對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分類時(shí)的最優(yōu)核參數(shù)。仿真實(shí)驗(yàn)在CPU速度為3.30 GHz，內(nèi)存為4 GB的計(jì)算機(jī)下進(jìn)行，基于Matlab平臺(tái)，對(duì)比研究了QGA優(yōu)化RVM(QGA-RVM)、GA優(yōu)化RVM(GA-RVM)以及PSO優(yōu)RVM(PSO-RVM)三種方法之間選取最優(yōu)核參數(shù)方面的性能，以此來分析QGA在核參數(shù)優(yōu)化中的優(yōu)越性。表1為UCI數(shù)據(jù)庫中4種測試數(shù)據(jù)的樣本個(gè)數(shù)、樣本維數(shù)及類數(shù)的說明。

表1　測試數(shù)據(jù)集說明

表2顯示了不同種方法優(yōu)化RVM后對(duì)UCI標(biāo)準(zhǔn)數(shù)據(jù)庫3種數(shù)據(jù)進(jìn)行50次分類實(shí)驗(yàn)后的平均分類準(zhǔn)確率以及平均運(yùn)行時(shí)間。從表2中的結(jié)果可見，對(duì)于Wine、Iris、Glass三種測試數(shù)據(jù)集，三種不同算法優(yōu)化后的RVM分類的準(zhǔn)確率效果相差不大。相比較而言，QGA-RVM的平均分類準(zhǔn)確率最高，GA-RVM和PSO-RVM對(duì)不同數(shù)據(jù)的測試結(jié)果互有高低。此結(jié)果說明了QGA算法在測試過程中每次都能自適應(yīng)的選取到最佳的核函數(shù)參數(shù)，而GA和PSO算法在測試過程中出現(xiàn)了核函數(shù)參數(shù)值陷入局部最優(yōu)、早熟收斂等問題，未收斂到最優(yōu)的核參數(shù)值。從測試的時(shí)間上看，QGA-RVM的平均運(yùn)行時(shí)間最低，體現(xiàn)了QGA算法在優(yōu)化RVM過程中并行計(jì)算的優(yōu)勢。

表2　不同算法對(duì)測試數(shù)據(jù)分類結(jié)果

3滾動(dòng)軸承故障診斷實(shí)例分析

3.1實(shí)驗(yàn)數(shù)據(jù)

滾動(dòng)軸承故障試驗(yàn)的數(shù)據(jù)來自于美國Case Western Reserve University軸承數(shù)據(jù)中心。試驗(yàn)采用的軸承型號(hào)為6205-2RS JEM SKF深溝球軸承，采樣頻率fs=12 000 Hz，軸承的轉(zhuǎn)速為1 748 r/min。實(shí)驗(yàn)中分別獲取不同載荷 (0 ph,1 ph,3 ph)、不同故障程度下三種軸承故障的振動(dòng)測試數(shù)據(jù)以及正常軸承數(shù)據(jù)。三種故障分別為外圈點(diǎn)蝕故障、內(nèi)圈點(diǎn)蝕故障及滾動(dòng)體點(diǎn)蝕故障。其中點(diǎn)蝕故障通過電火花加工生成，其故障程度為0.178～1.016cm。

3.2故障特征提取

EEMD(Ensemble Empirical Mode Decomposition)是一種基于EMD(Empirical Mode Decomposition，經(jīng)驗(yàn)?zāi)B(tài)分解)的改進(jìn)的信號(hào)自適應(yīng)分解方法，通過對(duì)待分解信號(hào)加入白噪聲后再進(jìn)行多次分解求平均，有效地減少了EMD分解中的模態(tài)混淆現(xiàn)象[10]。EEMD作為一種非線性、非平穩(wěn)信號(hào)有效的時(shí)頻分析方法，近年來，在旋轉(zhuǎn)機(jī)械故障診斷中得到了有效應(yīng)用[11-12]。原始的復(fù)雜信號(hào)可以被EEMD自適應(yīng)的分解為多個(gè)內(nèi)稟模態(tài)函數(shù)(IMF)之和，不同的IMF分量包含著不同的頻率信息成分。不同的機(jī)械故障發(fā)生時(shí)，振動(dòng)信號(hào)中就會(huì)產(chǎn)生不同的頻率成分，不同頻帶上振動(dòng)信號(hào)的能量就會(huì)不同，即各個(gè)IMF上的能量分布不同[13]。因此，提取滾動(dòng)軸承經(jīng)EEMD分解后的IMF能量作為故障特征。

設(shè)加入了白噪聲后的待分解滾動(dòng)軸承振動(dòng)信號(hào)為

xi(t)=x(t)+ni(t)

(13)

其中：x(t)為原始信號(hào)，ni(t)為第i次加入的白噪聲，i=1,2,…,M。其中M為初始的平均次數(shù)。

用原始的EMD算法將xi(t)自適應(yīng)分解為若干個(gè)IMF之和，即

(14)

其中：n=1,2,…,N，N為分解后IMF的個(gè)數(shù)，ci,n(t)為第i次加入白噪聲后分解得到的IMF分量，ri,N(t)為本次分解后剩余的殘差。

通過添加不同的白噪聲，重復(fù)上述兩個(gè)步驟M次后得到每次分解后IMF分量的集合為

[{c1,n(t)},c2,n(t),…,cM,n(t)]

(15)

對(duì)每次分解后的IMF分量進(jìn)行總體平均運(yùn)算來消除白噪聲對(duì)實(shí)際IMF的影響，得到最終的IMF

(16)

其中：cn(t)為EEMD分解后的第n個(gè)IMF。

滾動(dòng)軸承振動(dòng)信號(hào)x(t)經(jīng)過EEMD最終分解為

(17)

其中：r(t)為最終的剩余殘差。

EEMD分解滾動(dòng)軸承振動(dòng)信號(hào)形成的第n個(gè)IMF的能量可以表達(dá)為

(18)

假設(shè)r(t)忽略不計(jì)，信號(hào)整體的能量值E應(yīng)該等于IMF的能量之和，即

E=E1+E2+…+EN

(19)

能量值通常較大，為了便于故障診斷模型的最終診斷，對(duì)Ei進(jìn)行歸一化處理后構(gòu)建故障的特征向量

T=[P1,P2,…,PN]

(20)

3.3效果分析

選擇載荷3 ph、故障程度為1.63cm的實(shí)驗(yàn)數(shù)據(jù)為本研究的分析數(shù)據(jù)。每個(gè)樣本由2 048個(gè)連續(xù)采樣數(shù)據(jù)構(gòu)成，每種滾動(dòng)軸承狀態(tài)下有59個(gè)樣本。每種狀態(tài)下任意選取20個(gè)樣本作為RVM診斷模型的訓(xùn)練樣本(共80個(gè))，剩余的樣本(共156個(gè))構(gòu)成模型的測試樣本。

對(duì)滾動(dòng)軸承4種狀態(tài)下的振動(dòng)測試信號(hào)進(jìn)行EEMD分解，根據(jù)上述故障特征提取的步驟，分別提取4種狀態(tài)下的IMF能量作為滾動(dòng)軸承的故障特征。以外圈故障信號(hào)為例，設(shè)置M=100，白噪聲的幅值為0.2。圖2顯示了EEMD分解過程。見圖2，外圈故障信號(hào)的一個(gè)樣本被分解成為11個(gè)IMF:c1(t)，c2(t)，…，c11(t)和一個(gè)殘差r(t)。

圖2　外圈故障信號(hào)EEMD分解圖 Fig.2 EEMD decomposition results of outer race fault signals

經(jīng)過分析，提取前6個(gè)IMF的能量即可反映出各種狀態(tài)下能量的差別，以此構(gòu)成診斷模型的輸入特征向量。提取的IMF能量被分成訓(xùn)練樣本和測試樣本。圖3顯示了滾動(dòng)軸承4種狀態(tài)下的IMF能量分布圖。

圖3　不同狀態(tài)下EEMD能量分布圖 Fig.3 EEMD energy spectrum of different conditions

對(duì)于4種滾動(dòng)軸承狀態(tài)智能診斷的多分類問題，采用“一對(duì)一”策略，構(gòu)建6個(gè)RVM二類故障診斷器進(jìn)行診斷。將訓(xùn)練樣本分別輸入到QGA-RVM中，選用高斯徑向基核，利用QGA算法自適應(yīng)的選取各自最優(yōu)的核參數(shù)σ。將測試樣本分別輸入到各自訓(xùn)練好的QGA-RVM中，最終智能的輸出診斷結(jié)果。為了驗(yàn)證所提診斷模型的有效性和穩(wěn)定性，以10次故障診斷實(shí)驗(yàn)后的平均故障診斷率作為評(píng)價(jià)指標(biāo)，同時(shí)將滾動(dòng)軸承的IMF能量輸入到RVM、 SVM以及QGA優(yōu)化的SVM(QGA-SVM)中做對(duì)比分析研究。其中各個(gè)模型均采用高斯徑向基核函數(shù)，由于RVM及SVM沒有采用優(yōu)化算法，其核參數(shù)人為設(shè)置，最終的診斷實(shí)驗(yàn)結(jié)果見表3。

表3　不同種診斷模型的故障診斷實(shí)驗(yàn)結(jié)果

從表3中的結(jié)果可見，QGA-RVM診斷模型及QGA-SVM診斷模型的平均診斷率最高，達(dá)到100%，說明基于QGA的核參數(shù)選擇方法在RVM和SVM的核參數(shù)優(yōu)化中均達(dá)到了最優(yōu)值。而其它兩種診斷模型的平均診斷率相對(duì)較低，但也基本上能夠識(shí)別出主要故障，其主要原因是由于人為設(shè)置了核參數(shù)的值，未設(shè)置成最優(yōu)值。表3中同時(shí)也顯示了各個(gè)診斷模型的訓(xùn)練時(shí)間和測試時(shí)間。QGA-RVM、QGA-SVM和RVM的平均訓(xùn)練時(shí)間比SVM的平均訓(xùn)練時(shí)間要長，其中QGA-RVM的時(shí)間最長。由于RVM在分類訓(xùn)練過程中需要不斷的迭代，影響了其訓(xùn)練的效率。而QGA優(yōu)化RVM的過程需要再次的迭代，因而再次增加了模型的訓(xùn)練時(shí)間。但實(shí)際的故障診斷中，診斷模型都是在診斷之前已經(jīng)訓(xùn)練完成，因此不影響整體故障診斷過程的效率。對(duì)于診斷模型的測試時(shí)間，QGA-RVM的測試時(shí)間最少，RVM次之，QGA-SVM和SVM的測試時(shí)間最長。這是由于RVM需要的相關(guān)向量比SVM支持向量更少，即RVM的決策函數(shù)更加稀疏。雖然QGA-RVM和QGA-SVM都能準(zhǔn)確的診斷出故障，但QGA-RVM的測試時(shí)間更短。因此，基于以上分析，RVM更加適合處理小樣本問題，同時(shí)更加適合應(yīng)用于在線故障診斷。

4結(jié)論

提出了一種基于QGA優(yōu)化RVM核函數(shù)參數(shù)的方法，利用了QGA具有全局搜索尋優(yōu)能力強(qiáng)、種群多樣性、收斂速度快及并行計(jì)算時(shí)間短等特點(diǎn)，自適應(yīng)地選取最優(yōu)核參數(shù)，并將其應(yīng)用于滾動(dòng)軸承的智能故障診斷中。通過仿真，將GQA優(yōu)化的RVM與PSO以及GA優(yōu)化的RVM進(jìn)行數(shù)據(jù)分類實(shí)驗(yàn)對(duì)比，利用UCI 標(biāo)準(zhǔn)數(shù)據(jù)集中的三個(gè)數(shù)據(jù)作為被測數(shù)據(jù)，結(jié)果表明了QGA在核函數(shù)參數(shù)尋優(yōu)方面的優(yōu)勢。利用EEMD將滾動(dòng)軸承的故障信號(hào)自適應(yīng)的分解成若干個(gè)IMF，提取各個(gè)IMF能量作為故障特征，并將其輸入到QGA-RVM進(jìn)行最終的故障診斷。結(jié)果表明，該方法能夠快速、準(zhǔn)確地診斷出滾動(dòng)軸承的故障；同時(shí)，將EEMD能量特征輸入到SVM中進(jìn)行對(duì)比研究，同樣也顯示了RVM在故智能故障診斷中的優(yōu)越性。

參考文獻(xiàn)

[1]Tipping M E. Sparse Bayesian learning and the relevance vector machine[J]. Journal of Machine Learning Research, 2001,1(3):211-244.

[2]Tipping M E. Bayesian inference: An introduction to principles and practice in machine learning[J]. Advanced Lectures on Machine Learning, 2004,3176:41-62.

[3]周曉英, 李巍華, 丁康. 基于關(guān)聯(lián)向量機(jī)的齒輪故障檢測方法研究[J]. 振動(dòng)與沖擊, 2008,27(6):51-54.

ZHOU Xiao-ying, LI Wei-hua, DING Kang. Gear fault detection based on relevance vector machine[J]. Journal of Vibration and Shock, 2008,27(6):51-54.

[4]Widodo A, Yang B S, Kim E Y, et al. Fault diagnosis of low speed bearing based on acoustic emission signal and multi-class relevance vector machine[J]. Nondestructive Testing and Evaluation, 2009,24(4):313-328.

[5]Li N, Liu C, He C, et al. Gear fault detection based on adaptive wavelet packet feature extraction and relevance vector machine[J].Proceedings of the Institution of Mechanical Engineers part C-Journal of Mechanical Engineering Science, 2011,225(C11):2727-2738.

[6]He C, Liu C, Li Y, et al. Intelligent fault diagnosis of rotating machinery based on multiple relevance vector machines with variance radial basis function kernel[J]. Proceedings of the Institution of Mechanical Engineers Part C-Journal of Mechanical Engineering Science, 2011,225(C7):1718-1729.

[7]董超, 田聯(lián)房, 趙慧潔. 遺傳關(guān)聯(lián)向量機(jī)高光譜影像分類[J]. 上海交通大學(xué)學(xué)報(bào), 2011(10):1516-1520.

DONG Chao,TIAN Lian-fang, ZHAO Hui-jie. Hyperspectral image classification by genetic relevance vector machine[J]. Journal of Shanghai Jiaotong University, 2011(10):1516-1520.

[8]Du W L, Li A S, Ye P F, et al. Fault diagnosis of plunger pump in truck crane based on relevance vector machine with particle swarm optimization algorithm[J]. Shock and Vibration, 2013,20(4):781-792.

[9]楊俊安, 莊鎮(zhèn)泉, 史亮. 多宇宙并行量子遺傳算法[J]. 電子學(xué)報(bào), 2004(6):923-928.

YANG Jun-an, ZHUANG Zhen-quan, SHI Liang. Multi-universe parallel quantum genetic algorithm[J]. Acta Electronica Sinica, 2004(6):923-928.

[10]Wu Z, Guang N E. Ensemble empirical mode decomposition: A noise-assisted data analysis method[J]. Advances in Adaptive Data Analysis, 2009,1(1):1-41.

[11]沈長青, 謝偉達(dá), 朱忠奎, 等. 基于EEMD和改進(jìn)的形態(tài)濾波方法的軸承故障診斷研究[J]. 振動(dòng)與沖擊, 2013,32(2):39-43.

SHEN Chang-qing, Peter W.Tse, ZHU Zhong-kui, et al. Rolling element bearing fault diagnosis based on EEMD and improved morphological filtering method[J]. Journal of Vibration and Shock, 2013,32(2):39-43.

[12]周智, 朱永生, 張優(yōu)云, 等. 基于EEMD和共振解調(diào)的滾動(dòng)軸承自適應(yīng)故障診斷[J]. 振動(dòng)與沖擊, 2013,32(2):76-80.

ZHOU Zhi, ZHU Yong-sheng, ZHANG You-yun, et al. Adaptive fault diagnosis of rolling bearings based on EEMD and demodulated resonance[J]. Journal of Vibration and Shock, 2013,32(2):76-80.

[13]Yu Y, Yu D J, Cheng J S. A roller bearing fault diagnosis method based on EMD energy entropy and ANN[J]. Journal of Sound and Vibration, 2006,294(1-2):269-277.