王前東

（中國西南電子技術(shù)研究所，成都610036）

雙坐標(biāo)系多參量目標(biāo)的融合定位?

王前東

（中國西南電子技術(shù)研究所，成都610036）

雷達(dá)探測的距離、方位角與目標(biāo)應(yīng)答的高度屬于兩個(gè)不同坐標(biāo)系，無法使用坐標(biāo)變換公式進(jìn)行定位，為此提出了一種迭代求解方法。首先采用幾何法求出雷達(dá)探測坐標(biāo)系的俯仰角初始值，然后直接利用坐標(biāo)變換公式給出求解的迭代算法。Matlab仿真數(shù)值實(shí)例表明了經(jīng)度和緯度的定位誤差均小于0.000 02°。

艦載雷達(dá)；融合定位；雙坐標(biāo)系；坐標(biāo)變換；迭代原理

1 引言

一般對空艦載雷達(dá)是三維雷達(dá)，給出目標(biāo)在艦船球坐標(biāo)系下的距離、方位角、俯仰角，三參量屬于同一坐標(biāo)系，可以用坐標(biāo)變換公式對目標(biāo)進(jìn)行直接定位［1-3］。但實(shí)際應(yīng)用中，常常遇到的艦載雷達(dá)是二維雷達(dá)，可以提供目標(biāo)在艦船球坐標(biāo)系下的距離、方位角，不能獲得俯仰角；同時(shí)用其它方法可以獲得目標(biāo)在WGS-84大地坐標(biāo)系下的高度，不能獲得經(jīng)度與緯度。這樣在艦載雷達(dá)截獲目標(biāo)中構(gòu)成了由艦船球坐標(biāo)系下的距離、方位角、WGS-84大地坐標(biāo)系下的高度共3個(gè)坐標(biāo)組成的雙坐標(biāo)系，不能直接利用已有的坐標(biāo)變換公式對目標(biāo)進(jìn)行定位，公開文獻(xiàn)沒有查到對這種雙坐標(biāo)系進(jìn)行定位的方法。為此，本文提出了一種對雙坐標(biāo)系參量進(jìn)行融合定位的新方法。

2 坐標(biāo)變換的基礎(chǔ)知識(shí)［1-5］

2．1 坐標(biāo)系定義

要對目標(biāo)進(jìn)行定位，需要利用以下幾個(gè)常用坐標(biāo)系與相關(guān)參數(shù)。

（1）WGS-84直角坐標(biāo)系

即美國國防部1984年世界大地坐標(biāo)系（GPS采用WGS-84）。其坐標(biāo)原點(diǎn)是地球質(zhì)心，Z軸指向BIH1984.0定義的協(xié)議地極（CTP）方向，X軸指向BIH1984.0的零度子午面和CTP赤道的交點(diǎn)，Y軸和Z、X軸構(gòu)成右手坐標(biāo)系。

（2）WGS-84大地坐標(biāo)系

空間一點(diǎn)P的大地坐標(biāo)用大地緯度M、大地經(jīng)度L和大地高度H表示。其中P點(diǎn)的大地子午面和起始大地子午面所構(gòu)成的二面角L是P點(diǎn)的大地經(jīng)度；P點(diǎn)的法線和赤道面的夾角M是P點(diǎn)的大地緯度；P點(diǎn)沿法線到地球球面的距離H是大地高。

（3）艦船地理坐標(biāo)系

以艦船質(zhì)心位置為原點(diǎn)，X軸正方向指向正東，Y軸正方向指向正北，Z軸正方向垂直X-Y平面向上。

（4）艦船地理球坐標(biāo)系

方位角為目標(biāo)在水平面的投影與正北的夾角，順時(shí)針為正；俯仰角為目標(biāo)在水平面的投影與目標(biāo)到原點(diǎn)的連線間的夾角，向上為正；距離為目標(biāo)與艦船質(zhì)心之間的距離。

（5）艦船直角坐標(biāo)系

以艦船的質(zhì)心位置為坐標(biāo)原點(diǎn)，X-Y平面平行于艦主基準(zhǔn)面，Y軸正方向指向艦艏方向，X軸正方向指向右舷90°方向，Z軸正方向垂直向上。

（6）艦船甲板球坐標(biāo)系

方位角為目標(biāo)在艦主基準(zhǔn)面的投影離艦艏方向的夾角；俯仰角為目標(biāo)在艦主基準(zhǔn)面的投影與目標(biāo)到原點(diǎn)的連線間的夾角，向上為正。

（7）艏向角

艦船艏艉線在地理水平面內(nèi)的投影與真北的夾角。

（8）縱搖角

艦船艏艉線相對于地理坐標(biāo)系在其垂直面內(nèi)的搖擺角。

（9）橫搖角

艦船坐標(biāo)系相對于地理坐標(biāo)系在其肋骨面內(nèi)的搖擺角。

2．2 單坐標(biāo)系坐標(biāo)變換

單坐標(biāo)系坐標(biāo)變換步驟如下。

（1）從艦船甲板球坐標(biāo)距離、方位角、俯仰角組成的目標(biāo)位置S0=［rαθ］′，到艦船甲板直角坐標(biāo)位置S1=［x1y1z1］′的坐標(biāo)變換

（2）從艦船直角位置S1=［x1y1z1］′到艦船地理直角位置S2=［X2Y2Z2］′的坐標(biāo)變換

設(shè)艦船的姿態(tài)為橫搖角γ、艏向角β、縱搖角為ε，令

則有變換：

（3）目標(biāo)由艦船地理位置S2=［X2Y2Z2］′到地心坐標(biāo)位置S3=［X3Y3Z3］′的坐標(biāo)變換

設(shè)艦船在地心坐標(biāo)下的位置為S0=［X0Y0Z0］′，設(shè)艦船WGS-84大地坐標(biāo)經(jīng)度為l、緯度為m、高度為h，令

則有變換：

（4）目標(biāo)由WGS-84地心直角坐標(biāo)S3=［X3Y3Z3］′到WGS-84大地坐標(biāo)即經(jīng)度L、緯度M、高度H的變換

在上式中，求M和H的迭代方法是：迭代開始時(shí)，設(shè)：

隨后，每次迭代按下述公式進(jìn)行：

直至（Mi-Mi-1）和（Hi-Hi-1）小于某一所要求的限值為止。例如，在保證H的計(jì)算精度為0.001 m和M的計(jì)算精度為0.000 01″的情況下，一般僅迭代4次左右。

如果高度H已知，則坐標(biāo)變換變?yōu)?/p>

在上式中，求M迭代方法是：迭代開始時(shí)，設(shè)：

隨后，每次迭代按下述公式進(jìn)行：

（5）從WGS-84大地坐標(biāo)即經(jīng)度L、緯度M、高度H到WGS-84地心直角坐標(biāo)S3=［X3Y3Z3］′的變換

式中，e為橢球第一偏心率，a為橢球長半徑，N為橢球卯酉圈曲率半徑。

（6）目標(biāo)由地心坐標(biāo)位置S3=［X3Y3Z3］′到艦船地理位置S2=［X2Y2Z2］′的坐標(biāo)變換

設(shè)艦船在地心坐標(biāo)下的位置為S0=［X0Y0Z0］′，設(shè)艦船WGS-84大地坐標(biāo)經(jīng)度為l、緯度為m、高度為h，令

則有變換：

（7）從艦船地理直角位置S2=［X2Y2Z2］′到艦船直角位置S1=［x1y1z1］′的坐標(biāo)變換

設(shè)艦船的橫搖角、艏向角、縱搖角分別為γ、β、ε，令

則有變換：

（8）從艦船直角坐標(biāo)位置S1=［x1y1z1］′到艦船球坐標(biāo)位置S0=［rαθ］′的坐標(biāo)變換

3 基于迭代原理的雙坐標(biāo)系目標(biāo)融合定位

3．1 融合定位的迭代方法

設(shè)已知艦船的姿態(tài)為橫搖角γ，艏向角β，縱搖角為ε，設(shè)艦船在地心坐標(biāo)下的位置為S0=［x0y0z0］′。利用坐標(biāo)變換步驟（1）～（4），可以求出目標(biāo)由艦船甲板球坐標(biāo)位置S0=［rαθ］′到WGS-84大地坐標(biāo)的經(jīng)度L、緯度M的坐標(biāo)變換。設(shè)其中求解經(jīng)度L的復(fù)合函數(shù)為

設(shè)求解緯度M的復(fù)合函數(shù)為

同理，可以用坐標(biāo)變換步驟（5）～（8）求出目標(biāo)由WGS-84大地坐標(biāo)的經(jīng)度L、緯度M、高度H到艦船甲板球坐標(biāo)系下的俯仰角θ的坐標(biāo)變換。設(shè)求俯仰角θ的復(fù)合函數(shù)為

利用公式（17）～（19），可以求出俯仰角θ滿足如下方程：

此方程可以記為如下形式：

式中，gr，α，H為f1、f2、f3的復(fù)合函數(shù)，根據(jù)此方程容易知道求俯仰角θ的一般迭代方法：

3．2 幾何法求迭代初值

假定地球?yàn)橐磺蝮w，半徑為R，球心為O。目標(biāo)P的高度為H，目標(biāo)距離艦船M的距離為r，如圖1所示。

在三角形OPM中可以求出目標(biāo)在艦船地理球坐標(biāo)系下的俯仰角θ1滿足如下方程：

設(shè)艦船的橫搖角、縱搖角分別為γ、ε，從艦船甲板球坐標(biāo)系下的方位角為α、俯仰角為θ，則有如下方程［3］：

聯(lián)合公式（23）、（24）可求出俯仰角初值θ。

3．3 融合定位的求解步驟

迭代開始時(shí)，利用公式（23）、（24）求出θ的初值，代入距離r、方位角α，利用坐標(biāo)變換（1）～（4）與艦船位置和姿態(tài)參數(shù)，可以求出目標(biāo)的經(jīng)度L0、緯度M0；隨后，令i=i+1，每次迭代按下述公式進(jìn)行：

（1）代入經(jīng)度Li、緯度Mi、高度H，利用坐標(biāo)變換（5）～（8）可以求出俯仰角θi+1；

（2）代入距離r、方位角α、俯仰角θi+1，利用坐標(biāo)變換（1）～（4）可以求出經(jīng)度Li+1、緯度Mi+1；

令i=i+1，重復(fù)上述過程直到滿足所需精度。

4 Matlab軟件仿真試驗(yàn)結(jié)果

設(shè)艦船位置經(jīng)度為l=120，緯度為m=20，高度為h=0，艦船的橫搖角、艏向角、縱搖角分別為γ =10、β=10、ε=5，距離r=50 000，方位角α=330，仿真結(jié)果如表1所示。

從表1可以看出：真實(shí)的俯仰角從0°增加到27°時(shí)，本文的坐標(biāo)變換經(jīng)度與緯度定位誤差小于0.000 02°，俯仰角誤差小于0.005°，能滿足定位需求。

5 結(jié)論

本文針對二維艦載雷達(dá)目標(biāo)的坐標(biāo)定位問題進(jìn)行了探索，利用艦船球坐標(biāo)系下的距離、方位角、WGS-84大地坐標(biāo)系下的高度共3個(gè)坐標(biāo)組成的雙坐標(biāo)系，詳細(xì)推導(dǎo)了相應(yīng)的坐標(biāo)定位過程，并用Matlab軟件進(jìn)行了仿真驗(yàn)證。仿真結(jié)果表明：坐標(biāo)變換經(jīng)度與緯度定位誤差小于0.000 02°，俯仰角誤差小于0.005°。

本文采用的是迭代法，該方法可以推廣應(yīng)用到其它雙坐標(biāo)系中多參量的定位。實(shí)際應(yīng)用中，該方法計(jì)算量較大，應(yīng)根據(jù)實(shí)際需要進(jìn)行調(diào)整，如減少迭代次數(shù)以減少計(jì)算量，或優(yōu)化迭代初始值和增加迭代次數(shù)以提高定位精度等。

［1］曹正才.艦載雷達(dá)常用穩(wěn)定方式坐標(biāo)變換［J］.雷達(dá)與對抗，2010，30（1）：47-52.

CAO Zheng-cai.The commonly used coordinate transformation of stabilization modes for shipborne radars［J］.Radar＆ECM，2010，30（1）：47-52.（in Chinese）

［2］楊帆，李強(qiáng)，王衛(wèi)疆.坐標(biāo)轉(zhuǎn)換在機(jī)載預(yù)警雷達(dá)情報(bào)關(guān)聯(lián)中的應(yīng)用［J］.空軍雷達(dá)學(xué)院學(xué)報(bào)，2010（4）：28-31.

YANG Fan，LI Qiang，WANG Wei-jiang.Application of Coordinate Transformation too Intelligence Correlation of Airborne Early Warning Radar［J］.Journal of Air Force Radar A-cademy，2010（4）：28-31.（in Chinese）

［3］方成一，王振旺.艦載雷達(dá)天線電子穩(wěn)定方程的推導(dǎo)方法［J］.雷達(dá)與對抗，1999（2）：68-74.

FANG Cheng-yi，WANG Zhen-wang.Derivation of electronic steadiness equations of ship-borne radar［J］.Radar＆ECM，1999（2）：68-74.（in Chinese）

［4］常軍.機(jī)載雷達(dá)目標(biāo)的大地坐標(biāo)定位［J］.電訊技術(shù)，2003，43（2）：97-100.

CHANG Jun.Positioning of Airborne Radar Targets in 54 Coordinate［J］.Telecommunication Engineering，2003，43（2）：97 -100.（in Chinese）

［5］吳順君，梅曉春.雷達(dá)信號(hào)處理和數(shù)據(jù)處理技術(shù)［M］.北京：電子工業(yè)出版社，2008.

WU Shun-jun，MEI Xiao-chun.Radar signal processing and data processing［M］.Beijing：Publishing House of Electronic Industry，2008.（in Chinese）

WANG Qian-dong was born in Xichong，Sichuan Province，in 1977.He received the M.S.degree from Sichuan University in 2004.He is now an engineer.His research direction is data fusion in intelligence system.

Email：wangqiandong＠sohu.com

《電訊技術(shù)》征稿啟事

《電訊技術(shù)》（月刊）創(chuàng)刊于1958年，由中國西南電子技術(shù)研究所主辦，系國內(nèi)外公開發(fā)行的、理論與應(yīng)用相結(jié)合的綜合性電子專業(yè)科技刊物，為中文核心期刊。目前，已被英國IEE《科學(xué)文摘（SA）》INSPEC、美國《劍橋科學(xué)文摘（CSA）》、波蘭《哥白尼索引（IC）》等國外知名數(shù)據(jù)和國內(nèi)多個(gè)中文數(shù)據(jù)庫收錄。

本刊主要刊登涉及下列應(yīng)用方向和技術(shù)領(lǐng)域的述評(píng)、論文、新概念新技術(shù)新產(chǎn)品介紹：

·電子系統(tǒng)工程 ·通信 ·導(dǎo)航 ·識(shí)別

·飛行器測控 ·衛(wèi)星應(yīng)用 ·雷達(dá) ·信息戰(zhàn)

·共性技術(shù)（包括天線、射頻電路、信號(hào)處理、信息處理、監(jiān)視與控制、時(shí)間與頻率、先進(jìn)制造、電磁兼容等）。

本刊欄目有：系統(tǒng)總體技術(shù)、信號(hào)與信息處理技術(shù)、計(jì)算機(jī)、網(wǎng)絡(luò)及其應(yīng)用技術(shù)、信道技術(shù)、先進(jìn)制造技術(shù)、基礎(chǔ)技術(shù)、述評(píng)與展望。

歡迎業(yè)內(nèi)學(xué)者、專家及科技人員踴躍投稿。

來稿要求及注意事項(xiàng)：

（1）文稿務(wù)必主題明確，論述合理，邏輯嚴(yán)謹(jǐn)，數(shù)據(jù)可靠，敘述清楚，文字精煉。內(nèi)容應(yīng)保守國家機(jī)密，并提供所在單位的保密審查證明，引用他人作品應(yīng)給出來源。

（2）文稿一般不應(yīng)超過6 000字，綜述稿不超過8 000字。稿件應(yīng)附英文題名、作者名、單位名、摘要和關(guān)鍵詞，基金項(xiàng)目應(yīng)注明項(xiàng)目編號(hào)。中文題名一般不超過20個(gè)漢字，必要時(shí)可加副標(biāo)題。

（3）摘要應(yīng)包括目的、方法、結(jié)果和結(jié)論四要素，即用簡潔的語言說明文章要解決的問題，主要工作過程及所采用的技術(shù)手段和方法，研究所獲得的實(shí)驗(yàn)數(shù)據(jù)、結(jié)果及其意義。篇幅以200～300字為宜。

（4）關(guān)鍵詞以3～8個(gè)為宜。為便于文獻(xiàn)檢索，應(yīng)盡可能根據(jù)《中國圖書館分類法（第四版）》提供中圖分類號(hào)。

（5）文中涉及的物理量和計(jì)量單位應(yīng)符合國家有關(guān)標(biāo)準(zhǔn)。計(jì)量單位請用GB3100-3102-93《量和單位》規(guī)定的法定計(jì)量單位。注意區(qū)分各物理量符號(hào)的文種、大、小寫、正斜體（矢量和矩陣用黑斜體）、上、下角標(biāo)等。

（6）插圖和表格只給出必要的，且應(yīng)有圖題和表題。插圖最好采用計(jì)算機(jī)制作。照片以黑白為佳，也可采用掃描的電子文檔（精度高于400 dpi，tif、jpg、psd等格式均可）。

（7）文稿中引用他人的成果，務(wù)請寫明原作者姓名、題名、來源，一并在參考文獻(xiàn)中給出，并在正文中相應(yīng)位置進(jìn)行標(biāo)示，否則責(zé)任由來稿人自負(fù)。參考文獻(xiàn)只擇主要的，未公開發(fā)表的文獻(xiàn)請勿列入?yún)⒖嘉墨I(xiàn)。書寫格式請參見GB／T 7714-2005文后參考文獻(xiàn)著錄規(guī)則。

（8）投稿郵箱：dianxunjishu＠china.com，并盡量同時(shí)提供Word和pdf文檔，無需另寄打印稿。來稿請注明作者詳細(xì)通信地址、聯(lián)系電話和有效電子郵箱，并注明擬投欄目。

（9）本刊編輯部將在3個(gè)月之內(nèi)對來稿作出取舍，如逾期未收到處理意見或刊用通知，作者有權(quán)對稿件另行處理。稿件一經(jīng)刊用，本刊將酌情從優(yōu)支付稿酬并贈(zèng)送當(dāng)期樣刊，本刊支付的稿酬中已包含作者著作權(quán)使用費(fèi)。請勿一稿多投，否則后果自負(fù)。

Fusion Positioning of Multi-parameter Target in Double Coordinates

WANG Qian-dong
（Southwest China Institute of Electronic Technology，Chengdu 610036，China）

As the detection distance，azimuth and altitude are in different two coordinates，the positioning of the target can not be realized by coordinate transformation.So，an iterative method is presented for fusion positioning.Firstly，the initial values of elevation of spherical coordinate are solved by geometric arithmetic.Secondly，the iterative scheme is formed by coordinate transformation.A numerical example in Matlab software shows the positioning error of both the longitude and the latitude is less than 0.00002 degree.

shipborne radar；fusion positioning；double coordinates；coordinate transformation；iterative principle

TN957.52

A

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

王前東（1977—），男，四川西充人，2004年于四川大學(xué)獲碩士學(xué)位，現(xiàn)為工程師，主要從事情報(bào)數(shù)據(jù)融合處理研究。

1001-893X（2011）06-0066-05

2011-03-22；

2011-05-13