張東波，汪立新

（火箭軍工程大學(xué)，西安 710025）

隨著我國(guó)航空航天、武器裝備等領(lǐng)域的發(fā)展，火箭、導(dǎo)彈等飛行器裝備數(shù)量、使用頻次急劇增加。因此，如何在降低其生產(chǎn)、維護(hù)成本的同時(shí)提高其健康管理水平是需要重點(diǎn)解決的問題之一[1-3]。

慣導(dǎo)系統(tǒng)作為火箭、導(dǎo)彈等飛行器控制系統(tǒng)的關(guān)鍵部件，其旨在對(duì)飛行器的位置、姿態(tài)等信息進(jìn)行采集，是飛行器實(shí)現(xiàn)時(shí)序飛行、變軌突防的關(guān)鍵。多表冗余慣導(dǎo)系統(tǒng)通過用多個(gè)成本較低的陀螺儀、加速度計(jì)的交叉耦合，以提高慣導(dǎo)系統(tǒng)的精度和可靠度，同時(shí)降低了慣導(dǎo)系統(tǒng)的設(shè)計(jì)成本。

多表冗余慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)旨在根據(jù)多個(gè)陀螺儀、加速度計(jì)的測(cè)試信息對(duì)其開展在線故障檢測(cè)、診斷、定位，并根據(jù)故障診斷結(jié)果，利用剩余可用同類相關(guān)陀螺儀、加速度計(jì)開展故障陀螺儀、加速度計(jì)信息重構(gòu)，以保證慣導(dǎo)系統(tǒng)輸出信息的準(zhǔn)確性[4-8]。

目前多表冗余慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)主要面臨三個(gè)方面的問題：首先，飛行器慣導(dǎo)系統(tǒng)中所使用的陀螺儀、加速度計(jì)等傳感器均屬于軍品級(jí)，篩選過程嚴(yán)格，設(shè)計(jì)可靠性較高，進(jìn)而導(dǎo)致其在使用過程中出現(xiàn)故障的概率較低，長(zhǎng)周期通電狀態(tài)下所獲取的數(shù)據(jù)多數(shù)為慣導(dǎo)系統(tǒng)正常狀態(tài)下的數(shù)據(jù)，故障狀態(tài)下的數(shù)據(jù)較少[9]。其次，多表冗余設(shè)計(jì)下的慣導(dǎo)系統(tǒng)由多個(gè)陀螺儀、加速度計(jì)組成，受其飛行環(huán)境不確定性、不可預(yù)測(cè)性和多飛行狀態(tài)切換、部件安裝誤差等因素的影響，輸出結(jié)果具有一定的不確定性和隨機(jī)性[10]。最后，飛行器在飛行過程中速度較快，慣導(dǎo)系統(tǒng)一旦出現(xiàn)故障需要在最短時(shí)間內(nèi)檢測(cè)、診斷、定位故障，進(jìn)而快速進(jìn)行重構(gòu)容錯(cuò)，降低慣導(dǎo)系統(tǒng)故障給飛行器控制系統(tǒng)所帶來的影響，確保飛行器飛行過程的穩(wěn)定性。因此，在線故障診斷與容錯(cuò)重構(gòu)的實(shí)時(shí)性是需要重點(diǎn)關(guān)注的性能之一。

置信規(guī)則庫(kù)（Belief rule base,BRB）作為專家系統(tǒng)的一種，是基于證據(jù)推理、IF-THEN 規(guī)則、模糊理論等理論發(fā)展而來的。BRB 模型通過專家初始建立、數(shù)據(jù)學(xué)習(xí)訓(xùn)練的方式將知識(shí)與數(shù)據(jù)有效融合，同時(shí)通過證據(jù)推理融合克服知識(shí)不確定性所帶來的影響[11-15]。由于BRB 在小樣本建模中所具有的優(yōu)良性能，其被廣泛應(yīng)用于故障診斷、安全性評(píng)估、醫(yī)學(xué)決策等領(lǐng)域。

針對(duì)多表冗余設(shè)計(jì)方案下慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)過程中所面臨的三個(gè)難題：故障狀態(tài)下數(shù)據(jù)缺乏、專家知識(shí)不確定性和高實(shí)時(shí)性要求，本文基于BRB構(gòu)建了一種故障診斷與重構(gòu)框架，通過融合專家知識(shí)和測(cè)試數(shù)據(jù)最大化利用系統(tǒng)信息，并采用了故障傳感器輸出重構(gòu)策略對(duì)慣導(dǎo)系統(tǒng)內(nèi)部傳感器故障進(jìn)行容錯(cuò)；為了進(jìn)一步提高模型實(shí)時(shí)性，設(shè)計(jì)了一種新的基于輸出誤差影響度自適應(yīng)擴(kuò)散的在線更新信息調(diào)整策略，通過BRB 模型輸出誤差對(duì)故障診斷和重構(gòu)模型在線更新的信息量進(jìn)行反饋調(diào)節(jié)，在使其精度滿足應(yīng)用要求的同時(shí)提高在線建模的實(shí)時(shí)性，有效處理精度與實(shí)時(shí)性之間的沖突；為克服專家不確定性對(duì)模型精度的影響，構(gòu)建了考慮模型參數(shù)物理意義的優(yōu)化模型，通過少量不平衡數(shù)據(jù)對(duì)模型進(jìn)行學(xué)習(xí)，進(jìn)而達(dá)到數(shù)據(jù)與知識(shí)的混合。

1 基于BRB 的多表冗余慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)

本節(jié)對(duì)基于BRB 的慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)模型過程進(jìn)行了推理，如圖1 所示。首先由專家結(jié)合其自身經(jīng)驗(yàn)給定模型的初始參數(shù)，在慣導(dǎo)系統(tǒng)的監(jiān)測(cè)信息可用后，對(duì)慣導(dǎo)系統(tǒng)開展故障診斷，并根據(jù)故障診斷結(jié)果對(duì)其內(nèi)部故障的加速度計(jì)/陀螺儀開展重構(gòu)，進(jìn)而保證慣導(dǎo)系統(tǒng)輸出信息的可靠性。

圖1 慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)框架Fig.1 Fault diagnosis and reconstruction framework of inertial navigation system

1.1 慣導(dǎo)系統(tǒng)故障診斷模型推理

在基于BRB 的慣導(dǎo)系統(tǒng)故障診斷模型中，其輸入信息為慣導(dǎo)系統(tǒng)的測(cè)試信息、輸出為慣導(dǎo)系統(tǒng)中陀螺儀和加速度計(jì)等傳感器的故障情況。慣導(dǎo)系統(tǒng)故障診斷模型推理過程主要分為四步，即信息轉(zhuǎn)換、規(guī)則激活、證據(jù)推理（Evidential reasoning,ER）規(guī)則融合與輸出診斷。

Step 1：多元監(jiān)測(cè)信息轉(zhuǎn)換

由于慣導(dǎo)系統(tǒng)中的陀螺儀、加速度計(jì)等傳感器的監(jiān)測(cè)信息格式、量級(jí)各不相同，無法直接融合使用，需要將多元監(jiān)測(cè)信息轉(zhuǎn)換到統(tǒng)一框架下：

式中，x1(t),x2(t) …xM(t)表示慣導(dǎo)系統(tǒng)的M個(gè)監(jiān)測(cè)指標(biāo)的監(jiān)測(cè)信息，為故障診斷模型中第i個(gè)監(jiān)測(cè)指標(biāo)的輸入信息相對(duì)于第j條規(guī)則匹配度[7-10]。Aik和Ai(k+1)分別為兩條相鄰規(guī)則的參考等級(jí)，主要按照如下三個(gè)原則進(jìn)行確定：涵蓋數(shù)據(jù)變化范圍、間距合理、數(shù)量適中。L為模型中規(guī)則數(shù)量。首先，參考等級(jí)需要涵蓋監(jiān)測(cè)指標(biāo)的波動(dòng)范圍，以防止超出波動(dòng)范圍后數(shù)據(jù)無法轉(zhuǎn)換；其次，參考等級(jí)之間的距離決定了信息轉(zhuǎn)換的精度，距離太近會(huì)導(dǎo)致所需的參考等級(jí)數(shù)量過多，距離太遠(yuǎn)會(huì)導(dǎo)致轉(zhuǎn)換精度降低；最后，參考等級(jí)的數(shù)量直接決定了規(guī)則庫(kù)中的規(guī)則數(shù)量，尤其在監(jiān)測(cè)指標(biāo)多時(shí)，每個(gè)指標(biāo)的參考等級(jí)數(shù)量會(huì)影響模型的可實(shí)施性。

Step 2：規(guī)則激活

在BRB 模型中，不同規(guī)則中所包含的監(jiān)測(cè)指標(biāo)參考等級(jí)不同，在監(jiān)測(cè)信息輸入到故障診斷模型中時(shí)對(duì)于規(guī)則的激活程度也不同。結(jié)合上述信息轉(zhuǎn)換過程中所獲取的監(jiān)測(cè)信息匹配度，可以得到規(guī)則中所有監(jiān)測(cè)指標(biāo)信息相對(duì)于該條規(guī)則匹配度：

式中，αk(t)為t時(shí)刻模型指標(biāo)輸入信息相對(duì)于第k條規(guī)則的總體匹配度；分別為第i條監(jiān)測(cè)指標(biāo)的權(quán)重和相對(duì)權(quán)重；M為指標(biāo)的個(gè)數(shù)。

在得到指標(biāo)輸入信息相對(duì)于每條規(guī)則的匹配度后，可以進(jìn)一步計(jì)算得到該條規(guī)則的激活程度，即激活權(quán)重。

式中，wk(t)為第k條規(guī)則的激活權(quán)重；θk(t)表示第k條規(guī)則在規(guī)則庫(kù)中的權(quán)重，L為模型中規(guī)則的數(shù)量。

Step 3：ER 算法規(guī)則融合

在故障診斷模型中，不同的規(guī)則的激活權(quán)重不同，在整個(gè)規(guī)則庫(kù)中所發(fā)揮的作用也不同，并且每條規(guī)則輸出結(jié)果中的匹配度也不同。因此，需要通過ER 算法對(duì)模型中的規(guī)則進(jìn)行融合，以得到最終各個(gè)故障狀態(tài)的置信度。

ER 算法如下式所示：

式中，βn(t)為融合L條規(guī)則后第n個(gè)系統(tǒng)故障狀態(tài)的置信度；βj,k(t)表示規(guī)則輸出的置信度；μ、Δn(t)和μ(t)分別為ER 算法的中間系數(shù)，無物理意義；N表示輸出結(jié)果中的故障狀態(tài)或參考等級(jí)的個(gè)數(shù)。

Step 4：最鄰近距離故障診斷策略

ER 融合后的輸出結(jié)果為慣導(dǎo)系統(tǒng)各故障狀態(tài)的置信度，即β(t)=[β1(t),β2(t)…βN(t)]。為了進(jìn)一步確定當(dāng)前時(shí)刻慣導(dǎo)系統(tǒng)的故障狀態(tài)，本文采用了基于最鄰近距離的故障診斷策略：

式中，Ε為慣導(dǎo)系統(tǒng)標(biāo)準(zhǔn)故障向量，Eni為標(biāo)準(zhǔn)故障向量的參考值；J為范數(shù)的階次，其具體數(shù)值由專家給定。

根據(jù)故障診斷模型輸出向量與標(biāo)準(zhǔn)故障向量之間的距離，當(dāng)前時(shí)刻慣導(dǎo)系統(tǒng)的故障狀態(tài)為：

其中，En和Ee分別為第n個(gè)和第e個(gè)標(biāo)準(zhǔn)故障向量。

基于上述最臨近距離的計(jì)算公式，可以判斷當(dāng)前時(shí)刻慣導(dǎo)系統(tǒng)的故障狀態(tài)為第e種。本文針對(duì)的慣導(dǎo)系統(tǒng)故障為加速度計(jì)和陀螺儀等傳感器由于長(zhǎng)時(shí)間通電所導(dǎo)致的恒偏差、漂移、性能退化等非致命性故障，并假設(shè)慣導(dǎo)系統(tǒng)同一時(shí)刻僅有一個(gè)傳感器發(fā)生故障，即故障狀態(tài)為加速度計(jì)故障或陀螺儀故障。

1.2 故障傳感器輸出重構(gòu)模型推理

在故障傳感器輸出重構(gòu)模型中，其推理過程與上節(jié)故障診斷模型相似，區(qū)別在于故障傳感器輸出重構(gòu)模型中的輸入指標(biāo)為其相同類型可用傳感器的監(jiān)測(cè)信息，其對(duì)應(yīng)的指標(biāo)參考等級(jí)根據(jù)輸入指標(biāo)確定。輸出重構(gòu)模型規(guī)則輸出中的{H1,H2…HN}需要根據(jù)故障傳感器監(jiān)測(cè)信息確定。

在經(jīng)過ER 算法融合后，模型輸出的結(jié)果為各參考等級(jí)的置信度β(t)。故障傳感器輸出重構(gòu)屬于典型的擬合類輸出，需要進(jìn)一步將置信度向量轉(zhuǎn)換為重構(gòu)輸出：

式中，c f(t)為t時(shí)刻重構(gòu)后故障傳感器的監(jiān)測(cè)信息。u(b)為輸出重構(gòu)模型中第b個(gè)參考等級(jí)效用值，B為參考等級(jí)的個(gè)數(shù)。ηb(t)為重構(gòu)模型輸出第b個(gè)參考等級(jí)的置信度。

經(jīng)過重構(gòu)后的監(jiān)測(cè)信息可以替換故障傳感器的監(jiān)測(cè)信息，進(jìn)而能夠確保飛行器控制系統(tǒng)按照時(shí)序工作指令正常工作。

1.3 模型穩(wěn)定性分析

慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)模型穩(wěn)定表示監(jiān)測(cè)信息受環(huán)境干擾、人為誤操作等原因所導(dǎo)致的無規(guī)律、隨機(jī)波動(dòng)，使得監(jiān)測(cè)信息中存在一定的不確定性的情況下仍然能夠滿足系統(tǒng)要求。模型穩(wěn)定性分析旨在證明輸入信息存在不確定性時(shí)，BRB 模型參數(shù)設(shè)置合理情況下能夠保證輸出誤差趨向于0。假設(shè)慣導(dǎo)系統(tǒng)同一時(shí)刻僅有一個(gè)傳感器發(fā)生故障，僅考慮了加速度計(jì)性能下降這類非致命性故障下的診斷與重構(gòu)，暫未考慮加速度計(jì)、陀螺儀等傳感器本身安裝誤差、標(biāo)度因數(shù)誤差等造成的誤報(bào)。為此，本節(jié)對(duì)基于BRB 所構(gòu)建的慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)模型開展了穩(wěn)定性分析。

BRB 中所使用的ER 算法如式(5)-(7)所示，其輸出為β(t)。BRB 的穩(wěn)定性可以表示為：

因此，通過上述推導(dǎo)過程可以證明，在BRB 模型設(shè)置合理的情況下，基于BRB 所構(gòu)建的慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)框架能夠在輸入存在不確定性的情況下實(shí)現(xiàn)準(zhǔn)確故障診斷與信息重構(gòu)。

2 模型參數(shù)與更新信息量在線更新模型

2.1 模型在線更新信息調(diào)整策略

模型輸出誤差表示當(dāng)前時(shí)刻模型對(duì)于輸入信息與理想結(jié)果之間的非線性建模能力。模型輸出誤差增大表示當(dāng)前時(shí)刻模型參數(shù)不準(zhǔn)確，應(yīng)該進(jìn)一步學(xué)習(xí)調(diào)整；反之，當(dāng)模型誤差降低時(shí)表示當(dāng)前時(shí)刻模型參數(shù)準(zhǔn)確，后續(xù)可以適當(dāng)降低學(xué)習(xí)強(qiáng)度。模型在線更新信息量表示用于模型在線更新的信息量，增大更新信息量能夠提高建模精度，但同時(shí)會(huì)增加模型學(xué)習(xí)訓(xùn)練時(shí)間，降低其實(shí)時(shí)性；反之，模型建模精度降低，實(shí)時(shí)性增加。因此，做好在線更新過程中建模精度和實(shí)時(shí)性的平衡是保證所構(gòu)建的在線故障診斷與重構(gòu)框架性能的關(guān)鍵。基于此，本文設(shè)計(jì)了一種基于輸出誤差影響度自適應(yīng)擴(kuò)散的在線更新信息調(diào)整策略，通過模型輸出誤差影響度的自適應(yīng)擴(kuò)散實(shí)現(xiàn)對(duì)在線更新信息量的在線調(diào)整，平衡建模精度和實(shí)時(shí)性兩大性能指標(biāo)。

故障診斷模型與重構(gòu)模型的輸出結(jié)果可以表示為：

式中，ΛFDD(t)和ΛRE(t)分別為t時(shí)刻故障診斷與重構(gòu)模型的誤差；表示t-1時(shí)刻傳感器的真實(shí)輸出。

基于輸出誤差影響度自適應(yīng)擴(kuò)散的在線更新信息調(diào)整策略通過模型誤差調(diào)整不同階段歷史信息的使用比例，對(duì)在線更新信息量進(jìn)行動(dòng)態(tài)調(diào)整。調(diào)整策略可以表示為：

式中，O(t)為t時(shí)刻模型更新所使用的信息量；Η1,Η2…ΗQ表示按照時(shí)間窗口所劃分的Q個(gè)不同階段歷史信息，τ1,τ2…τQ為所對(duì)應(yīng)的誤差影響度自適應(yīng)因子，用來調(diào)整對(duì)應(yīng)階段的歷史信息使用比例，其初始值通過專家結(jié)合歷史數(shù)量量、模型復(fù)雜度、建模時(shí)實(shí)行要求等給定，進(jìn)而通過在線優(yōu)化進(jìn)行調(diào)整；Λ (t)表示模型輸出誤差。

2.2 模型參數(shù)優(yōu)化

由于BRB 屬于專家系統(tǒng)的一種，所構(gòu)建的故障診斷與重構(gòu)模型中的初始參數(shù)由專家結(jié)合其自身經(jīng)驗(yàn)知識(shí)給定。受專家認(rèn)知能力的影響，初始模型參數(shù)不一定準(zhǔn)確，存在一定的偏差，需要通過所獲得歷史信息進(jìn)行更新。因此，在確定在線更新信息量后，需要利用優(yōu)化算法對(duì)模型參數(shù)開展學(xué)習(xí)調(diào)整。

由于BRB 模型參數(shù)具有物理意義，因此在優(yōu)化過程中需要遵循一定的約束條件。故障診斷與重構(gòu)模型中的待優(yōu)化參數(shù)約束條件為：

式中，δi(t)和σi(t)分別為故障診斷模型和重構(gòu)模型中指標(biāo)的權(quán)重，θk(t)和νk(t)分別為對(duì)應(yīng)兩個(gè)模型中的規(guī)則權(quán)重。τ q、τmax、Q為建模過程中涉及的第q個(gè)專家的權(quán)重、專家權(quán)重的最大值和建模過程中使用專家的人數(shù)。

在優(yōu)化模型中，對(duì)慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)的優(yōu)化目標(biāo)為輸出誤差最小，即故障診斷率最高與重構(gòu)誤差（Mean Square Error,MSE）最小。診斷率與重構(gòu)誤差計(jì)算方法分別表示為：

對(duì)于慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)框架中兩個(gè)模型的優(yōu)化目標(biāo)為：

在本文中，優(yōu)化算法選用了考慮投影算子的協(xié)方差矩陣自適應(yīng)優(yōu)化策略（Projection Covariance Matrix Adaption Evolution Strategy,P-CMA-ES），對(duì)所構(gòu)建的故障診斷與重構(gòu)模型參數(shù)開展學(xué)習(xí)訓(xùn)練[16]。

3 實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證所提方法的有效性，本文基于某型號(hào)多表冗余慣導(dǎo)系統(tǒng)的多個(gè)加速度計(jì)故障診斷與重構(gòu)開展了實(shí)驗(yàn)校驗(yàn)。

3.1 實(shí)驗(yàn)設(shè)置

實(shí)驗(yàn)是在實(shí)驗(yàn)室所搭建的飛行器慣導(dǎo)系統(tǒng)健康管理平臺(tái)上開展。實(shí)驗(yàn)室所在地的初始對(duì)準(zhǔn)重力加速度為9.8014 2m/s2，地球自轉(zhuǎn)角速度設(shè)置為7.292115147 ×10-5rad/s，地球長(zhǎng)半徑設(shè)置為6378245m。慣導(dǎo)系統(tǒng)采樣周期為5ms。實(shí)驗(yàn)過程中加速度計(jì)故障設(shè)置為恒偏差故障和漂移故障。實(shí)驗(yàn)過程選用了某型號(hào)裝備慣導(dǎo)系統(tǒng)中三個(gè)加速度計(jì)，分別命名為A、B和C，實(shí)驗(yàn)過程中通過每9s 的累計(jì)脈沖信息變化量作為模型輸入信息。實(shí)驗(yàn)過程中，通過加速度計(jì)不同的累計(jì)使用時(shí)間來區(qū)分故障與否，即性能穩(wěn)定期的加速度計(jì)和長(zhǎng)時(shí)間使用后超差的加速度計(jì)。實(shí)驗(yàn)共設(shè)置了四種模式，即全部正常、加速度計(jì)A故障、加速度計(jì)B故障和加速度計(jì)C故障。

3.2 故障診斷與重構(gòu)模型構(gòu)建

在本節(jié)實(shí)驗(yàn)中，三個(gè)加速度計(jì)監(jiān)測(cè)數(shù)據(jù)的參考點(diǎn)各選擇了4 個(gè)，即高（H）、稍高（SH）、中（M）和低（L），根據(jù)歷史數(shù)據(jù)的分布可以進(jìn)一步確定三個(gè)加速度計(jì)的所對(duì)應(yīng)的參考等級(jí)，如表1 所示，表中的數(shù)據(jù)為累計(jì)脈沖量，無量綱。在故障診斷模型中，輸入信息為三個(gè)加速度計(jì)的監(jiān)測(cè)數(shù)據(jù)，輸出為加速度計(jì)的故障情況；在故障傳感器輸出重構(gòu)模型中，輸入信息為剩余兩個(gè)傳感器的監(jiān)測(cè)數(shù)據(jù)，輸出為重構(gòu)后的故障傳感器輸出。實(shí)驗(yàn)中，標(biāo)準(zhǔn)故障向量設(shè)置為：

表1 加速度計(jì)參考等級(jí)和參考值Tab.1 Reference points and value of the accelerometer

在本節(jié)實(shí)驗(yàn)中，J=2。

因此，根據(jù)規(guī)則的構(gòu)造方式，實(shí)驗(yàn)中的故障診斷模型和故障加速度計(jì)輸出重構(gòu)模型分別有64 條和16條規(guī)則。初始故障診斷模型和三個(gè)加速度計(jì)故障下的輸出重構(gòu)模型的初始參數(shù)由專家結(jié)合自身經(jīng)驗(yàn)給定，即根據(jù)規(guī)則中加速度計(jì)監(jiān)測(cè)數(shù)據(jù)參考等級(jí)的組合判斷當(dāng)前時(shí)刻加速度計(jì)的故障狀態(tài)。

3.3 模型訓(xùn)練與測(cè)試

實(shí)驗(yàn)過程中共獲取了800 組數(shù)據(jù)，其中隨機(jī)抽取400 組作為模型的訓(xùn)練數(shù)據(jù)，剩余的400 組數(shù)據(jù)作為模型的測(cè)試數(shù)據(jù)。在訓(xùn)練環(huán)節(jié)，全部的400 組數(shù)據(jù)用來對(duì)模型參數(shù)進(jìn)行訓(xùn)練；在測(cè)試環(huán)節(jié)，利用本文方法開展在線訓(xùn)練，并通過在線更新信息調(diào)整策略對(duì)在線更新信息量進(jìn)行調(diào)整。訓(xùn)練過程中，將400 組數(shù)據(jù)平均分成4 個(gè)時(shí)間窗口，其中誤差影響度自適應(yīng)因子初始設(shè)置為τ1=1，τ1=0.5，τ1=0.2和τ1=0.1。P-CAM-ES 優(yōu)化算法中訓(xùn)練迭代次數(shù)設(shè)置為200 次。優(yōu)化后的故障診斷模型和三個(gè)加速度計(jì)的輸出重構(gòu)模型分別如表2-5 所示。

表2 訓(xùn)練后的故障診斷模型Tab.2 The optimized fault diagnosis model

表3 加速度計(jì)A 輸出重構(gòu)模型Tab.3 Output reconstruction model for accelerometer A

表4 加速度計(jì)B 輸出重構(gòu)模型Tab.4 Output reconstruction model for accelerometer B

表5 加速度計(jì)C 輸出重構(gòu)模型Tab.5 Output reconstruction model for accelerometer C

訓(xùn)練后模型的故障診斷率基于式(32)計(jì)算為96.75%，輸出結(jié)果如圖2 所示。從圖中可以看出，經(jīng)過訓(xùn)練后的模型能夠準(zhǔn)確地對(duì)加速度計(jì)故障進(jìn)行診斷，誤診斷的點(diǎn)較少。實(shí)驗(yàn)重復(fù)20 次，每次實(shí)驗(yàn)過程中不同時(shí)刻的在線更新信息量和隨機(jī)4 次實(shí)驗(yàn)的全過程信息量變化如圖3 所示，其中訓(xùn)練數(shù)據(jù)量上限設(shè)置為180 組。

圖2 加速度計(jì)故障診斷結(jié)果Fig.2 Fault diagnosis output of accelerometers

圖3 在線更新信息量變化規(guī)律Fig.3 Online update information volume change law

從圖3(a)中可以看出，在線更新信息量基本上分為了5 個(gè)階段，不同階段占比不同，按照調(diào)節(jié)因子的擴(kuò)散程度進(jìn)行階段性調(diào)整。另外，縱向?qū)Ρ龋ㄈ鐖D3(b)）可以發(fā)現(xiàn)，訓(xùn)練初期對(duì)于數(shù)據(jù)量要求比較大，證明初期誤差變化較大；訓(xùn)練后期模型參數(shù)誤差進(jìn)一步降低，進(jìn)而對(duì)于在線更新數(shù)據(jù)量的需求有所降低。為了進(jìn)一步驗(yàn)證所提方法的有效性，與傳統(tǒng)的BRB 模型、模糊模型、神經(jīng)網(wǎng)絡(luò)模型開展了對(duì)比實(shí)驗(yàn)，其診斷率分別為97.00%、78.65%和70.20%?？梢钥闯?，相對(duì)于傳統(tǒng)BRB 模型，本文的建模方法精度與其基本一致，但在線更新所使用的信息量有了很大程度的降低；同時(shí)，對(duì)比其他兩種典型的知識(shí)驅(qū)動(dòng)和數(shù)據(jù)驅(qū)動(dòng)方法，建模精度分別提高了23.01%和37.82%，進(jìn)而驗(yàn)證了所提方法的有效性。

重構(gòu)模型的輸出如圖4 所示，從中可以看出，對(duì)于單個(gè)加速度計(jì)故障的情況下，所構(gòu)建的模型可以對(duì)其進(jìn)行準(zhǔn)確重構(gòu)在理想?yún)^(qū)域，誤差較大的點(diǎn)較少。

圖4 三個(gè)加速度計(jì)故障下的重構(gòu)輸出Fig.4 Reconstructed output of three failure accelerometers

為了進(jìn)一步驗(yàn)證重構(gòu)信息的準(zhǔn)確性，本文通過平臺(tái)開發(fā)的慣導(dǎo)系統(tǒng)導(dǎo)航精度計(jì)算軟件，利用重構(gòu)后的加速度計(jì)信息與故障狀態(tài)的加速度計(jì)信息分別計(jì)算了所對(duì)應(yīng)的導(dǎo)航精度，如圖5 所示。從圖中可以看出，經(jīng)過重構(gòu)后的加速度計(jì)監(jiān)測(cè)信息使飛行器的導(dǎo)航精度較正常偏大，但基本與正常精度一致，保證了加速度計(jì)故障下飛行器的飛行穩(wěn)定，驗(yàn)證了本文方法的有效性。

圖5 故障加速度計(jì)重構(gòu)后的導(dǎo)航精度Fig.5 Navigation accuracy after reconstruction of the failure accelerometers

4 結(jié)論

本文針對(duì)飛行器多表冗余慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)過程中所面臨的故障狀態(tài)樣本缺乏、專家知識(shí)不確定和高實(shí)時(shí)性要求等難題，建立了基于BRB 多表冗余慣導(dǎo)系統(tǒng)故障診斷與重構(gòu)框架，其中包含故障診斷模型和故障傳感器輸出重構(gòu)模型，提出了一種新的基于輸出誤差影響度自適應(yīng)擴(kuò)散的在線更新信息調(diào)整策略，通過模型輸出誤差實(shí)現(xiàn)在線更新信息量的自適應(yīng)動(dòng)態(tài)調(diào)整，能夠同時(shí)兼顧建模過程中的實(shí)時(shí)性與精度。

為驗(yàn)證本文方法的有效性，在某型號(hào)慣導(dǎo)系統(tǒng)加速度計(jì)中開展了應(yīng)用研究，實(shí)驗(yàn)結(jié)果表明故障診斷精度為96.75%，同時(shí)所使用的數(shù)據(jù)量能夠根據(jù)建模精度進(jìn)行實(shí)時(shí)調(diào)整。