程洪炳 倪世宏 黃國榮 劉華偉

(空軍工程大學(xué)工程學(xué)院，陜西 西安 710038)

0 引言

隨著導(dǎo)航系統(tǒng)種類的增加以及系統(tǒng)結(jié)構(gòu)的日益復(fù)雜化，導(dǎo)航傳感器故障率已成為影響高性能動態(tài)系統(tǒng)完好性的一個重要因素［1－2］。建立在組合導(dǎo)航濾波器之上的系統(tǒng)級檢測也因此成為保證系統(tǒng)輸出數(shù)據(jù)完好性的關(guān)鍵［3］。

在聯(lián)邦濾波器中，既可以基于子濾波器對每個傳感器單獨進(jìn)行故障檢測和隔離，也可以基于主濾波器進(jìn)行整體性能檢測［4－6］。靈活的結(jié)構(gòu)特性使得基于聯(lián)邦濾波器的故障檢測技術(shù)受到了許多學(xué)者的關(guān)注［7］。然而，聯(lián)邦濾波器結(jié)構(gòu)的特殊性卻使目前較為成熟的殘差檢驗算法存在一定風(fēng)險。針對這一問題，本文深入分析了殘差χ2檢驗應(yīng)用聯(lián)邦濾波器時存在的問題，并提出了改進(jìn)算法。

1 殘差故障統(tǒng)計檢驗原理

基于殘差的故障統(tǒng)計檢測在導(dǎo)航系統(tǒng)故障檢測算法中占有絕對優(yōu)勢。無論是觀測量故障還是狀態(tài)量故障，總會通過一定的奇偶關(guān)系反映在殘差變化上?；跉埐畹墓收辖y(tǒng)計檢驗原理就是通過檢驗殘差的變化，檢驗并隔離故障。

帶故障的離散系統(tǒng)模型如式(1)所示:

式中:X(k)為系統(tǒng)狀態(tài)，且 X(k)∈Rn;Φ(k，k－1)為一步狀態(tài)轉(zhuǎn)移陣，且 Φ(k，k－1)∈Rn×n;Γ(k－1)為系統(tǒng)的噪聲陣，且 Γ(k－1)∈Rn×r;H(k)為量測矩陣;w(k)、v(k)為相互獨立的高斯白噪聲序列，且w(k)∈Rr和 v(k)∈Rm;γf為隨機故障向量;f( k，t)f為分段函數(shù)。

f(k，tf)的表達(dá)式如下:

式中:tf為傳感器故障時刻。無故障時，f( k，tf)=0，卡爾曼濾波能夠給出最小方差意義上的最優(yōu)估計;反之，卡爾曼濾波變?yōu)榇蝺?yōu)，甚至失效。

基于殘差的導(dǎo)航系統(tǒng)故障檢測一般可描述如下。令r(k)為k時刻通過殘差生成器(最小二乘法或卡爾曼濾波器)得到的殘差序列，分布函數(shù)為F( k)，且有:

式中:μ描述了故障方向。故障檢驗算法就是要通過假設(shè)檢驗方法，檢驗殘差r(k)的均值由0到μ的變化。根據(jù)μ是否已知，檢測通常分為兩類:①若μ已知，稱為簡單故障檢測;②若μ未知，但已知μ∈Θ:{μTPr－1μ=b2}，稱為復(fù)合故障檢測。

導(dǎo)航系統(tǒng)故障檢測通常被認(rèn)為屬于復(fù)合故障檢測。因為對于一個多維殘差序列，很難事先確定方向和幅度都已知的μ。但簡單故障檢測算法是復(fù)合故障檢測算法的基礎(chǔ)。

2 殘差χ2檢驗算法描述

由于卡爾曼濾波器不受觀測量維數(shù)的限制，又不要求量測矩陣H(k)必須列滿秩。因此，多數(shù)導(dǎo)航系統(tǒng)內(nèi)部均采用卡爾曼濾波生成殘差［8］。預(yù)測殘差序列又稱新息序列，其表達(dá)式為:

根據(jù)卡爾曼濾波方程，新息的協(xié)方差陣為:

殘差χ2檢驗法實際上是一種基于預(yù)測殘差序列的 χ2檢驗法，其中預(yù)報值 X^(k，k－1)為:

對r(k)作如下二元假設(shè)。

H0:無故障 E{ r(k)}=0，E r(k) rT(){ k}=Pr(k)

H1:有故障E{ r(k)}=μ，E{［r(k)－μ］［r(k)－μ］T}=Pr(k)

故障檢測函數(shù)為:

式中:m為量測矩陣Z(k)的維數(shù)。

故障判定準(zhǔn)則為:

3 殘差χ2檢驗問題分析

將殘差χ2檢驗法應(yīng)用于聯(lián)邦濾波器傳感器故障單獨檢驗，局部估計預(yù)報值(k，k－1)一般通過以下兩種方法計算。

①使用局部估計進(jìn)行預(yù)測，即:

②使用全局融合估計進(jìn)行預(yù)測，即:

式中:下標(biāo)i和g分別表示局部濾波和全局濾波。

4 雙“狀態(tài)傳播器”的改進(jìn)

針對殘差χ2檢驗法對軟故障的不敏感性，提出了基于移動偽正常狀態(tài)的殘差χ2檢驗法和基于融合量測預(yù)報的故障檢測方法［9－10］。這兩種算法實際上都是基于“狀態(tài)傳播器”思想，能有效解決上述問題。為避免將一個已受污染的狀態(tài)傳播器作為參考系統(tǒng)，建議采用將雙“狀態(tài)傳播器”思想引入殘差χ2檢驗法，使用全局估計對其進(jìn)行交替重置，具體算法如下。

考慮子濾波器殘差為:

①狀態(tài)傳播器的重置

在系統(tǒng)正常工作狀態(tài)下，兩個狀態(tài)傳播器交替重置。設(shè)交替重置時間間隔為Δt，k0為重置時刻。(k)和PS(k)的計算過程如下。式中:上標(biāo)S、SA、SB分別表示最終狀態(tài)輸出、狀態(tài)傳播器A輸出、狀態(tài)傳播器B輸出;Δt應(yīng)該保證軟故障能夠積累到被檢測的程度，假設(shè)設(shè)置合適，則能夠保證重置時刻全局估計沒有受到污染，那就可以很好地解決軟故障跟蹤問題。

如果系統(tǒng)在tf時刻被檢測出故障，則精度會嚴(yán)重下降，因此，需要同時重置，即:

式中:j=k－Δt+1。

② 時間間隔Δt的選擇。

Δt的選取原則是不能太小也不能太大。

5 仿真分析

考慮如下三維系統(tǒng)［11］:

設(shè)仿真時間為1000s，誤檢率為0.1，查 χ2分布函數(shù)表，可得門限TDi=6.251(自由度為3)。假設(shè)500～800s之間發(fā)生如下形式的軟故障:Z1(k)=Z1( k －1)+0.05(k－500)。直接將殘差χ2檢驗法應(yīng)用于聯(lián)邦卡爾曼濾波器的傳統(tǒng)方法稱為算法一，本文提出改進(jìn)算法標(biāo)為算法二，這兩種算法的故障檢測函數(shù)及漏檢曲線如圖1和圖2所示。

圖1 故障檢測函數(shù)曲線Fig.1 The curves of fault detection function

從圖1可以看出，在故障發(fā)生期間，改進(jìn)算法的故障檢測函數(shù)值比傳統(tǒng)算法普遍要高，這在一定程度上克服了傳統(tǒng)算法對軟故障不敏感的缺陷。

圖2 漏檢曲線Fig.2 The curves of leak out of test

從圖2可以看出，改進(jìn)算法大約從548 s開始漏檢情況很少發(fā)生，而傳統(tǒng)算法漏檢頻發(fā)的狀況一直持續(xù)到656 s左右。這從漏檢風(fēng)險的角度再次證明改進(jìn)算法比傳統(tǒng)算法對軟故障檢測的靈敏度更高。

以上兩種算法的具體檢測風(fēng)險數(shù)據(jù)如表1所示。

表1 檢測風(fēng)險數(shù)據(jù)對比Tab.1 Inter-comparison of the detection risk

6 結(jié)束語

從上述研究可以看出，殘差χ2檢驗法在聯(lián)邦濾波器的系統(tǒng)級故障檢測應(yīng)用過程中，局部估計的次優(yōu)性和自適應(yīng)信息分配策略將對各子濾波器故障檢測靈敏度產(chǎn)生影響。基于雙“狀態(tài)傳播器”的改進(jìn)算法能夠有效改善局部估計次優(yōu)性和自適應(yīng)信息分配策略對故障檢測靈敏度的不一致性，并且仿真結(jié)果也驗證了該方法的有效性和正確性。

［1］Brumback B，Srinath M.A chi-square test for fault-detection in Kalman filters［J］.IEEE Transactions on Automatic Control，1987，32(6):552－554.

［2］Zhu Yanhua，Liu Jianye，Zeng Qinghua.Application of state chi—square test in damping Kalman filter of inertial attitude and heading reference system［J］.Chinese Journal of Scientific Instrument，2007，28(9):1569 －1576.

［3］Wong W E，Wei T T，Qi Y，et al.A crosstab-based statistical method for effective fault localization［C］∥20081stInternational Conference on Software Testing，Varification and Validation，2008:42 －51.

［4］Kehan G，Khoshgoftaar T M.A comprehensive empirical study of count models for software fault prediction［J］.IEEE Transactions on Reliability，2007，56(2):223 －236.

［5］Birsen Y.Statistical pattern analysis of partial discharge measurements for quality assessment of insulation systems in high-voltage electrical machinery［J］.IEEE Transactions on Industry Applications，2004，40(6):1579－1594.

［6］李金梁.基于聯(lián)邦結(jié)構(gòu)的多傳感器自主式組合導(dǎo)航技術(shù)研究［D］.西安:空軍工程大學(xué)，2007.

［7］Kerr T H，Brumback B D，Srinath M D.Comments on“a chi-square test for fault detection in Kalman Filters”［J］.IEEE Transactions on Automatic Control，1990，35(11):1277 －1278.

［8］秦永元，張洪鉞，汪叔華.卡爾曼濾波與組合導(dǎo)航原理［M］.西安:西北工業(yè)大學(xué)出版社，1998.

［9］王志勝.信息融合估計理論及其在航天器控制中的應(yīng)用研究［D］.西安:西北工業(yè)大學(xué)，2002.

［10］邱愷.基于聯(lián)邦結(jié)構(gòu)的組合導(dǎo)航多傳感器信息融合技術(shù)研究［D］.西安:空軍工程大學(xué)，2005.

［11］王仲生，何紅，陳錢.小波分析在發(fā)動機早期故障識別中的應(yīng)用研究［J］.西北工業(yè)大學(xué)學(xué)報:自然科學(xué)版，2006，24(1):68 －71.