姬新陽(yáng) 高 山 陳慶良 張海龍 宮福紅 范志鵬

(1.中國(guó)洛陽(yáng)電子裝備試驗(yàn)中心 洛陽(yáng) 471003;2.西安衛(wèi)星測(cè)控中心 西安 710043)

基于無(wú)人機(jī)GPS的測(cè)量雷達(dá)標(biāo)校方法研究

姬新陽(yáng)1高 山2陳慶良1張海龍1宮福紅1范志鵬1

(1.中國(guó)洛陽(yáng)電子裝備試驗(yàn)中心 洛陽(yáng) 471003;2.西安衛(wèi)星測(cè)控中心 西安 710043)

測(cè)量雷達(dá)任務(wù)前必須進(jìn)行標(biāo)校以消除系統(tǒng)誤差，通常使用方位標(biāo)、距離標(biāo)和光電綜合標(biāo)等標(biāo)校設(shè)施進(jìn)行標(biāo)校，但有些臨時(shí)陣地并不具備完整的標(biāo)校設(shè)施，此時(shí)采用常規(guī)標(biāo)校方法就存在一定的限制。文中通過(guò)對(duì)標(biāo)校原理進(jìn)行分析，設(shè)計(jì)了一種利用小型無(wú)人機(jī)加裝GPS設(shè)備對(duì)精測(cè)雷達(dá)進(jìn)行標(biāo)校的方法。通過(guò)模擬數(shù)據(jù)解算證明該方法有效、可行，且得到了該方法的適用航線。該標(biāo)校方法簡(jiǎn)便、快捷，具有不受場(chǎng)地和常規(guī)標(biāo)校條件限制的優(yōu)點(diǎn)。

標(biāo)校;無(wú)人機(jī);GPS;測(cè)量雷達(dá)

0 引言

通常，雷達(dá)測(cè)量精度的高低用其測(cè)量誤差的大小來(lái)表征和衡量，測(cè)量誤差是指測(cè)量值與真實(shí)值之間的偏差，測(cè)量誤差小即意味著測(cè)量精度高。精密跟蹤測(cè)量雷達(dá)的原始測(cè)量數(shù)據(jù)存在系統(tǒng)誤差和隨機(jī)誤差。為提高測(cè)量數(shù)據(jù)的準(zhǔn)確度，在數(shù)據(jù)處理中通過(guò)數(shù)據(jù)平滑濾波來(lái)抑制隨機(jī)噪聲;而系統(tǒng)誤差則通過(guò)標(biāo)校來(lái)修正。因此，標(biāo)校是精密跟蹤測(cè)雷達(dá)使用維護(hù)中的一項(xiàng)經(jīng)常性工作，準(zhǔn)確的標(biāo)校將有助于修正系統(tǒng)誤差，提高雷達(dá)測(cè)量精度。

目前雷達(dá)標(biāo)校方法主要有常規(guī)標(biāo)校、星體標(biāo)校和衛(wèi)星標(biāo)校三種。在不方便建立常規(guī)標(biāo)校設(shè)施，而又不具備星體標(biāo)校和衛(wèi)星標(biāo)校的條件時(shí)，如何進(jìn)行雷達(dá)標(biāo)校就成為一個(gè)亟待解決的問(wèn)題［1－6］。文中通過(guò)雷達(dá)跟蹤帶有GPS定位設(shè)備的小型無(wú)人機(jī)，由GPS設(shè)備提供的GPS數(shù)據(jù)作為真值對(duì)雷達(dá)的系統(tǒng)誤差進(jìn)行標(biāo)校，可有效解決該問(wèn)題，同時(shí)該方法對(duì)地形條件沒(méi)有要求且不需要常規(guī)標(biāo)校設(shè)施，使用范圍廣。

1 標(biāo)校原理

1.1 系統(tǒng)誤差模型

精測(cè)雷達(dá)在任務(wù)前通常通過(guò)標(biāo)校來(lái)減小系統(tǒng)誤差，提高雷達(dá)測(cè)量精度，常規(guī)標(biāo)校方法是通過(guò)對(duì)方位標(biāo)和光電綜合標(biāo)，進(jìn)行大盤水平測(cè)試等得到標(biāo)定參數(shù)后，通過(guò)系統(tǒng)誤差修正公式對(duì)測(cè)量數(shù)據(jù)進(jìn)行修正，消除系統(tǒng)誤差。分析時(shí)以方位角為例，俯仰角誤差與方位角誤差類似，不再分析。

雷達(dá)測(cè)量數(shù)據(jù)一次誤差信息中含有的主要誤差有:

(1)方位零值:A0。

(2)大盤不水平產(chǎn)生的方位角誤差:

ΔA= θMsin(A － AM)tanE= θMcosAMsinAtanE－θMsinAMcosAtanE，其中θM為最大不水平角度，AM為最大不水平方位角。

(3)方位軸與俯仰軸不正交產(chǎn)生的方位角誤差為:ΔA=δntanE，其中δM為方位軸與俯仰軸不正交度。

(4)光軸電軸不平行產(chǎn)生的方位角誤差為ΔA=ΔAGD·sec(E)，其中ΔAGD為方位上光電軸偏差。

(5)光軸機(jī)械軸不平行方位角產(chǎn)生誤差為ΔA=ΔAGJ·sec(E)，其中ΔAGJ為方位上光軸機(jī)械軸偏差。

根據(jù)以上誤差量可以得到方位角誤差模型如下:

其中:a0為方位零值;a1、a2為大盤不水平系數(shù);a3為方位軸、俯仰軸不正交系數(shù);a4為光電軸、光機(jī)軸和定向靈敏度總誤差。

1.2 解算方法

由誤差模型可知，如果誤差(ΔA，ΔE，ΔR)可以精確測(cè)量，就能得到一組數(shù)據(jù)值 ΔAi，，ΔEi，ΔRi，Ai，Ei，Ri(i=1，2，3，……，m)，將這些數(shù)據(jù)按照誤差模型進(jìn)行回歸分析，通過(guò)最小二乘法解算，即可求得各誤差系數(shù)［6－11］。用矩陣表述如下:

其中:a為誤差系數(shù)向量

根據(jù)求得的誤差系數(shù) a0、a1、a2、a3、a4，可以解算出方位零值、大盤不水平度、方位軸與俯仰軸不正交度、光電軸、光機(jī)軸和定向靈敏度等標(biāo)校參數(shù)，這樣就完成了雷達(dá)的標(biāo)校。

1.3 精度分析

文中標(biāo)校方法是利用GPS的測(cè)量數(shù)據(jù)提供真值，因此所得標(biāo)校參數(shù)的精度主要取決于GPS的精度和誤差模型，應(yīng)用差分GPS技術(shù)時(shí)的動(dòng)態(tài)GPS定位精度在0.1米以下，完全可以滿足雷達(dá)標(biāo)校要求。

2 標(biāo)校方法

將GPS定位設(shè)備放置在小型無(wú)人機(jī)上，無(wú)人機(jī)作為雷達(dá)校準(zhǔn)時(shí)的跟蹤目標(biāo)，GPS定位設(shè)備實(shí)時(shí)的發(fā)送GPS數(shù)據(jù)，地面接收設(shè)備實(shí)時(shí)記錄GPS數(shù)據(jù)。獲得GPS數(shù)據(jù)和雷達(dá)測(cè)量數(shù)據(jù)后，首先進(jìn)行坐標(biāo)轉(zhuǎn)換，然后進(jìn)行數(shù)據(jù)處理，就能得到雷達(dá)的系統(tǒng)誤差標(biāo)定參數(shù)。

2.1 坐標(biāo)轉(zhuǎn)換

由于GPS設(shè)備測(cè)量的是無(wú)人機(jī)的位置，并不是雷達(dá)的直接觀測(cè)量，在校準(zhǔn)前應(yīng)首先進(jìn)行坐標(biāo)變換，將GPS數(shù)據(jù)轉(zhuǎn)換到雷達(dá)站心極坐標(biāo)下，再進(jìn)行校準(zhǔn)。

轉(zhuǎn)換過(guò)程為:由地心大地坐標(biāo)系依次轉(zhuǎn)至地心空間直角坐標(biāo)系、站心空間直角坐標(biāo)系和站心極坐標(biāo)系［11－12］，地心空間直角坐標(biāo)系如圖1所示。

站心空間直角坐標(biāo)系的定義為:原點(diǎn)OR位于雷達(dá)站點(diǎn)T0，ZR軸與OR(T0)的橢球法線相重合，XR軸垂直于ZR軸指向橢球的短軸，而YR軸垂直于XRORZR平面，構(gòu)成左手坐標(biāo)系，目標(biāo)點(diǎn)TT的坐標(biāo)表示為(xT，yT，zT)。如圖2 所示。

站心極坐標(biāo)系定義為:以XRORZR平面為基準(zhǔn)，以O(shè)R為極點(diǎn)，以XR軸為極軸。目標(biāo)點(diǎn)TT的坐標(biāo)表示為(rT，αT，βT)，rT為目標(biāo)點(diǎn) TT距離觀測(cè)點(diǎn) OR的距離。αT為目標(biāo)點(diǎn)TT在站心極坐標(biāo)系的方位角，βT為目標(biāo)點(diǎn)TT在站心極坐標(biāo)系的俯仰角。雷達(dá)探測(cè)目標(biāo)得到的信息一般為站心極坐標(biāo)系下的坐標(biāo)。任一點(diǎn)T在地心大地坐標(biāo)系中的坐標(biāo)可表示為(B，L，H)，任一點(diǎn)T在地心空間直角坐標(biāo)系中的坐標(biāo)也可表示為(X，Y，Z)，這兩種坐標(biāo)相互可以轉(zhuǎn)換［13］。其轉(zhuǎn)換關(guān)系如下式。

其中，a為地球長(zhǎng)半徑:a=6378137m;e為橢球的第一偏心率:e2=0.00669437999013。

使用公式(8)可以將雷達(dá)站點(diǎn)T0的地心大地坐標(biāo)(B0，L0，H0)和目標(biāo)點(diǎn)TT的地心大地坐標(biāo)(BT，LT，HT)轉(zhuǎn)換為地心空間直角坐標(biāo)系的坐標(biāo)(X0，Y0，Z0) 和(XT，YT，ZT)。

使用公式(9)，根據(jù)雷達(dá)站點(diǎn)T0和目標(biāo)點(diǎn)TT的地心空間直角坐標(biāo)系的坐標(biāo)(X0，Y0，Z0) 和(XT，YT，ZT)可以求出目標(biāo)點(diǎn)TT在站心空間直角坐標(biāo)系中的坐標(biāo)(xT，yT，zT)。

由于雷達(dá)測(cè)得的目標(biāo)信息是站心極坐標(biāo)系下的坐標(biāo)，應(yīng)當(dāng)使用公式(11)將目標(biāo)點(diǎn)TT的站心空間直角坐標(biāo)系坐標(biāo)(XT，YT，ZT)轉(zhuǎn)換到站心極坐標(biāo)系下的坐標(biāo)(rT，αT，βT):

至此，已經(jīng)完成了地心大地坐標(biāo)系到站心極坐標(biāo)系的轉(zhuǎn)換。

2.2 測(cè)量數(shù)據(jù)處理

2.2.1 測(cè)量數(shù)據(jù)隨機(jī)誤差平滑

雷達(dá)測(cè)量數(shù)據(jù)通常帶有隨機(jī)誤差，需要對(duì)測(cè)量數(shù)據(jù)進(jìn)行平滑處理，以得到滿意的結(jié)果。平滑計(jì)算的目的即消除隨機(jī)誤差，在保證一定精度的條件下，壓縮測(cè)量數(shù)據(jù)量，本文采用多項(xiàng)式平滑方法［14］:

設(shè)一組測(cè)量數(shù)據(jù){ti，xi}(i=1，2…，n)(ti為采樣時(shí)間，xi為測(cè)量數(shù)據(jù))，設(shè)多項(xiàng)式為:

通過(guò)最小二乘法原理，使偏差最小，即:

根據(jù)求極值方法，將(12)帶入(13)，并求偏導(dǎo)數(shù)，得出:

一般情況下多項(xiàng)式次數(shù)取3次。

2.2.2 異常值的檢測(cè)與剔除

本文采用外推擬和法進(jìn)行野值的識(shí)別和檢驗(yàn)。外推擬和法通常取前4個(gè)或5個(gè)連續(xù)正常的測(cè)量數(shù)據(jù)，采用一階或二階多項(xiàng)式進(jìn)行計(jì)算。δ為相應(yīng)測(cè)量量測(cè)量誤差的均方差。假設(shè)連續(xù)5個(gè)測(cè)量數(shù)據(jù)為:xi－5、xi－4、xi－3、xi－2、xi－1，由最小二乘線形估計(jì)外推獲

2.2.3 拉格朗日插值方法

由于雷達(dá)測(cè)量數(shù)據(jù)與GPS數(shù)據(jù)的采樣時(shí)刻和采樣頻率并不一致，因此，需要將雷達(dá)測(cè)量數(shù)據(jù)插值到GPS數(shù)據(jù)對(duì)應(yīng)時(shí)刻，擬采用10階拉格朗日插值法進(jìn)行計(jì)算，計(jì)算公式如下:

當(dāng)選取的點(diǎn)數(shù)少于10個(gè)時(shí)，舍去該待插值點(diǎn)。

2.3 標(biāo)定參數(shù)解算

標(biāo)定參數(shù)解算采用最小二乘法進(jìn)行，通過(guò)Matlab編程實(shí)現(xiàn)。解算時(shí)輸入雷達(dá)測(cè)量數(shù)據(jù)和GPS數(shù)據(jù)，首先進(jìn)行野值識(shí)別和處理、數(shù)據(jù)平滑和插值，其次對(duì)GPS數(shù)據(jù)進(jìn)行坐標(biāo)轉(zhuǎn)換，最后進(jìn)行標(biāo)定參數(shù)解算，并利用GPS真值對(duì)得到的標(biāo)定參數(shù)進(jìn)行檢驗(yàn)，滿足指標(biāo)要求后，標(biāo)校過(guò)程結(jié)束。

3 模擬數(shù)據(jù)解算結(jié)果

由于各種條件限制，缺少以GPS為真值的測(cè)量數(shù)據(jù)，文中采用模擬數(shù)據(jù)對(duì)標(biāo)校方法進(jìn)行了檢驗(yàn)。模擬數(shù)據(jù)沒(méi)有考慮隨機(jī)誤差，通過(guò)解算得到了在理想情況下該標(biāo)校方法的適用范圍。

根據(jù)對(duì)誤差模型的分析，參與解算的數(shù)據(jù)對(duì)解算結(jié)果影響非常大，一是方位角必須在0°～360°四個(gè)象限內(nèi)均有數(shù)據(jù)，才能解算出大盤水平誤差;二是由于最后一項(xiàng)誤差為Z·sec(E)，當(dāng)仰角E為3°時(shí)，sec(E)為1.00137，這樣仰角E在3°以下時(shí)，最后一項(xiàng)誤差Z·sec(E)與零值a0存在合并到一起的可能，也將會(huì)導(dǎo)致結(jié)算結(jié)果不準(zhǔn)確。因此，分兩種情況進(jìn)行模擬解算，一是保持仰角數(shù)據(jù)線性分布在3°～40°不變，改變方位角數(shù)據(jù)的范圍，分析解算結(jié)果是否準(zhǔn)確;二是保持方位角數(shù)據(jù)線性分布在0°～360°不變，改變仰角數(shù)據(jù)的范圍，分析解算結(jié)果是否準(zhǔn)確。模擬解算結(jié)果如表1和表2所示。

通過(guò)對(duì)模擬數(shù)據(jù)的解算結(jié)果進(jìn)行分析，可以得到該標(biāo)校方法的適用范圍，主要有以下三點(diǎn):

(1)方位角數(shù)據(jù)必須均勻分布在0°～360°四個(gè)象限內(nèi)，每個(gè)象限內(nèi)不少于兩組數(shù)據(jù)，每組數(shù)據(jù)在40個(gè)點(diǎn)以上;

(2)仰角數(shù)據(jù)必須是變化的，最小變化范圍為 17°;

(3)仰角最低為6°;

根據(jù)以上三點(diǎn)就可以設(shè)計(jì)出該標(biāo)校方法的航線范圍，一是仰角數(shù)據(jù)必須是變化的，同時(shí)仰角不能低于6°，范圍至少為6°～23°;二是方位角數(shù)據(jù)可以為分段的，但是最少為八段，而且均勻分布在四個(gè)象限，每段數(shù)據(jù)方位角變化范圍最小為10度，在條件允許的情況下方位角數(shù)據(jù)范圍越大解算結(jié)果越準(zhǔn)確。

在實(shí)際標(biāo)校時(shí)可以采用小型無(wú)人機(jī)以雷達(dá)為中心，做圓周飛行，同時(shí)保持飛行高度連續(xù)變化，飛完一周后即可利用獲得的數(shù)據(jù)計(jì)算出雷達(dá)的系統(tǒng)誤差標(biāo)定參數(shù)，完成雷達(dá)的標(biāo)校工作。

表1 保持仰角數(shù)據(jù)線性分布在3°～40°時(shí)的解算結(jié)果

表2 保持方位角數(shù)據(jù)線性分布在四個(gè)象限時(shí)的解算結(jié)果

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

通過(guò)利用小型無(wú)人機(jī)攜帶GPS設(shè)備對(duì)雷達(dá)進(jìn)行標(biāo)校，具有不受場(chǎng)地和常規(guī)標(biāo)校條件限制的優(yōu)點(diǎn)，隨著小型無(wú)人機(jī)技術(shù)的發(fā)展和成熟，該標(biāo)校方法對(duì)于滿足特定場(chǎng)地條件下雷達(dá)的標(biāo)校具有一定的應(yīng)用前景。

［1］朱起悅.應(yīng)用差分GPS技術(shù)進(jìn)行雷達(dá)標(biāo)校［J］. 電訊技術(shù)，2006，46(1):108 －110.

［2］吳三元，侯志楠，王瑞恒.基于GPS的雷達(dá)標(biāo)校方法［J］.雷達(dá)與對(duì)抗，2010，30(1):4－6.

［3］杜娟，朱華邦.基于 GPS技術(shù)的艦載探測(cè)系統(tǒng)標(biāo)校方法研究［J］.火控雷達(dá)技術(shù)，2013.，42(1):17－20.

［4］姚景順，楊世興.艦載雷達(dá)的動(dòng)態(tài)標(biāo)校［J］.火力與指揮控制，2008，33(3):128－130.

［5］潘紹仁，察豪.基于AIS的艦載雷達(dá)標(biāo)校方法研究［J］.艦船科學(xué)技術(shù)，2009，31(6):101－104，127.

［6］張舒彥，鮑慶龍，林財(cái)永，陳曾平.一種非合作雙基地雷達(dá)目標(biāo)跟蹤方法科學(xué)技術(shù)與工程［J］.2015，15(29):153 －157.

［7］豐少偉，黃雙華，楊云生.基于差分GPS技術(shù)的艦載雷達(dá)通用標(biāo)校方法［J］.現(xiàn)代雷達(dá)，2012，34(5):14 －16.

［8］趙罄，王世峰，佟首峰等.艦載雷達(dá)零位標(biāo)校新法［J］.火力與指揮控制，2010，35(2):169－171.

［9］羅軍，商允力，曾浩，等.應(yīng)用AIS信息進(jìn)行岸基對(duì)海雷達(dá)標(biāo)校［J］.電訊技術(shù)，2009，49(9):87－89.

［10］姚兆寧.艦載精密測(cè)量雷達(dá)星體標(biāo)校方法及應(yīng)用［J］.現(xiàn)代雷達(dá)，1999，(4):8－12.

［11］苑文亮，唐小明，朱洪偉等.基于 ADS－B數(shù)據(jù)的雷達(dá)標(biāo)校新方法［J］.艦船電子工程，2010，30(3):147 －150.

［12］邊少鋒，李文魁.衛(wèi)星導(dǎo)航系統(tǒng)概論［M］.北京:電子工業(yè)出版社，2005，02.

［13］寇艷紅.GPS原理與應(yīng)用(第二版)［M］.北京:電子工業(yè)出版社，2007，07.

［14］郭軍海.基于最優(yōu)節(jié)點(diǎn)樣條逼近的觀測(cè)數(shù)據(jù)平滑方法［J］.中國(guó)空間科學(xué)技術(shù)，2000，20(3):43－48.

Calibration Method of Measurement Radar Based on UAV with GPS

Ji Xinyang1，Gao Shan2，Chen Qingliang1，Zhang Hailong1，Gong Fuhong1，F(xiàn)an Zhipeng1
(1.Luoyang Electronic Equipment Test Center of China，Luoyang 471003;2.Xi'an Satellite Control Center，Xi'an 710043)

Measurement radar must be calibrated to eliminate system error before execute the mission，and azimuth，range and optoelectronic integrated calibration facilities are generally used for calibration.But to some temporary positions without complete calibration facilities，calibration by using conventional calibration method is restricted.Through analysis of calibration principle，a calibration method for measurement radar is designed by using UAV with GPS.Simulation result proved that the method is effective and feasible，and suitable route can be obtained.The calibration method is simple and convenient，and is not restricted by site and conventional calibration conditions.

calibration;UAV;GPS;measurement radar

TN953+.5

A

1008-8652(2017)01-073-06

2016-12-11

姬新陽(yáng)(1989－)，男，工程師。研究方向?yàn)闇y(cè)量雷達(dá)技術(shù)。