朱立新，孟 ，王 江，馬春來(lái)

（1.解放軍電子工程學(xué)院，合肥 230037；2.電子制約技術(shù)安徽省重點(diǎn)實(shí)驗(yàn)室，合肥 230037）

局部狀態(tài)檢測(cè)的GPS/INS組合導(dǎo)航轉(zhuǎn)發(fā)式干擾檢驗(yàn)方法

朱立新1，2，孟1，2，王 江1，馬春來(lái)1，2

（1.解放軍電子工程學(xué)院，合肥 230037；2.電子制約技術(shù)安徽省重點(diǎn)實(shí)驗(yàn)室，合肥 230037）

GPS/INS組合導(dǎo)航，轉(zhuǎn)發(fā)式干擾，局部狀態(tài)，χ2檢驗(yàn)

引言

GPS/INS（Global Positioning System/Inertial Naviga-tion System）組合導(dǎo)航系統(tǒng)克服了各自缺點(diǎn)，使組合后的導(dǎo)航精度高于兩個(gè)系統(tǒng)單獨(dú)工作的精度［1］。GPS轉(zhuǎn)發(fā)式干擾是一種通過(guò)轉(zhuǎn)發(fā)虛假時(shí)延信號(hào)，欺騙接收機(jī)定位的手段［2］。一旦GPS接收機(jī)鎖定虛假信號(hào)的載波及碼環(huán)，接收機(jī)在“無(wú)意識(shí)”中對(duì)錯(cuò)誤的信號(hào)進(jìn)行解算，從而產(chǎn)生錯(cuò)誤的定位信息。對(duì)于GPS/INS組合導(dǎo)航系統(tǒng)而言，GPS接收機(jī)部分一旦被成功欺騙，整個(gè)系統(tǒng)的導(dǎo)航數(shù)據(jù)將被污染，并輸出錯(cuò)誤的定位信息。

GPS轉(zhuǎn)發(fā)式干擾是非合作式的，干擾帶來(lái)的噪聲難以用建模的方法進(jìn)行濾除。因此，可考慮采用故障檢測(cè)的方法判斷系統(tǒng)是否受到誘偏干擾，并剔除受到干擾的數(shù)據(jù)。

本文在介紹轉(zhuǎn)發(fā)式干擾基本原理的前提下，提出局部狀態(tài)χ2檢驗(yàn)的方法，驗(yàn)證了該方法的可行性，并且通過(guò)仿真實(shí)驗(yàn)進(jìn)行了驗(yàn)證，結(jié)果表明該方法能夠有效檢測(cè)GPS轉(zhuǎn)發(fā)式干擾。

1 GPS/INS組合導(dǎo)航系統(tǒng)非線性模型

1.1 狀態(tài)方程

以導(dǎo)航參數(shù)誤差量為狀態(tài)變量，建立GPS/INS組合導(dǎo)航系統(tǒng)狀態(tài)方程，誤差較小，變化率較低，使每一步迭代中，濾波器能夠在更長(zhǎng)的時(shí)間內(nèi)進(jìn)行狀態(tài)估計(jì)，易于滿足導(dǎo)航中對(duì)實(shí)時(shí)性的要求。其狀態(tài)方程如下：

其中φE、φN、φU分別為平臺(tái)東、北、天誤差角，δvE、δvU、δvN為系統(tǒng)東、北、天向速度誤差，δL、δλ、δh為當(dāng)?shù)刈鴺?biāo)系下地位置誤差，εbx、εby、εbz為陀螺常值漂移，εrx、εry、εrz為陀螺隨機(jī)漂移，Δx、Δy、Δz為加速度計(jì)零偏，δtu為GPS鐘差等效距離，δtru為時(shí)鐘頻差等效距離變化率，F(xiàn)（t）是狀態(tài)轉(zhuǎn)移矩陣，G（t）為噪聲系數(shù)矩陣，W（t）為噪聲向量，由于過(guò)于復(fù)雜不再展開(kāi)，可參見(jiàn)文獻(xiàn)［3］。

1.2 量測(cè)方程

以偽距作為觀測(cè)量，結(jié)合星歷數(shù)據(jù)輸出量測(cè)更新值，建立量測(cè)方程。GPS接收機(jī)解算得到衛(wèi)星j之間的偽距可表示為［3］：

可得式（3）：

其中λI、LI、hI分別為慣導(dǎo)輸出的載體在地球大地坐標(biāo)系中的經(jīng)度、緯度及高度。量測(cè)方程如下：

2 GPS轉(zhuǎn)發(fā)式干擾基本原理

GPS轉(zhuǎn)發(fā)式干擾是通過(guò)轉(zhuǎn)發(fā)信號(hào)的時(shí)延控制及同步算法，使得被干擾目標(biāo)偏離原有軌跡，從而實(shí)現(xiàn)對(duì)GPS接收機(jī)及其載體的運(yùn)行軌跡控制的技術(shù)手段［4-5］。圖1所示即為典型的單站轉(zhuǎn)發(fā)多信號(hào)的基本原理。其原理是：在A點(diǎn)部署干擾機(jī)，通過(guò)轉(zhuǎn)發(fā)較大功率的帶有一定時(shí)延的GPS信號(hào)，將干擾機(jī)有效區(qū)域內(nèi)的任一C點(diǎn)的接收機(jī)誘騙至B點(diǎn)，使該區(qū)域內(nèi)的接收機(jī)在“無(wú)意識(shí)”狀態(tài)下解算出錯(cuò)誤的位置與授時(shí)信息［6］。

圖1 GPS轉(zhuǎn)發(fā)式干擾原理

設(shè)干擾機(jī)所在A點(diǎn)坐標(biāo)為（xA，yA，zA），解算得4顆衛(wèi)星S1、S2、S3、S4的偽距為ρ1A、ρ2A、ρ3A、ρ4A。A點(diǎn)干擾機(jī)的授時(shí)誤差為A，則A點(diǎn)的GPS位置方程可寫(xiě)為：

GPS轉(zhuǎn)發(fā)式干擾誘偏下，目標(biāo)接收機(jī)同時(shí)收到來(lái)自衛(wèi)星及干擾機(jī)的信號(hào)。由于干擾機(jī)轉(zhuǎn)發(fā)的干擾信號(hào)強(qiáng)度比GPS衛(wèi)星的信號(hào)強(qiáng)度強(qiáng)得多，加之前期進(jìn)行的大功率白噪聲壓制干擾，信號(hào)很容易解鎖。信號(hào)解鎖后，目標(biāo)接收機(jī)則接收到干擾機(jī)發(fā)送的事先加入時(shí)延的GPS衛(wèi)星信號(hào)［7］。設(shè)ρ1B、ρ2B、ρ3B、ρ4B為衛(wèi)星至誘偏后映射的目標(biāo)點(diǎn)B點(diǎn)的偽距，式（6）則為B點(diǎn)定位方程。

式（5）和式（6）做差可得時(shí)延控制方程：

3 基于局部狀態(tài)的χ2檢驗(yàn)方法

χ2檢驗(yàn)方法是組合導(dǎo)航系統(tǒng)中一類應(yīng)用最為廣泛的故障檢測(cè)方法［8］。它通過(guò)檢驗(yàn)所構(gòu)造的n維隨機(jī)變量的均值及方差是否滿足高斯分布來(lái)判斷系統(tǒng)是否發(fā)生故障。按所構(gòu)造的隨機(jī)變量不同χ2檢驗(yàn)可分為殘差χ2檢驗(yàn)法、狀態(tài)χ2檢驗(yàn)法、雙狀態(tài)χ2檢驗(yàn)法。

其中殘差χ2檢驗(yàn)法利用當(dāng)前量測(cè)值與量測(cè)值的一步預(yù)測(cè)之差建立故障檢測(cè)函數(shù)，通過(guò)對(duì)殘差的統(tǒng)計(jì)特性進(jìn)行分析來(lái)確定系統(tǒng)是否發(fā)生故障［9］。這種檢驗(yàn)方法對(duì)于突變的硬故障的檢測(cè)效果較好，而對(duì)于緩慢變化的軟故障難以及時(shí)檢測(cè)出來(lái)。

狀態(tài)χ2檢驗(yàn)的方法能夠較靈敏地檢測(cè)出緩慢變化的軟故障，適于檢測(cè)GPS轉(zhuǎn)發(fā)式干擾。但需要注意的是：若GPS干擾機(jī)有效實(shí)施了對(duì)組合導(dǎo)航系統(tǒng)中GPS接收機(jī)的欺騙干擾，在無(wú)虛警的情況下，一定能夠判斷出系統(tǒng)有故障；若判斷出系統(tǒng)存在故障，并不能確定這種故障是由GPS轉(zhuǎn)發(fā)式干擾造成的。換言之，GPS接收機(jī)受到欺騙干擾是判斷系統(tǒng)發(fā)生故障的充分而非必要條件。另外，狀態(tài)檢驗(yàn)的方法不僅需要濾波器提供量測(cè)更新信息，還需要單獨(dú)的“狀態(tài)遞推器”，這無(wú)疑增加了計(jì)算量，對(duì)于20維的狀態(tài)進(jìn)行χ2檢驗(yàn)更是如此。

鑒于此，本文在分析故障變化特性與故障診斷函數(shù)特性之間的關(guān)系的基礎(chǔ)上，提出了局部狀態(tài)χ2檢驗(yàn)。所謂局部狀態(tài)χ2檢驗(yàn)，是指僅采用對(duì)局部狀態(tài)即xχ=［δtu，δtru］進(jìn)行χ2檢驗(yàn)。這種方法可以實(shí)現(xiàn)對(duì)轉(zhuǎn)發(fā)式干擾的有效檢測(cè)。下面將詳細(xì)論證采用該方法的可行性。

3.1 觀測(cè)值的變化

假設(shè)GPS接收機(jī)被欺騙后，鎖定干擾機(jī)的碼環(huán)，在C點(diǎn)提取出偽距：

3.2 系統(tǒng)觀測(cè)模型的變化

在未受到轉(zhuǎn)發(fā)式干擾的情況下，GPS/INS組合導(dǎo)航系統(tǒng)觀測(cè)方程：

若受到加入偽距噪聲的干擾，則模型在理論上應(yīng)為：

3.3 χ2檢驗(yàn)中的變化

前面已經(jīng)闡述了轉(zhuǎn)發(fā)式干擾加入的時(shí)延噪聲與等效時(shí)鐘誤差校正距離的關(guān)系，下面將說(shuō)明這種時(shí)延噪聲對(duì)χ2檢驗(yàn)的影響。為了方便描述，給出一般的量測(cè)更新方程：

式中，Zk為量測(cè)值且Zk=Pc，若受到干擾可寫(xiě)作：

因此，未受到干擾的更新方程可寫(xiě)作：

若受到干擾則方程變?yōu)椋?/p>

4 仿真實(shí)驗(yàn)與結(jié)果分析

為了驗(yàn)證上節(jié)提出的局部狀態(tài)χ2檢驗(yàn)法，實(shí)驗(yàn)中模擬多站轉(zhuǎn)發(fā)時(shí)延干擾，假設(shè)干擾機(jī)所在坐標(biāo)分別為A1：（30.8°N，117.4°E，6 500 m），A2：（30.8° N，117.8°E，6 500 m），A3：（30.4°N，117.6°E，6 500 m），A4：（31.2°N，117.6°E，6 500 m）。載體真實(shí)軌跡如圖2粗實(shí)線所示，載體在200 s～300 s進(jìn)入干擾區(qū)域，按預(yù)想計(jì)劃，GPS接收機(jī)受到誘偏如圖2細(xì)虛線所示。

圖2 真實(shí)軌跡與預(yù)想誘偏效果

在不加入干擾的情況下，單獨(dú)進(jìn)行一次實(shí)驗(yàn)，狀態(tài)檢驗(yàn)如圖3。

圖3 未加入干擾時(shí)延的χ2狀態(tài)檢測(cè)

在200 s～300 s，加入干擾時(shí)延后的狀態(tài)檢驗(yàn)如圖4。

由圖可知，200 s～300 s的故障檢測(cè)函數(shù)值明顯比其他時(shí)間段要高很多。t=200 s加入時(shí)延干擾后，t=203 s即判斷出偽距故障的發(fā)生。局部狀態(tài)χ2檢驗(yàn)?zāi)軌蛞暂^高的靈敏度和較高的準(zhǔn)確性判斷轉(zhuǎn)發(fā)式干擾的發(fā)生。這就避免了GPS接收機(jī)將干擾數(shù)據(jù)引入INS反饋校正的發(fā)生。

圖4 加入干擾時(shí)延的χ2狀態(tài)檢測(cè)

5 結(jié)束語(yǔ)

為了檢測(cè)GPS轉(zhuǎn)發(fā)式干擾，根據(jù)故障診斷的啟示，采用局部狀態(tài)χ2檢驗(yàn)的方法，根據(jù)對(duì)鐘差等效距離及頻差等效距離變化率的故障檢測(cè)的結(jié)果，對(duì)INS的輸出進(jìn)行選擇性校正。仿真實(shí)驗(yàn)?zāi)M了典型的多站GPS轉(zhuǎn)發(fā)式干擾環(huán)境，結(jié)果表明了基于鐘差χ2狀態(tài)檢驗(yàn)方法能夠有效判斷干擾，但正如前文所述，GPS接收機(jī)受到轉(zhuǎn)發(fā)式干擾是判斷系統(tǒng)發(fā)生故障的充分而非必要條件。

［1］劉建業(yè)，曾慶化，趙 偉.導(dǎo)航系統(tǒng)理論與應(yīng)用［M］.西安：西北工業(yè)大學(xué)出版社，2010.

［2］張 頌，楊景曙，苗 苗.誘偏暨導(dǎo)航一體化平臺(tái)位置誤差補(bǔ)償算法［J］.現(xiàn)代雷達(dá)，2012，34（2）：18-19.

［3］王惠南.GPS導(dǎo)航原理與應(yīng)用［M］.北京：科學(xué)出版社，2003.

［4］焦 遜，陳永光，孫云輝.對(duì)GPS接收機(jī)實(shí)施欺騙式干擾的定位誤差研究［J］.電子對(duì)抗技術(shù)，2003，15（6）：21-24.

［5］王妍潔，董緒榮.轉(zhuǎn)發(fā)式GPS干擾仿真研究［J］.全球定位系統(tǒng)，2004，29（6）：16-18.

［6］劉延斌，蘇五星，閆抒升.轉(zhuǎn)發(fā)式欺騙信號(hào)干擾GPS接收機(jī)的效能分析［J］.空軍雷達(dá)學(xué)院學(xué)報(bào)，2004，18（4）：4-6.

［7］王 偉，陶業(yè)榮，王國(guó)玉，等.GPS欺騙干擾原理研究與建模仿真［J］.火力與指揮控制，2009，34（6）：115-118.

［8］Tallarida R J.Murray R B.Chi-Square Test［J］.Manual of Pharmacologic Calculations，1986：140-142.

［9］王仲濤，王 龍，王 楠.一種組合導(dǎo)航低成本地面驗(yàn)證技術(shù)［J］.四川兵工學(xué)報(bào)，2012，33（3）：59-61.

［10］錢華明.故障診斷與容錯(cuò)技術(shù)及其在組合導(dǎo)航系統(tǒng)中的應(yīng)用研究［D］.哈爾濱：哈爾濱工程大學(xué)，2004.

A Test Method for GPS/INS Integrated Navigation Repeater Jamming Based on Partial States Test

ZHU Li-xin1，2，MENG Yan1，2，Wang Jiang1，MA Chun-lai1，2
（1.Electronic Engineering Institute of PLA，Hefei 230037，China；
2.Province Key Laboratory of Electronic Restriction，Hefei 230037，China）

For the problem of GPS/INS Integrated Navigation Repeater Jamming，a method based on partial states chi-square test is presented.This method is based on the relationship between analysis of the characteristics of fault changes and the characteristics of fault diagnosis function.According to Demonstrate the feasibility of the method via changes in the observed value，changes in the observed model，and state value changes in the chi-square test，it is concluded that Repeater Jamming is a unnecessary and sufficient condition for GPS/INS Integrated Navigation's system fault.At last，by adding cajolery to simulation trajectory to simulate the Repeater Jamming，the simulation results show that this method can judge the occurrence of Repeater Jamming in high sensitivity and high accuracy.

GPS/INS Integrated Navigation，Repeater Jamming，partial states，chi-square test

TN967.2

A

1002-0640（2014）10-0146-04

2013-08-09

2013-10-17

朱立新（1963- ），女，安徽省無(wú)為人，副教授。研究方向：信息處理、目標(biāo)跟蹤、衛(wèi)星導(dǎo)航等。