楊華東, 燕永敦, 楊 侃

(1.海軍工程大學(xué)兵器學(xué)院, 湖北 武漢 430033; 2. 海軍研究院, 北京 100161)

0 引 言

“飽和攻擊”是反艦導(dǎo)彈作戰(zhàn)的一般樣式[1-2]，即同時從多方向、多層次向目標(biāo)發(fā)射多枚導(dǎo)彈,使敵反導(dǎo)能力處于飽和狀態(tài),以提高反艦導(dǎo)彈的突防概率[3-4]?！巴瑫r到達(dá)”是“飽和攻擊”戰(zhàn)術(shù)實施的必要條件,如何才能實現(xiàn)“同時到達(dá)”,主要有以下幾種技術(shù)途徑:

一是控制導(dǎo)彈的發(fā)射順序和發(fā)射間隔時間。不同于無人飛行器、無人潛航器等低速、長航時航行器,可以通過調(diào)節(jié)飛行速度來控制飛行時間,以達(dá)到“同時到達(dá)”的目的[5-11]。反艦導(dǎo)彈的飛行速度一般不進(jìn)行精確控制,全程的飛行時間通常以導(dǎo)彈平均速度計算。因此,可以通過控制導(dǎo)彈的發(fā)射順序和發(fā)射間隔時間,使得不同航程的導(dǎo)彈近似同時到達(dá)目標(biāo)。文獻(xiàn)[12]基于這種思想研究了不同攻擊角情況下,航路的飛行時間差。這種方法優(yōu)點在于實現(xiàn)簡單,但其缺陷在于:發(fā)射平臺從進(jìn)入攻擊陣位發(fā)射第一枚導(dǎo)彈開始,需要一直保持陣位,直到同一攻擊波次的導(dǎo)彈全部發(fā)射完畢才能撤離。

二是通過航路規(guī)劃控制導(dǎo)彈的飛行航程。即在認(rèn)為導(dǎo)彈平均速度為定值的前提下,通過航路規(guī)劃技術(shù)調(diào)節(jié)各枚導(dǎo)彈的航程,使得導(dǎo)彈的飛行航路既能滿足導(dǎo)彈的飛行約束和突防要求,又能保證同一攻擊波次所有導(dǎo)彈的飛行航程相等,從而使得同一攻擊波次的各枚導(dǎo)彈“同時到達(dá)”攻擊目標(biāo)。文獻(xiàn)[13]利用“橢圓上任意一點到兩個焦點距離之和相等”的幾何特性,將攻擊艦和目標(biāo)艦視為橢圓的兩個焦點,提出了只有一個航路點的情況下航程相等的規(guī)劃方法。文獻(xiàn)[14]基于類似思想,研究了航路點解空間的分布。但這種方法受導(dǎo)彈飛行限制條件和戰(zhàn)場環(huán)境約束,應(yīng)用場合十分受限。本文針對上述兩種技術(shù)途徑當(dāng)前存在缺陷,進(jìn)行了反艦導(dǎo)彈同時臨空控制策略探索研究。

1 海洋環(huán)境因素對導(dǎo)彈飛行時間修正

導(dǎo)彈飛行速度受大氣溫度、風(fēng)速風(fēng)向的影響存在一定的波動,對導(dǎo)彈飛行時間預(yù)測誤差較大,會影響打擊目標(biāo)時臨空時間的協(xié)同[15]。本文通過獲取導(dǎo)彈飛行航路上各航路點的溫度、風(fēng)速、風(fēng)向等海洋環(huán)境參數(shù),進(jìn)行導(dǎo)彈預(yù)估速度修正,分段計算各航路的準(zhǔn)確飛行時間,從而獲得總的飛行時間,用于導(dǎo)彈臨空時間的協(xié)同計算,通過發(fā)射時序控制,精確實現(xiàn)同一波次攻擊的導(dǎo)彈同時臨空。

首先,通過引入海洋環(huán)境參數(shù),查閱各個航路點所在地理位置點對應(yīng)的環(huán)境溫度,按照溫度修正式(1),計算同一條航路中各航路點導(dǎo)彈修正速度Vi(i=1,2,…,n),n為該條航路的航路點總數(shù)。

(1)

式中,K0為溫度修正常數(shù);V(δZf)陣風(fēng)影響速度修正常數(shù);Ki為航路第i個航路點所在地理位置點的環(huán)境溫度。

然后,依次計算發(fā)射點至第i個(i=1,2,…,n)航路點各段航路的飛行時間為

(2)

最后,按式(3)計算每條航路導(dǎo)彈的飛行時間。

(3)

2 基于發(fā)射間隔控制實現(xiàn)同時到達(dá)的方法

基于臨空時間協(xié)同實現(xiàn)同時到達(dá)的方法實現(xiàn)簡單,適用于單艘艦艇和編隊對?；鹆ν粨?。

2.1 打擊單艦?zāi)繕?biāo)發(fā)射間隔時間控制

基于發(fā)射間隔時間控制,實現(xiàn)打擊單艦?zāi)繕?biāo)協(xié)同方法,通俗地講即為:同一攻擊波次中航程大的導(dǎo)彈先發(fā)射,航程小的導(dǎo)彈后發(fā)射,技術(shù)原理如下:

設(shè)同一攻擊波次中,攻擊同一目標(biāo)的最大航路數(shù)量為N條,設(shè)此N條航路的航程為Li,i=1,2,…,N。

按照第1節(jié)所述的方法,可求取考慮海洋環(huán)境參數(shù)化修正后的導(dǎo)彈飛行時間Tzi,i=1,2,…,N。

航路中飛行時間最大者為

Tmax=max(Tzi),i=1,2,…,N

(4)

于是,各航路發(fā)射間隔時間為

ΔTi=Tmax-Tzi,i=1,2,…,N

(5)

考慮到導(dǎo)彈最小發(fā)射間隔時間安全性約束,ΔTi取為

ΔTi=max{ΔTi,Tmin,Tc},i=1,2,…,N

(6)

Tmin為在保證武器發(fā)射安全性前提下,武器系統(tǒng)允許的導(dǎo)彈最小發(fā)射間隔時間,由武器系統(tǒng)戰(zhàn)技性能決定。Tc為考慮多枚導(dǎo)彈連射時戰(zhàn)術(shù)決策要求的最小發(fā)射間隔時間,由上級指揮控制系統(tǒng)效率決定。

設(shè)打擊目標(biāo)要求的臨空時刻為T0,則首枚導(dǎo)彈的發(fā)射時刻為

(7)

同一攻擊波次中,其他各枚導(dǎo)彈的發(fā)射時刻為

(8)

按式(8)計算的發(fā)射時刻控制導(dǎo)彈發(fā)射,即可保證同一攻擊波次中各枚導(dǎo)彈同時到達(dá)目標(biāo)。

現(xiàn)代艦艇超低空反導(dǎo)(硬抗擊)和電子對抗(軟抗擊)手段日臻完善,反艦導(dǎo)彈突防難度很大,單一波次多彈同時到達(dá)攻擊,可能也難以達(dá)成作戰(zhàn)目的。因此,本文還考慮了多波次協(xié)同攻擊的問題,各波次之間的攻擊時間間隔ΔTbj(j=1,2,…,K)還應(yīng)當(dāng)滿足以下條件:

ΔTbj=min{Tqh,Tzh,Tzt},j=1,2,…,K

(9)

式中,K為打擊同一目標(biāo)的最大攻擊波次數(shù)量;Tqh為艦載防空導(dǎo)彈火力通道切換時間;Tzh為艦載末端反導(dǎo)火炮轉(zhuǎn)火時間;Tzt為艦載電子戰(zhàn)武器再裝填時間。

2.2 打擊編隊目標(biāo)同時臨空時間控制

在臨空時間協(xié)同方面,打擊編隊目標(biāo)不同于打擊單目標(biāo)。針對艦艇編隊目標(biāo),在實現(xiàn)多方向同時臨空飽和攻擊戰(zhàn)術(shù)時,不能按照單艦?zāi)繕?biāo)的協(xié)同方法進(jìn)行計算[16],而必須考慮外圍艦艇與中心艦艇的協(xié)同存在距離差和臨空時間差問題。如果采用按照單一目標(biāo)的協(xié)同臨空時間,則打擊中心目標(biāo)的導(dǎo)彈需提前發(fā)射,先穿過外圍艦艇構(gòu)成的火力防御體系,有可能被外圍艦艇攔截,打擊中心艦艇使用導(dǎo)彈的突防概率將下降,這不符合飽和攻擊要求,也不符合對編隊目標(biāo)協(xié)同攻擊的需求[17-18]。

如打擊編隊內(nèi)部艦艇的導(dǎo)彈,為了避免被編隊外圍艦艇攔截,應(yīng)根據(jù)不同打擊目標(biāo)設(shè)置不同的臨空時間和臨空時間差值。通常,打擊編隊中心目標(biāo)和外圍目標(biāo)的臨空時間要有時間差,打擊近距目標(biāo)和遠(yuǎn)距目標(biāo)也要有時間差,最終實現(xiàn)對編隊的同時臨空協(xié)同攻擊。

基于臨空時間協(xié)同的打擊編隊目標(biāo)實施方法,本文按照以下方法實現(xiàn):

步驟1先計算出來的編隊中心,以此中心點為參照系,計算編隊各目標(biāo)到中心點的距離,以此距離為突防臨空的距離差值。

設(shè)編隊中心點經(jīng)緯度為(λ0、φ0),編隊中心點取編隊內(nèi)各艦艇位置的經(jīng)緯度(λi,φi)均值點具體計算公式為

(10)

(11)

步驟2將每個目標(biāo)到編隊中心的距離差值ΔLi(i=1,2,…,M)與攻擊該目標(biāo)的導(dǎo)彈航路航程Li(i=1,2,…,N)相加,獲得各導(dǎo)彈航路到編隊中心的航程值。

步驟3以各航路到編隊中心點的航程為計算條件,計算按照式(1)～式(3)計算考慮海洋環(huán)境參數(shù)化修正后的各航路對應(yīng)的導(dǎo)彈飛行時間Tzi,i=1,2,…,N,按照第2.1節(jié)所述打擊單艦?zāi)繕?biāo)協(xié)同方法,選取最長飛行時間Tmax對應(yīng)的航路為首發(fā)導(dǎo)彈航路。以Tmax為基準(zhǔn),按式(5)～式(8),可依次計算同一攻擊波次中其他導(dǎo)彈的發(fā)射時間間隔ΔTi,并進(jìn)一步可獲得每條航路對應(yīng)的導(dǎo)彈發(fā)射時間Ti,形成多發(fā)導(dǎo)彈同時臨空攻擊編隊目標(biāo)的導(dǎo)彈發(fā)射時間和發(fā)射順序。

3 基于航程相等實現(xiàn)同時到達(dá)的方法

同一攻擊波次各枚導(dǎo)彈航程相等是實現(xiàn)導(dǎo)彈同時到達(dá)的有效手段。與發(fā)射間隔控制方法相比,其最大的優(yōu)勢在于:同一攻擊波次的導(dǎo)彈可同時發(fā)射,導(dǎo)彈發(fā)射完畢后,發(fā)射平臺可迅速撤離,有助于保證發(fā)射平臺安全。特別地,對飛行中有操控要求的導(dǎo)彈而言,航程相等可保證齊射的多枚導(dǎo)彈同時抵達(dá)操控區(qū),以滿足操作手可同時操控多枚導(dǎo)彈的要求。

基于航程相等實現(xiàn)同時到達(dá)的航路規(guī)劃計算實現(xiàn)方法有兩種:一是規(guī)劃的航路沿彈目連線嚴(yán)格對稱,二是采用動態(tài)橢圓的幾何方法。

3.1 航路對稱法

如圖1所示,考慮導(dǎo)彈航路有3個航路點的情況,A點和A′點為末段進(jìn)入方向上的點,滿足|OA|≥L0,|OA′|≥L0(L0為最后一個航路點到目標(biāo)的最近距離約束),B點為2條航路共同航路點,OB連線為各個末段進(jìn)入方向的中間方向,A點和A′點關(guān)于OB連線對稱,G為A點和A′連線的中點。

圖1 三航路點對稱航路示意圖Fig.1 Route symmetry diagram of three navigating points

算法實現(xiàn)步驟設(shè)計如下。

步驟1令2條航路的末段進(jìn)入方向分別為α1和α2, 則中心對稱線方位為

(12)

步驟2求取2個末段進(jìn)入方向的夾角的半角為

(13)

步驟3求取∠ABG的最大允許角度為

(14)

式中,θmin為導(dǎo)彈航路最小轉(zhuǎn)角,為常數(shù)。

步驟4AA′中點G到第一航路點B的最小距離為

(15)

于是,可求得第一航路點B。

步驟5判斷條件|AB|>L1(L1為相鄰兩個航路點之間的最小距離約束,為常數(shù))是否成立,若不成立則伸長B點,伸長步長為e1(常數(shù)),直到條件|AB|>L1滿足或總航程S超出限制為止。

步驟6若步驟5計算成功且航路AB段或AB′段碰障礙,則伸長A點,伸長步長為e2(常數(shù)),重新進(jìn)行上述步驟1～步驟6,直到碰障規(guī)避成功或總航程S超出限制為止。

若上述三航路點對稱計算不成功,則進(jìn)行4個航路點對稱計算,如圖2所示。

圖2 四航路點對稱航路示意圖Fig.2 Route symmetry diagram of four navigating points

4個航路點算法是在3個航路點算法基礎(chǔ)上的改進(jìn),即B點繞A點向偏離中心對稱軸的方向轉(zhuǎn)動得到新的航路點B點和B′點,然后再按第3.1節(jié)描述的算法計算2條航路的共同航路點C,算法設(shè)計如下:

步驟6.1判斷A點位于中心對稱軸的左側(cè)還是右側(cè),若A點位于左側(cè),則順時針轉(zhuǎn)動B點,否則逆時針轉(zhuǎn)動B′點。

步驟6.2按無障礙物情況規(guī)劃出初始航路OAB,計算B點相對A點和A′點的方位AngAB和AngA′B,此兩方位即分別為B點和B′點偏離中心對稱軸轉(zhuǎn)動的初始方位。

步驟6.3進(jìn)入B點轉(zhuǎn)動循環(huán),轉(zhuǎn)動B點和B′點,轉(zhuǎn)動步長ε1(常數(shù))。

設(shè)轉(zhuǎn)動后B點和B′點的中點為G,則G點經(jīng)緯度坐標(biāo)按如下式計算，即

(16)

步驟6.4當(dāng)B點和B′點分別轉(zhuǎn)過OA和OA′連線方向時,進(jìn)入B點伸長循環(huán),按B點和B′點轉(zhuǎn)動后的方位伸長B點和B′點,伸長步長為e1。然后重新進(jìn)入步驟6.1—B點轉(zhuǎn)動循環(huán)。B點伸長循環(huán)中,如果BC段航路或B′C′段航路碰障礙,則說明B點伸長后碰障,伸長也無意義,則退出B點伸長級循環(huán),進(jìn)入B點轉(zhuǎn)動級循環(huán)。

步驟6.5若步驟6.3或步驟6.4任意一步中規(guī)劃的航路規(guī)避障礙物成功,則退出B點伸長和轉(zhuǎn)動循環(huán),并輸出航路規(guī)劃結(jié)果,規(guī)劃計算結(jié)束。

3.2 動態(tài)橢圓法

航路嚴(yán)格對稱方法實現(xiàn)簡單,但規(guī)避障礙物能力有限,當(dāng)作戰(zhàn)海域地理環(huán)境較復(fù)雜或需要規(guī)避的障礙物較多時,航路對稱法難以規(guī)劃成功[19-20]。因此,本文還提出了基于動態(tài)橢圓的幾何計算方法,并可保證在規(guī)避復(fù)雜障礙物條件下,航路的航程精確相等。 動態(tài)橢圓幾何方法描述如下:

橢圓的幾何特性:橢圓上的點到其兩個焦點的距離之和相等。如圖3所示,若視發(fā)射點和目標(biāo)點為橢圓的兩個焦點,則當(dāng)導(dǎo)彈航路只有一個航路點且位于由發(fā)射點和目標(biāo)點確定的橢圓軌跡上時,則任意兩條航路的航程嚴(yán)格相等。

圖3 橢圓幾何法示意圖Fig.3 Ellipse geometric method diagram

但考慮到導(dǎo)彈實際性能約束條件,純橢圓幾何方法的工程應(yīng)用價值有限:受導(dǎo)彈發(fā)射扇面、導(dǎo)彈機動性能以及末段搜捕目標(biāo)性能約束等因素影響,當(dāng)航路只有1個航路點時,該航路點到目標(biāo)的距離(P1B或P2B)、航路點到發(fā)射點的距離(P1A或P2A)以及航路的轉(zhuǎn)彎角(∠AP1B或∠AP2B),在設(shè)計上均有嚴(yán)格的約束,那么橢圓上航路點P1、P2的選擇空間就極其有限。尤其當(dāng)航路的攻擊進(jìn)入角度(∠P1BA、∠P2BA)有特定要求且航路上有地理障礙物時,基于純橢圓幾何的航路規(guī)劃方法則完全失效。

針對上述橢圓幾何方法的缺陷,提出了動態(tài)橢圓航路規(guī)劃算法,如圖4所示。考慮飛機平臺協(xié)同對海突擊的情況,F點為飛機起飛位置,A點和A′點為末段進(jìn)入方向上的點,B點、B′點和C點均為航路點。

圖4 動態(tài)橢圓法規(guī)劃航路示意圖Fig.4 Dynamic ellipse method diagram

算法實現(xiàn)步驟設(shè)計如下:

步驟1求解F點相對于2個末段進(jìn)入點A點和A′點的方位,分別記為α1和α2。

步驟2進(jìn)入B點和B′點伸長循環(huán)。求解F點—末段進(jìn)入點A—目標(biāo)點O構(gòu)成的夾角∠FAO,并確定B點和B′點轉(zhuǎn)動的初始方位β1和β2,順時針轉(zhuǎn)動時,β1和β2值為正。

(18)

然后,按B點和B′點轉(zhuǎn)動的初始方位β1和β2,及|AB|=|AB′|=L1的條件計算得到新的B點和B′點,再按第3.1節(jié)中給出的算法計算得到B點和B′點的中點G點。

設(shè)兩航路共同點為C。求出B點相對于G的方位為AngGC,則C點相對于G點的方位滿足:

(19)

取|GC=nL1|,n為L1倍數(shù),為常數(shù)，則可求得C點經(jīng)緯度。

接著,判斷C點相對于G點方位的合法性。分別求取C點到B點的距離|CB|、C點到A點距離|CA|,如果|BA|≥|CA|成立,則C點相對于G點的方位初值不合法,需要重新計算C點，即

(20)

再按|GC=n*L1|,重新計算C點經(jīng)緯度,重新計算的C點必是合法的,如圖4所示。

步驟4進(jìn)行動態(tài)橢圓二次優(yōu)化:取兩條航路中航程大的那條,剔除其航路點B或B′,為保持兩條航路航程相等,以A′點和C點為橢圓焦點,建立動態(tài)橢圓,B/B′點即為橢圓上的點,滿足:

|A′B′|+|B′C′|=|AC|

(21)

步驟5如果航路不碰障礙且航路二次優(yōu)化成功,輸出航路規(guī)劃結(jié)果,規(guī)劃計算結(jié)束。

步驟6當(dāng)以下條件任意一個滿足時,退出B點伸長循環(huán):

①航路計算成功;

②航路碰障;

③航路總航程S超出限制。

步驟7滿足下列條件之一時,進(jìn)入B點轉(zhuǎn)動循環(huán):

①航路碰障礙;

②航路二次優(yōu)化不成功。

B點轉(zhuǎn)動循環(huán)角度步長ε2,B點轉(zhuǎn)動后的方位為

(22)

(23)

步驟8B點轉(zhuǎn)動后,若規(guī)劃成功,則說明航路二次優(yōu)化規(guī)劃成功,規(guī)劃計算結(jié)束。

若下列條件之一被滿足,說明航路二次優(yōu)化規(guī)劃不成功,則無解:

②航路總航程S超出限制。

動態(tài)橢圓法適用性強,當(dāng)?shù)乩碚系K比較復(fù)雜、禁飛區(qū)數(shù)量比較多時,航路對稱算法可能無法規(guī)劃,動態(tài)橢圓航路規(guī)劃算法由于其靈活性強,在復(fù)雜障礙物情況下,仍可以規(guī)劃出經(jīng)二次優(yōu)化后的導(dǎo)彈可飛的航路并可保證攻擊同一目標(biāo)的任意2條航路的航程精確相等。

3.3 航程精確計算方法

精確計算航路的航程是實現(xiàn)各條航路航程相等的基礎(chǔ),本文采用幾何解析方法計算多航路點航路的航程,計算方法描述如下。

考慮多個航路點情況,如圖5所示。

圖5 航程計算示意圖Fig.5 Schematic diagram of voyage calculation

N1、N2、N3、N4為規(guī)劃的航路點,M為目標(biāo)點,O1、O2、O3為航路轉(zhuǎn)彎圓心,R為航路轉(zhuǎn)彎半徑,θ1、θ2、θ3為航路轉(zhuǎn)彎圓心角。

則圖5所示的4航路點航程L4滿足:

(24)

由圖5可導(dǎo)出4個航路點巡航段航程L4計算公式為

L4=|N1N2|+|N2N3|+|N3N4|+|N4M|+

(25)

進(jìn)一步,可推廣到n個航路點航程Ln精確計算公式為

Ln=|N1N2|+|N2N3|+|Nn-2Nn-1|+…+|Nn-1M|+

(26)

4 仿真驗證

首先,在50～550 km范圍內(nèi),考察了海洋環(huán)境因素對導(dǎo)彈飛行時間修正的效果,數(shù)據(jù)統(tǒng)計結(jié)果如表1所示。

表1 飛行時間誤差數(shù)據(jù)統(tǒng)計

表中,Di表示飛行距離;ΔTw表示無修正預(yù)估飛行時間與理論飛行時間偏差;ΔTy表示有修正預(yù)估飛行時間與理論飛行時間偏差。飛行時間統(tǒng)計分析表明:在考慮海洋環(huán)境因素對導(dǎo)彈飛行時間修正條件下,按導(dǎo)彈飛行時間修正式(3)計算得到的導(dǎo)彈飛行時間計算值與導(dǎo)彈飛行時間的平均偏差下降到11 s,較無海洋環(huán)境因素修正的44 s平均誤差有了大幅度的提高,平均臨空時間差值減小了75%,按照對海洋環(huán)境參數(shù)化修正后的導(dǎo)彈飛行時間控制精度得到了大幅提高,有效地提高了多彈同時臨空的協(xié)同攻擊效果。

然后,計算了同一平臺發(fā)射的多枚導(dǎo)彈分別沿不同的航路攻擊同一目標(biāo)的發(fā)射時間間隔。如表2所示。

表2 發(fā)射間隔時間計算結(jié)果

由表2統(tǒng)計結(jié)果可見,對擁有較大動力航程的亞聲速反艦導(dǎo)彈(如航程大于150 km),同一攻擊波次中首枚導(dǎo)彈發(fā)射到末枚導(dǎo)彈發(fā)射之間的間隔可長達(dá)數(shù)分鐘,這段時間內(nèi),被攻擊艦艇有充分的時間進(jìn)行反導(dǎo)抗擊準(zhǔn)備,逐一抗擊來襲導(dǎo)彈。若控制發(fā)射間隔時間,使得導(dǎo)彈同時到達(dá)目標(biāo),則有助于實現(xiàn)飽和攻擊,提高突防效果。

最后,在考慮障礙物規(guī)避條件下(障礙物以多邊形描述),按三航路點對稱算法、四航路點對稱算法和動態(tài)橢圓法分別計算了導(dǎo)彈突防航路,其中任意2條航路航程近似相等(誤差小于10 m),以實現(xiàn)導(dǎo)彈同時臨空,計算結(jié)果如表3所示。采用航路對稱法和動態(tài)橢圓法,均具有較高的計算精度,任意2條航路的總航程差小于10 m。此外,由于本文所提出的航程相等的航路規(guī)劃方法,均是基于解析幾何的數(shù)學(xué)方法,有唯一解,且無迭代計算,具有支持快速計算的特點,非常適于各類平臺火控實時航跡解算,也可適用于彈載或機載綜控機飛行中在線航跡規(guī)劃計算。

表3 航程相等實現(xiàn)同時臨空計算結(jié)果

5 結(jié) 論

反艦導(dǎo)彈飽和攻擊,同時到達(dá)是必要條件。本文基于導(dǎo)彈的發(fā)射順序和發(fā)射間隔時間控制和通過航路規(guī)劃實現(xiàn)導(dǎo)彈的航程相等兩種思路,在考慮海洋環(huán)境因素對導(dǎo)彈飛行時間影響條件下,提出了反艦導(dǎo)彈“同時到達(dá)”打擊單艦和編隊目標(biāo)的臨空時間協(xié)同方法和航路規(guī)劃方法,可為指揮員組織實施反艦導(dǎo)彈協(xié)同飽和攻擊提供決策參考,提高導(dǎo)彈攻擊和兵力行動的效率。

參考文獻(xiàn)：

[1] 張金春, 張晶, 張書宇. 反艦導(dǎo)彈對海目標(biāo)飽和攻擊輔助決策系統(tǒng)研究[J]. 戰(zhàn)術(shù)導(dǎo)彈技術(shù), 2014 (2): 21-27.

ZHANG J C, ZHANG J, ZHANG S Y, et.al. Research on assistant decision-making syetem about saturation attacks of anti-ship missile to sea targets[J].Tactical Missile Technology,2014(2):21-27.

[2] 劉樹衍,劉暢平,王航宇,等.偽人工勢場在多彈飽和攻擊航路規(guī)劃中的應(yīng)用研究[J].海軍工程大學(xué)學(xué)報,2015,27(5):108-112.

LIU S K, LIU C P, WANG H Y, et al. Coordinative route planning of anti-ship missile saturation attack based on pseudo artificial potential field[J]. Journal of Naval University of Engineering, 2015, 27(5): 108-112.

[3] 曾家有, 王國衛(wèi), 鐘建林, 等. 多平臺艦艦導(dǎo)彈飽和攻擊幾個協(xié)同問題與模型[J]. 兵工學(xué)報, 2014, 35(2):256-261.

ZENG J Y, WANG G W, ZHONG J L, et al. Research on cooperative saturation attack problems and models of ship-to-ship missiles from multi-ship platforms[J]. Acta ArmamentarⅡ, 2014, 35(2):256-261.

[4] WANG G H, SUN X F, ZHANG L P, et al. Saturation attack based route planning and threat avoidance algorithm for cruise missiles[J]. Journal of Systems Engineering and Electronics, 2011, 22(6):948-953.

[5] YAO P, WANG H L, SU Z K. Real-time path planning of unmanned aerial vehicle for target tracking and obstacle avoidance in complex dynamic environment[J]. Aerospace Scienceand Technology, 2015, 47(4):269-279.

[6] BABEL L. Flight path planning for unmanned aerial vehicles with landmark-based visual navigation[J]. Robotics and Autonomous Systems, 2014,62(2):142-150.

[7] XU C F, DUAN H B, LIU F. Chaotic artificial bee colony approach to uninhabited combat air vehicle (UCAV) path planning[J]. Aerospace Science & Technology, 2010,14(8):535-541.

[8] ZHANG X Y, DUAN H B. An improved constrained differential evolution algorithm for unmanned aerial vehicle global route planning[J]. Applied Soft Computing, 2015,26(C):270-284.

[9] LUA W C, LEE M T, WANG M W. Route planning for light-sport aircraft in constrained airspace[C]∥Proc.of the 7th Asian-Pacific Conference on Aerospace Technology and Science,2013:140-146.

[10] ZHANG M, SU C, LIU Y, et al. Unmanned aerial vehicle route planning in the presence of a threat envoronment based on a virtual globe platform[J]. ISPRS International Journal of GEO-Infomation, 2016, 184(5): 1-30.

[11] 袁利平, 陳宗基, 周銳, 等. 多無人機同時到達(dá)的分散化控制方法[J]. 航空學(xué)報, 2010, 31(4): 797-805.

YUAN L P, CHEN Z J, ZHOU R, et al. Decentralized control for simultaneous arrival of multiple uavs[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(4): 797-805.

[12] 曾家有, 宋友凱, 車志宇. 基于航路規(guī)劃的反艦導(dǎo)彈發(fā)射順序和間隔研究[J]. 航天控制, 2009, 27(2):22-25.

ZENG J Y, SONG Y K, CHE Z Y. Research on launch sequence and interval onf anti-ship missile based on route planning[J]. Aerospace Control, 2009, 27(2):22-25.

[13] 王建國, 徐曉剛, 高軍亮, 等. 航路規(guī)劃技術(shù)在反艦導(dǎo)彈遠(yuǎn)程打擊中的應(yīng)用[J]. 飛航導(dǎo)彈, 2006 (9): 23-25.

WANG J G, XU X G, GAO J L, et al. The application of route planning technology in long range attack of anti-ship missile[J]. Aerodynamic Missile Journal, 2006(9):23-25.

[14] LIU G, LAO S Y, HOU L L, at el. OARPER-MAFO algorithm for anti-ship missile path planning[J]. Aerospace Science and Technology, 2015,47:135-145.

[15] LI F, HAO B, ZHAO J H, et al. Cruise missile multiple route planning based on hybrid particl swarm optimization[J].Journal of Beijing Institute onf Technology,2013,22(3):354-360.

[16] LIU G, LAO S Y, HOU L L, et al. A simulation system of anti-ship missile path planning oriented ship formation[J]. Advanced Materials Research, 2012, 532(6): 645-649.

[17] DESARAJU V R, HOW J P. Decentralized path planning for multi-agent teams with complex constraints[J]. Autonomous Robots, 2012, 32(4): 385-403.

[18] TAM C K, BUCKNALL R. Cooperative path planning algorithm for marine surface vessels[J]. Ocean Engineering, 2013, 57(2): 25-33.

[19] WANG Q J, PENG J, NI B H, et al. Route planning for military aircrafts based on waypoint searching method[J]. Computer Engineering, 2014, 40(2):280-283.

[20] LIU Y.Design and implementation of obstacle assessment in flight route[J].Command and Information System and Technology,2012,3(6):25-29.