柏如龍，霍立寰*，陳 兵，王海清，江 漫，楊 光，許瑞杰，廖桂生

(1.中國電子科技集團(tuán)公司第五十四研究所，河北 石家莊 050081；2.河北省電磁頻譜認(rèn)知與管控重點實驗室，河北 石家莊 050081；3.西安電子科技大學(xué) 雷達(dá)信號處理國家重點實驗室，陜西 西安 710071)

0 引言

無源干擾源定位技術(shù)能夠被動接收電磁波信號對干擾源進(jìn)行測向定位，是現(xiàn)代電子對抗的重要手段[1-3]。以美國海軍“NOSS”系列的海洋監(jiān)視衛(wèi)星為典型代表的基于衛(wèi)星平臺的時差定位體制，具有覆蓋范圍廣、作用距離遠(yuǎn)、接收隱蔽、生存能力強的優(yōu)點，能夠迅速、準(zhǔn)確地發(fā)現(xiàn)干擾源目標(biāo)[4-6]，在各領(lǐng)域都有著重要作用和廣泛的應(yīng)用前景。

目前的衛(wèi)星時頻差定位主要包括時頻差聯(lián)合定位和多時差定位[7-12]，利用了不同衛(wèi)星接收信號間的到達(dá)時間差(TDOA)和到達(dá)頻率差(FDOA)。相比于時差測量，頻差測量需要更多的參考站進(jìn)行標(biāo)校。因此，僅僅基于時差的三星時差定位體制成本更低、定位結(jié)果更為可靠。三星時差定位基于兩條時差線進(jìn)行交會定位，受衛(wèi)星位置和干擾源位置的影響，兩條時差線在某些時刻存在兩個交點，難以確認(rèn)干擾源的真實位置。文獻(xiàn)[13-18]的研究主要集中在時差測量優(yōu)化、定位方程解算等方面，包括不同體制信號的時差測量和校正方法，以及在不同定位場景下的定位方法。對于定位過程中出現(xiàn)的虛假點目標(biāo)，僅基于目標(biāo)位置的先驗知識進(jìn)行區(qū)分，當(dāng)先驗知識缺失時則無法判斷虛假點。因此，亟需提出有效、穩(wěn)健的虛假點消除方法。

針對上述問題，本文提出了基于多時刻測量的三星時差定位虛假點消除方法。首先建立了基于衛(wèi)星多時刻測量時差結(jié)果的三星時差定位模型，然后基于WGS-84地球橢球模型建立多時刻定位方程組，最后推導(dǎo)了牛頓迭代的求解方法。通過多時刻的測量數(shù)據(jù)直接獲取干擾源真實位置附近的定位結(jié)果，以該定位結(jié)果為初始值進(jìn)行逐時刻精確定位。該方法能夠在去除虛假點的基礎(chǔ)上，得到更高精度的定位結(jié)果；同時在先驗初始值偏差和時差測量誤差較大的情況下仍具有良好的穩(wěn)健性。

1 多時刻時差定位模型

圖1為建立大地坐標(biāo)系模型，包括地球橢球、指向0°經(jīng)度方向的X軸，指向北極方向的Z軸，以及根據(jù)右手定則得到的Y軸。主星接收到干擾源信號，兩顆輔星能夠收到干擾源的旁瓣信號。每個時刻能夠測量得到主星和輔星一、主星和輔星二的兩個到達(dá)時間差，形成如圖所示的時差面，在假設(shè)目標(biāo)零高程的條件下，兩個時差面和地球橢球面的交點即為三星定位的結(jié)果。由圖1可以發(fā)現(xiàn)，兩個時差面相交于一條直線，該直線與地球橢球面相交于兩個點，一個為干擾源目標(biāo)真實位置，另一個為虛假位置。在沒有先驗信息的情況下，僅依靠單時刻的時差難以區(qū)分真實位置和虛假位置。

圖1 三星干擾源定位示意Fig.1 Configuration of the tri-satellite interference source location

圖2給出了時差線存在虛假點時的時差線相交情況，可以發(fā)現(xiàn)，多個時刻的時差線在真實位置處逐漸收斂于一點，而虛假點位置的時差線相交較為發(fā)散。因此，可以利用多個時刻的時差測量結(jié)果進(jìn)行定位，從而判斷出干擾源的真實位置。

圖2 三星時差線相交情況示意Fig.2 Intersection of tri-satellite time difference line

假設(shè)干擾源的坐標(biāo)為(x,y,z)，第p個時刻時主星坐標(biāo)為(xp0,yp0,zp0)，輔星一坐標(biāo)為(xp1,yp1,zp1)，輔星二坐標(biāo)為(xp2,yp2,zp2)，干擾源與輔星一和輔星二在第p個時刻的時差測量值分別為τp1和τp2，滿足：

(1)

式中，c為光速，rp0，rp1和rp2分別為干擾源與主星、輔星一和輔星二在第p個時刻的距離。

根據(jù)距離的計算方式，以及地球橢球面模型，構(gòu)造如下定位方程組。

(2)

2 迭代求解算法

針對多時刻的測量結(jié)果構(gòu)成的定位方程組，直接求解較為困難，考慮采用牛頓迭代方法求解。將以上方程組進(jìn)行如下表示：

(3)

F(x,y,z)=0。

(4)

相應(yīng)的Hesse矩陣為：

(5)

若第k次迭代中干擾源坐標(biāo)為Pk=(xk,yk,zk)，則根據(jù)下式進(jìn)行迭代。

Pk+1=Pk-[F′(xk,yk,zk)]-1F(xk,yk,zk)。

(6)

基于多時刻的三星定位方法步驟如下：

步驟1：根據(jù)P個時刻的時差τp1和τp2，計算對應(yīng)的主輔星與干擾源的距離差(rp1-rp0)和(rp2-rp0)；

步驟2：初始化，設(shè)迭代次數(shù)k=0，干擾源目標(biāo)初始位置P0；

步驟3：計算該次迭代過程中的Hesse矩陣F′(xk,yk,zk)，并根據(jù)式(6)得到下一次迭代中干擾源位置的估計結(jié)果Pk+1；

步驟5：判斷目標(biāo)運動屬性，若為運動目標(biāo)，則輸出的定位結(jié)果作為初始值，對每個時刻的測量結(jié)果利用球形迭代方法得到定位結(jié)果的精確值。

算法流程如圖3所示。

圖3 基于多時刻測量的三星定位方法Fig.3 Tri-satellite localization method based on multi-time TDOA

3 仿真實驗

在下面的仿真試驗中，單時刻定位中定位時刻個數(shù)P=1，多時刻定位中定位時刻個數(shù)P>1，所提方法為利用多時刻定位去除虛假點后，將得到的定位結(jié)果作為單時刻定位的初始迭代值的方法。

3.1 多時刻定位影響因素分析實驗

多時刻定位的性能主要與選取時刻個數(shù)、時間總長度有關(guān)，通過仿真實驗探究其變化關(guān)系,如圖4所示。

(a) 定位誤差隨時刻個數(shù)變化情況

(b) 定位誤差隨時長變化圖4 多時刻定位方法性能影響試驗Fig.4 Performance experiment of the multi-time location method

干擾源位置為125°E，30°N，時差測量誤差為0.5 μs。首先仿真了時間總時長為20 min，多時刻定位和所提方法的定位誤差隨采用時刻個數(shù)的變化關(guān)系，可以發(fā)現(xiàn)：多時刻定位性能隨時刻個數(shù)增加而提高，并不斷趨近于所提方法，所提方法的定位誤差基本不變。然后仿真了采用時刻個數(shù)為10，多時刻定位和所提方法的定位誤差隨總時長的變化關(guān)系，可以發(fā)現(xiàn)：多時刻定位性能隨總時長增加而提高，并不斷趨近于所提方法，所提方法的定位誤差同樣保持穩(wěn)定，能夠?qū)崿F(xiàn)高精度地快速定位。

3.2 迭代初始值影響實驗

通過仿真實驗，對比迭代初始值對算法的性能影響。干擾源目標(biāo)位置為125°E，30°N，時差測量誤差在±0.5 μs內(nèi)均勻分布，蒙特卡洛次數(shù)為500，采用3個連續(xù)時刻測量得到的時差值進(jìn)行去除虛假點的目標(biāo)定位，時間間隔為5 min，迭代初始值誤差為初次迭代的先驗經(jīng)緯度偏差的均方根誤差。圖5給出了迭代初始值誤差在0°～20°情況下，單時刻定位、多時刻定位和所提方法定位誤差隨迭代初始值誤差的變化關(guān)系?？梢园l(fā)現(xiàn)：當(dāng)?shù)跏贾抵饾u增大時，單時刻定位的定位誤差惡化明顯，多時刻的定位方法基本不受迭代初始值的影響，但是在迭代初始誤差很較小時，定位效果相比單時刻定位較差。而所提方法在多時刻定位的基礎(chǔ)上，能夠獲得同時優(yōu)于單時刻和多時刻的定位精度。

圖5 不同方法的定位誤差隨迭代初始值誤差的變化關(guān)系Fig.5 Location error of different methods versus the initial value error of the iteration

3.3 測量誤差影響實驗

對比3種方法隨時差測量誤差的變化關(guān)系如圖6所示。干擾源目標(biāo)位置為125°E，30°N，多時刻定位采用5個連續(xù)時刻測量得到的時差值以去除虛假點，時間間隔為4 min，迭代初始值偏差的均方根為1°，蒙特卡洛次數(shù)為500。可以發(fā)現(xiàn)，隨著時差測量誤差增加，3種方法的定位誤差增大。其中，單時刻定位在時差測量誤差較小時，主要受迭代初始值影響，受時差測量誤差影響不明顯。多時刻定位方法由于采用解方程組的方式，受時差測量誤差影響較大。而基于多時刻測量的單時刻定位方法，在時差測量誤差增大時具有良好的穩(wěn)健性。

圖6 不同方法的定位誤差隨時差測量誤差的變化關(guān)系Fig.6 Location error of different methods versus the TDOA measurement error

4 實測數(shù)據(jù)驗證

圖7為某次實測數(shù)據(jù)的處理情況。

(a) 代價函數(shù)變化情況

(b) 時差線相交情況圖7 某次實測數(shù)據(jù)處理結(jié)果Fig.7 Results of the real data processing

實測數(shù)據(jù)處理后，多時刻定位方法能夠定位誤差為23.11 km，基于多時刻測量的方法得到定位誤差為9.33 km。從迭代過程中代價函數(shù)δ的變化情況可以看出，進(jìn)行迭代2～3次即可得到滿足要求的結(jié)果。從時差線相交情況可以看出，由于時差測量存在誤差，真實點位置多時刻的時差線也未能完全相交于一點。相比于真實點位置，虛假點處的時差線更為發(fā)散，因此能夠采用基于多時刻測量的方法去除虛假點。多時刻定位算法能夠定位至輻射源目標(biāo)真實位置附近，然而精度受時差測量誤差影響，難以得到較高的定位精度。而所提方法通過后續(xù)的逐點定位，在去除虛假點的基礎(chǔ)上，能夠大幅提高多時刻的定位精度。

5 結(jié)束語

三星時差定位通過兩條時差線的交點確定干擾源位置，然而，若兩條時差線存在兩個交點時，則難以確定真實點位置。本文提出基于多時刻測量的虛假點消除方法。該方法利用多個時刻的時差測量數(shù)據(jù)判斷干擾源真實位置所在區(qū)域，然后將該結(jié)果作為單時刻定位方程牛頓迭代方法求解的初始值，最終得到去除虛假點后高精度的定位結(jié)果。仿真實驗分析了對定位性能的影響因素，驗證了較大測量誤差下和先驗初始值偏差情況下的穩(wěn)健性，實測數(shù)據(jù)處理結(jié)果驗證了所提方法能夠有效、快速地得到較高精度的定位結(jié)果。