摘 要：

針對滑翔制導(dǎo)炮彈多炮齊射打擊固定目標(biāo)的協(xié)同攻擊場景，為充分發(fā)揮整體控制能力，增強(qiáng)打擊效果，研究著眼于彈群的協(xié)同方案彈道規(guī)劃問題。首先，考慮到攻擊不同距離目標(biāo)時飛行方案的側(cè)重點與彈道特性差異，為滑翔制導(dǎo)炮彈進(jìn)行戰(zhàn)區(qū)劃分，并在此基礎(chǔ)上建立尋優(yōu)指標(biāo)因戰(zhàn)區(qū)而異的全彈道協(xié)同規(guī)劃模型。其次，針對實戰(zhàn)中的戰(zhàn)斗隊形配置問題，提出一種簡單通用的協(xié)作范圍確定方法。隨后，針對現(xiàn)有協(xié)同方式在處理多階段規(guī)劃問題時存在的適用性或最優(yōu)性方面的局限，提出一種雙層自適應(yīng)協(xié)同策略。最后，分析特殊作戰(zhàn)任務(wù)引入的額外約束條件對戰(zhàn)區(qū)劃分與協(xié)作范圍的影響。仿真結(jié)果表明，所提方法可以簡捷有效地確定可行工況，輔助戰(zhàn)斗隊形的配置;同時，相比于傳統(tǒng)的分布式策略，在所提的協(xié)同策略作用下，不同戰(zhàn)區(qū)內(nèi)的方案彈道均具有優(yōu)越性。

關(guān)鍵詞：

彈道規(guī)劃; 多炮齊射; 同時彈著; 戰(zhàn)區(qū)劃分; 協(xié)同策略; 協(xié)作范圍

中圖分類號：

TJ 413.+6

文獻(xiàn)標(biāo)志碼： A""" DOI：10.12305/j.issn.1001-506X.2024.09.26

Cooperative trajectory programming and cooperation range determination method for

gliding-guided projectiles in multi-artillery-salvo-shooting scenario

YIN Qiulin， CHEN Qi*， WANG Zhongyuan， WANG Qinghai

（School of Energy and Power Engineering， Nanjing University of Science and Technology，

Nanjing 210094， China）

Abstract：

Aiming at the scenario of coordinated attack by gliding-guided projectiles in multi-artillery-salvo-shooting to strike fixed targets， the trajectory planning problem for collaboration schemes focusing on projectile groups is investigated to fully leverage the overall control capacity and amplify the strike outcomes. Firstly， considering the differences in focal points and trajectory characteristics of flighting scheme when engaging targets at various distances， operational zone for gliding-guided projectiles is divided， and on this basis， a full trajectory cooperative programming model with optimization indicators varying by zone is established. Secondly， in response to the challenges posed by combat formation configurations in practical warfare， a universally applicable and straightforward method for delineating the scope of collaboration is introduced. Then， to overcome the restrictions pertaining to applicability or optimality encountered by existing methods of collaboration within multi-stage planning contexts， a bi-level adaptive cooperation strategy is proposed. Finally， the influence of additional constraints introduced by specific operational tasks on both the allocation of operational zone and the extent of cooperation is analyzed. Simulation results show that the proposed method can swiftly and effectively ascertain viable operating scenarios， thus aiding in combat formation arrangements; additionally， compared with conventional distributed strategies， it is evident that the proposed collaborative strategy yields superior trajectories within various operational zones.

Keywords：

trajectory programming; multi-artillery-salvo-shooting; simultaneous attack; operational zone division; collaborative strategy; cooperation range

0 引 言

滑翔制導(dǎo)炮彈射程遠(yuǎn)、制導(dǎo)精度高、毀傷效果好且成本低廉，打擊固定或低速機(jī)動目標(biāo)時具有顯著優(yōu)勢，是現(xiàn)代戰(zhàn)爭中兵器技術(shù)發(fā)展的重點方向之一。但其體積小、舵控能力有限，通常采用方案制導(dǎo)的方式，方案彈道規(guī)劃在滑翔制導(dǎo)炮彈總體設(shè)計中發(fā)揮著關(guān)鍵作用，方案彈道設(shè)計的好壞將直接影響滑翔制導(dǎo)炮彈的作戰(zhàn)效果［1］。當(dāng)前戰(zhàn)場環(huán)境愈發(fā)復(fù)雜，攻防對抗日趨嚴(yán)密，多樣化的作戰(zhàn)任務(wù)很難依靠單一彈丸獨立完成，多彈協(xié)同的攻擊模式便應(yīng)運(yùn)而生。由于其突防能力強(qiáng)、打擊范圍廣、抗電子干擾能力強(qiáng)且效費(fèi)比高，逐漸成為近年來的研究熱點［2-5］，同時彈著也作為新的約束條件融入了方案彈道的規(guī)劃過程中［6-8］。不同炮位的多發(fā)彈丸從不同的距離和方向同時發(fā)射，可以分散防御系統(tǒng)的火力，提高彈丸生存能力;多彈同時命中，可以提高打擊的強(qiáng)度，確保毀傷效果。結(jié)合工程實踐中常見的滑翔制導(dǎo)炮彈發(fā)射條件與作戰(zhàn)場景，研究一種在多炮齊射場景下，可以充分發(fā)揮彈群控制潛力，并實現(xiàn)同時彈著的協(xié)同方案，具有重要的現(xiàn)實意義。

滑翔制導(dǎo)炮彈射程遠(yuǎn)，作戰(zhàn)空域大。在彈型確定的情況下，針對不同距離處的攻擊目標(biāo)，其飛行方案具有不同的側(cè)重。例如，目標(biāo)較近時，重點保證存速以增加毀傷效果;目標(biāo)較遠(yuǎn)時，則更加重視降低控制能量的損耗以確保射程;目標(biāo)位置居中時，則兩者兼而有之。因此，有必要對滑翔制導(dǎo)炮彈進(jìn)行戰(zhàn)區(qū)劃分，并且根據(jù)不同戰(zhàn)區(qū)內(nèi)的飛行方案的側(cè)重建立評價指標(biāo)不同的彈道規(guī)劃模型。

為實現(xiàn)同時彈著的目標(biāo)，協(xié)調(diào)彈群末端攻擊時刻顯得尤為關(guān)鍵。文獻(xiàn)［9］采用人為設(shè)定末端攻擊時刻的開環(huán)方式規(guī)劃能量最優(yōu)協(xié)同彈道，這種方式對經(jīng)驗依賴較大，若設(shè)定時間不合理將導(dǎo)致規(guī)劃效果不佳甚至失敗，而且其本質(zhì)上是各彈的獨立規(guī)劃，并非真正意義上的協(xié)同。文獻(xiàn)［10］采用領(lǐng)彈-從彈策略，選取最長的飛行時間作為標(biāo)準(zhǔn)，其他各彈以此進(jìn)行協(xié)調(diào)。這種方式對于從彈的控制能力要求較高，常用于導(dǎo)彈的末制導(dǎo)段，不適用于制導(dǎo)炮彈。文獻(xiàn)［11］采用分布式策略對導(dǎo)彈進(jìn)行獨立規(guī)劃，計算各自飛行時間，并求取平均值作為依據(jù)進(jìn)行協(xié)調(diào)。這種方式允許各發(fā)彈丸擁有獨立的規(guī)劃目標(biāo)，靈活性較強(qiáng)。但多數(shù)情況下，平均值并不是最優(yōu)解，即便是加權(quán)平均，權(quán)重的設(shè)計也依賴于經(jīng)驗，合理性值得商榷。針對上述問題，文獻(xiàn)［12］將所有動力學(xué)方程聯(lián)立于統(tǒng)一時間變量，避免飛行時間的協(xié)調(diào)，且末端時刻由算法生成，保證了最優(yōu)性。這種統(tǒng)一自變量的方式，局限性在于其僅適用于單一階段的協(xié)同問題，在處理并非全程有動力飛行的多階段滑翔制導(dǎo)炮彈彈道規(guī)劃時，無法確保各發(fā)彈丸自身各階段的獨立性。但是，這種集中式的策略與本文彈群整體規(guī)劃的思想非常契合，如果能夠把分布式策略的靈活性與集中式策略的最優(yōu)性結(jié)合起來，設(shè)計一種可以處理多階段問題的協(xié)同策略，將對滑翔制導(dǎo)炮彈的方案設(shè)計有所裨益。

針對多炮齊射場景下的協(xié)同規(guī)劃問題，現(xiàn)有文獻(xiàn)大多深耕于規(guī)劃算法的改進(jìn)與協(xié)同策略的研究，以某些可行的工況作為算例進(jìn)行仿真來驗證所提方法的有效性，但對于不可行工況試錯的過程與選擇可行工況的方法鮮有提及。倘若戰(zhàn)場環(huán)境發(fā)生變化，規(guī)劃結(jié)果可能會不理想甚至規(guī)劃失敗。針對具體的固定目標(biāo)，確定出特定工況下的最優(yōu)協(xié)同方案固然重要，但不同炮口位置的成員具體可以與多大范圍內(nèi)的其他組員進(jìn)行協(xié)作，從而確定戰(zhàn)斗隊形，也是工程應(yīng)用所關(guān)心的實際問題。同時，不同戰(zhàn)區(qū)內(nèi)彈丸協(xié)調(diào)飛行時間的方式與能力大小，協(xié)同彈道所呈現(xiàn)的特征，以及不同的作戰(zhàn)任務(wù)需求對戰(zhàn)區(qū)劃分與協(xié)作范圍的影響也值得深入探討。截至目前，針對此類問題的研究仍顯不足。

根據(jù)上述分析，本文以某型滑翔制導(dǎo)炮彈為研究對象，對多炮齊射場景下的協(xié)同彈道規(guī)劃問題展開了研究。首先，根據(jù)打擊距離進(jìn)行了戰(zhàn)區(qū)劃分，并建立著眼于彈群整體的協(xié)同規(guī)劃模型;隨后，提出一種簡單通用的協(xié)作范圍確定方法以及一種可以處理多階段問題的協(xié)同策略;最后，研究分析不同攻擊距離協(xié)同彈道的特征以及協(xié)作范圍，并初步探討特殊任務(wù)需求（如灌頂打擊）對戰(zhàn)區(qū)劃分與協(xié)作范圍的影響，為后續(xù)相關(guān)研究提供了思路。通過數(shù)值仿真驗證了所提方法的有效性與優(yōu)越性。

1 多炮齊射彈道規(guī)劃問題描述

1.1 滑翔制導(dǎo)炮彈全彈道飛行運(yùn)動模型

在初步設(shè)計階段，暫時不考慮末制導(dǎo)過程。彈丸從發(fā)射到命中目標(biāo)，飛行方案主要分為4個階段：起飛段、助推段、爬升段和滑翔段。各階段及其劃分節(jié)點如圖1所示。

在設(shè)計方案彈道時，為了把握住問題的本質(zhì)，以便研究飛行彈道的特性，可將制導(dǎo)炮彈視作一個可操縱的質(zhì)點。假設(shè)彈丸在飛行過程中滿足瞬時平衡條件、氣象環(huán)境標(biāo)準(zhǔn)、控制系統(tǒng)理想工作，且無外界因素的干擾［13］。其質(zhì)心運(yùn)動方程組為

3 協(xié)作范圍確定

當(dāng)目標(biāo)位置固定時，為了確定當(dāng)前炮口處的成員具體可以與多大范圍內(nèi)的其他組員進(jìn)行協(xié)作，從而輔助戰(zhàn)斗隊形的配置;并進(jìn)一步研究特殊作戰(zhàn)任務(wù)需求（如灌頂打擊）引入的額外約束條件對該范圍的影響，本文提出一種簡單通用的協(xié)作范圍確定方法。

3.1 攻擊時間范圍確定

為提高打擊效果，假設(shè)所有炮口在發(fā)射階段均對準(zhǔn)目標(biāo)，不考慮故意瞄歪再繞行的情況。此時，各發(fā)彈丸的飛行軌跡均近似在各自炮目連線所構(gòu)成的彈道面內(nèi)。故本節(jié)將各發(fā)制導(dǎo)炮彈的運(yùn)動模型定義在各自的坐標(biāo)系中，每個坐標(biāo)系的原點Oi取為各自的炮口位置，Oixi軸指向目標(biāo)方向為正，Oyi與Ozi軸的定義參照地面坐標(biāo)系。各發(fā)彈丸的射程定義為各自x方向的位移而非地面坐標(biāo)系中的x坐標(biāo)，后續(xù)展示協(xié)同飛行軌跡時會通過坐標(biāo)轉(zhuǎn)換的方式將所有彈丸的軌跡向地面坐標(biāo)系進(jìn)行統(tǒng)一，特此說明。已知本文研究的制導(dǎo)炮彈最小射程Rmin為20 km，記炮口與目標(biāo)之間的距離為R，根據(jù)第1.2節(jié)的描述，在不同戰(zhàn)區(qū)的約束條件下，以射程最大為目標(biāo)進(jìn)行方案彈道規(guī)劃：

minJ=-R（13）

可以得到，當(dāng)前作戰(zhàn)任務(wù)下的制導(dǎo)炮彈戰(zhàn)區(qū)范圍［Rs，min，Rs，max］，s∈{near，mid，far}。選取合適的步長將戰(zhàn)區(qū)分段，在每個節(jié)點N處，為約束條件增加當(dāng)前炮目間距RN，分別以最小、最大飛行時間為目標(biāo)進(jìn)行彈道規(guī)劃：

minJ1=tf（14）

minJ2=-tf（15）

可以得到當(dāng)前射程RN下的飛行時間范圍［tmin，RN，tmax，RN］。在仿真過程中觀察到，算法在規(guī)劃飛行時間最大時會出現(xiàn)控制指令高頻跳變的情況，且射程越小，這種情況越易發(fā)生。因此，在規(guī)劃過程中，以式（16）所示的動態(tài)系數(shù)kRN（本文取區(qū)間［0.618，1］對炮目間距RN∈［Rnear，min，Rfar，max］的線性映射）對式（3）中的幅值進(jìn)行修正，以留下部分裕量。得到RN對應(yīng)的tmax，RN后，將其作為約束條件加入規(guī)劃模型中進(jìn)行驗證，并以控制量光滑、不飽和為原則對不合理的數(shù)據(jù)進(jìn)行調(diào)整，如此可以從理論上確?；柚茖?dǎo)炮彈的控制能力足以覆蓋此飛行時間范圍。

kRN=（1－0.618）（RN－Rnear，min）Rfar，max－Rnear，min+0.618

αmax，RN=kRNαmax

βmax，RN=kRNβmax（16）

3.2 協(xié)作范圍確定

將第3.1節(jié)中所有節(jié)點的計算結(jié)果編制成表格并進(jìn)行插值計算（當(dāng)節(jié)點數(shù)量足夠多時，插值結(jié)果可近似作為真值處理），可以快速得到射程范圍內(nèi)任意炮位的打擊時間范圍?？紤]n發(fā)齊射彈丸，只要滿足關(guān)系：

max（tmin，Ri）≤min（tmax，Ri）， i=1，2，…，n（17）

即存在公共的攻擊時間范圍（此范圍已留有一定裕量，并非極限情況），便一定可以實現(xiàn)齊射協(xié)同攻擊。

同理，對于特定射程R*的炮位，只需在制導(dǎo)炮彈射程范圍之內(nèi)（超出范圍則取邊界值）向兩邊求取R1與R2，以滿足如下關(guān)系：

tmin，R*=tmax，R1

tmax，R*=tmin，R2（18）

便可以得到制導(dǎo)炮彈在R*處的具體協(xié)作范圍［R1，R2］，通過這樣的方式，可以快速確定出可行的齊射工況，避免了不斷的試錯過程。對于不同的作戰(zhàn)任務(wù)需求，只需修改相應(yīng)的約束條件，進(jìn)行相似的處理即可，具有一定的通用性。對仿真結(jié)果進(jìn)行觀察分析，可以得到不同約束條件對協(xié)作范圍的影響，為工程應(yīng)用提供一定的參考。

4 仿真驗證

為驗證本文所提方法的有效性與優(yōu)越性，并對確定協(xié)作范圍的過程進(jìn)行直觀展示，本節(jié)以某型滑翔制導(dǎo)炮彈為對象，模擬不同戰(zhàn)場環(huán)境進(jìn)行多炮齊射協(xié)同規(guī)劃仿真。彈道規(guī)劃算法采用hp自適應(yīng)偽譜法［32-33］，仿真參數(shù)如表1所示。

4.1 常規(guī)任務(wù)下不同戰(zhàn)區(qū)內(nèi)齊射協(xié)同規(guī)劃

以無額外約束條件的情況作為常規(guī)作戰(zhàn)任務(wù)場景，根據(jù)第3節(jié)所提的方法，首先代入相關(guān)約束條件與仿真參數(shù)，得到常規(guī)作戰(zhàn)任務(wù)要求下，滑翔制導(dǎo)炮彈的射程范圍為20～80 km。其中，死區(qū)20～26 km，近區(qū)26～44 km，中區(qū)44～72 km，遠(yuǎn)區(qū)72～80 km。不同射程對應(yīng)的飛行時間范圍如圖3所示。

從圖3結(jié)果可知，滑翔制導(dǎo)炮彈的飛行時間下界隨射程近似線性變化，而飛行時間上界存在一定范圍的浮動。近區(qū)與中區(qū)之內(nèi)浮動較為明顯，這是由于飛行方案的側(cè)重并非增程所致;而在遠(yuǎn)區(qū)內(nèi)則近似線性變化。整體上看，飛行時間上下界之間的差異隨射程先增后減，近似呈現(xiàn)兩邊低中間高的特征。說明在常規(guī)任務(wù)下，中部戰(zhàn)區(qū)發(fā)射的彈丸協(xié)調(diào)飛行時間能力更強(qiáng)，協(xié)作范圍更廣。不同射程對應(yīng)的具體協(xié)作范圍如圖4所示。

圖4結(jié)果印證了對圖3結(jié)果的分析。同時可以觀察到，死區(qū)之內(nèi)（20～26 km）彈丸協(xié)作范圍較小且變化劇烈，并不適合協(xié)同作戰(zhàn)。在工程實踐中，可根據(jù)圖4結(jié)果快速確定可行工況，輔助戰(zhàn)斗隊形的配置。

4.1.1 近區(qū)齊射協(xié)同規(guī)劃

本節(jié)模擬戰(zhàn)場環(huán)境，以圖4為參考，在圖2所示制導(dǎo)炮彈近區(qū)120°扇形攻擊區(qū)域內(nèi)，以40°為間隔（目的是從理論上均勻覆蓋整個作戰(zhàn)區(qū)域，實際中可以靈活變化，中、遠(yuǎn)區(qū)的處理相同），選取距離目標(biāo)31 km、33 km、35 km、37 km的4個位置組成戰(zhàn)斗隊形，進(jìn)行齊射協(xié)同規(guī)劃，仿真結(jié)果如圖5所示。

從圖5結(jié)果可知，近區(qū)協(xié)作范圍內(nèi)的彈丸成功實現(xiàn)了齊射協(xié)同攻擊，證明了本文所提方法的有效性。由于沒有增程需要，飛行方案側(cè)重于提高末端存速，近區(qū)內(nèi)的協(xié)同彈丸均采用最大射角以充分利用重力勢能，并利用負(fù)的舵偏使彈道下壓，協(xié)同彈道呈拋物線狀。由于同時彈著的約束，炮目間距越大的彈丸彈道曲率越大，落角也越大，消耗的控制能量也越多。

為進(jìn)一步體現(xiàn)本文方法的優(yōu)越性，與現(xiàn)有文獻(xiàn)中最為常見的分布式協(xié)同策略［28］（distributed cooperation strategy， DCS）進(jìn)行對比（其他策略更加不適合本文研究的多階段齊射攻擊場景，篇幅所限僅對比分布式策略），仿真結(jié)果如表2和表3所示。

從表2和表3的對比結(jié)果可以看出，DCS以平均值方式協(xié)調(diào)出的飛行時間并非最優(yōu)解，飛行時間減少了7.97 s，射角留有余量，點火時間推遲（開舵時間提前是因為全程飛行時間變短所致），彈道曲率較BACS小，并未充分發(fā)揮彈群的控制能力。因此，彈群的末端存速相比BACS降低了34.02 m/s。在本算例中，BACS的規(guī)劃目標(biāo)函數(shù)值相比DCS縮減了1.97%。

4.1.2 中區(qū)齊射協(xié)同規(guī)劃

本節(jié)模擬戰(zhàn)場環(huán)境，在中區(qū)內(nèi)選取距離目標(biāo)53 km、57 km、61 km、65 km的4個位置組成戰(zhàn)斗隊形，進(jìn)行齊射協(xié)同規(guī)劃，仿真結(jié)果如圖6所示。

從圖6結(jié)果可知，中區(qū)協(xié)作范圍內(nèi)的彈丸成功實現(xiàn)了齊射協(xié)同攻擊，證明了本文所提方法的有效性。一定的增程需求使彈道出現(xiàn)滑翔特征，彈丸落角相比近區(qū)降低。由于飛行方案兼顧存速與控制能量，中區(qū)內(nèi)的協(xié)同彈丸調(diào)節(jié)飛行時間的方式相比近區(qū)更加多樣，末發(fā)彈丸射角留有余量，彈丸存速均高于所設(shè)計的閾值。炮目間距越大的彈丸末端彈道曲率越大，落角也越大，消耗的控制能量也越多。類似地，為體現(xiàn)本文方法的優(yōu)越性，與分布式協(xié)同策略的對比結(jié)果如表4和表5所示。

從表4和表5的對比結(jié)果可以看出，DCS協(xié)調(diào)出的飛行時間并非最優(yōu)解，增加了4.01 s。這4 s左右的差異對炮目間距最小的首發(fā)彈丸產(chǎn)生了較大影響（觀察圖3可知，195 s的飛行時間已逼近此處滑翔制導(dǎo)炮彈的能力上限），迫使其晚點火、降低彈道高，犧牲了大量存速并消耗大量控制能量以完成時間的協(xié)同。在本算例中，BACS的規(guī)劃目標(biāo)函數(shù)值相比DCS縮減了84.72%，充分證明了以平均值方式協(xié)調(diào)飛行時間的劣勢所在。

4.1.3 遠(yuǎn)區(qū)齊射協(xié)同規(guī)劃

本節(jié)模擬戰(zhàn)場環(huán)境，在遠(yuǎn)區(qū)內(nèi)選取距離目標(biāo)73 km、75 km、77 km、79 km的4個位置組成戰(zhàn)斗隊形，進(jìn)行齊射協(xié)同規(guī)劃，仿真結(jié)果如圖7所示。

從圖7結(jié)果可知，遠(yuǎn)區(qū)協(xié)作范圍內(nèi)的彈丸成功實現(xiàn)了齊射協(xié)同攻擊，證明了本文所提方法的有效性。彈丸全部采用最大射角，不再是為了提高存速，而是為了滿足迫切的增程需求，彈道滑翔特征明顯，落角相比中區(qū)進(jìn)一步降低。由于飛行方案側(cè)重于降低控制能量的損耗，彈丸末端存速均為所設(shè)計的閾值。類似地，為體現(xiàn)BACS的優(yōu)越性，與DCS的對比結(jié)果如表6和表7所示。

從表6和表7的對比結(jié)果可以看出，DCS協(xié)調(diào)出的飛行時間為274 s（已逼近首發(fā)彈丸能力上限），相比BACS增加了12.12 s。此差異對前3發(fā)彈丸均產(chǎn)生了一定影響，并由更多的控制能量損耗來彌補(bǔ)。在本算例中，BACS的規(guī)劃目標(biāo)函數(shù)值相比DCS縮減了25.30%。

4.2 額外約束對作戰(zhàn)區(qū)域與協(xié)作范圍的影響

為進(jìn)一步研究額外約束條件對齊射協(xié)同作戰(zhàn)區(qū)域與協(xié)作范圍的影響，本節(jié)增加考慮針對坦克等頂部薄弱目標(biāo)進(jìn)行大落角打擊的任務(wù)需求，為協(xié)同彈丸增加落角約束|θf|min=85°進(jìn)行仿真，過程與第4.1節(jié)相仿。得到灌頂打擊任務(wù)下，滑翔制導(dǎo)炮彈的戰(zhàn)區(qū)劃分如下：近區(qū)26～44 km，中區(qū)44～67 km，遠(yuǎn)區(qū)67～80 km。增加落角約束的影響如圖8所示。

從圖8結(jié)果可知，在兩種不同的任務(wù)需求下，彈丸飛行時間下界隨射程變化的趨勢相近，增加大落角約束會增加末端彈道的曲率，使飛行時間延長，且隨著射程增大，影響程度增加;而飛行時間上界的變化趨勢出現(xiàn)了明顯差異。近區(qū)之內(nèi)，大落角約束對飛行時間上界并無影響，這是因為該區(qū)域內(nèi)彈丸延長飛行時間與增加存速所采用的方式均為增加彈道曲率，大落角是必然的結(jié)果，無需額外約束;中區(qū)之內(nèi)，飛行方案同時兼顧存速與控制能量，增加彈道曲率不再是最佳的方式，彈丸會犧牲一定存速裕量，延長飛行時間以滿足落角約束;而遠(yuǎn)區(qū)之內(nèi)，彈丸有限的舵控能力主要服務(wù)于增程，此時增加末端彈道的曲率反而會減小飛行時間上界。上下界之間的差值表示了當(dāng)前射程彈丸的協(xié)調(diào)能力，并在協(xié)作范圍的變化上得到體現(xiàn)。從整體上看，增加大落角約束會使彈丸的協(xié)作范圍有所降低。

5 結(jié) 論

本文針對多炮齊射場景下滑翔制導(dǎo)炮彈的多階段全彈道協(xié)同規(guī)劃問題進(jìn)行研究。為充分挖掘彈群的控制潛力，以整體為對象建立了協(xié)同彈道規(guī)劃模型，并在此基礎(chǔ)上提出一種雙層自適應(yīng)協(xié)同策略。在該策略下，彈群以最優(yōu)末端時刻進(jìn)行打擊的同時，保留了自身各階段的獨立性。此外，根據(jù)不同射程范圍內(nèi)協(xié)同彈道的側(cè)重與特性差異，為制導(dǎo)炮彈進(jìn)行了戰(zhàn)區(qū)劃分，提出一種通用的確定戰(zhàn)區(qū)與協(xié)作范圍的方法，并初步研究了不同作戰(zhàn)任務(wù)需求引入的額外約束對范圍的影響。仿真結(jié)果表明：

（1） 本文所提出的方法可以簡捷有效地獲取滑翔制導(dǎo)炮彈的戰(zhàn)區(qū)分布與任意炮口位置的齊射協(xié)作范圍，從而輔助戰(zhàn)斗隊形的配置;不同戰(zhàn)區(qū)內(nèi)的飛行方案側(cè)重不同，齊射彈丸協(xié)調(diào)飛行時間的主要方式與彈道特性也不同。

（2） 在所提出的雙層自適應(yīng)協(xié)同策略下，齊射彈丸可以滿足全部約束條件，成功實現(xiàn)協(xié)同打擊。在本文算例中，近、中、遠(yuǎn)3區(qū)內(nèi)彈群的整體規(guī)劃目標(biāo)函數(shù)值相較于現(xiàn)有的分布式策略分別減小1.97%、84.72%、25.30%，具有一定的優(yōu)越性。同時，本文算例也證明以平均值方式協(xié)調(diào)飛行時間的分布式策略易造成規(guī)劃結(jié)果不理想甚至失敗的局限所在。

（3） 不同的作戰(zhàn)任務(wù)引入的額外約束條件主要影響彈丸的最大飛行時間，進(jìn)而影響戰(zhàn)區(qū)分布與協(xié)作范圍。灌頂打擊任務(wù)下，滑翔制導(dǎo)炮彈的整體協(xié)作范圍較常規(guī)任務(wù)有所縮減。本文提供一種通用的處理思路，不同場景下的影響程度應(yīng)視具體情況而定。

參考文獻(xiàn)

［1］ 徐秋坪. 滑翔制導(dǎo)炮彈彈道規(guī)劃及其自抗擾控制系統(tǒng)研究［D］. 南京： 南京理工大學(xué)， 2018.

XU Q P. Trajectory optimization and design of active-disturbance-rejection control system for gliding guided projectiles［D］. Nanjing： Nanjing University of Science and Technology， 2018.

［2］ 魯濟(jì)帥. 多彈協(xié)同攻擊策略關(guān)鍵技術(shù)研究［D］. 西安： 西安工業(yè)大學(xué)， 2023.

LU J S. Research on key technologies of multi-bomb coordinated attack strategy［D］. Xi’an： Xi’an Technological University， 2023.

［3］ 王嬌嬌. 智能彈藥多彈協(xié)同航跡規(guī)劃研究［D］. 西安： 西安工業(yè)大學(xué)， 2022.

WANG J J. Research on multi missile cooperative track planning of intelligent ammunition［D］. Xi’an： Xi’an Technological University， 2022.

［4］ 姚禹正， 余文斌， 楊立軍， 等. 多導(dǎo)彈協(xié)同制導(dǎo)技術(shù)綜述［J］. 飛航導(dǎo)彈， 2021（6）： 112-121.

YAO Y Z， YU W B， YANG L J， et al. A review of multi-missile cooperative guidance technology［J］. Aerodynamic Missile Journal， 2021（6）： 112-121.

［5］ 趙建博， 楊樹興. 多導(dǎo)彈協(xié)同制導(dǎo)研究綜述［J］. 航空學(xué)報， 2017， 38（1）： 22-34.

ZHAO J B， YANG S X. Review of multi-missile cooperative guidance［J］. Acta Aeronautica et Astronautica Sinica， 2017， 38（1）： 22-34.

［6］ 崔正達(dá)， 魏明英， 李運(yùn)遷. 基于速度預(yù)測的防空導(dǎo)彈中制導(dǎo)末段協(xié)同彈道規(guī)劃方法［J］. 系統(tǒng)工程與電子技術(shù)， 2023， 45（9）： 2912-2921.

CUI Z D， WEI M Y， LI Y Q. Cooperate trajectory planning method in later part of midcourse based on velocity estimation［J］. Systems Engineering and Electronics， 2023， 45（9）： 2912-2921.

［7］ 陳子豪， 李娟， 劉暢， 等. 時間約束下對地面目標(biāo)協(xié)同攻擊任務(wù)規(guī)劃方法［J］. 系統(tǒng)工程與電子技術(shù)， 2023， 45（8）： 2353-2360.

CHEN Z H， LI J， LIU C， et al. Task planning method for coordinated attacks on ground targets under time constraints［J］. Systems Engineering and Electronics， 2023， 45（8）： 2353-2360.

［8］ 尹秋霖， 陳琦， 王中原， 等. 考慮時空協(xié)同的滑翔制導(dǎo)炮彈單炮多發(fā)快速彈道規(guī)劃［J］. 兵工學(xué)報， 2023， 45（3）： 798-809.

YIN Q L， CHEN Q， WANG Z Y， et al. Rapid trajectory planning for glide-guided projectiles in single-gun multishot scenarios considering time and space cooperation［J］. Acta Armamentarii， 2023， 45（3）： 798-809.

［9］ CHEN Q， WANG Z Y. Optimal trajectory for time-on-target of a guided projectile using direct collocation method［C］∥Proc.of the International Conference on Mechatronic Sciences， Electric Engineering and Computer， 2013： 2803-2806.

［10］ 許強(qiáng)強(qiáng)， 葛健全， 楊濤， 等. 面向突防的多導(dǎo)彈協(xié)同彈道規(guī)劃方法［J］. 中國慣性技術(shù)學(xué)報， 2018， 26（4）： 524-530.

XU Q Q， GE J Q， YANG T， et al. Cooperative trajectory planning for penetration of multiple missiles［J］. Journal of Chinese Inertial Technology， 2018， 26（4）： 524-530.

［11］ JIANG H， AN Z， YU Y N， et al. Cooperative guidance with multiple constraints using convex optimization［J］. Aerospace Science and Technology， 2018， 79： 426-440.

［12］ LIU S Y， LI S， LIU Y J. Research on anti-ship missile cooperative trajectory planning based on gauss pseudospectral method［C］∥Proc.of the 5th Chinese Conference on Swarm Intelligence and Cooperative Control， 2022： 1163-1176.

［13］ 錢杏芳， 林瑞雄， 趙亞男. 導(dǎo)彈飛行力學(xué)［M］. 北京： 北京理工大學(xué)出版社， 2006.

QIAN X F， LIN R X， ZHAO Y N. Missile flight mechanics［M］. Beijing： Press of Beijing Institute of Technology， 2006.

［14］ 陳琦. 滑翔增程制導(dǎo)炮彈彈道優(yōu)化及制導(dǎo)控制方法研究［D］. 南京： 南京理工大學(xué)， 2017.

CHEN Q. Investigation on trajectory optimization and guidance amp; control scheme for a type of gliding guided projectiles［D］. Nanjing： Nanjing University of Science and Technology， 2017.

［15］ 韓子鵬. 彈箭外彈道學(xué)［M］. 北京： 北京理工大學(xué)出版社， 2008.

HAN Z P. External ballistics of projectiles［Z］. Beijing： Press of Beijing Institute of Technology， 2008.

［16］ ZHANG J， SHENG H L， CHEN Q， et al. A four-dimensional space-time automatic obstacle avoidance trajectory planning method for multi-UAV cooperative formation flight［J］. Drones， 2022， 6（8）： 192.

［17］ YI K， SHE S H， ZHANG S， et al. Research on cooperative strategy of multiple rockets［J］. MATEC Web of Conferences， 2020， 316： 4009.

［18］ DONG W， WEN Q Q， XIA Q L， et al. Multiple-constraint cooperative guidance based on two-stage sequential convex programming［J］. Chinese Journal of Aeronautics， 2020， 33（1）： 296-307.

［19］ ZHOU B， HUANG B， SU Y M， et al. Two-layer leader-follower optimal affine formation maneuver control for networked unmanned surface vessels with input saturations［J］. International Journal of Robust and Nonlinear Control， 2024， 34（5）： 3631-3655.

［20］ LIU J Z， FAN C Y， PENG Y C， et al. Emergent leader-follower relationship in networked multiagent systems［J］. SCIENCE CHINA Information Sciences， 2023， 66（12）： 229201.

［21］ WU X Q， MAO B， WU X Q， et al. Dynamic event-triggered leader-follower consensus control for multiagent systems［J］. SIAM Journal on Control and Optimization， 2022， 60（1）： 189-209.

［22］ 孫世巖， 陳俊丞， 梁偉閣， 等. 制導(dǎo)炮彈單炮多發(fā)同時彈著實現(xiàn)方法研究［J］. 海軍工程大學(xué)學(xué)報， 2021， 33（1）： 84-88.

SUN S Y， CHEN J C， LIANG W G， et al. Research on implementation method of multiple round simultaneous impact of guided projectiles［J］. Journal of Naval University of Engineering， 2021， 33（1）： 84-88.

［23］ ZHOU T， ZHENG F F， LIU X B， et al. A two-stage trajectory planning model for cooperative truck platooning on freeways［J］. Journal of Intelligent Transportation Systems， 2023， 27（2）： 217-237.

［24］ XU G T， LONG T， WANG Z， et al. Trust-region filtered sequential convex programming for multi-UAV trajectory planning and collision avoidance［J］. ISA Transactions， 2022， 128： 664-676.

［25］ TAJALLI M， NIROUMAND R， HAJBABAIE A. Distributed cooperative trajectory and lane changing optimization of connected automated vehicles： freeway segments with lane drop［J］. Transportation research part C： Emerging Technologies， 2022， 143： 103761.

［26］ ZHANG X C， ZHAO H T， WEI J B， et al. Cooperative trajectory design of multiple UAV base stations with heterogeneous graph neural networks［J］. IEEE Trans.on Wireless Communications， 2022， 22（3）： 1495-1509.

［27］ 劉超越， 張成. 基于分布式并行偽譜-神經(jīng)網(wǎng)絡(luò)脈沖導(dǎo)彈多階段協(xié)同軌跡優(yōu)化［J］. 兵工學(xué)報， 2020， 41（10）： 1988-2000.

LIU C Y， ZHANG C. Multi-stage cooperative trajectory optimization of dual-pulse missile based on decentralized parallel pseudospectral-neural network algorithm［J］. Acta Armamentarii， 2020， 41（10）： 1988-2000.

［28］ LIU C Y， ZHANG C， XIONG F F. Multistage cooperative trajectory planning for multimissile formation via bi-level sequential convex programming［J］. IEEE Access， 2020， 8： 22834-22853.

［29］ XU L， CAO X B， DU W B， et al. Cooperative path planning optimization for multiple UAVs with communication constraints［J］. Knowledge-Based Systems， 2023， 260： 110164.

［30］ OUYANG Y K， LI B， ZHANG Y M， et al. Fast and optimal trajectory planning for multiple vehicles in a nonconvex and cluttered environment： benchmarks， methodology， and experiments［C］∥Proc.of the IEEE International Conference on Robotics and Automation， 2022： 10746-10752.

［31］ ZONG F， HE Z B， ZENG M， et al. Dynamic lane changing trajectory planning for CAV： a multi-agent model with path preplanning［J］. Transportmetrica B： Transport Dynamics， 2022， 10（1）： 266-292.

［32］ 李征， 彭博， 陳海東， 等. 可重復(fù)使用航天器時間協(xié)同飛行軌跡優(yōu)化［J］. 計算機(jī)仿真， 2020， 37（1）： 40-45.

LI Z， PENG B， CHEN H D， et al. Time-coordination reentry trajectory design for reusable launch vehicle［J］. Computer Simulation， 2020， 37（1）： 40-45.

［33］ PATTERSON M A， RAO A V. GPOPS-II： a MATLAB software for solving multiple-phase optimal control problems using HP-adaptive Gaussian quadrature collocation methods and sparse nonlinear programming［J］. ACM Transactions on Mathematical Software， 2014， 41（1）： 1-37.

作者簡介

尹秋霖（1993—），男，博士研究生，主要研究方向為滑翔制導(dǎo)炮彈彈道規(guī)劃與導(dǎo)引控制。

陳 琦（1989—），男，副研究員，博士，主要研究方向為彈箭外彈道設(shè)計及優(yōu)化、飛行動力學(xué)建模與仿真、制導(dǎo)與控制。

王中原（1958—），男，研究員，博士，主要研究方向為彈箭飛行控制理論與技術(shù)、外彈道仿真與設(shè)計。

王慶海（1994—），男，博士研究生，主要研究方向為制導(dǎo)炮彈彈道規(guī)劃、制導(dǎo)控制一體化。