摘 要：煙氣輪機(jī)機(jī)組是利用余熱發(fā)電的原理回收高溫?zé)崮茉偕娔?，但是由于多方面的原因，煙氣輪機(jī)機(jī)組時(shí)常發(fā)生故障，對生產(chǎn)的經(jīng)濟(jì)性和安全性造成巨大損失。為了保障設(shè)備正常運(yùn)行，對設(shè)備進(jìn)行實(shí)時(shí)監(jiān)測，及時(shí)判斷故障發(fā)生的可疑點(diǎn)，減少故障發(fā)生的頻率，同時(shí)保護(hù)環(huán)境，采用對原始信號的3層小波分解，提取信號的特征向量，再采用CSVC和V參數(shù)支持向量機(jī)二種方法進(jìn)行分類實(shí)驗(yàn)，得到V參數(shù)的支持向量機(jī)分類效果最佳。



關(guān)鍵詞：煙氣輪機(jī); 支持向量機(jī); 故障診斷; 分類

中圖分類號：

TN9834; TP277

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

文章編號：1004373X(2012)05

0043

05



Application of V parameter of SVM for fault diagnosis in gas turbine



WANG Shufang

(Information Technology College， Beijing Institute of Petrochemical Technology， Beijing 102600， China)



Abstract：

The gas turbine unit usually uses the principle of hightemperature heat recovery renewable energy， but due to various reasons， the failure of gas turbine unit causes huge losses for economical efficiency and security of production. To ensure normal operation， a realtime monitoring was performed， the suspicious points were determined to reduce the frequency of failure and protect the environment. The signal feature vector was extracted through decomposing three layers wavelet of the original signal， a classified experiment was performed with CSVC and V parameters， the best V parameters of support vector machine classification was obtained.



Keywords： gas turbine; support vector machine; fault diagnosis; classification





收稿日期：20110928

0 引 言

煉油廠中的煙氣輪機(jī)機(jī)組是利用原油催化裂化裝置在生產(chǎn)過程中產(chǎn)生的高溫余熱驅(qū)動煙氣輪機(jī)的空氣壓縮機(jī)做功或給發(fā)電機(jī)機(jī)組提供輸入能量，發(fā)電機(jī)組進(jìn)而輸出電能。這種方法是目前國內(nèi)外催化裂化工藝環(huán)節(jié)回收能量最有效的方法。以燕山石化公司煉油廠YL18000A型煙氣輪機(jī)直接發(fā)電機(jī)組為例，該機(jī)組平均每天再生發(fā)電10 368 kWh，按當(dāng)?shù)仉妰r(jià)每千瓦時(shí)0.50元計(jì)算，每年節(jié)約人民幣5 184萬元。煙氣輪機(jī)機(jī)組的投入運(yùn)行不僅節(jié)約了大量的能源，為石化企業(yè)帶來了巨大的經(jīng)濟(jì)效益，而且通過回收利用催化裂化裝置產(chǎn)生的煙氣，在保護(hù)環(huán)境方面也發(fā)揮了很好的作用。但是由于機(jī)械設(shè)備質(zhì)量、工藝過程參數(shù)配置、電氣電子設(shè)備控制方法等各方面的原因，煙氣輪機(jī)機(jī)組振動故障率一致比較高，因此做好煙機(jī)機(jī)組的狀態(tài)監(jiān)測和故障診斷，采取相應(yīng)控制措施，保證機(jī)組安全、穩(wěn)定、長期運(yùn)行具有十分重要的現(xiàn)實(shí)意義。

目前，國外有些公司正在從事煙氣輪機(jī)發(fā)電系統(tǒng)的故障診斷研究，但其成果引進(jìn)價(jià)格昂貴，更重要的是其產(chǎn)品難以與現(xiàn)有設(shè)備配套，引進(jìn)就要對現(xiàn)有設(shè)備進(jìn)行配套改造，投資很大。國內(nèi)也有一些科研院所在做這方面的初步研究。課題研究的具體方法是對從煉油廠采集的煙氣輪機(jī)故障發(fā)生時(shí)段的數(shù)據(jù)信息進(jìn)行篩選，對故障信號用小波包分析方法提取故障發(fā)生時(shí)刻的特征信息，最后將特征樣本通過支持向量機(jī)訓(xùn)練器進(jìn)行訓(xùn)練，并獲得支持向量與分類模型，用于對測試樣本進(jìn)行分類和函數(shù)回歸，通過對相關(guān)已知信號數(shù)據(jù)的分析得到煙起輪機(jī)的瞬時(shí)運(yùn)行狀態(tài)。

在此，首先使用Matlab平臺的小波分析工具對信號進(jìn)行奇異點(diǎn)分析，完成信號特征提取工作。然后在此基礎(chǔ)上進(jìn)行故障診斷研究，使用的分類及回歸程序是基于林智仁老師的一個(gè)支持向量機(jī)工具的Java實(shí)現(xiàn)的，該支持向量機(jī)工具支持CSVC支持向量機(jī)、V參數(shù)的支持向量機(jī)，在控制臺環(huán)境下，實(shí)現(xiàn)對樣本進(jìn)行二分類、多分類、支持向量回歸等功能，圖形化的實(shí)驗(yàn)結(jié)果及分析能直觀地反映診斷的效果。 

1 信號的采集

煙氣輪機(jī)在工作過程中，由于發(fā)電機(jī)葉片引起的軸振動、煙氣輪機(jī)設(shè)備的溫度、入口溫度、出口溫度和煙機(jī)內(nèi)的壓力等多個(gè)參數(shù)都是體現(xiàn)故障的指標(biāo)。實(shí)踐證明，故障發(fā)生時(shí)煙機(jī)機(jī)組設(shè)備包括發(fā)電機(jī)的橫向、縱向、軸向振動值，還有煙機(jī)溫度、壓力都會發(fā)生較大變動。而振動值的變化最能體現(xiàn)故障發(fā)生的特征，故采用振動值數(shù)據(jù)經(jīng)過小波處理之后提取的特征向量作為樣本訓(xùn)練支持向量機(jī)。作為故障診斷的依據(jù)，信號數(shù)據(jù)采集點(diǎn)情況如圖1所示。圖1中煙氣輪機(jī)工作狀態(tài)下的實(shí)時(shí)監(jiān)測點(diǎn)分布圖包括兩個(gè)與煙機(jī)相連的發(fā)電機(jī)的各個(gè)檢測點(diǎn)振動、溫度、壓力差參數(shù)值。

圖2是從煉油廠控制中心工作站采集到的一張數(shù)據(jù)波形圖(樣本)。

圖2 煙氣輪機(jī)振動量及過程量波形

圖中橫軸代表時(shí)間量，可以具體到秒，縱軸代表的是振動量和過程量，即不同測試點(diǎn)的具體振動值，單位是μm。圖2中過程量一共4條曲線，其中不同顏色的曲線代表不同測試點(diǎn)的過程量參數(shù)，藍(lán)色曲線代表煙機(jī)設(shè)備在Yt7701A測試點(diǎn)通過振動檢測儀獲取的振動數(shù)據(jù)信息。從圖2中可以看出，在煙氣輪機(jī)設(shè)備上具體監(jiān)測位置和旋轉(zhuǎn)軸上的振動數(shù)據(jù)；以此類推，粉色曲線代表Yt7701A的振動波形；綠色曲線代表Yt7702A的振動波形；灰色曲線代表Yt7702B的振動波形，每個(gè)測試點(diǎn)有三個(gè)方向的監(jiān)測振動數(shù)據(jù)，A代表水平方向的振動，B代表垂直方向，而軸向的即是過程量波形顯示的數(shù)據(jù)，在煙機(jī)實(shí)體圖中以Zt表示。圖2顯示的波形是在Zt7701A檢測點(diǎn)獲取的振動數(shù)據(jù)信息。通過實(shí)時(shí)監(jiān)測，可以對煙機(jī)的工作狀態(tài)進(jìn)行直觀的判斷，圖2選取的是燕山石化煉油廠2004年11月24日—2005年11月29日一年的振動數(shù)據(jù)?？梢灾庇^地看出，這一年中在不同時(shí)期，煙機(jī)中不同位置的振動變化，其中波形發(fā)生劇烈變化的地方就可能是故障發(fā)生點(diǎn)，從煉油廠的歷史記錄中可以查詢到這一年中的幾次故障數(shù)據(jù)、故障相關(guān)圖片信息。

1.1 信號特征的提取

對故障信號的預(yù)分析方法是在Matlab平臺載入原始信號，截取故障部分信號，利用小波包將頻帶進(jìn)行多層劃分，對多分辨分析沒有細(xì)分的高頻部分進(jìn)一步分解，并能夠根據(jù)被分析信號的特征，自適應(yīng)地選擇相應(yīng)頻帶，使之與信號頻譜相匹配，從而提高時(shí)頻分辨率。

1.2 故障信號層次分解與重構(gòu)

采用Daubechies小波包對信號依次進(jìn)行了3，4，6，10層分解，對信號的高頻部分做更加細(xì)致的刻畫，現(xiàn)以6層分解情況為例進(jìn)行說明。利用Matlab平臺的小波工具箱的信號單尺度一維離散小波分解函數(shù)wpdec對信號進(jìn)行6層分解，生成小波分解二叉樹，然后使用Matlab小波工具箱中的wrpcoef對第6層分解的細(xì)節(jié)系數(shù)進(jìn)行重構(gòu)，可以獲得各個(gè)頻帶子信號的重構(gòu)系數(shù)向量。圖3是某個(gè)小波分解故障信號的具體情況，圖4是重構(gòu)后的各層細(xì)節(jié)信號。

其中，圖3(a)是分解的完整二叉樹結(jié)構(gòu)，6層分解一共能夠獲取26 =64個(gè)頻段的節(jié)點(diǎn)信息，樹的葉節(jié)點(diǎn)分別為(6，0)，(6，1)，(6，2)等，以此類推到(6，63)；圖3(b)是節(jié)點(diǎn)(6，1)的波形圖，通過觀察不同的節(jié)點(diǎn)信息，可發(fā)現(xiàn)分解層數(shù)越高，不同頻段的波形信息差別越明顯。

1.3 特征向量

特征向量的選擇在智能化故障診斷中是一個(gè)非常重要的環(huán)節(jié)，因?yàn)樘卣魇求w現(xiàn)故障的直接標(biāo)準(zhǔn)。選擇節(jié)點(diǎn)系數(shù)平方和歸一化之后的信息作為該頻段的能量指標(biāo)，因?yàn)楣收闲盘柵c正常信號在同一頻段的能量是不同的，所以可以將此作為機(jī)器學(xué)習(xí)的訓(xùn)練樣本。

分別計(jì)算分解后信號第6層頻段的系數(shù)節(jié)點(diǎn)系數(shù)平方和，如果xij(i=0，1，…，7;j=0，1，…，n)代表各頻段重構(gòu)信號的系數(shù)信息，則∑nj=1|xij|2即體現(xiàn)各層分解系數(shù)序列的能量，然后對能量信息做歸一化處理(用各個(gè)系數(shù)能量除以總能量)，并寫入數(shù)據(jù)文件中。研究中對于小波基及分解層次進(jìn)行了多次試驗(yàn)，對比計(jì)算量及分解效果，確定對原始信號的基于DB2小波的三層分解取得故障信號特征，作為支持向量機(jī)訓(xùn)練器的輸入樣本進(jìn)行訓(xùn)練作為實(shí)驗(yàn)方案。

圖5是實(shí)驗(yàn)中對原始信號進(jìn)行基于小波包三層分解的系數(shù)計(jì)算出的特征向量，三層分解并重構(gòu)信號在各層的系數(shù)之后，以平方和相加的方式計(jì)算第三層的所有系數(shù)的能量值，并進(jìn)行歸一化處理，最后輸出到文本格式的文件中，以每個(gè)頻段的能量指標(biāo)作為特征能夠體現(xiàn)故障的特點(diǎn)，作為支持向量機(jī)訓(xùn)練器的輸入樣本，能夠達(dá)到故障診斷實(shí)驗(yàn)的目的。

2 支持向量機(jī)簡介

支持向量機(jī)方法在解決小樣本、非線性等模式識別問題中表現(xiàn)出的優(yōu)勢能很好地應(yīng)用到機(jī)器學(xué)習(xí)領(lǐng)域，并能推廣到函數(shù)擬合等其他機(jī)器學(xué)習(xí)問題中［1］。

2.1 基于支持向量機(jī)的二分類方法

煙氣輪機(jī)的故障狀態(tài)與正常狀態(tài)的區(qū)分可以用支持向量機(jī)二分類方法實(shí)現(xiàn)。首先對原始故障數(shù)據(jù)的預(yù)處理，再分別選取正常狀態(tài)與故障狀態(tài)的特征向量作為支持向量機(jī)的訓(xùn)練數(shù)據(jù)，通過訓(xùn)練器進(jìn)行訓(xùn)練，進(jìn)而實(shí)現(xiàn)分類工作。支持向量機(jī)從線性可分的最優(yōu)分類面發(fā)展而來，基本思想可用圖6的二維情況說明。圖中，實(shí)心點(diǎn)和空心點(diǎn)代表兩類樣本，H為分類線，H1和H2分別為過各類中離分類線最近的樣本，且平行于分類線的直線，它們之間的距離叫作分類間隔。能將兩類正確分開，而且使分類間隔最大的分類線為最優(yōu)分類線，其方程為x·w+b=0，對它進(jìn)行歸一化，使得對于線性可分的樣本集(xi，yi)，i=1，…，n， x∈Rd，y∈{+1，－1}，滿足：



yi［(w·xi)+b］－1≥0，i=1，…，n

(1)



一種最大間隔分類器是這樣的：固定函數(shù)間隔為1(函數(shù)間隔為1的超平面有時(shí)稱為正則超平面)，如果w是權(quán)重向量，要在正點(diǎn)x+和負(fù)點(diǎn)x-上實(shí)現(xiàn)函數(shù)間隔為1，可以如下計(jì)算函數(shù)間隔：



+b=+1

(2)

+b=－1

(3)



同時(shí)，為計(jì)算幾何間隔，必須歸一化w。幾何間隔γ是所得分類器的函數(shù)間隔：



γ=1/2(－)

=1/‖w‖2

(4)

此時(shí)分類間隔等于2/‖w‖，使間隔最大等價(jià)于使‖w‖2最小。滿足式(5)且使1/2‖w‖2最小的分類面就叫做最優(yōu)分類面，H1，H2上的訓(xùn)練樣本點(diǎn)就稱為支持向量［24］。

2.2 核函數(shù)

核函數(shù)是SVM方法的關(guān)鍵所在。核函數(shù)K的選取需要滿足Mercer條件。選擇不同的核函數(shù)可以產(chǎn)生不同的支持向量機(jī)。目前沒有固定的核函數(shù)，一般對于不同的問題選擇不同的核函數(shù)。常用核函數(shù)有線性核函數(shù)、多項(xiàng)式核函數(shù)、徑向基函數(shù)(RBF)、 多層感知機(jī)核函數(shù)。

2.3 參數(shù)選擇 

通常參數(shù)的不同選擇會對支持向量機(jī)泛化性性能有較大影響，目前對于泛化性的優(yōu)化集中在支持向量機(jī)模型的選擇上。

2.3.1 支持向量機(jī)參數(shù)

目前常用的幾種參數(shù)有C參數(shù)的支持向量機(jī)CSVC，即C參數(shù)的支持向量機(jī)分類算法及V參數(shù)的支持向量機(jī)算法。

(1) C參數(shù)

對于支持向量機(jī)中的優(yōu)化問題，會有一個(gè)C值對應(yīng)著權(quán)值二范數(shù)的最優(yōu)選擇，這個(gè)C值也給出最優(yōu)界，從而在給定權(quán)值二范數(shù)下找到松弛變量范數(shù)的最小值。目前還沒有一個(gè)統(tǒng)一的方法來決定C的最佳取值［5］，所以現(xiàn)在廣泛采用的方法是對C在0到正無窮范圍內(nèi)的不同值進(jìn)行實(shí)驗(yàn)，對C先賦一個(gè)固定的初值，通過實(shí)驗(yàn)結(jié)果精度不斷調(diào)整參數(shù)值，直到達(dá)到滿意效果，所以出現(xiàn)了一種V參數(shù)的支持向量機(jī)算法。

(2) V參數(shù)

V在一個(gè)很小的范圍 ［0，1］ 內(nèi)變化，其結(jié)果與C在0到無窮之間變化是一致的。V參數(shù)在物理意義上實(shí)際是給出了拉格朗日參數(shù)的加和的下界，從目標(biāo)函數(shù)中除去了線性項(xiàng)，同時(shí)V也是間隔誤差的訓(xùn)練集比例的上界，也是全部支持向量數(shù)目與全部樣例數(shù)目的比例的下界。因此，V是一個(gè)更加透明的參數(shù)，與特征空間的尺度無關(guān)。所以實(shí)驗(yàn)中采用V參數(shù)支持向量機(jī)是一種較好的選擇，不僅降低計(jì)算復(fù)雜度，還不影響精度。

2.3.2 核函數(shù)參數(shù)

Vapnik等人在研究中發(fā)現(xiàn)，不同的核函數(shù)對SVM性能的影響不大，反而核函數(shù)的參數(shù)對SVM性能的影響較大，核函數(shù)本質(zhì)上是提高了學(xué)習(xí)器的分線性處理能力，因此選擇合適的核函數(shù)參數(shù)，對學(xué)習(xí)機(jī)器的性能至關(guān)重要。常用的幾種核函數(shù)參數(shù)有：多項(xiàng)式核的參數(shù)d、高斯核的參數(shù)δ、多層感知機(jī)核的參數(shù)V，C。

目前關(guān)于如何確定SVM參數(shù)的方法具有不確定性，各種方法都有優(yōu)缺點(diǎn)，所以在實(shí)際情況下，對不同的參數(shù)對SVM性能的影響進(jìn)行優(yōu)化是一個(gè)值得深入研究的課題［6］。

3 V參數(shù)支持向量機(jī)在故障診斷中的應(yīng)用

3.1 研究方案

研究使用的一種支持向量機(jī)分類器用Java代碼基于MyEclipse環(huán)境下開發(fā)。其中核心算法采用了林智仁老師的Libsvm軟件包中的支持向量機(jī)算法。研究采用廣泛使用的擇優(yōu)選擇的方法，分別采用CSVM、nuSVM、徑向基和多項(xiàng)式形式作為核函數(shù)，并采用了LSSVM的核心算法通過Matlab平臺進(jìn)行模型訓(xùn)練及測試分析，同時(shí)也對幾種常見核函數(shù)進(jìn)行運(yùn)算結(jié)果的對比，比較其分類及預(yù)測效果。

分類流程圖如圖7所示。

在CSVC與V參數(shù)的分類實(shí)驗(yàn)中對原始信號進(jìn)行小波包3層分解，選取前各個(gè)頻帶的特征向量作為訓(xùn)練樣本的輸入/輸出選擇一個(gè)節(jié)點(diǎn)，兩個(gè)分類標(biāo)識為:-1表示故障發(fā)生點(diǎn);+1表示正常工作狀態(tài)。在最小二乘支持向量機(jī)故障診斷實(shí)驗(yàn)中，采用函數(shù)估計(jì)的方法對模型進(jìn)行預(yù)測，并針對不同的支持向量機(jī)參數(shù)進(jìn)行實(shí)驗(yàn)結(jié)果對比，比較實(shí)驗(yàn)效果。

3.2 訓(xùn)練樣本與測試樣本

訓(xùn)練數(shù)據(jù)集選取煙氣輪機(jī)工作正常狀態(tài)與故障狀態(tài)兩組數(shù)據(jù)，各22個(gè)樣本，總共44個(gè)訓(xùn)練樣本，作為訓(xùn)練數(shù)據(jù)，選取的樣本數(shù)據(jù)分別是煙氣輪機(jī)正常工作狀態(tài)與故障狀態(tài)下對原始信號進(jìn)行小波包分析處理過的特征向量，訓(xùn)練數(shù)據(jù)集用來訓(xùn)練支持向量機(jī)，進(jìn)而建立二分類故障分類器，測試數(shù)據(jù)集用于對支持向量機(jī)訓(xùn)練完成后檢驗(yàn)訓(xùn)練的結(jié)果，即支持向量機(jī)泛化能力的好壞。針對煙氣輪機(jī)某一個(gè)測試點(diǎn)如(Yt7701A)三個(gè)方向的振動數(shù)據(jù)，包含水平振動、垂直振動和軸向振動，進(jìn)行數(shù)據(jù)采集工作，每個(gè)方向的振動都有可能造成煙氣輪機(jī)的故障，樣本選擇多個(gè)測試點(diǎn)歷史數(shù)據(jù)中1年內(nèi)某個(gè)月或者幾個(gè)月以天為單位的數(shù)據(jù)，包括三個(gè)方向的振動值。其中，水平方向有222個(gè)數(shù)據(jù)，垂直方向上有308個(gè)數(shù)據(jù)，軸向有308個(gè)數(shù)據(jù)，采集的時(shí)刻為零點(diǎn)，負(fù)值代表的是規(guī)定的負(fù)方向。

測試樣本的選擇是正常工作狀態(tài)下的9組數(shù)據(jù)，故障狀態(tài)下的6組數(shù)據(jù)，一共15組振動數(shù)據(jù)。利用小波工具箱進(jìn)行樣本特征提取，將振動信號的能量信息作為特征向量，也作為支持向量機(jī)訓(xùn)練器的訓(xùn)練樣本。

3.3 SVM故障診斷

將SVM應(yīng)用到特定的實(shí)際問題中需要解決大量的設(shè)計(jì)問題，如算法的設(shè)計(jì)、核函數(shù)的設(shè)計(jì)等等。首先是為給定的應(yīng)用選擇適當(dāng)?shù)暮?。研究中使用了線性核函數(shù)、多項(xiàng)式核函數(shù)、徑向基核函數(shù)(RBF)和Sigmoid核函數(shù)進(jìn)行實(shí)驗(yàn)對比。其次就是實(shí)現(xiàn)哪種SVM，如果數(shù)據(jù)需要分類，就要決定是否使用最大間隔，或者某種軟間隔方法，這里的關(guān)鍵因素是數(shù)據(jù)的維數(shù)。一旦核與優(yōu)化條件確定，將進(jìn)行一系列的實(shí)驗(yàn)，調(diào)整參數(shù)不斷變化直到達(dá)到滿意的性能。

3.3.1 基于CSVC的故障診斷

標(biāo)準(zhǔn)的CSVC分類的算法形式是：

式中：C就是懲罰因子，C參數(shù)的值設(shè)置得越大表示將對錯(cuò)誤分類的懲罰越大;ξi是引入的松弛變量，它允許一定程度上違反間隔約束。實(shí)踐中，C是一個(gè)變化范圍很大的數(shù)，優(yōu)化性能的評價(jià)是通過使用獨(dú)立的驗(yàn)證集或者交叉驗(yàn)證的方法，對C在一定范圍內(nèi)變化，‖w‖會有相應(yīng)的連續(xù)變化。通過實(shí)驗(yàn)發(fā)現(xiàn)，對于不同范圍的C值，針對線性核的支持向量機(jī)，以及針對高斯徑向基核、多項(xiàng)式核和S核的支持向量機(jī)的不同參數(shù)等多次實(shí)驗(yàn)發(fā)現(xiàn)，分類效果的精確率都達(dá)到98.0%。

3.3.2 V參數(shù)SVM的故障診斷

軟間隔技術(shù)CSVC支持向量機(jī)有一個(gè)問題即參數(shù)C的選擇。首先范圍可以確定，典型的方法是在這個(gè)范圍內(nèi)試驗(yàn)，直到找到對特定訓(xùn)練集最好的選擇。同時(shí)特征空間也會進(jìn)一步影響參數(shù)的尺度。因此引入?yún)?shù)V提出如下的優(yōu)化問題［7］：

在這個(gè)參數(shù)化的過程中，l是樣本總數(shù)，V給出了αi的下界，它從目標(biāo)函數(shù)中去除了線性項(xiàng)?？梢钥闯?，v是間隔誤差的訓(xùn)練集比例的上界，同時(shí)v又是支持向量全部數(shù)目與全部樣例數(shù)比例的下界。因此，v給出了問題的一個(gè)更加透明的參數(shù)，它與空間的尺度無關(guān)，而僅與數(shù)據(jù)的噪聲程度有關(guān)?？梢奦參數(shù)可以控制支持向量的數(shù)量和誤差，所以這里采用V參數(shù)SVM方法處理故障診斷分類過程。

由于每個(gè)支持向量機(jī)只進(jìn)行二分類工作，一共只有兩種狀態(tài)，一種是正常狀態(tài)，設(shè)定為＋1，另一種是故障狀態(tài)，設(shè)定為－1。執(zhí)行訓(xùn)練器主程序svm_train.java，對提取的特征訓(xùn)練數(shù)據(jù)集進(jìn)行訓(xùn)練，選擇支持向量機(jī)所需的參數(shù)，不同支持向量機(jī)參數(shù)有差異，訓(xùn)練時(shí)在一定范圍內(nèi)對參數(shù)設(shè)置進(jìn)行實(shí)驗(yàn)，以達(dá)到最佳實(shí)驗(yàn)效果。研究中采用Libsvm中交叉實(shí)驗(yàn)法選取核函數(shù)的參數(shù)，最后將訓(xùn)練結(jié)果保存到訓(xùn)練生成的格式為.mdl模板文件中。以NuSVM訓(xùn)練為例，采用NuSVM訓(xùn)練后的結(jié)果保存在Nu_svm.mdl中，結(jié)果如圖8所示。

從圖8的訓(xùn)練結(jié)果中可看出，支持向量機(jī)類型為nu_svc，而參數(shù)V值設(shè)置當(dāng)前值為0.6，核函數(shù)當(dāng)前選取的類型為RBF(徑向基核函數(shù))，r值為0.125，采用二分類支持向量機(jī)。一共參與訓(xùn)練的44個(gè)樣本數(shù)據(jù)中，有總共28個(gè)支持向量。隨后可看出正負(fù)支持向量數(shù)目及支持向量的具體節(jié)點(diǎn)值。按照這樣的方式，再進(jìn)行各個(gè)參數(shù)值范圍測試，訓(xùn)練結(jié)束后，已經(jīng)成功建立支持向量機(jī)模型，選取最優(yōu)的支持向量機(jī)模型。采用選取好的20個(gè)測試樣本進(jìn)行分類實(shí)驗(yàn)，測試結(jié)果如圖9所示。

按照這樣的流程，將多種方法進(jìn)行多次測試，最終以算術(shù)平均值作為結(jié)果，對多種方法的比較結(jié)果如表1所示。

支持向量機(jī)平均訓(xùn)練時(shí)間 /s分類平均準(zhǔn)確率 /%

CSVC1.292.0

NuSVM1.395.0

基于多項(xiàng)式核函數(shù)2.285.0

基于RBF核函數(shù)4.298.5



從表1可以看出，采用NuSVM V參數(shù)支持向量機(jī)分類器無論是在訓(xùn)練時(shí)間還是在正判率上，都比采用核函數(shù)的分類器有優(yōu)勢，CSVM分類器適用于經(jīng)過歸一化后的煙氣輪機(jī)原始數(shù)據(jù)。

4 結(jié) 語

本文借鑒和學(xué)習(xí)了國內(nèi)外成熟的相關(guān)理論及實(shí)驗(yàn)方法，通過對小樣本支持向量機(jī)的二分類應(yīng)用于煙氣輪機(jī)振動故障診斷分類的大量實(shí)驗(yàn)。主要包括對CSVM

(下轉(zhuǎn)第51頁)