彭華峰，夏暢雄，曹金坤

（西南電子電信技術(shù)研究所，成都610041）

基于最小二乘融合估計的雙星時頻差定位?

彭華峰，夏暢雄，曹金坤

（西南電子電信技術(shù)研究所，成都610041）

針對單參考站條件下具有測量誤差和星歷誤差時定位精度不高的特點，提出了一種基于多次觀測數(shù)據(jù)的最小二乘融合估計定位算法，該算法無需增加觀測條件即可有效提高輻射源定位精度。分析了測量誤差、星歷誤差對單參考站單次定位及融合定位精度的影響，推導(dǎo)了測量誤差、星歷誤差對定位誤差的傳遞公式，提出了含星歷誤差影響的最小二乘融合估計加權(quán)算法。通過Monte-Carlo仿真驗證了誤差分析結(jié)果和定位算法，并比較了加權(quán)最小二乘估計定位和單次定位的性能。仿真試驗表明：在相同觀測精度條件下，加權(quán)最小二乘融合定位可極大地提高輻射源定位精度，最大提高10倍以上。

雙星；時差頻差；定位；估計；數(shù)據(jù)融合

1 引言

測量信號的到達(dá)時頻差對目標(biāo)進(jìn)行定位是最重要的定位技術(shù)之一［1-4］。利用兩顆衛(wèi)星測量上行信號的到達(dá)時頻差可以實現(xiàn)對地面輻射源定位，在導(dǎo)航、遙感及干擾源定位等方面有廣泛的應(yīng)用［3-4］。

利用同步軌道雙星對地面進(jìn)行定位時，由于衛(wèi)星與輻射源距離遠(yuǎn)、相對運動較小、信號信噪比低等因素影響，定位精度比較低，在一些應(yīng)用中通常難以同時找到滿足條件的多顆衛(wèi)星進(jìn)行測量以進(jìn)一步提高定位精度。工程實踐中，為了克服星歷誤差及轉(zhuǎn)發(fā)器時延等的影響，通常需要利用多個參考站校正誤差影響［3-4］，但要滿足同時多參考站條件要求相對比較苛刻，影響了系統(tǒng)實際效能的發(fā)揮。如何在現(xiàn)有雙星定位條件下，通過融合多次觀測數(shù)據(jù)進(jìn)行融合估計實現(xiàn)單參考站條件下的輻射源定位，充分發(fā)揮單參考站條件下的定位能力是目前定位系統(tǒng)所需要解決的重要難點。

本文通過詳細(xì)推導(dǎo)誤差傳遞方程，將星歷誤差傳遞到測量誤差后，然后取誤差協(xié)方差的倒數(shù)為權(quán)值進(jìn)行加權(quán)最小二乘融合估計，提高單參考站條件下定位精度，并通過計算機(jī)仿真比較了加權(quán)融合估計和非融合估計的定位精度，驗證加權(quán)最小二乘融合估計算法的有效性。

2 雙星定位原理

2.1 定位原理

雙星定位原理為：利用兩顆衛(wèi)星接收到達(dá)信號的時頻差（TDOA／FDOA），聯(lián)立地球橢球方程可以求解得到目標(biāo)的位置向量。

假設(shè)t時刻目標(biāo)位置為Rs，參考站位置為Rr，主星位置為Rm，鄰星的位置為Rn，則測量方程為

式中，rij=Ri-Rj，vij=Vi·j（Ri-Rj）／Ri-Rj，（i=m，n；j=r，s），Vij=Vi-Vj為相對速度向量；vij為相對速率；c是電磁波在空間的傳播速度，ri（i =m，n）為衛(wèi)星i與地面輻射源的距離。

此外，地球球面為橢球面，其曲面方程為

式中，a=6 378.14 km、b=6 356.755 km分別為地球的長短半徑。

2.2 單次定位算法

基于單個時刻的兩個觀測量求解地面目標(biāo)的位置時，已經(jīng)提出了多種算法，如最大似然估計法、牛頓迭代算法等。本文給出牛頓迭代算法做單次定位仿真。

令

3 融合估計定位原理［5］

3.1 誤差傳遞方程

對測量方程兩邊全微分可得

上式即為絕對星歷誤差表示時的誤差方程。

本文采用相對星歷誤差進(jìn)行分析。令相對星歷誤差為

式中，G=C+D。式（9）即為相對誤差表示時的誤差方程。

3.2 加權(quán)最小二乘融合定位算法

假設(shè)測量誤差、主星星歷誤差和鄰星星歷誤差的協(xié)方差矩陣分別為Cov（ΔZ）、Cov（ΔXm）和Cov（ΔδX），則誤差方程右邊等價為一個新的觀測量ΔY，其協(xié)方差為

由加權(quán)最小二乘估計理論可知，取觀測量的測量誤差的協(xié)方差的倒數(shù)為其權(quán)值時估計效能最好，因此上式表明當(dāng)有星歷誤差時，需要將星歷誤差傳遞到測量誤差后，取傳遞后的測量誤差協(xié)方差的倒數(shù)為權(quán)值才能實現(xiàn)最優(yōu)估計。設(shè)W為Cov（ΔY）的對角線元素的倒數(shù)，則估計結(jié)果為

4 仿真分析

為了分析雙星時頻差對地面輻射源的定位，利用STK仿真產(chǎn)生了兩顆同步軌道衛(wèi)星的星歷并隨機(jī)生成了一個地面目標(biāo)。模擬生成真實的測量參數(shù)，并分別對測量數(shù)據(jù)和衛(wèi)星星歷加正態(tài)高斯分布的隨機(jī)噪聲。假設(shè)衛(wèi)星星歷為正態(tài)分布的隨機(jī)噪聲是合理的，這是因為衛(wèi)星星歷可以通過實時測量和幾何解算得到（如三站測距定位等），通過幾何解算得到的星歷符合正態(tài)高斯分布。測量數(shù)據(jù)和星歷數(shù)據(jù)的采樣間隔為1min。本文所有研究都是在該仿真條件下進(jìn)行的。

4.1 單次定位仿真

采用Monte-Carlo分別仿真了在無星歷誤差、無相對星歷誤差和有相對星歷誤差條件下雙星定位系統(tǒng)對輻射源的定位誤差，如圖1～3所示。設(shè)定的參數(shù)中，測量精度為300m、0.1 Hz，主星星歷誤差為20 km和2m／s，相對星歷誤差為5 km和0.5m／s。

圖1 無星歷誤差時單次定位Monte-Carlo仿真Fig1Monte-Carlo simulation of positioning withoutephemeris error

圖2 無相對星歷誤差時Monte-Carlo仿真Fig.2 Monte-Carlo simulation of positioning without relative ephemeris error

圖3 有相對星歷誤差時的Monte-Carlo仿真Fig.3 Monte-Carlo simulation of positioning with relative ephemeris error

從仿真結(jié)果可以看出，主星星歷誤差對定位結(jié)果影響很大，相對星歷誤差對定位結(jié)果影響更大。當(dāng)無星歷誤差時，單次定位誤差在100 km以內(nèi)；無相對星歷誤差時，單次定位誤差在250 km以內(nèi)；存在相對星歷誤差時，單次定位誤差達(dá)到600 km。

4.2 融合定位仿真分析

采用60組數(shù)據(jù)進(jìn)行融合估計，圖4～6對應(yīng)無星歷誤差、無相對星歷誤差和有相對星歷誤差時的Monte-Carlo仿真結(jié)果。圖中虛線圈對應(yīng)定位誤差橢圓。從Monte-Carlo仿真結(jié)果可以看出，利用該加權(quán)估計算法得到的定位誤差仍然滿足高斯分布，驗證了該加權(quán)估計方法的正確性。

對比圖1～3可知，通過融合估計后，定位精度得到極大提高。當(dāng)無星歷誤差時，定位精度由100 km提升至10 km以內(nèi)；無相對星歷誤差時，定位精度由250 km提升至15 km以內(nèi)；有星歷誤差時定位精度由600 km提升至100 km以內(nèi)。

圖4 無星歷誤差時融合定位Monte-Carlo仿真（60組數(shù)據(jù)融合）Fig.4 Monte-Carlo simulation of fusion positioning withoutephemeris error（60 observations fused）

圖5 無相對星歷誤差時融合定位Monte-Carlo仿真（60組數(shù)據(jù)融合）Fig.5Monte-Carlo simulation of fusion positioning without relative ephemeris error（60 observations fused）

圖6 有相對星歷誤差時融合定位Monte-Carlo仿真（60組數(shù)據(jù)融合）Fig.6Monte-Carlo simulation of fusion positioning with relative ephemeris error（60 observations fused）

5 結(jié)論

通過推導(dǎo)單參考站條件下融合定位的誤差傳遞方程，證明當(dāng)存在星歷誤差時，需要將星歷誤差傳遞到測量誤差，取傳遞后的測量誤差方差倒數(shù)進(jìn)行加權(quán)。理論分析和Monte-Carlo定位仿真結(jié)果均表明，基于該加權(quán)方法的最小二乘估計算法在無星歷誤差時，定位精度由100 km提高至10 km；無相對星歷誤差時，定位精度由250 km以內(nèi)提高到15 km以內(nèi)；有相對星歷誤差時，定位精度由600 km以內(nèi)提高到100 km以內(nèi)。

本文從理論上分析和計算機(jī)仿真研究了具有隨機(jī)星歷噪聲條件下觀測數(shù)據(jù)加權(quán)融合估計算法，但利用網(wǎng)上公開星歷時的最優(yōu)加權(quán)估計算法是下一步需要研究的問題。

［1］Ho K C，Xu Wenwei.An Accurate Algebraic Solution for Moving Source Location Using TDOA and FDOA Measurements［J］.IEEE Transactions on Signal Processing，2004，52（9）：2453-2463.

［2］Ho K C，Lu Xiaoning，Kovavisaruch L.Source Localization Using TDOA and FDOA Measurements int the Presence of Receiver Location Errors：Analysis and Solution［J］.IEEE Transactions on Signal Processing，2007，55（2）：684-696.

［3］Ho K C，Chan Y T.Geolocation of a known altitude object from TDOA and FDOA Measurements［J］.IEEE Transactions on Aerospace and Electronic Systems，1997，33（3）：770-783.

［4］Haworth D P，Smith N G，Bardelli R.Interferrence localization for EUTELSAT Satellites-The First European Transmitter Location System［J］.International Journal of SatelliteCommunications，1997，15（2）：155-183.

［5］韓崇昭，朱洪艷，段戰(zhàn)勝，等.多源信息融合［M］.北京：清華大學(xué)出版社，2006. HAN Chong-zhao，ZHU Hong-yan，DUAN Zhan-sheng，etal.Multisource information fusion［M］.Beijing：Tsinghua U-niversity Press，2006.（in Chinese）

PENG Hua-feng was born in Xinhua，Hunan Province，in 1979.He received the M.S.degree from Sichuan University in 2006.He is now an engineer and currently working toward the Ph.D.degree.His research concerns target location，tracking and satellite orbitmeasurement.

Email：hualong913＠sina.com

夏暢雄（1979—），男，湖南岳陽人，2008年于西南電子電信技術(shù)研究所獲博士學(xué)位，現(xiàn)為高級工程師，主要研究方向為目標(biāo)定位等；

XIA Chang-xiong was born in Yueyang，Hunan Province，in 1979.He received the Ph.D.degree from Southwest Electronics and Telecommunication Technology Research Institute in 2008.He is now a senior engineer.His research concerns target location.

Email：xiaxiong＠sohu.com

曹金坤（1980—），男，安徽阜陽人，2006年于信息工程大學(xué)獲碩士學(xué)位，現(xiàn)為博士研究生，主要研究方向為信號處理、目標(biāo)定位等。

CAO Jin-kunwasborn in Fuyang，AnhuiProvince，in 1980. He received the M.S.degree from Information Engineering University in 2006.He is currently working toward the Ph.D.degree.His research interests include target location，tracking and satellite orbit measurement.

Email：kingquencho＠gmail.com

TN 電子技術(shù)、通信技術(shù)類核心期刊表

Geo-location Using Dual-satellite′s TDOA and FDOA Data by Least Square Fusion Estimation

PENGHua-feng，XIAChang-xiong，CAO Jin-kun
（Southwest Electronics and Telecommunication Technology Research Institute，Chengdu 610041，China）

Aiming at the problem ofgeo-location with only one reference station，which has feature of low position precision when there aremeasurement error and ephemeris error，a Least Square（LS）estimationmethod is proposed here by fusingmulti-time′s observation，bywhich the location accuracy can be greatly improved without increasing the observation conditions.In condition of one reference station，when measurement error and ephemeris error are added，the one observation geo-location method and multi-observation geo-location method are studied.The influence of observation error and ephemeris error is analysed.The equation of error transfer from observation error and ephemeris error to geo-location error is derived.A weighted algorithm for LS estimation geo-location is proposed and compared with the one observation geo-location algorithm.By Monte -Carlo simulation，the analysis conclusion and weighted LSestimation geo-location algorithm are proved，and also the performance ofone-observation andmulti-observation algorithm is compared.Simulation result shows that theweighted LSmulti-observation geo-location estimator can greatly improve the location precision.The best result shows thatmulti-observation algorithm can be 10 times better than one-observation geo-location algorithm.

dual-satellite；time and frequency difference of arrival；location；estimation；data fusion

序號刊名序號刊名1電子學(xué)報23半導(dǎo)體技術(shù)2光電子、激光24激光雜志3液晶與顯示25光通信技術(shù)4紅外與激光工程26電路與系統(tǒng)學(xué)報5電波科學(xué)學(xué)報27電子元件與材料6紅外與毫米波學(xué)報28電子科技大學(xué)學(xué)報7電子與信息學(xué)報29應(yīng)用光學(xué)8通信學(xué)報30應(yīng)用激光9北京郵電大學(xué)學(xué)報31數(shù)據(jù)采集與處理10激光與紅外32光電子技術(shù)11西安電子科技大學(xué)學(xué)報33光通信研究12系統(tǒng)工程與電子技術(shù)34電子器件13現(xiàn)代雷達(dá)35電信科學(xué)14紅外技術(shù)36電訊技術(shù)15微電子學(xué)37電子技術(shù)應(yīng)用16半導(dǎo)體光電38電視技術(shù)17光電工程39壓電與聲光18微波學(xué)報40重慶郵電大學(xué)學(xué)報.自然科學(xué)版19激光技術(shù)41功能材料與器件學(xué)報20信號處理42南京郵電大學(xué)學(xué)報.自然科學(xué)版21激光與光電子學(xué)進(jìn)展43微納電子技術(shù)22固體電子學(xué)研究與進(jìn)展

TN971

A

10.3969／j.issn.1001-893x.2012.04.003

彭華峰（1979—），男，湖南新化人，2006年于四川大學(xué)獲碩士學(xué)位，現(xiàn)為博士研究生、工程師，主要研究方向為目標(biāo)定位、跟蹤和測軌等；

1001-893X（2012）04-0435-05

2011-12-21；

2012-02-17