秦耀璐 楊淑萍 束鋒,,3 孫琳琳 陸錦輝 許正文 桂林卿

(1.南京理工大學(xué)電子工程與光電技術(shù)學(xué)院,南京 210094; 2.中國電波傳播研究所 電波環(huán)境特性及?；夹g(shù)重點實驗室,青島 266107; 3.東南大學(xué) 移動通信國家重點實驗室,南京 210096)

引 言

無源定位是指在不發(fā)射電磁波的條件下獲得目標(biāo)的位置. 利用小衛(wèi)星平臺,能進(jìn)行大范圍的搜索和電子偵察,基于衛(wèi)星的高精度無源定位技術(shù)在導(dǎo)航、跟蹤和電子對抗等軍事和民用領(lǐng)域有著廣泛的應(yīng)用. 到達(dá)時間差(time difference of arrival, TDOA)定位技術(shù)利用由目標(biāo)輻射源信號到達(dá)不同觀測站的時間差確定的多個雙曲面相交進(jìn)行定位,到達(dá)頻率差(frequency difference of arrival, FDOA)定位技術(shù)則利用目標(biāo)輻射源與觀測站相對運動產(chǎn)生的多普勒頻移差進(jìn)行定位. 相較于其他無源定位技術(shù),無源TDOA定位精度較高. 若目標(biāo)輻射源與觀測站之間存在相對運動,利用FDOA信息可以進(jìn)一步提高定位性能,TDOA/FDOA聯(lián)合定位能很大程度地提高定位精度,因此聯(lián)合定位技術(shù)有廣闊的應(yīng)用前景[1]. Ho在文獻(xiàn)[2]中提出了利用TDOA/FDOA對運動目標(biāo)進(jìn)行定位的經(jīng)典兩階加權(quán)的最小二乘(weighted least squares, WLS) 算法. 近年來,定位領(lǐng)域的國內(nèi)外學(xué)者提出了許多利用TDOA和FDOA測量信息進(jìn)行定位的算法[3-9]. 文獻(xiàn)[3]提出了基于加權(quán)代價函數(shù)的TDOA/FDOA聯(lián)合定位算法,避免了加權(quán)WLS算法存在的只能在一定的噪聲環(huán)境中使得定位誤差達(dá)到克拉美·羅界(Cramer-Rao lower bound, CRLB)的問題,使得定位對于噪聲的穩(wěn)定性大幅提高. 文獻(xiàn)[4]提出了利用最大比合并技術(shù)基于最大似然原理的TDOA/FDOA定位算法,該算法的定位性能在有無地球約束的條件下均優(yōu)于單獨定位性能,并且逼近CRLB,但最大比合并的聯(lián)合搜索定位算法實現(xiàn)復(fù)雜度較高. 文獻(xiàn)[5]提出了基于網(wǎng)格搜索的加權(quán)最大似然代價函數(shù)定位算法,該算法與最大比合并的定位性能近似相等,但在復(fù)雜度上,節(jié)省了一半左右的計算量. 由于最大似然估計問題的非凸性質(zhì),如果沒有良好的初步估計,難以獲得全局最優(yōu)解,因此,文獻(xiàn)[6]將定位問題重新定義為加權(quán)最小二乘問題,將定位問題最終化成半正定松弛問題求解．TDOA/FDOA聯(lián)合定位方法通常不具有直接的解決方案,并且需要數(shù)值方法來確定發(fā)射器的位置,文獻(xiàn)[7]提出了解決非線性優(yōu)化問題基于粒子群算法的TDOA/FDOA聯(lián)合定位方法. 文獻(xiàn)[8]研究了采用最陡下降法和最大似然估計法的分布式被動傳感網(wǎng)絡(luò)中移動目標(biāo)源的TDOA/FDOA定位算法. 在最新的研究中,文獻(xiàn)[9]將目標(biāo)位置估計問題轉(zhuǎn)化為二次約束二次優(yōu)化問題,使用拉格朗日乘數(shù)法求得了近似解析解. TDOA/FDOA定位算法層出不窮,因而對評判定位算法性能的標(biāo)準(zhǔn)提出了迫切要求. 研究聯(lián)合定位的CRLB[10],可以衡量各種定位算法的性能優(yōu)劣,分析影響定位精度的各種因素.

Ho和Chan于1997年推導(dǎo)了多星TDOA和FDOA聯(lián)合定位的誤差CRLB[11],由此確定了無偏估計器中最優(yōu)無偏估計器所能達(dá)到的最好定位性能. 孫仲康在文獻(xiàn)[12]中推導(dǎo)了TDOA測量誤差較小情況下的多星TDOA定位方差. Pattison和Chou在文獻(xiàn)[13]中分析了雙星聯(lián)合定位的幾何精度因子(geometric dilution of precision, GDOP)和方差,發(fā)現(xiàn)TDOA定位曲面、FDOA定位曲面與地球表面夾角為90°時定位精度最好. 文獻(xiàn)[14]和文獻(xiàn)[15]對雙星聯(lián)合定位算法和誤差分析進(jìn)行了更深入的研究. 由于Ho推導(dǎo)的聯(lián)合定位方差公式CRLB并不適合于TDOA/FDOA三星單獨定位,故文獻(xiàn)[16-17]討論了三星定位誤差分析. 文獻(xiàn)[18]分析了多星TDOA/FDOA聯(lián)合定位中時差誤差和頻差誤差對定位精度的影響. Ho等[2]于2007年在考慮接收機位置和速度存在測量誤差情況下,推導(dǎo)了運動目標(biāo)多星聯(lián)合定位的方差下界,即CRLB,但是未能給出顯式表達(dá)式.

本文將應(yīng)用Taylor級數(shù)、矩陣?yán)碚摵徒y(tǒng)計理論等多種手段進(jìn)行簡化處理獲得靜態(tài)目標(biāo)的方差CRLB顯式表達(dá)式,有效地簡化定位精度分析. 同時考慮到三星單獨定位不能直接應(yīng)用該式,我們將雙星的聯(lián)合定位誤差推導(dǎo)思路應(yīng)用于三星單獨定位情況,并獲得其定位誤差公式. 進(jìn)一步,本文分析和討論了多星衛(wèi)星構(gòu)型對定位性能的影響. 仿真結(jié)果表明,FDOA較容易產(chǎn)生定位盲區(qū),且該盲區(qū)無法通過衛(wèi)星構(gòu)型設(shè)計來避免,而TDOA的定位盲區(qū)可以通過設(shè)計良好的衛(wèi)星構(gòu)型完全消除. 此外,TDOA/FDOA聯(lián)合定位時,能有效避免TDOA或FDOA帶來的定位盲區(qū),且能顯著提高定位精度.

1 系統(tǒng)模型

三維空間的多星無源定位系統(tǒng)由M個衛(wèi)星組成,位置向量為ri=(xi,yi,zi)T(i=1,2,…,M),未知目標(biāo)源位置坐標(biāo)向量為r=(x,y,z)T. 從未知目標(biāo)源到衛(wèi)星i的單位位置向量可表示為

(1)

目標(biāo)源到任何兩對衛(wèi)星j和i的TDOA 和FDOA可表示為

(2)

(3)

式中:c表示光速;fc為目標(biāo)源載波頻率;vi和v分別為衛(wèi)星i和目標(biāo)源的運行速度. 考慮衛(wèi)星位置誤差與速度誤差,對式(2)和式(3)進(jìn)行Taylor展開,得到如下表達(dá)式:

Δτji=τji(r+Δr,ri+Δri,rj+Δrj)-τji(r,ri,rj)

(4)

Δfji=fji(r+Δr,ri+Δri,rj+Δrj,vi+Δvi,

vj+Δvj)-fji(r,ri,rj,vi,vj)

(5)

(6)

(7)

式中,βi和pij表示為

pij=ri-rj．

(8)

選擇衛(wèi)星1為主星,收集所有的時差和頻差構(gòu)成TDOA和FDOA向量,共有N=M-1個獨立時差和頻差對,即

Δτ=(Δτ21,Δτ31,…,ΔτM1)T,

Δf=(Δf21,Δf31,…,ΔfM1)T.

(9)

將時差向量和頻差向量合成大的向量,定義如下向量矩陣表達(dá)式

(10)

考慮TDOA和FDOA測量誤差、衛(wèi)星自定位誤差,以及衛(wèi)星速度測量誤差的影響,利用式(4)～(7),定位誤差方程可近似表示為向量矩陣形式

Y(r,rs+Δrs,vs+Δvs)-Y(r,rs,vs).

(11)

式中:r表示目標(biāo)位置;rs+Δrs表示測量獲得的衛(wèi)星位置,rs為衛(wèi)星實際位置,Δrs為衛(wèi)星位置測量誤差;vs+Δvs表示測量獲得衛(wèi)星的速度向量,vs為衛(wèi)星實際速度向量,Δvs為衛(wèi)星速度向量測量誤差;而

Y(r,rs+Δrs,vs+Δvs)-Y(r,rs,vs)=Δrs+Δvs．

(12)

令

(13)

于是式(11)更改為

(14)

由式(14)有

(15)

式中:左邊表示定位誤差;右邊第一項表示TDOA/FDOA測量誤差,第二項表示衛(wèi)星自定位誤差,最后一項表示衛(wèi)星速度測量誤差,均服從高斯分布．式(15)中各矩陣向量可表示如下:

T=(-e,IM-1)T.

(16)

式中:

Δri=Δr1+Δri1,i∈{2,3,…,M};

e=(1,1,…,1)T;

(17)

2 多星定位誤差分析

2.1 多星TDOA/FDOA聯(lián)合定位誤差

基于式(14)和(15)得到接收的時差頻差聯(lián)合向量概率密度如下[10]:

(18)

我們得到Fisher信息矩陣[10]

J=E{[lnp(Y;r)][lnp(Y;r)]T}

(19)

式中,

Gr=(Gτ,rGf,r),

(20)

Gτ,r與Gf,r均為3×N(N≥3)的矩陣．定義

u=(u1,u2,…,uM),β=(β1,β2,…,βM),

T=(-e,IM-1)T,

(21)

利用矩陣T,結(jié)合式(8)和式(21),得到:

(22)

(23)

以及

Σ=E{ΔYΔYH}

(24)

GDOP2(x,y,z) =tr{J-1}

(25)

當(dāng)考慮地球約束[11]時,

(26)

式中,re為地球半徑．則受地球約束的多星TDOA/FDOA聯(lián)合定位的GDOP為

GDOP2(x,y,z)=tr{J-1-J-1F(FTJ-1F)-1FTJ-1}.

(27)

式中,F=(2x,2y,2z)T. 上面推導(dǎo)的多星GDOP,也可稱為CRLB聯(lián)合公式,適用于四星及以上單獨定位系統(tǒng)和三星及以上聯(lián)合定位系統(tǒng). 因此,我們將在下節(jié)推導(dǎo)考慮地球約束后的三星TDOA單獨定位誤差公式.

2.2 三星單獨定位誤差

上節(jié)推導(dǎo)的多星聯(lián)合公式不適用于三星單獨時差或頻差定位,主要原因是僅有兩個測量參數(shù),導(dǎo)致Fisher信息矩陣為奇異矩陣,無法求逆. 本節(jié)將以TDOA為例,專門推導(dǎo)三星單獨TDOA定位的GDOP.

對于三星單獨定位,兩測量參數(shù)加上地球約束三個方程恰好可聯(lián)立求解目標(biāo)的位置向量,其思路類似于雙星TDOA/FDOA. 類似于文獻(xiàn)[19],考慮衛(wèi)星位置測量誤差后,如式(2) 所示的TDOA定位方程以及式(26)地球約束依據(jù)Taylor級數(shù)可展開為:

(28)

(29)

0=(2x,2y,2z)T|r=r0Δr

(30)

式中:Δr為定位目標(biāo)的誤差項;Δr1為衛(wèi)星1也就是主星的絕對測量誤差,Δri1為衛(wèi)星i相對于主星的相對測量位置誤差,即Δri=Δr1+Δri1,i∈{2,3}. 為表達(dá)方便,下面推導(dǎo)中忽略r=r0. 將式(28)～(30)組合成矩陣向量形式,有

Δd=HΔr+H1Δr1+ΔH2．

(31)

式中:

H=(-u21,-u31,2r)T;

H1=(u2+u1,u3+u1,03)T;

(32)

其中ui和uji見式(1)和式(6),03表示全零的三維列向量. 由式(31)得到

Δr=H-1(Δd-H1Δr1-ΔH2)．

(33)

根據(jù)GDOP的定義,我們有

GDOP2(x,y,z) =tr{E(ΔrΔrT)}

(34)

3 衛(wèi)星構(gòu)型設(shè)計討論

地球半徑re=6 378.137 km,以五星為例分析多星構(gòu)型,圖1是五星的定位系統(tǒng)圖. 地心為O,其中主星為A,副星為B、C、D和E. 星下點定義為主星和地心的連線與地球表面的交點,星下點為S. 主星高程為h,假定目標(biāo)位于P處,θ表示主星-副星連線與主星-星下點連線夾角.AB=AC=AD=AE為主副星之間的距離. 給定主星的位置,四顆副 星的經(jīng)緯度坐標(biāo)及高程隨θ變化而變化. 由于衛(wèi)星間距一般在10 km到200 km之間變化,此值同衛(wèi)星高程與地球半徑相比可忽略. 考慮到五星編隊飛行,于是五顆衛(wèi)星速度向量可近似為相等.

圖1 五星定位系統(tǒng)框圖Fig.1 Block diagram of five-satellite localization system

定理1 目標(biāo)位于地球表面且距離星下點距離為Rkm范圍無定位盲區(qū)充分條件:主星-副星連線與主星-星下點連線夾角θ滿足φ<θ<180°-φ.

證明:

當(dāng)目標(biāo)P位于距離星下點Rkm圓周上時,對應(yīng)的∠AOP為

(35)

利用此角度可獲得線段AP長度等于

(36)

于是,可得∠PAO=φ如下:

(37)

若φ>θ,顯而易見,在AC延長線與地球表面交點至P點一段存在定位盲區(qū). 為了保證在星下點Rkm范圍內(nèi)不出現(xiàn)定位盲區(qū),令φ≤θ,考慮目標(biāo)P距離星下點Rkm,且位于主星-副星連線AC的延長線與地球表面交點附近,此時φ=θ,即主星A、副星C、目標(biāo)P在同一條直線上. 令主星A為衛(wèi)星1,副星C為衛(wèi)星3,由ui的定義式,且由于衛(wèi)星間距變化與同衛(wèi)星高程和地球半徑相比可忽略,得到u3=u1,即

u31=u3-u1=0．

(38)

僅考慮單獨TDOA定位時差測量誤差,式(25)簡化為

GDOP2(x,y,z) =tr{J-1}

(39)

(40)

根據(jù)Cauchy-Schwarz不等式可得

=9．

(41)

由式(41),有

(42)

(43)

目標(biāo)P位于主星-副星連線的延長線與地球表面的交點附近,利用式(38),得到‖u31‖2=0. 由式(43)可知,相較于其他區(qū)域,目標(biāo)位于主星-副星連線的延長線與地球表面的交點附近時,定位性能下降,形成定位盲區(qū).

同理,當(dāng)θ>180°-φ,即θ+φ>180°,AE的反向延長線與地球表面的交點至目標(biāo)P一段存在定位盲區(qū).

令θ=180°-φ,即目標(biāo)P位于AE的反向延長線與地球表面交點附近,可以得到,此時定位性能下降,該處形成定位盲區(qū),證明同上.

綜上,為了保證在星下點Rkm范圍內(nèi)不出現(xiàn)定位盲區(qū),四顆副星與主星連線同主星到星下點連線的夾角θ必須滿足下面不等式:

φ<θ<180°-φ．

(44)

證畢.

(45)

(46)

通過分析式(42)和(46)可知,增大αi可提高定位性能,當(dāng)保持主副星之間的間距不變時,αi隨相鄰主星-副星連線間夾角βi的增大而增大. 如圖1所示具有良好對稱結(jié)構(gòu)的衛(wèi)星構(gòu)形βi均為β,即∠BAC=∠BAE=∠CAD=∠DAC=β,利用余弦定理可得β與θ關(guān)系為

cosβ=cosθ2．

(47)

當(dāng)θ=90°時,β取得最大值90°.因而主星-副星連線與主星-星下點連線夾角θ趨于90°時,則主星-目標(biāo)連線與副星-目標(biāo)連線的夾角αi越大,定位性能越好.

4 仿真分析

設(shè)置系統(tǒng)仿真參數(shù)如下:TDOA測量誤差為50 ns,FDOA測量誤差5 Hz. 采用五星定位系統(tǒng),主星位于上方頂點,高程為h=1 000 km,東經(jīng)130°,北緯0°,其他四顆副星位于主星下方構(gòu)成正方形,主星到副星間距保持150 km不變,主星到其他四副星連線同主星到星下點連線夾角均為θ.

4.1 衛(wèi)星構(gòu)型對GDOP的影響

圖2、圖3和圖4分別繪出了θ=40°,θ=70°時,TDOA單獨定位的GDOP曲線、FDOA單獨定位的GDOP曲線以及TDOA/FDOA聯(lián)合定位的GDOP曲線.

(a) θ=40°

(b) θ=70°圖2 TDOA單獨定位的GDOP曲線Fig.2 GDOP of TDOA localization

(a) θ=40°

(b) θ=70°圖3 FDOA單獨定位的GDOP曲線Fig.3 GDOP of FDOA localization

(a) θ=40°

(b) θ=70°圖4 TDOA/FDOA聯(lián)合定位的GDOP曲線Fig.4 GDOP of joint TDOA and FDOA localization

由圖2(a)可知,當(dāng)θ=40°時,TDOA單獨定位存在四個定位盲區(qū),它們由主星-副星連線與地球交點附近區(qū)域形成盲區(qū). 在設(shè)計衛(wèi)星構(gòu)型時盡量保證所有主星-副星連線與地球不相交,則可避免這些定位盲區(qū). 圖2(b)示出了θ=70°時定位盲區(qū)已經(jīng)消除. 比較圖2(a)和圖2(b)發(fā)現(xiàn)TDOA在星下點的單獨定位精度由2 km提升到0.8 km,遠(yuǎn)離星下點精度也有顯著提升,因此增加θ角可提高TDOA定位精度.

仔細(xì)觀察圖3會發(fā)現(xiàn):有一條通過星下點的定位盲區(qū),此盲區(qū)內(nèi)FDOA的定位誤差較大．該盲區(qū)實際上是垂直于衛(wèi)星速度方向通過主星的平面同地球相交圓周形成,由于衛(wèi)星繞地球飛行,此平面與地球表面肯定相交,因此定位盲區(qū)是沒有地球約束FDOA定位系統(tǒng)本身所固有的. 產(chǎn)生物理機制如下:此圓周上任意一點到衛(wèi)星連線近似垂直于衛(wèi)星速度方向,因此速度在此連線上沒有分量,產(chǎn)生的Doppler頻率為零,不起任何作用,也就是不提供任何定位信息. 因此角度變化對FDOA定位盲區(qū)消除沒有影響.

結(jié)合圖4(a)和圖4(b)可知θ角度變化對聯(lián)合定位精度基本沒有影響．比較圖2、圖3和圖4可知,TDOA/FDOA聯(lián)合定位時,不存在定位盲區(qū),且相較于TDOA單獨定位和FDOA單獨定位,定位精度明顯提高.

4.2 衛(wèi)星位置誤差、衛(wèi)星速度誤差對GDOP的影響

系統(tǒng)仿真參數(shù)設(shè)置如下:主星到其他四副星連線同主星到星下點連線夾角均為θ=70°.

圖5繪出了衛(wèi)星位置絕對誤差和相對誤差為1 000 m、100 m,及100 m、10 m兩種情況下TDOA/FDOA聯(lián)合定位精度隨離星下點距離的變化曲線. 由圖5可知:當(dāng)衛(wèi)星位置絕對誤差為100 m,相對誤差為10 m時,衛(wèi)星位置誤差對定位精度近乎沒有影響;在衛(wèi)星位置絕對誤差為1 000 m,相對位置誤差為100 m時,衛(wèi)星位置誤差較大,衛(wèi)星位置對定位精度的影響已經(jīng)凸顯,相較于衛(wèi)星位置絕對誤差為100 m相對誤差為10 m,定位精度下降約1.2 km. 圖6繪出了衛(wèi)星速度誤差分別為2 m/s,10 m/s時,TDOA/FDOA聯(lián)合定位精度隨離星下點距離的變化曲線. 由圖6可知,衛(wèi)星速度誤差對定位精度影響較大. 相較于無衛(wèi)星速度誤差的TDOA/FDOA聯(lián)合定位,對距星下點2 000 km左右的區(qū)域,衛(wèi)星誤差為2 m/s時,定位精度下降約3 km;衛(wèi)星誤差為10 m/s時,定位精度下降約6 km. 比較圖5與圖6可知,衛(wèi)星位置誤差提升一個數(shù)量級,定位精度約下降1.2 km,而衛(wèi)星速度誤差由2 m/s增大至10 m/s,定位精度平均下降約2.7 km. 速度測量精度對定位精度的影響較大,而測距精度對目標(biāo)定位精度的影響較小. 圖7繪出了同時考慮衛(wèi)星位置誤差和速度誤差,TDOA/FDOA聯(lián)合定位精度隨離星下點距離的變化曲線. 比較圖5和圖6可知,當(dāng)衛(wèi)星位置誤差較小時,速度誤差為影響定位精度的主要因素. 同時圖7示出綜合考慮位置誤差與速度誤差,當(dāng)位置誤差達(dá)到一定量級時,衛(wèi)星位置誤差將會使定位精度產(chǎn)生大幅度下降.

圖5 衛(wèi)星位置測量誤差對定位精度的影響Fig.5 The effect of satellite position measurement errors on location accuracy

圖6 衛(wèi)星速度測量誤差對定位精度的影響Fig.6 The effect of satellite velocity measurement errors on location accuracy

圖7 衛(wèi)星位置和速度測量誤差對定位精度的影響Fig.7 The effect of satellite position measurement and velocity measurement errors on location accuracy

5 結(jié) 論

本文推導(dǎo)了多星TDOA/FDOA聯(lián)合定位的誤差公式,該公式充分考慮了TDOA測量誤差、FDOA測量誤差、衛(wèi)星自定位誤差和速度測量誤差. 并且對衛(wèi)星自定位誤差與速度誤差進(jìn)行仿真分析,發(fā)現(xiàn)衛(wèi)星速度誤差是影響定位精度的關(guān)鍵因素. 最后對衛(wèi)星構(gòu)型進(jìn)行了簡單的討論. 由此可見,良好的衛(wèi)星構(gòu)形、高精度的衛(wèi)星速度測量與TDOA/FDOA聯(lián)合定位相結(jié)合可以保證相當(dāng)高的定位精度．