專利名稱：一種機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種機(jī)械故障預(yù)測方法，特別是關(guān)于一種基于信息熵加權(quán)和時間因素 加權(quán)的機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法。
背景技術(shù)：
面對機(jī)電設(shè)備故障預(yù)測的非線性、非平穩(wěn)動態(tài)問題，傳統(tǒng)的線性化處理方法效果 欠佳，而一些本質(zhì)是非線性的動態(tài)系統(tǒng)分析處理方法在故障預(yù)測方面具有應(yīng)用前景。神經(jīng) 網(wǎng)絡(luò)預(yù)測方法具有自學(xué)習(xí)功能，以及非線性、非局域性、非定常性等特點(diǎn)，可以通過恰當(dāng)選 擇網(wǎng)絡(luò)層次和隱層單元數(shù)，能夠以任意精度逼近任意連續(xù)非線性函數(shù)及其各階導(dǎo)數(shù)的特 性，因而在故障預(yù)測中被廣泛。目前，采用神經(jīng)網(wǎng)絡(luò)進(jìn)行故障預(yù)測的方法大體是將設(shè)備運(yùn) 行狀態(tài)的時間序列數(shù)據(jù)依次輸入到神經(jīng)網(wǎng)絡(luò)的輸入層，然后采用通用神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練和 預(yù)測。在這種方法中，網(wǎng)絡(luò)輸入值對網(wǎng)絡(luò)輸出預(yù)測值的貢獻(xiàn)程度在概率上基本上是等同的， 而且所用的神經(jīng)網(wǎng)絡(luò)大多是靜態(tài)網(wǎng)絡(luò)，不適合動態(tài)系統(tǒng)的實(shí)時辨識。在實(shí)際應(yīng)用中，網(wǎng)絡(luò)結(jié) 構(gòu)中隱含層節(jié)點(diǎn)數(shù)一般是依靠試算的辦法，但是該方法計算量較大，并且不容易確定所得 預(yù)測模型結(jié)構(gòu)的優(yōu)劣。

發(fā)明內(nèi)容
針對上述問題，本發(fā)明的目的是提供一種故障預(yù)測可靠性較高、計算量較小且能 適用于動態(tài)系統(tǒng)實(shí)時辨識的機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法。為實(shí)現(xiàn)上述目的，本發(fā)明采取以下技術(shù)方案一種機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù) 測方法，其步驟如下(1)通過遠(yuǎn)程監(jiān)測診斷中心獲取對故障敏感的某測點(diǎn)傳感器輸出的 一段連續(xù)振動信號；(2)利用3o方法和插值方法對獲取的振動數(shù)據(jù)分別進(jìn)行異常值剔除 和缺失數(shù)據(jù)填補(bǔ)，得到振動數(shù)據(jù)序列{X”…，Xn} ； (4)對振動數(shù)據(jù)序列{X”…，XJ進(jìn)行歸 一化處理后得到序列{Xl，-,xj ； (5)根據(jù)歸一化處理后的序列{Xl，…，xn}，計算基于信 息熵加權(quán)系數(shù)wSi，進(jìn)而計算出經(jīng)熵加權(quán)后的振動數(shù)據(jù)序列{yi，…，yn} ； (6)由于時間因素 影響，利用時間加權(quán)對熵加權(quán)后的振動數(shù)據(jù)序列{y”…，yj進(jìn)行時間加權(quán)計算后，得到數(shù) 據(jù)序列Iz”…，zj ； (7)利用數(shù)據(jù)序列{Zl，…，zj構(gòu)建非線性的動態(tài)遞歸神經(jīng)網(wǎng)絡(luò)預(yù)測 模型，并采用黃金分割法確定隱含層最優(yōu)節(jié)點(diǎn)數(shù)，進(jìn)而確定神經(jīng)網(wǎng)絡(luò)最佳預(yù)測模型結(jié)構(gòu)，進(jìn) 行故障趨勢預(yù)測，得到預(yù)測結(jié)果{+，…，uj，m為預(yù)測點(diǎn)的數(shù)目；(8)對趨勢預(yù)測結(jié)果進(jìn)行 反歸一化處理，得到實(shí)際預(yù)測結(jié)果。所述步驟(5)中，所述振動序列經(jīng)所述熵加權(quán)后得到序列{yi，…，yj為yi
=xiWsi i = 1,2,…n，式中，信息熵加權(quán)系數(shù)其中Ei為反映振動數(shù)據(jù)攜
帶信息量的信息熵，E, = -Pilog2Pi(i = 1,2,…，n), Pi為各振動數(shù)據(jù)提供信息的概率， Pi=^L~ 二 1,2, .., )。
所述步驟(6)中，所述振動序列經(jīng)所述時間加權(quán)后得到數(shù)據(jù)序列{Zl，…，zj為<formula>formula see original document page 5</formula>
zj = yiwNi i = 1,2,…n，式中，wNi為時間加權(quán)系數(shù)，其為:WNi =‘
，其中<formula>formula see original document page 5</formula>所述步驟(7)中，所述隱含層最優(yōu)節(jié)點(diǎn)數(shù)的確定方法如下①根據(jù)公式ni = log2n 和 二^/^^ +僅分別確定隱含層節(jié)點(diǎn)數(shù)所在的區(qū)間的最小值n_min和最大值n_maX，其 中，n為輸入神經(jīng)元數(shù)，m為輸出神經(jīng)元數(shù)，a為[1，10]之間的常數(shù)；②計算隱含層節(jié)點(diǎn) 數(shù)所在區(qū)間端點(diǎn)處的誤差均方E(n_min)和E(n_maX)；③比較端點(diǎn)處誤差的大小，當(dāng)E(n_ min) ^E(n_max)時，從隱含層節(jié)點(diǎn)數(shù)所在區(qū)間最小值的右側(cè)進(jìn)行黃金搜索，搜索點(diǎn)i = n_ min+0.618(n_max-n_min)，則設(shè)置njnax = i ；反之，從隱含層節(jié)點(diǎn)數(shù)所在區(qū)間最大值的左 側(cè)進(jìn)行黃金搜索，搜索點(diǎn)i = n_max-0. 618 (n_max-n_min)，則設(shè)置n_min = i ；④根據(jù)已設(shè) 定的均方誤差下限值判斷黃金點(diǎn)是否滿足誤差要求，當(dāng)滿足誤差要求時，最佳隱含層節(jié)點(diǎn) 數(shù)即為搜索點(diǎn)i ；否則，進(jìn)入步驟②循環(huán)。本發(fā)明由于采取以上技術(shù)方案，其具有以下優(yōu)點(diǎn)1、本發(fā)明由于采用動態(tài)神經(jīng)網(wǎng) 絡(luò)模型進(jìn)行預(yù)測，利用信息熵對所輸入的表征設(shè)備運(yùn)行的振動信號，進(jìn)行信息融合以獲得 對機(jī)電設(shè)備運(yùn)行狀態(tài)的一致性描述，建立基于信息熵的權(quán)矩陣，同時考慮時間因素對網(wǎng)絡(luò) 輸入的影響，建立新息加權(quán)的動態(tài)神經(jīng)網(wǎng)絡(luò)預(yù)測模型，在動態(tài)神經(jīng)網(wǎng)絡(luò)的預(yù)測中，采用黃金 分割法確定隱含層節(jié)點(diǎn)數(shù)的數(shù)目，可以高效地預(yù)測設(shè)備運(yùn)行狀態(tài)，因此提高了故障預(yù)測可 靠性。2、本發(fā)明由于對來自傳感器的原始信息進(jìn)行信息熵融合，得到的信息熵加權(quán)系數(shù)，并 按照信息在信息熵上信息量的貢獻(xiàn)程度進(jìn)行加權(quán)，因此可以充分有效地利用傳感器信息進(jìn) 行預(yù)測。3、本發(fā)明由于考慮時間因素影響，進(jìn)而建立新息加權(quán)矩陣，因此確定了網(wǎng)絡(luò)輸入值 對網(wǎng)絡(luò)輸出預(yù)測值的貢獻(xiàn)程度。4、本發(fā)明由于采用黃金分割法搜尋確定隱含層節(jié)點(diǎn)數(shù)，因 此大大簡化了計算量，并且容易確定所得預(yù)測模型結(jié)構(gòu)的優(yōu)劣。5、本發(fā)明由于利用本質(zhì)為 非線性拓?fù)浣Y(jié)構(gòu)的動態(tài)神經(jīng)網(wǎng)絡(luò)預(yù)測模型，因此使得預(yù)測模型具有動態(tài)自適應(yīng)特點(diǎn)，能夠 適應(yīng)工況條件和環(huán)境變化等，實(shí)現(xiàn)了非線性故障預(yù)測。本發(fā)明可以廣泛應(yīng)用于各種機(jī)電設(shè) 備故障預(yù)測分析中。


圖1是本發(fā)明的整體流程示意圖，圖2是本發(fā)明的基于黃金分割法確定隱含層最優(yōu)節(jié)點(diǎn)數(shù)流程示意圖。
具體實(shí)施例方式下面結(jié)合附圖和實(shí)施例對本發(fā)明進(jìn)行詳細(xì)的描述。如圖1所示，本發(fā)明采用動態(tài)神經(jīng)網(wǎng)絡(luò)模型進(jìn)行預(yù)測，利用信息熵加權(quán)對輸入的 代表設(shè)備運(yùn)行的振動信號進(jìn)行信息融合，進(jìn)而獲得對機(jī)電設(shè)備運(yùn)行狀態(tài)的一致性描述。建立基于信息熵加權(quán)的同時，考慮時間因素對網(wǎng)絡(luò)輸入的影響，建立新息加權(quán)的動態(tài)神經(jīng)網(wǎng) 絡(luò)預(yù)測模型，在動態(tài)神經(jīng)網(wǎng)絡(luò)的預(yù)測中，采用黃金分割法確定隱含層最優(yōu)節(jié)點(diǎn)數(shù)，以實(shí)現(xiàn)高 效的預(yù)測設(shè)備運(yùn)行狀態(tài)。其具體步驟如下(1)通過現(xiàn)有遠(yuǎn)程監(jiān)測診斷中心獲取對故障敏感的某測點(diǎn)傳感器輸出的一段連續(xù) 振動信號數(shù)據(jù)；(2)利用3 σ方法和插值方法對獲取的振動數(shù)據(jù)分別進(jìn)行異常值剔除和缺失數(shù)據(jù) 填補(bǔ)，得到振動數(shù)據(jù)序列{X” -,XJ ；(4)對振動數(shù)據(jù)序列{Χ1;…，XJ進(jìn)行歸一化處理后得到序列Ix1,…，Xj，即
<formula>formula see original document page 6</formula>(5)根據(jù)歸一化處理后的序列Ix1,…，xn}，計算基于信息熵加權(quán)系數(shù)Wsi，經(jīng)過熵 加權(quán)后得到振動數(shù)據(jù)序列…，yj為<formula>formula see original document page 6</formula>上式中，信息熵加權(quán)系數(shù)=，其中Ei為反映振動數(shù)據(jù)攜帶信息量的信
息熵，Ei = -PiIog2Pi (i = 1,2,…，η)，Pi為各振動數(shù)據(jù)提供信息的概率， <formula>formula see original document page 6</formula>
(6)由于時間序列數(shù)據(jù)隨振動數(shù)據(jù)的新舊程度不同，會對輸出預(yù)測值的貢獻(xiàn)大 小也不相同，因此用時間加權(quán)對熵加權(quán)后的振動數(shù)據(jù)序列…，yn}進(jìn)行時間加權(quán)計 算，進(jìn)而能體現(xiàn)出樣本數(shù)據(jù)中新舊數(shù)據(jù)對預(yù)測值的貢獻(xiàn)大小，經(jīng)時間加權(quán)后得到數(shù)據(jù)序列 Iz1,…，zn}為<formula>formula see original document page 6</formula>η (3)上式中，wNi為時間加權(quán)系數(shù)，其為
<formula>formula see original document page 6</formula>Wm 吐,，其中 kn hi
<formula>formula see original document page 6</formula>(7)利用數(shù)據(jù)序列Iz1,…，zn}構(gòu)建非線性的動態(tài)遞歸神經(jīng)網(wǎng)絡(luò)預(yù)測模型，并在動 態(tài)神經(jīng)網(wǎng)絡(luò)預(yù)測模型中，采用黃金分割法確定隱含層最優(yōu)節(jié)點(diǎn)數(shù)，進(jìn)而確定神經(jīng)網(wǎng)絡(luò)最佳 預(yù)測模型結(jié)構(gòu)，進(jìn)行故障趨勢預(yù)測，得到預(yù)測結(jié)果Iu1,…，Uj，m為預(yù)測點(diǎn)的數(shù)目；(8)對趨勢預(yù)測結(jié)果進(jìn)行反歸一化處理，得到實(shí)際預(yù)測結(jié)果Iv1,…，vm}，預(yù)測結(jié)果 Vi為(4)如圖2所示，上述步驟(7)中，隱含層最優(yōu)節(jié)點(diǎn)數(shù)的確定方法如下①根據(jù)經(jīng)驗(yàn)公式叫=Iog2Ii確定隱含層節(jié)點(diǎn)數(shù)所在的區(qū)間的最小值n_min，根據(jù) 經(jīng)驗(yàn)公式％ = V^ + α確定隱含層節(jié)點(diǎn)數(shù)所在的區(qū)間的最大值n_maX，進(jìn)而確定了隱含層 節(jié)點(diǎn)數(shù)所在區(qū)間為[n_min，n_maX]；其中，η為輸入神經(jīng)元數(shù)，m為輸出神經(jīng)元數(shù)，α為[1，10]之間的常數(shù)；②計算隱含層節(jié)點(diǎn)數(shù)所在區(qū)間端點(diǎn)處的均方誤差E (n_min)和E (n_max)；③比較端點(diǎn)處誤差的大小，當(dāng)E (n_min)彡E (n_max)時，從隱含層節(jié)點(diǎn)數(shù)所在區(qū)間 最小值的右側(cè)進(jìn)行黃金搜索，搜索點(diǎn)i = n_min+0. 618 (n_max-n_min)，則設(shè)置njnax = i ； 否則，從隱含層節(jié)點(diǎn)數(shù)所在區(qū)間最大值的左側(cè)進(jìn)行黃金搜索，搜索點(diǎn)i = n_max-0.618(n_ max-n_min)，貝U設(shè)置 n—min = i ；④根據(jù)已設(shè)定的均方誤差下限值判斷黃金點(diǎn)是否滿足誤差要求，當(dāng)滿足誤差要求 時，最佳隱含層節(jié)點(diǎn)數(shù)即為搜索點(diǎn)i ；否則，進(jìn)入步驟②循環(huán)。上述各實(shí)施例僅是本發(fā)明的優(yōu)選實(shí)施方式，在本技術(shù)領(lǐng)域內(nèi)，凡是基于本發(fā)明技 術(shù)方案上的變化和改進(jìn)，不應(yīng)排除在本發(fā)明的保護(hù)范圍之外。
權(quán)利要求
一種機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法，其步驟如下(1)通過遠(yuǎn)程監(jiān)測診斷中心獲取對故障敏感的某測點(diǎn)傳感器輸出的一段連續(xù)振動信號；(2)利用3σ方法和插值方法對獲取的振動數(shù)據(jù)分別進(jìn)行異常值剔除和缺失數(shù)據(jù)填補(bǔ)，得到振動數(shù)據(jù)序列{X1，…，Xn}；(4)對振動數(shù)據(jù)序列{X1，…，Xn}進(jìn)行歸一化處理后得到序列{x1，…，xn}；(5)根據(jù)歸一化處理后的序列{x1，…，xn}，計算基于信息熵加權(quán)系數(shù)wSi，進(jìn)而計算出經(jīng)熵加權(quán)后的振動數(shù)據(jù)序列{y1，…，yn}；(6)由于時間因素影響，利用時間加權(quán)對熵加權(quán)后的振動數(shù)據(jù)序列{y1，…，yn}進(jìn)行時間加權(quán)計算后，得到數(shù)據(jù)序列{z1，…，zn}；(7)利用數(shù)據(jù)序列{z1，…，zn}構(gòu)建非線性的動態(tài)遞歸神經(jīng)網(wǎng)絡(luò)預(yù)測模型，并采用黃金分割法確定隱含層最優(yōu)節(jié)點(diǎn)數(shù)，進(jìn)而確定神經(jīng)網(wǎng)絡(luò)最佳預(yù)測模型結(jié)構(gòu)，進(jìn)行故障趨勢預(yù)測，得到預(yù)測結(jié)果{u1，…，um}，m為預(yù)測點(diǎn)的數(shù)目；(8)對趨勢預(yù)測結(jié)果進(jìn)行反歸一化處理，得到實(shí)際預(yù)測結(jié)果。
2.如權(quán)利要求1所述的一種機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法，其特征在于所述 步驟(5)中，所述振動序列經(jīng)所述熵加權(quán)后得到序列{y1;…，yn}為Yi = XiWsi i = 1,2, ".η,式中，信息熵加權(quán)系數(shù)二 }，其中Ei為反映振動數(shù)據(jù)攜帶信息量的信息熵，Ei="PiIog2Pi (i = 1,2,…，η)，Pi為各振動數(shù)據(jù)提供信息的概率， (ζ. = 1，2，."，《)0/=1
3.如權(quán)利要求1所述的一種機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法，其特征在于所述 步驟(6)中，所述振動序列經(jīng)所述時間加權(quán)后得到數(shù)據(jù)序列Iz1,…，ζη}為Zi = Y1Wni i = 1,2,…η, 式中，wNi為時間加權(quán)系數(shù)，其為<formula>formula see original document page 2</formula>，其中 <formula>formula see original document page 2</formula>
4.如權(quán)利要求1所述的一種機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法，其特征在于所述 步驟(7)中，所述隱含層最優(yōu)節(jié)點(diǎn)數(shù)的確定方法如下①根據(jù)公式Il1= Iog2Ii和巧=V^ + 分別確定隱含層節(jié)點(diǎn)數(shù)所在的區(qū)間的最小值 n_min和最大值n_maX，其中，η為輸入神經(jīng)元數(shù)，m為輸出神經(jīng)元數(shù)，α為[1，10]之間的常 數(shù)；②計算隱含層節(jié)點(diǎn)數(shù)所在區(qū)間端點(diǎn)處的誤差均方E(n_min)和E(n_maX)；③比較端點(diǎn)處誤差的大小，當(dāng)E(n_min)^E(n_max)時，從隱含層節(jié)點(diǎn)數(shù)所在區(qū)間最小 值的右側(cè)進(jìn)行黃金搜索，搜索點(diǎn)i = n_min+0. 618 (n_max-n_min)，則設(shè)置njnax = i ；反之，從隱含層節(jié)點(diǎn)數(shù)所在區(qū)間最大值的左側(cè)進(jìn)行黃金搜索，搜索點(diǎn)i = n_max-0. 618(n_max-n_ min),貝設(shè)置 n_min = i ；④根據(jù)已設(shè)定的均方誤差下限值判斷黃金點(diǎn)是否滿足誤差要求，當(dāng)滿足誤差要求時， 最佳隱含層節(jié)點(diǎn)數(shù)即為搜索點(diǎn)i ；否則，進(jìn)入步驟②循環(huán)。
全文摘要
本發(fā)明涉及一種機(jī)電設(shè)備神經(jīng)網(wǎng)絡(luò)故障趨勢預(yù)測方法，(1)獲取對故障敏感的某測點(diǎn)傳感器輸出的一段連續(xù)振動信號；(2)利用3σ方法和插值方法對振動數(shù)據(jù)分別進(jìn)行異常值剔除和缺失數(shù)據(jù)填補(bǔ)；(3)對振動數(shù)據(jù)序列進(jìn)行歸一化處理；(4)根據(jù)歸一化處理后的序列，計算出經(jīng)熵加權(quán)后的振動數(shù)據(jù)序列；(5)由于時間因素影響，利用時間加權(quán)對熵加權(quán)后的振動數(shù)據(jù)序列進(jìn)行時間加權(quán)計算；(6)利用由步驟(6)得到的數(shù)據(jù)序列構(gòu)建非線性的動態(tài)遞歸神經(jīng)網(wǎng)絡(luò)預(yù)測模型，并采用黃金分割法確定隱含層最優(yōu)節(jié)點(diǎn)數(shù)；(7)對趨勢預(yù)測結(jié)果進(jìn)行反歸一化處理，得到實(shí)際預(yù)測結(jié)果。本發(fā)明由于采用動態(tài)神經(jīng)網(wǎng)絡(luò)模型進(jìn)行預(yù)測，提高了故障預(yù)測可靠性。本發(fā)明可以廣泛應(yīng)用于各種機(jī)電設(shè)備故障預(yù)測分析中。
文檔編號G01M99/00GK101799367SQ201010101328
公開日2010年8月11日 申請日期2010年1月27日 優(yōu)先權(quán)日2010年1月27日
發(fā)明者徐小力, 王少紅, 王紅軍, 陳濤 申請人:北京信息科技大學(xué)