李 敬，張 琛

(中國(guó)船舶重工集團(tuán)公司第七二三研究所， 江蘇 揚(yáng)州 225101)

0 引 言

相控陣?yán)走_(dá)是雷達(dá)信號(hào)處理技術(shù)、新型器件與計(jì)算機(jī)技術(shù)結(jié)合后發(fā)展到高階段的產(chǎn)物，是隨著電子計(jì)算機(jī)和微波移相技術(shù)的發(fā)展而誕生的。雷達(dá)探測(cè)的重要任務(wù)之一就是探測(cè)目標(biāo)位置，目標(biāo)測(cè)向是目標(biāo)位置參數(shù)估計(jì)中的重要研究?jī)?nèi)容，因相控陣天線具有波束掃描的快速性、靈活性與波束形狀的捷變能力，其在波束形成上是多種多樣的，測(cè)向方法也不盡相同。相控陣?yán)走_(dá)一般采用同時(shí)多波束進(jìn)行目標(biāo)探測(cè)，為了測(cè)量目標(biāo)的精確方位，必須采用相應(yīng)的方法利用這些多波束進(jìn)行方位估計(jì)。本文針對(duì)三波束同時(shí)接收的工作方式介紹一種高精度測(cè)角方法，即利用拋物線模型擬合實(shí)現(xiàn)三波束比幅測(cè)角。

1 天線方向圖

比幅法測(cè)角是利用天線收到的回波信號(hào)幅度值來(lái)做角度測(cè)量的，該幅度值的變化規(guī)律取決于天線方向圖以及天線收發(fā)方式。本文所介紹的工作方式是：發(fā)射一個(gè)寬波束，采用3個(gè)相同且彼此部分重疊的波束同時(shí)接收的方式[1]，接收示意圖見(jiàn)圖1。

圖1 三波束接收示意圖

圖1中左、中、右3個(gè)波束同時(shí)接收目標(biāo)回波，Bw為半功率波束寬度，三波束彼此重疊的角度為Bw/2。圖2是1路發(fā)射波束3路接收波束的方向圖。圖中黑實(shí)線表示指向?yàn)榉ň€方向(0°)的發(fā)射波束，虛線依次從左到右分別表示3 dB波束寬度為Bw的接收波束的左(波束1)、中(波束2)、右(波束3)三波束。

圖2 發(fā)射波束和3路接收波束方向圖

2 三波束比幅測(cè)角

三波束同時(shí)接收的目標(biāo)回波，在幅度上有一定的對(duì)應(yīng)關(guān)系，即目標(biāo)靠近哪個(gè)接收波束，哪個(gè)波束的回波幅度就大，且回波幅度最大的位置就是目標(biāo)所在位置，為此可以利用開(kāi)口向下的拋物線函數(shù)y=ax2+c進(jìn)行三波束目標(biāo)回波幅度的函數(shù)擬合，求出目標(biāo)所在的具體方位。如圖3所示，該拋物線以y軸為對(duì)稱軸，x軸表示中心波束指向偏離目標(biāo)的角度偏差，y軸表示回波幅度。

設(shè)左、中、右三波束測(cè)得的目標(biāo)參數(shù)為(xlyl)、(xcyc)和(xryr)。xl，xc和xr分別為左、中、右三波束的方位指向，yl，yc和yr分別為三波束測(cè)得的目標(biāo)回波幅度。如果目標(biāo)位于中心波束，不考慮副瓣影響，理論上左右兩波束幅度值相等，分別位于拋物線的兩邊。此時(shí)xc=0表示中心波束指向目標(biāo)，目標(biāo)幅度最大值c=yc落在y軸上。當(dāng)目標(biāo)偏向左波束，位于中心波束和左波束之間時(shí)，回波幅度的最大值已不是中心波束上測(cè)得的回波幅度值了，則中心波束指向xc≠0。如圖3所示，xc代表了目標(biāo)偏離中心波束的角度偏差。此時(shí)目標(biāo)指向的回波幅度仍然是最大，仍然位于拋物線最高點(diǎn)，最大值仍然是c，即目標(biāo)最大幅度永遠(yuǎn)位于y軸上。由此設(shè)目標(biāo)在拋物線上的參數(shù)為(xpyp)，xp>0,yc

圖3 拋物線擬合實(shí)現(xiàn)三波束比幅測(cè)角

將左、中、右三波束測(cè)得的目標(biāo)參數(shù)(xlyl)、(xcyc)和(xryr)代入拋物線模型得：

(1)

若三波束之間兩兩夾角為波束寬度的一半，即Bw/2，則有：

(2)

將公式(2)代入公式(1)求得：

(3)

從公式(3)可知yl，yc和yr分別為三波束測(cè)得的目標(biāo)回波幅度，為已知值，Bw為接收波束寬度，也是已知的，由此可求得擬合參數(shù)a和c，該拋物線就確定了。同時(shí)還可求得目標(biāo)角度修正值xc。相控陣?yán)走_(dá)在工作的時(shí)候，波束指向是確定的，將波束指向用xc值修正后就可得出目標(biāo)的真方位和回波幅度。

3 數(shù)據(jù)驗(yàn)證與結(jié)論

拋物線模型擬合實(shí)現(xiàn)三波束比幅測(cè)角的驗(yàn)證是通過(guò)在設(shè)備上測(cè)量民航機(jī)來(lái)驗(yàn)證的，見(jiàn)圖4。圖4中的點(diǎn)是民航機(jī)在不同波位上的視頻回波，為了方便顯示，在每個(gè)波位上只繪出中心波束的回波視頻，框內(nèi)是民航機(jī)的回波視頻。通過(guò)采用拋物線模型擬合實(shí)現(xiàn)求每組波束目標(biāo)的位置后，再進(jìn)行點(diǎn)跡凝聚和航跡處理，實(shí)現(xiàn)對(duì)目標(biāo)的穩(wěn)定跟蹤，與ADS-B數(shù)據(jù)對(duì)比，試驗(yàn)測(cè)得方位跟蹤精度均方根小于0.3°。

圖4 民航機(jī)回波視頻

圖5是對(duì)回波數(shù)據(jù)進(jìn)行比幅測(cè)角，求目標(biāo)方位和俯仰、視頻凝聚后的點(diǎn)跡數(shù)據(jù)，框內(nèi)是民航機(jī)的點(diǎn)跡數(shù)據(jù)。從圖5可以看出，民航機(jī)在整個(gè)飛行過(guò)程是均勻飛行的。表1是對(duì)應(yīng)圖4框內(nèi)視頻回波的數(shù)值，圖4中顯示的是地理坐標(biāo)系下目標(biāo)的位置，表1中對(duì)應(yīng)的是陣面坐標(biāo)系下目標(biāo)的位置，陣面架設(shè)參數(shù)方位-45°，傾斜角20°，比幅測(cè)角數(shù)據(jù)最終要從陣面坐標(biāo)轉(zhuǎn)換到地理坐標(biāo)系下[2]。

圖5 民航機(jī)點(diǎn)跡與航跡

由于數(shù)據(jù)量比較大，只給出了3個(gè)掃描周期的視頻數(shù)據(jù)及對(duì)應(yīng)的參數(shù)，每個(gè)粗框內(nèi)對(duì)應(yīng)著一個(gè)掃描周期的視頻回波，每個(gè)周期大約有3～4個(gè)波位的立體角內(nèi)可以探測(cè)到目標(biāo)。每個(gè)粗框內(nèi)的數(shù)據(jù)可以凝聚成圖5中的一個(gè)點(diǎn)。每個(gè)波位對(duì)應(yīng)大約2～3個(gè)距離單元，每個(gè)距離單元對(duì)應(yīng)目標(biāo)左、中、右(如圖1所示的)3個(gè)回波波束的幅度值，目標(biāo)幅度和方位分別是對(duì)左、中、右三波束進(jìn)行拋物線擬合實(shí)現(xiàn)三波束比幅測(cè)角得到的目標(biāo)方位值和該方位對(duì)應(yīng)的目標(biāo)幅度值。

表1 民航機(jī)三波束視頻回波

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

文獻(xiàn)[3]中的映射內(nèi)插法雖然也采用拋物線模型，但是用了y=ax2+bx+c模型，而工程中已知的只有3個(gè)波束的幅值和波束間隔，在解算目標(biāo)方位時(shí)，少一個(gè)已知量，因此添加了一個(gè)約束條件α=-b/2a，這就導(dǎo)致函數(shù)的不靈活性，而本文巧妙地將目標(biāo)位置移至y軸上，這樣就變成了y=ax2+c模型，既滿足了已知參數(shù)，又不加任何約束，可以根據(jù)已知參數(shù)，擬合最佳拋物線系數(shù)，對(duì)于多波束(大于3個(gè)接收波束)也可采用此方法。

利用拋物線擬合實(shí)現(xiàn)三波束比幅測(cè)向不需要進(jìn)行繁瑣的定標(biāo)，處理簡(jiǎn)單，且由相控陣?yán)走_(dá)原理知，隨著陣列掃描方位的變化，發(fā)射波束和接收波束會(huì)相應(yīng)展寬，即Bw隨著天線掃描角度變化而展寬，因此隨著角度變化拋物線形狀是變化的，該方法可以充分利用3對(duì)測(cè)量數(shù)據(jù)進(jìn)行拋物線擬合求得a和c，從而求出目標(biāo)指向。同時(shí)該算法還可以求得目標(biāo)的幅度參數(shù)，如果設(shè)備的穩(wěn)定性好，幅度參數(shù)作為后續(xù)數(shù)據(jù)處理會(huì)是一個(gè)很重要的參數(shù)，可以為后續(xù)雜波圖積累和航跡關(guān)聯(lián)提供依據(jù)。

必須注意的是比幅測(cè)角要求同一時(shí)刻測(cè)得的目標(biāo)才能進(jìn)行比幅，如果對(duì)某個(gè)運(yùn)動(dòng)目標(biāo)掃描時(shí)不連續(xù)，是無(wú)法進(jìn)行比幅測(cè)角的，目標(biāo)可能跨多個(gè)距離單元。