吳其華，劉進(jìn)，艾小鋒，王雪松，肖順平

(國防科技大學(xué)a.電子信息系統(tǒng)復(fù)雜電磁環(huán)境效應(yīng)國家重點(diǎn)實(shí)驗(yàn)室；b.理學(xué)院，湖南 長沙　410073)

?

探測跟蹤技術(shù)

反導(dǎo)任務(wù)下一種相控陣?yán)走_(dá)搜索參數(shù)優(yōu)化方法*

吳其華a，劉進(jìn)a，艾小鋒a，王雪松，肖順平a

(國防科技大學(xué)a.電子信息系統(tǒng)復(fù)雜電磁環(huán)境效應(yīng)國家重點(diǎn)實(shí)驗(yàn)室；b.理學(xué)院，湖南 長沙410073)

摘要：搜索參數(shù)優(yōu)化設(shè)計是相控陣?yán)走_(dá)探測性能提升的重要途徑。在相控陣?yán)走_(dá)反導(dǎo)任務(wù)背景下，針對導(dǎo)彈類目標(biāo)建立了以最大期望發(fā)現(xiàn)距離為目標(biāo)函數(shù)的搜索參數(shù)優(yōu)化模型，在目標(biāo)勻速運(yùn)動和加速運(yùn)動的條件下分別給出了雷達(dá)搜索波束最佳駐留時間以及相應(yīng)目標(biāo)期望發(fā)現(xiàn)距離的優(yōu)化表達(dá)式，仿真證明了算法的有效性。

關(guān)鍵詞：搜索參數(shù)優(yōu)化；期望發(fā)現(xiàn)距離；反導(dǎo)；相控陣?yán)走_(dá)；最佳駐留時間；雷達(dá)探測

0引言

相控陣?yán)走_(dá)由于其實(shí)時多任務(wù)、靈活的資源調(diào)度方式等優(yōu)點(diǎn)，在戰(zhàn)略防御、精確打擊等現(xiàn)代戰(zhàn)爭的核心環(huán)節(jié)中發(fā)揮著至關(guān)重要的作用。從相控陣?yán)走_(dá)內(nèi)部任務(wù)角度，搜索任務(wù)在其中占有重要位置，它是通過對興趣區(qū)域建立波位表，根據(jù)設(shè)定的相應(yīng)參數(shù)進(jìn)行搜索操作，這項(xiàng)任務(wù)也是后續(xù)建立目標(biāo)跟蹤的前提與基礎(chǔ)。因此，搜索性能優(yōu)化對于提升相控陣?yán)走_(dá)性能有著重要的意義。

相控陣?yán)走_(dá)搜索性能優(yōu)化一個重要途徑是對搜索參數(shù)的優(yōu)化設(shè)計。文獻(xiàn)[1-4]研究了單個空域下搜索參數(shù)優(yōu)化問題。其中，文獻(xiàn)[1]針對美海軍UESA雷達(dá)通過調(diào)整搜索幀周期與雷達(dá)靈敏度的方法，實(shí)現(xiàn)對于雷達(dá)資源的有效利用。文獻(xiàn)[2]研究了搜索資源受限情況下搜索參數(shù)設(shè)計，通過調(diào)整雷達(dá)的探測距離與搜索幀周期在滿足多幀積累檢測概率要求條件下使得目標(biāo)初始跟蹤距離最大。文獻(xiàn)[3]建立了遠(yuǎn)程相控陣?yán)走_(dá)搜索參數(shù)設(shè)計模型，通過調(diào)整搜索幀周期和單次檢測概率保證目標(biāo)最小捕獲時間。文獻(xiàn)[4]則加入虛警概率因素，聯(lián)合搜索波束駐留時間調(diào)整實(shí)現(xiàn)優(yōu)化。文獻(xiàn)[5-11]則從單空域向多區(qū)域優(yōu)化搜索發(fā)展，關(guān)注了多空域情況下最優(yōu)搜索資源分配的問題，相應(yīng)地提出了最佳搜索負(fù)載[5]、最優(yōu)搜索數(shù)據(jù)率[6]等模型。然而，現(xiàn)有研究假定目標(biāo)服從SwerlingⅠ型分布的前提，此模型符合由若干個相似的散射體組成的目標(biāo)，沒有考慮到實(shí)際戰(zhàn)情中復(fù)雜多變的目標(biāo)特性，當(dāng)目標(biāo)起伏模型失配時所得到的結(jié)果則不是最優(yōu)的，因此具有一定的局限性。本文從相控陣?yán)走_(dá)擔(dān)負(fù)的反導(dǎo)任務(wù)出發(fā)，針對所關(guān)心的導(dǎo)彈類目標(biāo)以期望發(fā)現(xiàn)距離為目標(biāo)函數(shù)建立搜索參數(shù)優(yōu)化模型，給出相應(yīng)的搜索波束最佳駐留時間與最大目標(biāo)期望發(fā)現(xiàn)距離，并通過仿真分析其優(yōu)化性能。

1相控陣?yán)走_(dá)反導(dǎo)任務(wù)相關(guān)模型

1.1導(dǎo)彈目標(biāo)起伏統(tǒng)計模型及檢測概率

現(xiàn)代戰(zhàn)爭中，來自于空中的軍事打擊已成為國防安全的主要威脅，其中導(dǎo)彈類目標(biāo)機(jī)動靈活、破壞性大，是戰(zhàn)略防御的主要對象，反導(dǎo)也成為先進(jìn)相控陣?yán)走_(dá)的重要作戰(zhàn)任務(wù)。

根據(jù)目標(biāo)檢測理論，檢測概率與起伏類型有關(guān)，導(dǎo)彈類目標(biāo)主要由一個大RCS的散射體與很多小的、RCS相等的散射體組成，服從Swerling Ⅲ模型[3]。Swerling Ⅲ模型滿足在任何一次掃描時從目標(biāo)收到的回波脈沖幅度在這一次掃描的整個過程中恒定不變，且從掃描到掃描不相關(guān)[12]。其目標(biāo)檢測概率計算公式[13]如下：

(1)

(2)

式中：np為脈沖積累個數(shù)；VT為檢測門限；SNR為信噪比。

不使用脈沖積累(np=1)時，假設(shè)雷達(dá)虛警概率為Pfa，有Pfa=e-VT，得檢測概率：

(3)

通常情況SNR?1，式(3)簡化為

(4)

式(4)即為反導(dǎo)任務(wù)中導(dǎo)彈類目標(biāo)檢測概率與信噪比關(guān)系的最終表達(dá)式。

1.2目標(biāo)期望發(fā)現(xiàn)距離

相控陣?yán)走_(dá)搜索性能優(yōu)化準(zhǔn)則有多種表達(dá)形式，如最大目標(biāo)跟蹤起始距離、最小發(fā)現(xiàn)時間等。本文則從統(tǒng)計意義上關(guān)注搜索優(yōu)化問題，將搜索性能表述為：搜索過程中雷達(dá)發(fā)現(xiàn)逼近目標(biāo)的期望距離最大，定義為目標(biāo)期望發(fā)現(xiàn)距離，指突防目標(biāo)進(jìn)入雷達(dá)探測威力范圍到被最終探測到時距離雷達(dá)的平均距離，如圖 1所示。

圖1　目標(biāo)期望發(fā)現(xiàn)距離示意圖Fig.1　Sketch map of expected discovering range

假設(shè)雷達(dá)搜索幀周期為tf，記雷達(dá)發(fā)現(xiàn)目標(biāo)的時間為t，對于采用均勻波位掃描策略的雷達(dá)來說，易得目標(biāo)進(jìn)入雷達(dá)的探測范圍后到被雷達(dá)發(fā)現(xiàn)的時間服從均勻分布，有t∈U(0,tf)，則期望發(fā)現(xiàn)時間tE為

(5)

目標(biāo)最終被發(fā)現(xiàn)時的平均運(yùn)動距離是關(guān)于tE的函數(shù)，記為R(tE)，則有R(tE)=R(0.5tf)。假設(shè)雷達(dá)的探測距離為RD(對應(yīng)檢測概率為Pd)，目標(biāo)期望發(fā)現(xiàn)距離可以表述為

RDE=RD-R(tE)=RD-R(0.5tf).

(6)

搜索參數(shù)優(yōu)化目標(biāo)即是使式(6)中所表示的目標(biāo)期望發(fā)現(xiàn)距離最大。

2基于期望發(fā)現(xiàn)距離的搜索參數(shù)優(yōu)化

分析式(6)，目標(biāo)期望發(fā)現(xiàn)距離與雷達(dá)的探測距離以及目標(biāo)在進(jìn)入雷達(dá)威力范圍后二分之一個搜索幀周期內(nèi)的運(yùn)動距離有關(guān)。

首先分析雷達(dá)探測距離的影響。根據(jù)雷達(dá)方程[14]：

(7)

式中：Pav為平均發(fā)射功率(W)；τs為搜索波束駐留時間(s)；Gt為發(fā)射天線增益；Gr為接收天線增益；λ為雷達(dá)波長(m)；σ為目標(biāo)RCS(m2)；k為玻爾茲曼常數(shù)；T0為接收機(jī)噪聲溫度，常溫下為290 K；Fn為接收機(jī)噪聲系數(shù)；L為雷達(dá)系統(tǒng)損耗；SNR為回波信噪比。

由式(7)得到：

(8)

(9)

對于導(dǎo)彈類目標(biāo)，將式(9)代入式(4)得到檢測概率與雷達(dá)探測距離的關(guān)系為

(10)

由式(10)得，雷達(dá)探測距離是關(guān)于搜索波束駐留時間的函數(shù)。將式(9)變形得

(11)

其次分析目標(biāo)在進(jìn)入雷達(dá)威力范圍后二分之一個搜索幀周期內(nèi)的運(yùn)動距離。經(jīng)過分析，R(tE)除了與tf有關(guān)，還與目標(biāo)的運(yùn)動模型有關(guān)。下面就目標(biāo)勻速運(yùn)動與加速運(yùn)動情況分別進(jìn)行研究。

2.1勻速運(yùn)動模型下搜索參數(shù)優(yōu)化

首先分析較為簡單的目標(biāo)勻速運(yùn)動的情況，假設(shè)突防彈以徑向速度v靠近雷達(dá)站，則R(tE)可寫為

R(tE)=R(0.5tf)=0.5vtf.

(12)

假定相控陣?yán)走_(dá)的搜索波束駐留時間為τs，搜索區(qū)域內(nèi)安排波位數(shù)為Ns，搜索任務(wù)所占雷達(dá)資源的比例為SR，則有

(13)

即在搜索波位與搜索資源比例一定時，tf與τs滿足線性關(guān)系。因此，R(tE)是關(guān)于τs的函數(shù)。

將式(11)，(12)，(13)代入式(6)，得到優(yōu)化表達(dá)式：

(14)

式(14)表明目標(biāo)期望發(fā)現(xiàn)距離RDE是關(guān)于搜索駐留時間τs的函數(shù)，基于期望發(fā)現(xiàn)距離的優(yōu)化即是對搜索波束駐留時間τs進(jìn)行優(yōu)化。

假定一組雷達(dá)參數(shù)：虛警概率Pfa=10-6，檢測信噪比SNRD=22.7 dB，系統(tǒng)常數(shù)Ω0=1.93×1022，搜索占雷達(dá)資源比例SR=100%，搜索波位數(shù)Ns=1 000，目標(biāo)速度v=4 km/s，得到RDE與τs的關(guān)系如圖 2所示。

圖2　期望發(fā)現(xiàn)距離與搜索波束駐留時間關(guān)系圖Fig.2　Relation between expected discovering range and beam dwell time

從圖2可以看出，當(dāng)搜索波束駐留時間τs增大時，增加了雷達(dá)的探測距離，目標(biāo)期望發(fā)現(xiàn)距離也變大；但如果搜索波束駐留時間τs過大，則導(dǎo)致目標(biāo)在雷達(dá)一個搜索幀周期內(nèi)運(yùn)動距離變大，最終將抵消掉由于探測距離的增加所帶來的期望發(fā)現(xiàn)距離的改善。因此存在一個最佳的搜索波束駐留時間τs，使得目標(biāo)期望發(fā)現(xiàn)距離RDE最大。

對式(14)求導(dǎo)，得到：

(15)

(16)

將式(16)代入式(14)得到最大期望發(fā)現(xiàn)距離表達(dá)式：

(17)

式(16)即為突防彈勻速運(yùn)動條件下相控陣?yán)走_(dá)搜索波束最佳駐留時間表達(dá)式，式(17)為對應(yīng)最大目標(biāo)期望發(fā)現(xiàn)距離。

2.2勻加速運(yùn)動模型下搜索參數(shù)優(yōu)化

考慮到對于實(shí)際的彈道類導(dǎo)彈目標(biāo)，其運(yùn)動特性符合自由空間的“二體軌道模型”[15]，因此上述討論中的勻速運(yùn)動模型并不能完全準(zhǔn)確描述導(dǎo)彈目標(biāo)的運(yùn)動狀態(tài)?？紤]到雷達(dá)探測過程所花費(fèi)時間較短，導(dǎo)彈運(yùn)動參數(shù)變化較為緩慢，在主動關(guān)機(jī)點(diǎn)以后，導(dǎo)彈只受到地球引力的作用，因此在雷達(dá)探測的這段較短時間內(nèi)導(dǎo)彈的運(yùn)動狀態(tài)可以用勻加速運(yùn)動模型來進(jìn)行描述，并假設(shè)其加速度為一個g。

假設(shè)導(dǎo)彈進(jìn)入雷達(dá)探測距離處的速度為v0，根據(jù)運(yùn)動學(xué)方程，在期望發(fā)現(xiàn)時間內(nèi)目標(biāo)的運(yùn)動距離可以表述為

(18)

同勻速運(yùn)動模型下的推理過程，可以得到對應(yīng)期望發(fā)現(xiàn)距離的表達(dá)式：

(19)

(20)

3仿真分析

上述優(yōu)化方法進(jìn)行仿真分析驗(yàn)證，假定相控陣?yán)走_(dá)只有搜索與跟蹤兩種任務(wù)形式，且跟蹤任務(wù)搶占搜索資源，搜索資源SR=1-TR(TR為跟蹤資源比例)，仿真參數(shù)設(shè)置如下：某脈沖體制S波段相控陣?yán)走_(dá)波長λ=9 cm，發(fā)射功率Pt=150 kW，最高占空比為5%，發(fā)射天線增益Gt=40 dB，接收天線增益Gr=40 dB，目標(biāo)RCSσ=1 m2，接收機(jī)噪聲系數(shù)Fn=3 dB，系統(tǒng)損耗L=5 dB；虛警概率Pfa=10-6，檢測概率Pd=0.99，搜索波位數(shù)Ns=1 000，戰(zhàn)情為導(dǎo)彈突防，進(jìn)入雷達(dá)探測距離處時速度v=4 km/s。分別在勻速運(yùn)動模型與勻加速運(yùn)動模型假設(shè)下將本文方法優(yōu)化結(jié)果與模型失配時優(yōu)化結(jié)果以及固定搜索波束駐留時間模型進(jìn)行比較，其中固定搜索波束駐留時間取τs=10 ms，所得結(jié)果如圖 3所示。

圖3　不同跟蹤資源下期望發(fā)現(xiàn)距離圖Fig.3　Expected discovering range with different tracking resources

可見，無論是在勻速運(yùn)動還是在勻加速運(yùn)動的條件下，隨著搜索資源的減少，3種方法對應(yīng)的目標(biāo)期望發(fā)現(xiàn)距離都呈下降的趨勢，雷達(dá)搜索性能降低；將本文優(yōu)化方法與固定駐留時間模型以及Swerling Ⅰ模型建模的期望發(fā)現(xiàn)距離結(jié)果比較，看出由于對導(dǎo)彈目標(biāo)進(jìn)行了更為準(zhǔn)確的建模，本文優(yōu)化方法在各種搜索資源比例下都獲得了更大的目標(biāo)期望發(fā)現(xiàn)距離，實(shí)現(xiàn)了更優(yōu)的搜索性能。另外由于用Swerling Ⅰ模型建模時目標(biāo)模型的失配，在一定范圍內(nèi)其優(yōu)化性能甚至不如固定搜索波束駐留時間的結(jié)果。另一方面，對圖 3a)， 3b)進(jìn)行比較發(fā)現(xiàn)，目標(biāo)勻速運(yùn)動和勻加速運(yùn)動模型對應(yīng)曲線的走勢非常相似，這是由于在突防彈高速運(yùn)動的情況下，在雷達(dá)探測的較短時間內(nèi)其速度等運(yùn)動參數(shù)相對變化很小，導(dǎo)致采用勻速運(yùn)動模型和勻加速運(yùn)動模型的區(qū)別并不是很明顯。

圖4　不同搜索資源下搜索性能改善因子圖Fig.4　Factor of search performance improvement with different search resources

4結(jié)束語

導(dǎo)彈目標(biāo)是國土安全的重要威脅，本文針對相控陣?yán)走_(dá)反導(dǎo)任務(wù)，對導(dǎo)彈類目標(biāo)以最大期望發(fā)現(xiàn)距離為目標(biāo)函數(shù)建立搜索參數(shù)優(yōu)化的模型，并針對突防彈勻速運(yùn)動與勻加速運(yùn)動分別給出搜索波束最佳駐留時間與相應(yīng)期望發(fā)現(xiàn)距離的表達(dá)式。仿真結(jié)果表明本文方法由于對目標(biāo)更為準(zhǔn)確的建模獲得了相比固定搜索波束駐留時間以及未考慮目標(biāo)特性的Swerling Ⅰ模型建模更優(yōu)的搜索性能。

參考文獻(xiàn)：

[1]ZATMAN M.Radar Resource Management for UESA[C]∥The 2002 IEEE Radar Conference，Massachusettes，IEEE,2002:73-76.

[2]盧建斌,胡衛(wèi)東,郁文賢.相控陣?yán)走_(dá)資源受限時最優(yōu)搜索性能研究[J].系統(tǒng)與控制工程,2004,26(10):1388-1390.

LU Jian-bin, HU Wei-dong, YU Wen-xian. Research on Optimal Search Performance of Phased Array Radars with Limited Resources[J].Systems Engineering and Electronics, 2004,26(10):1388-1390.

[3]鄧貴福，劉華林，胥雷.遠(yuǎn)程相控陣?yán)走_(dá)搜索參數(shù)優(yōu)化設(shè)計[J].雷達(dá)科學(xué)與技術(shù),2012,10 (1):32-36.

DENG Gui-fu,LIU Hua-lin,XU Lei. Optimization of Search Parameters of Long Range Phased Array Radar[J].Radar Science and Technology,2012,10(1):32-36.

[4]張華睿，楊宏文，郁文賢.相控陣?yán)走_(dá)最優(yōu)搜索參數(shù)設(shè)計研究[J].兵工學(xué)報,2012, 33(9):1062-1065.

ZHANG Hua-rui,YANG Hong-wen,YU Wen-xian. Design of Optimal Search Operation Parameters for Phased Array Radar[J].Acta Armamentarii,Acta Armamentarii 2012,33(9):1062-1065.

[5]張永杰，李少洪，朱海冰.多功能相控陣?yán)走_(dá)優(yōu)化搜索建模與仿真研究[J].系統(tǒng)仿真學(xué)報, 2008,20(16):4248-4251.

ZHANG Yong-jie,LI Shao-hong,ZHU Hai-bing. Research on Modeling and Simulation of Optimal Search of Multifunction Phased Array Radar[J].Journal of System Simulation, 2008,20(16):4248-4251.

[6]程小楓，吳少鵬，涂剛毅.基于權(quán)重的分區(qū)最優(yōu)搜索數(shù)據(jù)率研究[J].現(xiàn)代雷達(dá),2013, 35(12):27-30.

CHENG Xiao-feng, WU Shao-peng, TU Gang-yi. A Study on Regional Optimal Search Data Rate Based on Weight[J].Modern Radar, 2013, 35(12):27-30.

[7]徐斌，楊晨陽，李少洪，等.相控陣?yán)走_(dá)中的最優(yōu)分區(qū)搜索算法[J].電子學(xué)報,2000,28(12):69-73.

XU Bin, YANG Chen-yang, LI Shao-hong,et al. Optimal Region Search Strategy in Phased Array Radar[J]. Acta Electronica Sinaca, 2000,28(12):69-73.

[8]徐斌，楊晨陽，李少洪，等.相控陣?yán)走_(dá)中的自適應(yīng)搜索研究[J].電子學(xué)報,2001,29(12):1719-1722.

XU Bin, YANG Chen-yang, LI Shao-hong,et al. A Study on Adaptive Search Strategy in Phased Array Radar[J]. Acta Electronica Sinaca, 2001,29(12):1719-1722.

[9]周穎，王雪松，王國玉，等.相控陣?yán)走_(dá)最優(yōu)搜索隨機(jī)規(guī)劃研究[J].現(xiàn)代雷達(dá),2005,27(4):60-63.

ZHOU Ying, WANG Xue-song, WANG Guo-yu, et al. Probabilistic Programming of Optimal Search Strategy for Phased.Array Radar[J].Mordern Radar, 2005,27(4):60-63.

[10]ZHANG Yong-jie, LI Shao-hong. Study on the Optimal and Adaptive Search of Multifunction Phased Array Radar[C]∥IEEE 8th International Conference on Signal Processing,2006:2755-2758.

[11]Payam Haftbaradaran, Mahmoud Kamarei, Reza Fatemi Mofrad. The Optimal Search for Multifunction Phased Array Radar[C]∥Antennas & Propagation Conference, Loughborough, UK,2009:609-612.

[12]Merrill I Skoinik.雷達(dá)系統(tǒng)導(dǎo)論 [M].3版.左群聲，徐國良，馬林，等，譯.北京：電子工業(yè)出版社,2010:50-52.

Merrill I Skoinik. Introduction to Radar Systems [M].3rd Edition. ZUO Qun-sheng, XU Guo-liang, MA Lin,et al，Translated.Beijing:Publishing House of Electronics Industry,2010:50-52.

[13]Bassem R Mahafza，Atef Z Elsherbeni.雷達(dá)系統(tǒng)設(shè)計MATLAB仿真[M]. 朱國富，黃曉濤，黎向陽，等,譯.北京：電子工業(yè)出版社,2009:70-71.

Bassem R Mahafza，Atef Z Elsherbeni. MATLAB Simulations of Radar Systems Design[M].ZHU Guo-fu,HUANG Xiao-tao, LI Xiang-yang, et al,Translated.Beijing: Publishing House of Electronics Industry,2009:70-71.

[14]丁鷺飛,耿富錄，陳建春.雷達(dá)原理[M].北京：電子工業(yè)出版社，2009：3-4.

DING Lu-fei, GENG Fu-lu, CHEN Jian-chun. Radar Principle[M].Beijing: Publishing House of Electronics Industry,2009：3-4.

[15]王雪松，肖順平，馮德軍，等.現(xiàn)代雷達(dá)電子戰(zhàn)系統(tǒng)建模與仿真[M].北京：電子工業(yè)出版社，2010:27-32.

WANG Xue-song, XIAO Shun-ping, FENG De-jun, et al.Modeling and Simulation of Modern Radar and Electronic Warfare Systems[M].Beijing: Publishing House of Electronics Industry,2010:27-32.

Method of Search Parameters Optimization of Phased Array Radar for Antimissile Mission

WU Qi-huaa, LIU Jina, AI Xiao-fenga, WANG Xue-song, XIAO Shun-pinga

(National University of Defense Technology,a.State Key Laboratory of Complex Electromagnetic Environment Effects on Electronics and Information System; b.School of Science,Hunan Changsha 410073, China)

Abstract:Optimization of search parameters is an important approach for phased array radar search performance improvement. For antimissile mission, an optimization model of search parameters is built based on analysis of missile features, then the expression of radar’s best search beam dwell time and corresponding expected discovering range are provided under in the condition of both uniform motion and uniformly accelerated motion. The simulation proves the validity of this method.

Key words:optimization of search parameters; expected discovering range; antimissile; phased array radar; optimal dwell time; radar detection

*收稿日期：2015-03-29；修回日期：2015-07-07

基金項(xiàng)目：國家自然科學(xué)基金(61101180,61401491)

作者簡介：吳其華(1990-)，男，江蘇鹽城人。碩士生，主要研究方向?yàn)榭臻g信息獲取與處理、相控陣?yán)走_(dá)技術(shù)。 E-mail：wuqihualan@163.com

通信地址：410073湖南省長沙市開福區(qū)德雅路109號國防科技大學(xué)電子科學(xué)與工程學(xué)院CEMEE實(shí)驗(yàn)室

doi:10.3969/j.issn.1009-086x.2016.02.026

中圖分類號：TN958.92

文獻(xiàn)標(biāo)志碼：A

文章編號：1009-086X(2016)-02-0165-06