陳 磊，方 青

(中國電子科技集團(tuán)公司第38研究所,合肥 230031)

?

一種基于波束增益差的兩坐標(biāo)雷達(dá)估高方法

陳 磊，方 青

(中國電子科技集團(tuán)公司第38研究所,合肥 230031)

兩坐標(biāo)雷達(dá)只能提供探測目標(biāo)的距離、方位信息，不能獲取目標(biāo)的高度或俯仰信息。為了估算目標(biāo)的高度，提出了一種基于波束增益差的兩坐標(biāo)雷達(dá)估高方法，根據(jù)雷達(dá)天線的“仰角-增益”測試數(shù)據(jù)，利用目標(biāo)點跡屬性中的高低波束增益差，求出目標(biāo)的仰角，進(jìn)而估算目標(biāo)的高度信息；同時采用加權(quán)平均值的方法平滑目標(biāo)的估高值。該方法簡單實用，方便工程實現(xiàn)，實驗中通過與目標(biāo)二次雷達(dá)的高度信息的對比，以驗證該方法的有效性。

兩坐標(biāo)雷達(dá)；波束增益差；估高算法；高度平滑

0 引 言

兩坐標(biāo)雷達(dá)，顧名思義，只能提供探測目標(biāo)2個維度的坐標(biāo)信息(距離和方位信息)。要獲取目標(biāo)的高度或俯仰信息，可采用以下幾種方式：(1)多部兩坐標(biāo)雷達(dá)組網(wǎng)測高的方式[1-3]；(2)兩坐標(biāo)雷達(dá)搭配測高雷達(dá)[4];(3)三坐標(biāo)雷達(dá)直接測高[5-7]。雷達(dá)組網(wǎng)測高利用多部雷達(dá)融合信息(例如距離和角度)來估計目標(biāo)的高度，進(jìn)行目標(biāo)三維定位，但具體實現(xiàn)時，一方面需要多部雷達(dá)的組網(wǎng)融合，另一方面系統(tǒng)的工作依賴于各雷達(dá)站之間的網(wǎng)絡(luò)通信情況，可靠性不高。兩坐標(biāo)搭配測高雷達(dá)的工作方式，測量過程較為復(fù)雜，實時性不高，且需要額外的測高雷達(dá)設(shè)備的支持。三坐標(biāo)雷達(dá)通常采用的是多波束比幅測高的方法，雖然可以實時、精確地測算目標(biāo)高度，但三坐標(biāo)雷達(dá)的生產(chǎn)和維護(hù)成本更高，技術(shù)更為復(fù)雜。因此，單部兩坐標(biāo)雷達(dá)測高的問題開始吸引了相關(guān)學(xué)者的研究，并出現(xiàn)了一些研究方法[8-9]。矯志寧等人[8]提出了一種兩坐標(biāo)雷達(dá)的測高方法研究，利用目標(biāo)幅度和距離變化信息進(jìn)行目標(biāo)估高，但該方法僅適用于飛行高度比較穩(wěn)定的目標(biāo)，而且需要長時間的積累，不具備普適性; 嚴(yán)俊坤等人[9]提出的基于單雷達(dá)的兩坐標(biāo)三維定位方法，是在假設(shè)目標(biāo)勻速運動的前提下，利用目標(biāo)的動態(tài)模型，結(jié)合兩坐標(biāo)雷達(dá)關(guān)于目標(biāo)多個時刻的測量信息對目標(biāo)高度進(jìn)行最大似然估計。該方法也具有一定的局限性，且運算較復(fù)雜，又僅是仿真實驗的驗證。

針對以上問題，本文依托于某XX波段航空管制一次雷達(dá)項目，提出了一種基于波束增益差的兩坐標(biāo)雷達(dá)估高算法，根據(jù)雷達(dá)天線的“仰角-增益”測試數(shù)據(jù)，利用目標(biāo)點跡高低波束增益差信息，求出目標(biāo)的仰角，進(jìn)而估算目標(biāo)的高度，并采用加權(quán)平均值的方法平滑目標(biāo)的估高值。

1 估高原理

目標(biāo)高度并非一個直接測量值，而是一個導(dǎo)出值，根據(jù)目標(biāo)的斜距和仰角計算可得[5-6]。當(dāng)大氣隨高度的變化梯度為-0.039×10-8m時，ae=(4/3)a=8 490km，a=6 370km，為地球曲率半徑。大氣折射使雷達(dá)電波傳播路徑發(fā)生彎曲，采用等效半徑后，可認(rèn)為電波仍沿直線傳播，則：

(1)

式中：ha為雷達(dá)天線安裝高度；R為目標(biāo)斜距；β為目標(biāo)仰角。

由式(1)可知， 只要測出目標(biāo)仰角β與斜距R，就可以導(dǎo)出目標(biāo)高度。

傳統(tǒng)的兩坐標(biāo)雷達(dá)收發(fā)共用波束，沒有其他多余的信息可用，無法估算目標(biāo)仰角。本文依托的某XX波段航空管制一次雷達(dá)項目，天線發(fā)射采用1個波束，接收采用高低2個波束的方式。根據(jù)天線每種波束的“仰角-波束增益”測試數(shù)據(jù)，可以得出“仰角-波束增益差”的數(shù)據(jù)，如圖1所示。

圖1 仰角-波束增益差圖

由圖1可知，在0°～6°仰角范圍內(nèi),波束增益差呈現(xiàn)單調(diào)減的趨勢，在仰角大于6°時，波束增益差則呈現(xiàn)曲線震蕩的形式，所以為了避免仰角-波束增益差的二值性，這里只選取0°～6°仰角作為估高區(qū)，并將這部分的數(shù)據(jù)作為配置文件存儲。

2 一種基于波束增益差的估高方法

在雷達(dá)探測范圍內(nèi)，探測到的民航飛機(jī)目標(biāo)一般的最大高度為12 000 m左右(由飛機(jī)自帶的GPS信息或二次雷達(dá)信息可得)，在0～6°仰角范圍內(nèi)，對應(yīng)的估高距離區(qū)間：最小距離為12 000/sin6°=114 801m，最大距離為對應(yīng)的雷達(dá)最大探測距離。在雷達(dá)探測目標(biāo)的過程中，由信號處理模塊過濾和凝聚點跡，并送給數(shù)據(jù)處理模塊處理，點跡所帶的屬性參數(shù)中有高波束增益和低波束增益的信息。數(shù)據(jù)處理模塊經(jīng)過點跡預(yù)處理、點航相關(guān)、自動起始、濾波跟蹤等過程實現(xiàn)對目標(biāo)的發(fā)現(xiàn)、起始和跟蹤，在形成的目標(biāo)航跡參數(shù)中有每幀的點跡屬性信息，根據(jù)點跡屬性可提取出每幀的波束增益差，再通過查“仰角-增益差”對照表可以獲得對應(yīng)的仰角，進(jìn)而根據(jù)公式(1)估算出當(dāng)前幀的目標(biāo)高度。

在算法實現(xiàn)中，為了避免奇異值的出現(xiàn)，需要剔除估算的奇異值，奇異值的剔除一方面需要考慮目標(biāo)可能出現(xiàn)的最大高度，另一方面需要結(jié)合估高的歷史值，以2幀平穩(wěn)為原則：假定目標(biāo)單幀最大的高度變化范圍為Δh，若第3幀目標(biāo)估高相比第1幀超過2Δh，則第3幀高度將作為奇異值剔除，并用合理值代替[10]。估高過程中還采取加權(quán)平均法對估高值進(jìn)行平滑，具體的算法流程如下：

STEP1：讀取“仰角-波束增益差”配置文件，得到目標(biāo)對應(yīng)的“波束增益差”和“仰角”對照表，其中仰角以0.1°步進(jìn)初始化目標(biāo)估高歷史值數(shù)組，保留5幀歷史值。

STEP2：待目標(biāo)穩(wěn)定跟蹤后，即目標(biāo)跟蹤歷史幀超過10幀，且目標(biāo)點落入估高的距離范圍時，可對目標(biāo)估高值進(jìn)行更新，并根據(jù)目標(biāo)當(dāng)前幀點跡屬性中的高波束增益和低波束增益求出當(dāng)前幀的波束增益差。

STEP3：由波束增益差查表求出對應(yīng)的仰角，分為以下幾種情況：

(1) 若實際增益差大于對照表中0°仰角對應(yīng)的增益，則仰角就取0°；

(2) 若實際增益差小于對照表中6°仰角對應(yīng)的增益，則仰角就取6°；

(3) 若實際增益差在對照表的增益差范圍內(nèi)，因為對照表的單調(diào)性，可采取“二分查找”法找尋仰角，若表中恰好有該增益值，則選取對應(yīng)的仰角；若沒有，則選取與該增益值距離最近的2個仰角，采用“線性插值”法計算獲取對應(yīng)的仰角。

STEP4：根據(jù)公式(1)，求出目標(biāo)當(dāng)前幀該仰角對應(yīng)的估高值：

(1) 若估高值大于12 000m，則認(rèn)為是奇異值，將其剔除，并用12 000m作為當(dāng)前幀的估高值；

(2) 若當(dāng)前幀的估高值與上2幀的估高值的差大于2×Δh，則認(rèn)為是奇異值，將其剔除，并以前兩幀估高值為基點的當(dāng)前幀線性預(yù)測值作為當(dāng)前幀的估高值。

STEP5：對STEP4中計算的當(dāng)前估高值，進(jìn)行平滑處理：

(1) 若估高值落入估高范圍內(nèi)的目標(biāo)首點，則直接采用該值作為當(dāng)前估高值；

(2) 非首點情況，需根據(jù)估高的歷史值，結(jié)合當(dāng)前值進(jìn)行加權(quán)平滑處理，本算法保留最大歷史幀為5。每幀平滑系數(shù)的選取如表 1所示。

表1 每幀平滑系數(shù)

最后以平滑的估高值作為當(dāng)前幀的估高值。

STEP6：若在估高距離范圍內(nèi)，目標(biāo)當(dāng)前幀沒有跟蹤點，即目標(biāo)點跟蹤丟失，則同樣以前兩幀估高值為基點預(yù)測當(dāng)前幀的估高值；若目標(biāo)連續(xù)丟失三點，則將估高值清零。

3 實驗驗證

該估高方法目前已應(yīng)用于最新研制的某XX波段航空管制一次雷達(dá)數(shù)據(jù)處理系統(tǒng)模塊。該航空管制雷達(dá)同時搭配了本地同軸的二次雷達(dá)，數(shù)據(jù)處理模塊在進(jìn)行點航相關(guān)時，與同一批目標(biāo)航跡相關(guān)的一次雷達(dá)點跡、二次雷達(dá)點跡會進(jìn)行融合處理，融合后的點跡帶有一次點跡和二次點跡的屬性信息，之后對融合后的點跡再進(jìn)行濾波跟蹤。本節(jié)利用目標(biāo)的二次雷達(dá)高度作為目標(biāo)估高算法性能分析的對比對象，二次雷達(dá)高度是目標(biāo)飛機(jī)上應(yīng)答器與二次雷達(dá)之間的應(yīng)答信號，一般與目標(biāo)的真實飛行高度比較接近，可作為目標(biāo)真實高度上報航空管制中心。

這里選取觀測范圍包含估高區(qū)間的目標(biāo)作為實驗對象，在估高距離區(qū)間[150km，350km]內(nèi)，目標(biāo)共觀測到76幀，觀察目標(biāo)估高值與目標(biāo)二次雷達(dá)高度的差值隨距離變化的情況如圖 2所示。目標(biāo)在距離范圍[240km,320km]內(nèi)，估高值與高度真值最為接近，且在[260km,280km]距離內(nèi)，目標(biāo)真實飛行高度發(fā)生變化時，估高值也隨之變化。在整個估高區(qū)間內(nèi)，高度差絕對值的均值為524m，估高值的相對均方根誤差為591m。之后再選取8批估高距離區(qū)間內(nèi)的目標(biāo)進(jìn)行觀察，估高值與二次雷達(dá)高度差值的情況如表 2所示。由表 2可知，這些目標(biāo)的估高值與高度真值差值的絕對均值大約在580m左右，估高值的相對均方根誤差大約在640m左右，具有比較高的估高準(zhǔn)確度。多批目標(biāo)估高的觀察結(jié)果也說明了該算法的有效性和普適性。

圖2 目標(biāo)估高結(jié)果

表2 目標(biāo)估高結(jié)果分析

4 結(jié)束語

本文提出了一種基于波束增益差的兩坐標(biāo)雷達(dá)估高方法，通過目標(biāo)點跡高低波束增益差，求出目標(biāo)的仰角，估算目標(biāo)高度，剔除奇異值，并采用加權(quán)平均值的方法對其進(jìn)行平滑。該方法簡單實用，已在新型的兩坐標(biāo)一次雷達(dá)數(shù)據(jù)處理中得到了應(yīng)用。最后，通過與本地同軸的二次雷達(dá)的目標(biāo)高度數(shù)據(jù)的對比，驗證了該方法的有效性。

[1] 羅志勇,何佳洲.2D組網(wǎng)雷達(dá)目標(biāo)真實高度的最小二乘估計[J].指揮控制與仿真，2007，29(3)：41-44.

[2] 孫景明,劉兆磊,王國宏.2D組網(wǎng)雷達(dá)目標(biāo)高度的最大似然估計[J].現(xiàn)代雷達(dá)，2005,27(9):7-9.

[3] 熊偉,潘旭東,彭應(yīng)寧,等.分布式2D 雷達(dá)網(wǎng)的高度估計技術(shù)[J].信息與控制，2010，39(4):408-412.

[4] 龐志兵．測高雷達(dá)與二維雷達(dá)數(shù)據(jù)融合問題研究[J]．指揮控制與仿真，2006,28(4):10-12.

[5] 劉玉杰，刁丹丹，張濤,等.一種改進(jìn)的三坐標(biāo)雷達(dá)測高方法[J].現(xiàn)代電子技術(shù)，2014，37(5):16-17.

[6] 邵余峰.高機(jī)動三坐標(biāo)雷達(dá)比幅測高誤差分析及修正[J].現(xiàn)代雷達(dá)，2007，29(4):79-81.

[7] 徐綿起,王斌,徐瀚智.三坐標(biāo)地面情報雷達(dá)測高誤差修正方法[J].雷達(dá)科學(xué)與技術(shù),2013，11(1):21-25.

[8] 矯志寧，陶晶.兩坐標(biāo)雷達(dá)測高方法研究[J].現(xiàn)代雷達(dá)，2011，33(12):18-21.

[9] 嚴(yán)俊坤，夏雙志，戴奉周，等.基于單部兩坐標(biāo)雷達(dá)的目標(biāo)三維定位算法[J].電波科學(xué)學(xué)報,2013,28(3):405-411.

[10]林曉斌，張承志，王寶,等.雷達(dá)數(shù)據(jù)處理軟件中的高度平滑算法[J].計算機(jī)光盤軟件與應(yīng)用，2013(1):228-229.

A Height Estimation Approach of 2D Radar Based on Beam Gain Difference

CHEN Lei,FANG Qing

(No.38 Research Institute of CETC,Hefei 230031,China)

2D radars can only provide the distance and azimuth information of target,not provide the height or elevation information of target.For estimating the height of target,this paper presents a height estimation method for 2D radar based on beam gain difference,according to the "elevation gain" test data of radar antenna,calculates the target elevation by using high and low beam gain difference of target point properties and estimates the height information of targets,at the same time uses the weighted average method to smooth the height estimation value of target.The algorithm is simple and applicable,and easy to be realized in engineering.The experiment validates the validity of the method through comparing the height estimation with the second survellance radar (SSR) height information of target.

2D radar;beam gain difference;height estimation approach;height smoothing

2015-04-28

國家科技支撐計劃課題，機(jī)場場面監(jiān)視雷達(dá)系統(tǒng)，項目編號：2011BAH24B05

TN957.51

A

CN32-1413(2015)04-0033-03

10.16426/j.cnki.jcdzdk.2015.04.009