徐 曌，張 斌

(桂林電子科技大學(xué) 計算機與信息安全學(xué)院，廣西 桂林 541000)

0 引 言

工業(yè)4.0環(huán)境下，智能設(shè)備通過多種傳感器采集實時參數(shù)產(chǎn)生大量的數(shù)據(jù)。利用數(shù)據(jù)挖掘進行故障診斷就是將狀態(tài)參數(shù)進行關(guān)聯(lián)、分類、推理，得到?jīng)Q策樹[1]，并利用決策樹對智能設(shè)備進行故障診斷。相關(guān)研究表明：工業(yè)機械遇到故障時，技術(shù)人員需要花70%～90%的時間來尋找故障原因和部位，而故障維修只占10%～30%的時間；在維修成本方面，預(yù)測維修成本只占故障后維修成本的40%[2]。

機械設(shè)備結(jié)構(gòu)功能的逐漸增強，使得現(xiàn)階段故障診斷的難度不斷提高。文獻[3]是基于知識庫邏輯理論對規(guī)則進行處理，但沒有給出規(guī)則間的相互關(guān)系，診斷效率不高。文獻[4]是基于知識庫邏輯理論，利用粗糙集直接從故障樣本中得到診斷規(guī)則，效率不高；并且在缺乏屬性間的關(guān)系時，必須對整個規(guī)則庫遍歷搜索，而且當(dāng)缺失某些關(guān)鍵信息時診斷結(jié)果會受到影響。文獻[5]結(jié)合了故障樹與模糊推理用于動態(tài)系統(tǒng)的故障診斷；文獻[6]利用粗糙集研究了故障樹基本事件的排序問題。當(dāng)系統(tǒng)比較復(fù)雜時，這類方法的搜索過程也變得非常復(fù)雜，診斷結(jié)果的準(zhǔn)確性也不高。

為解決上述問題，文中結(jié)合可辨識矩陣與決策樹的特點，提出基于可辨識矩陣約簡和C4.5決策樹的故障診斷方法，并對其進行了驗證。

1 約簡決策樹故障診斷方法

首先確定故障診斷的訓(xùn)練樣本集，接著對訓(xùn)練樣本集里的連續(xù)屬性進行離散化處理，得到離散樣本集。然后用可辨識矩陣約簡算法對離散樣本集進行屬性約簡，并根據(jù)聚類比選擇出最佳約簡矩陣，得到最佳約簡樣本集。再用C4.5算法對最佳約簡樣本集進行分析，生成最終決策樹。最后通過測試樣本集驗證決策樹的故障診斷準(zhǔn)確性。

訓(xùn)練樣本需要全面反映設(shè)備的情況，這樣的樣本經(jīng)過訓(xùn)練后生成的決策才能對新的數(shù)據(jù)做出正確的判斷。以轉(zhuǎn)子磨粉機為例，有轉(zhuǎn)子不平衡、轉(zhuǎn)子不對中、油膜振蕩三種常見機械故障。由歷史振動波形測量數(shù)據(jù)建立如表1所示的樣本訓(xùn)練集T'=。其中，U，C和D是非空有限集。U是故障狀態(tài)域；C=(C1,C2,…,C7)，表示7種故障條件屬性集，C1表示0.01～0.40f、C2表示0.41～0.50f、C3表示0.51～0.99f、C4表示1f、C5表示2f、C6表示3～5f、C7表示大于等于5f；D=(F1，F(xiàn)2，F(xiàn)3)是故障決策屬性，其中F1為不平衡、F2為不對中、F3為油膜振蕩，且C∩D=?；X=(X1,X2,…,X9)，表示9個故障樣本。

表1 訓(xùn)練樣本集T'

2 連續(xù)屬性離散化算法

連續(xù)屬性的離散化是非常必要的，因為粗糙集理論以及多種離散型數(shù)據(jù)挖掘算法只能解決離散數(shù)據(jù)，因此需要對訓(xùn)練樣本集進行離散化處理。離散化方法主要包括兩大類：無監(jiān)督的離散化方法和有監(jiān)督的離散化方法?；陟氐碾x散化方法就是在處理數(shù)據(jù)時考慮了條件屬性和決策屬性之間的依賴關(guān)系。

假設(shè)樣本集有N個故障樣本Xj(j=1,2,…,N)，M個故障條件屬性，K個故障類型。應(yīng)用一種基于最小熵原理的離散算法將樣本集中的連續(xù)屬性Ci進行離散化處理[7]。步驟如下：

步驟2：離散分類數(shù)r=3。

(1)

(2)

(3)

(4)

(5)

步驟5：根據(jù)最大隸屬度原則，將樣本集中的連續(xù)屬性值按照式(6)進行離散化處理，處理過后的條件屬性值記為CdiN(N=1,2…):

(6)

其中，P1i和P2i是連續(xù)屬性離散化過程中的兩個斷點值，其中P1i=(SEC1i+PRIi)/2，P2i=(SEC2i+PRIi)/2，j=1,2,…,N。

將樣本訓(xùn)練集T'中連續(xù)屬性進行離散化處理后得到離散樣本集T''。

3 屬性約簡及最佳約簡優(yōu)選

粗糙集理論[8]可以在缺乏先驗知識的前提下，以分類能力為基礎(chǔ)很好地進行模糊數(shù)據(jù)的分類和處理。粗糙集的關(guān)鍵內(nèi)容是屬性約簡，能夠有效地刪除樣本集中的冗余信息，形成精簡的規(guī)則庫，提高規(guī)則庫的讀取使用效率[9]。文中使用可辨別矩陣作為樣本集屬性的約簡算法。

假設(shè)T''=(U',C',D),其中包括N個經(jīng)過離散化處理的樣本Yj(j=1,2,…,N)，M個條件屬性和K個故障類型?；诳杀鎰e矩陣約簡[10]和最大集群比[11]原則，文中采用如下算法對T''進行約簡并選擇出聚類率最高的屬性組合。

步驟1：按照式(7)計算出T''的可辨識矩陣C=(cij)N×N。

C(i,j)=

(7)

步驟4：輸出屬性約簡結(jié)果。

步驟5：按照式(8)計算步驟5中得到的所有約簡集的聚類比。最佳約簡集即為聚類比Ra最大的一組。

Ra=(N-NR)/N-1

(8)

其中，N為規(guī)則表中離散化的樣本數(shù)目；NR為T''約簡過后的論域個數(shù)。

用步驟1對離散樣本集T''進行處理，得到對應(yīng)的可辨識矩陣，如表2所示。

表2 對應(yīng)于T''的可辨識矩陣

根據(jù)步驟2～4進一步對T''各個屬性進行約簡，可得：{C3,C4,C7}，{C2,C4,C6,C7}，{C1,C2,C4,C7}，{C1,C5,C6,C7}，{C4,C5,C7}，{C1,C2,C4,C5}，{C2,C3,C4,C5,C7}，其中{C3,C4,C7}的聚類率最大，所以把{C3,C4,C7}作為最佳約簡集，從而得到最優(yōu)約簡診斷決策表T'''，如表3所示。

4 C4.5故障決策樹算法

4.1 C4.5算法

C4.5算法[12]是目前最具影響的決策樹算法之一，以信息增益率作為確定分支屬性的標(biāo)準(zhǔn)，更好地處理算法在使用信息增益選擇分支屬性時偏向于取值較多的屬性的缺陷[13]。核心思想是通過分析訓(xùn)練集的數(shù)據(jù)，在數(shù)據(jù)集上遞歸地建立一個決策樹，該決策樹由一組規(guī)劃來表示，尤其是被識別的模式有連續(xù)特征值屬性時，甚至是在缺失某些重要特征屬性時，C4.5算法也是有效的。該算法用已知故障類別的樣本得到順序測試特征，直到所有特征被正確歸類。在得到?jīng)Q策樹之后，只需要判斷部分的條件屬性，即可得到正確的故障類型[14]。C4.5算法使用后剪枝方法簡化決策樹，從而提高多特征下的辨識效果。過程如下：

表3最優(yōu)約簡診斷決策表T'''

UC3C4C7D-X1331F1X2221F1X3331F1X4222F2X5223F2X6322F2X7213F3X8112F3X9312F3

步驟1：計算所有屬性的信息增益率，選擇具有最大信息增益率的屬性作為根節(jié)點的分支屬性。

(1)設(shè)數(shù)據(jù)集為T，|T|作為數(shù)據(jù)集T的樣本數(shù)，共有k個類。|Ci|是T中屬于類別Ci的樣本數(shù)，先驗概率可以表示為Pi=|Ci|/|T|，并且對T分類的信息熵為：

(9)

(10)

(11)

并由式(9)～(11)可以得到按照Ak劃分?jǐn)?shù)據(jù)集的信息增益為：

(12)

其中Ak的信息熵為：

(13)

最終由式(12)、(13)得到Ak的信息增益率為：

(14)

步驟2：不斷利用上述方法建立樹的分支，將信息增益率最大的屬性作為子節(jié)點，直到完成決策樹的創(chuàng)建為止。

4.2 生成故障診斷決策樹

利用C4.5算法對T'''進行分析，計算可得C4的信息增益率為0.8，同理可得C3和C7的信息增益率分別為0.18和0.6。因為C4的信息增益率最大，所以將C4作為根節(jié)點建立決策樹，并且在每個節(jié)點使用C4.5算法進行分析，最后生成的決策樹如圖1所示。使用C4.5算法直接對T'''進行分析，得到的決策樹如圖2所示。并且用離散化測試樣本對最終決策樹和C4.5算法生成的決策樹進行效果驗證。

圖1 約簡的決策樹(文中算法)

圖2 C4.5算法決策樹

5 實驗結(jié)果及分析

(1)決策樹構(gòu)建速度的比較。

5次使用相同訓(xùn)練數(shù)據(jù)集得到?jīng)Q策樹的時間，可以得出文中方法生成決策樹的平均時間比直接使用C4.5算法快約31%。

(2)決策樹復(fù)雜度的比較。

首先，直接使用C4.5算法生成的決策樹(圖2)所需要的故障特征向量為4維，而文中方法(圖1)僅需故障特征向量為2維。證明文中方法在保證的診斷準(zhǔn)確性的同時，有效地減少了故障特征數(shù)據(jù)獲取的工作量，并且降低了決策樹的復(fù)雜度。

(3)故障診斷準(zhǔn)確率的比較。

使用兩個決策樹對相同故障數(shù)據(jù)進行準(zhǔn)確性測試，如圖3所示。其中對于F1故障類型的診斷準(zhǔn)確率基本相同，達到了95.3%。C4.5算法對于F2故障類型診斷正確率為95.3%，而文中方法的診斷準(zhǔn)確率達到了96%。C4.5算法對于F3故障類型的診斷正確率為96.5%，而文中方法的診斷準(zhǔn)確率達到了97.8%。兩種方法的準(zhǔn)確率出現(xiàn)了較為明顯的差異。其次，文中方法可以準(zhǔn)確識別故障類型，說明將粗糙集和決策樹結(jié)合的方法不僅減少了規(guī)則庫搜索的步驟，同時具有很強的容錯能力。

圖3 故障診斷準(zhǔn)確率

綜上可以證明，文中方法不僅在降低最終決策樹復(fù)雜度的基礎(chǔ)上提高了決策樹的構(gòu)建速率，還提高了對部分故障類型的診斷準(zhǔn)確率。

6 結(jié)束語

應(yīng)用可辨識矩陣約簡算法和決策樹結(jié)合的方法能有效地對條件屬性進行約簡，消除冗余信息，并減少決策樹的節(jié)點數(shù)，降低決策樹復(fù)雜度，最終高效率地構(gòu)建決策樹，實現(xiàn)對機械快速準(zhǔn)確的故障識別。然而，數(shù)據(jù)挖掘技術(shù)是建立在海量的數(shù)據(jù)樣本的基礎(chǔ)上[15]，文中的故障樣本數(shù)據(jù)是比較典型的故障類別數(shù)據(jù)，不排除存在具有兩種甚至更多的故障特征數(shù)據(jù)。當(dāng)樣本數(shù)據(jù)更加充足時，決策規(guī)則庫也將有更強的工程實用性。

[1] 陳紹煒,王 聰,趙 帥.決策樹算法在電路故障診斷中的應(yīng)用[J].計算機工程與應(yīng)用,2013,49(12):233-236.

[2] GILLUM M N,ARMIJO C B.Optimizing the frequency of the rotary knife on a roller gin stand[J].Hydrobiologia,2000,43(4):809-817.

[3] 王 慶,巴德純,王曉冬.智能故障診斷的粗糙決策模型[J].東北大學(xué)學(xué)報:自然科學(xué)版,2005,26(1):80-83.

[4] 束洪春,孫向飛,司大軍.基于粗糙集理論的配電網(wǎng)故障診斷研究[J].中國電機工程學(xué)報,2001,21(10):73-77.

[5] CHANG S Y,LIN C R,CHANG C T.A fuzzy diagnosis approach using dynamic fault trees[J].Chemical Engineering Science,2002,57(15):2971-2985.

[6] KHOO L P,TOR S B,LI J R.A rough set approach to the ordering of basic events in a fault tree for fault diagnosis[J].International Journal of Advanced Manufacturing Technology,2001,17(10):769-774.

[7] 賈智皓,劉 方.服務(wù)于粗糙集信息處理的數(shù)據(jù)離散化技術(shù)[J].數(shù)字技術(shù)與應(yīng)用,2014(11):77-79.

[8] 馬文萍,黃媛媛,李 豪,等.基于粗糙集與差分免疫模糊聚類算法的圖像分割[J].軟件學(xué)報,2014,25(11):2675-2689.

[9] 李明祥.基于粗糙集理論的數(shù)據(jù)挖掘方法的研究[D].青島:山東科技大學(xué),2003.

[10] 武志峰,吉根林.一種基于決策矩陣的屬性約簡及規(guī)則提取算法[J].計算機應(yīng)用,2005,25(3):639-642.

[11] 滕書華,魯 敏,楊阿鋒,等.基于一般二元關(guān)系的粗糙集加權(quán)不確定性度量[J].計算機學(xué)報,2014,37(3):649-665.

[12] 徐 鵬,林 森.基于C4.5決策樹的流量分類方法[J].軟件學(xué)報,2009,20(10):2692-2704.

[13] 李學(xué)明,李海瑞,薛 亮,等.基于信息增益與信息熵的TFIDF算法[J].計算機工程,2012,38(8):37-40.

[14] 黃愛輝.決策樹C4.5算法的改進及應(yīng)用[J].科學(xué)技術(shù)與工程,2009,9(1):34-36.

[15] 云玉屏.基于C4.5算法的數(shù)據(jù)挖掘應(yīng)用研究[D].哈爾濱:哈爾濱理工大學(xué),2008.