金 天，丁東方，叢 麗

（北京航空航天大學(xué) 電子信息工程學(xué)院，北京 100191）

隨著導(dǎo)航技術(shù)的發(fā)展，系統(tǒng)的可靠性成為普遍關(guān)注的研究熱點，為了保證系統(tǒng)的可靠性，必須對系統(tǒng)進行實時的故障檢測和隔離（FDI）。70年代開始，很多學(xué)者對FDI進行了研究，其中應(yīng)用較為廣泛的是基于模型的FDI方法[1-4]，該方法通過設(shè)置一定的觀測器來檢測系統(tǒng)是否存在故障，比較有代表意義的是3δ檢驗和殘差χ2檢驗，然而兩種算法對一些大的突變故障具有良好的檢測效果，但對一些小于10δ的小幅值突變故障檢測效果不是太理想[5-6]。為此，鑒于兩種檢測存在的問題，文中提出了一種利用AR模型參數(shù)變化來進行故障診斷的算法，其模型簡單，計算量小，工程中也便于實現(xiàn)，同時它不僅對大的突變故障具有良好的檢測效果，對一些小于10δ的小幅值突變故障也有較好的檢測效果。

1 AR模型參數(shù)估計原理

對于時間序列{Xt}t=1，2，…N，AR（n）模型的表達式為

參數(shù)估計方法就是按照一定的方法估計出 φ1，φ2，…φn，這n個參數(shù)，其中n為AR模型的階數(shù)。at為均值為0，方差為σ2a的白噪聲。AR模型參數(shù)估計一般采用最小二乘算法進行精估計[7-8]。

由于GPS偽距新息數(shù)據(jù)是連續(xù)的時間序列，可以進行AR建模，當(dāng)其連續(xù)的時間序列中某些數(shù)據(jù)發(fā)生故障時，將改變系統(tǒng)的動態(tài)特性從而影響偽距新息序列的特性，使得其不再具有白噪聲特性，均值不再為零，另外其方差陣也將偏離濾波理論值，從而導(dǎo)致AR模型參數(shù)發(fā)生變化，因此可以利用AR模型參數(shù)的變化進行這類系統(tǒng)的故障檢測，即這類系統(tǒng)的故障檢測問題就轉(zhuǎn)化為如何盡可能快地檢測AR模型參數(shù)的改變。

2 基于AR模型參數(shù)變化的導(dǎo)航系統(tǒng)故障檢測模型

基于AR模型參數(shù)變化的導(dǎo)航系統(tǒng)故障檢測模型如圖2所示。

圖1 仿真系統(tǒng)模型框圖Fig.1 The simulation diagram of the system model

其基本思路為：確定AR模型階數(shù)為n，時間序列個數(shù)為m，可得到 n個參數(shù)分別為 φ1，φ2，…φn，。 利用 N 組 GPS 正常工作時輸出的新息序列來訓(xùn)練AR模型，則可得到N組φ1，φ2，…φn。由于N組參數(shù)相互獨立則可通過求出其標(biāo)準(zhǔn)差δ和μ均值，將其整理成標(biāo)準(zhǔn)正態(tài)分布，因此可以用卡方分布檢驗法來進行故障判別。根據(jù)奈曼-皮爾遜準(zhǔn)則，由給定的虛警概率 Pf，就可由 χ2分布求出門限值 TD，然后利用公式：λk=φk-μ2，k=1，…，n 是否大于 TD進行判斷，當(dāng) λ（k）>TD時，系統(tǒng)故障發(fā)生，當(dāng)λ（k）≤TD時，系統(tǒng)工作正常。

3 仿真分析驗證

為了驗證算法性能，利用軟件構(gòu)建了整個仿真系統(tǒng)，著重比較其在實際應(yīng)用中的故障檢測效果。

3.1 AR模型參數(shù)變化

設(shè)定AR模型階數(shù)為5，每組時間序列個數(shù)為60，選取100組GPS偽距新息時間序列利用MATLAB進行仿真，并分別在每組新息數(shù)據(jù)最后2個、10個、30個數(shù)中加入5δ（δ=6 m，為偽距噪聲的標(biāo)準(zhǔn)差）的故障，進行仿真驗證，這里我們分別比較5個參數(shù)的變化。

仿真結(jié)果分別如圖2～圖6所示。

圖2 參數(shù)1變化Fig.2 Parameter 1 change

圖3 參數(shù)2變化Fig.3 Parameter 2 change

從圖2中可以看出當(dāng)每組GPS新息加入故障時，參數(shù)1都有明顯的突變，而其余參數(shù)變化不大，因此可以利用參數(shù)1的變化來進行故障檢測，而且當(dāng)每組新息數(shù)據(jù)只有后2個有故障時，其參數(shù)1就會有明顯的突變，可以看出其對故障相當(dāng)敏感。

圖4 參數(shù)3變化Fig.4 Parameter 3 change

圖5 參數(shù)4變化Fig.5 Parameter 4 change

圖6 參數(shù)5變化Fig.6 Parameter 5 change

3.2 基于AR模型參數(shù)變化的導(dǎo)航系統(tǒng)故障檢測模型仿真

首先選取900組GPS正常輸出的新息序列 （階數(shù)為5，時間序列個數(shù)為60）來訓(xùn)練AR模型，可以獲得900個φ值，從而可以得到其均值u為-0.309 4，標(biāo)準(zhǔn)差δ為0.153 6，設(shè)其虛警概率為10-4，自由度為1，則可求出其門限T0為15.136 7。

本文構(gòu)建的仿真系統(tǒng)利用了GPS/INS組合導(dǎo)航系統(tǒng)中一顆衛(wèi)星輸出的1 000 s GPS偽距新息數(shù)據(jù)進行仿真驗證，并從第61 s開始加入不同幅度的故障。其基本思路為：先利用前60 s數(shù)據(jù)求出參數(shù)1，然后剔除第一個新息加入第61 s輸出的新息數(shù)據(jù)計算一個新的參數(shù)1，以此類推，可以實時得到參數(shù)1，實時得到的參數(shù)1都代入公式求出此刻的 λ1。

仿真結(jié)果分析：由圖7中可以看出當(dāng)在61秒加故障時，在62秒時已經(jīng)可以檢測出故障，但當(dāng)故障延續(xù)57秒后它的參數(shù)會恢復(fù)正常，這時已經(jīng)檢測不到故障。為了避免這種情況，我們需要對其進行改進來確保數(shù)據(jù)的實時檢測和持續(xù)報警。

圖7 故障檢測仿真示意圖Fig.7 The fault detection simulation diagram

4 帶重置的AR模型參數(shù)故障檢測算法

新的設(shè)計實現(xiàn)原理為首先預(yù)存59個新息數(shù)據(jù)，其中58個為GPS正常工作時的新息數(shù)據(jù)，第59個數(shù)據(jù)為故障數(shù)據(jù)，當(dāng)GPS傳感器開始工作時將其實時新息數(shù)據(jù)送入AR模型中進行實時計算，此時時間序列個數(shù)為60，可以得到其參數(shù)，并通過公式進行判斷。如果連續(xù)出現(xiàn)十次故障，則將其數(shù)據(jù)剔除，返回并再次啟用原來預(yù)存的新息數(shù)據(jù)進行實時計算，這樣如果GPS新息出現(xiàn)故障時，它可以實現(xiàn)連續(xù)報警，如果GPS新息恢復(fù)正常，它可以正常工作。其仿真結(jié)果如圖8所示。

圖8 帶重置的AR模型參數(shù)故障檢測仿真示意圖Fig.8 With reset AR model parameters of the fault detection simulation diagram

由圖4中可以看出帶重置的AR模型參數(shù)故障檢測算法能夠在5δ以上檢測到故障，且能夠?qū)崿F(xiàn)實時檢測和持續(xù)報警。

5 結(jié)論

本文提出了一種基于AR模型參數(shù)的故障檢測算法，利用AR模型中參數(shù)1的變化來進行導(dǎo)航系統(tǒng)故障檢測，并利用實際衛(wèi)星偽距新息數(shù)據(jù)進行了驗證。驗證結(jié)果表明其能夠快速檢測到小幅值突變故障，而這些小幅值突變故障利用傳統(tǒng)的殘差卡方檢測是檢測不到的，因此方法在實際應(yīng)用中具有一定意義。

由于仿真驗證是在較理想的條件下實現(xiàn)的，并且對于小于5δ的突變故障其檢驗效果也不夠理想。因此下一步工作要解決如何對小于5δ的故障進行識別和判斷。

[1]DEYSTJ J，DECKERTJ C．Maximum Likelihood Failure Detection Techniques Applied to the Shuttle RCSJet[J]．Journal of Spacecraft and Rockets，1976，13（2）：65-74．

[2]ISERMANN R．Process Fault Detection Based on Modeling and Estimation Methods A Survey[J]．Automatica，1984，20（4）：387-404．

[3]WILSON E．Experiments in Neural Network Control of a Free-flying Space Robot[D]．Stanford：Stanford University，1995．

[4]PATTON R J，F(xiàn)RANK P M，CLARK R N．Issues of Fault Diagnosis for Dynamic System[M]．Berlin：Springer Verlag，2000．

[5]DA R，LIN C F．Failure Detection of Dynamical Systems with the State Chi-square Test[J]．Journal of Guidance，Control and Dynamics，1994，17（2）：271-277．

[6]DA R，LINCF．Sensitivity Analysis Algorithm for the State Chi-squareTest [J]． Journal of Guidance．Control and Dynamics，1996，19（1）：219-222．

[7]董言治，劉松濤.基于Matlab的時間序列分析和動態(tài)數(shù)據(jù)建模[J].計算機工程，2003，29（12）：170～172 DONG Yan-zhi，LIU Song-tao.Time series analysis and dynamic data modeling based on matlab[J].Computer Engineering，2003，29（12）：170-172.

[8]朱增鋒，張鵬杰.AR模型的整體最小二乘解算及其應(yīng)用[J].測繪科學(xué)2012，10（2）；15-21 ZHU Zeng-feng，ZHANG Peng-jie.AR model of least square solution and its application as a whole[J].Science of Surveying and Mapping，2012，10（2）：15-21.

[9]王朋輝，范勝林，劉建業(yè).GPS/SINS超緊組合導(dǎo)航的性能分析[J].導(dǎo)航與控制，2010，9（2）:1-6.WANG Peng-hui，F(xiàn)AN Wei-lin，LIUJian-ye.The performance analysis of GPS/SINS super tight combination navigation[J].Navigation and Control，2010，9（2）:1-6.

[10]尤敏．導(dǎo)航系統(tǒng)地故障診斷技術(shù)研究[D]．北京：北京航空航天大學(xué)，2002．

[11]李春科.改進的狀態(tài)χ2故障檢驗方法及其應(yīng)用[D].西安：西北工業(yè)大學(xué)，2008.