郭正玉　韓治國(guó)

摘 要：針對(duì)空空導(dǎo)彈協(xié)同攻擊高速大機(jī)動(dòng)目標(biāo)問(wèn)題，提出了一種帶有期望攻擊角約束的多彈協(xié)同制導(dǎo)律。首先，在縱向平面內(nèi)建立彈-目幾何關(guān)系，建立含有攻擊角約束的導(dǎo)彈視線方向和視線法向方向的多彈協(xié)同制導(dǎo)模型;其次，針對(duì)視線方向制導(dǎo)模型，基于代數(shù)圖論和有限時(shí)間一致性理論設(shè)計(jì)了視線方向上的協(xié)同制導(dǎo)律，保證3枚導(dǎo)彈能夠?qū)崿F(xiàn)對(duì)目標(biāo)的協(xié)同攻擊，并利用現(xiàn)有觀測(cè)器對(duì)目標(biāo)的機(jī)動(dòng)能力進(jìn)行估計(jì);再次，基于快速非奇異終端滑?？刂评碚撛O(shè)計(jì)了視線法向方向上的制導(dǎo)律，保證3枚導(dǎo)彈均能夠精確命中目標(biāo)，同時(shí)保證彈-目的視線角收斂到期望的終端視線角，視線角速率收斂到0;最后，仿真驗(yàn)證了所設(shè)計(jì)制導(dǎo)律的有效性。

關(guān)鍵詞： 協(xié)同制導(dǎo);攻擊角約束;快速非奇異終端滑模;有限時(shí)間

中圖分類號(hào)：TJ765.3 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)： 1673-5048（2020）03-0062-05

0 引言

空空導(dǎo)彈經(jīng)過(guò)70多年的發(fā)展，由最初的無(wú)制導(dǎo)火箭彈發(fā)展到現(xiàn)在的制導(dǎo)方式多樣化，遠(yuǎn)、中、近距系列化和海、陸、空三軍通用化的空空導(dǎo)彈家族，成為空中對(duì)抗的主要武器和各軍事強(qiáng)國(guó)優(yōu)先發(fā)展的武器裝備[1-2]。

隨著作戰(zhàn)環(huán)境的日趨復(fù)雜，作戰(zhàn)目標(biāo)的不斷出現(xiàn)，影響了空空導(dǎo)彈制空作戰(zhàn)的效能。采用多枚空空導(dǎo)彈協(xié)同作戰(zhàn)，使導(dǎo)彈之間能夠直接進(jìn)行信息交互，有效對(duì)抗復(fù)雜環(huán)境，提高對(duì)目標(biāo)的探測(cè)和毀傷，特別是針對(duì)大機(jī)動(dòng)目標(biāo)，能夠從多個(gè)角度、多個(gè)方位協(xié)同攻擊，從而極大地提高作戰(zhàn)效能。因此，研究空空導(dǎo)彈的多彈協(xié)同制導(dǎo)方法有著重要的理論和現(xiàn)實(shí)意義。

協(xié)同制導(dǎo)的一個(gè)典型應(yīng)用是導(dǎo)彈同時(shí)發(fā)射對(duì)目標(biāo)進(jìn)行攻擊。一般地，導(dǎo)彈同時(shí)發(fā)射的制導(dǎo)律設(shè)計(jì)，可以根據(jù)攻擊目標(biāo)的不同分為兩類[3]，發(fā)射前預(yù)設(shè)攻擊時(shí)間的制導(dǎo)律和在飛行時(shí)的實(shí)時(shí)攻擊制導(dǎo)律。第一類制導(dǎo)律的開(kāi)創(chuàng)者Jeon等人，通過(guò)在比例導(dǎo)引的基礎(chǔ)上疊加一個(gè)對(duì)攻擊時(shí)間的誤差反饋，使得導(dǎo)彈的實(shí)際飛行時(shí)間接近預(yù)計(jì)的攻擊時(shí)間。進(jìn)一步地，Lee等人[4]以及Harl等人[5]分別通過(guò)最優(yōu)控制和滑?？刂评碚摪l(fā)展了該類制導(dǎo)律。然而，第一類制導(dǎo)律有一個(gè)明顯的缺點(diǎn)，就是很難對(duì)機(jī)動(dòng)目標(biāo)，尤其是對(duì)不確定其運(yùn)動(dòng)規(guī)律的機(jī)動(dòng)目標(biāo)設(shè)計(jì)預(yù)計(jì)的攻擊時(shí)間，因此第二類制導(dǎo)律應(yīng)運(yùn)而生。為了使攻擊時(shí)間動(dòng)態(tài)同步，研究者們提出了領(lǐng)-從式協(xié)同制導(dǎo)策略和集中式協(xié)同制導(dǎo)策略[6-7]。此外，為提高制導(dǎo)系統(tǒng)的適應(yīng)性和可擴(kuò)展性，Zhao等人[8]提出了一種通過(guò)控制攻擊時(shí)間的分布式協(xié)同制導(dǎo)算法。李強(qiáng)[9]針對(duì)協(xié)同制導(dǎo)問(wèn)題，在視線方向，基于一致性理論，設(shè)計(jì)了控制指令，實(shí)現(xiàn)了彈間協(xié)同制導(dǎo);在視線法向，設(shè)計(jì)了有限時(shí)

間收斂的滑模制導(dǎo)律。Zhai等人[10]通過(guò)考慮對(duì)帶有誘餌的目標(biāo)，運(yùn)用多個(gè)攔截導(dǎo)彈，提出了一種基于區(qū)域覆蓋的攔截算法，最大程度地提高了聯(lián)合攔截的成功率。Shaferman等人 [11]提出了一種領(lǐng)彈和從彈之間的協(xié)同策略，通過(guò)導(dǎo)彈和目標(biāo)的信息共享提高了從彈的攔截性能。 為減少傳遞誤差，Wang等人[12]提出了一種協(xié)同中段制導(dǎo)律，使得中段的導(dǎo)彈可以運(yùn)用末端導(dǎo)彈所收集的目標(biāo)信息。

本文在上述研究成果的基礎(chǔ)上，針對(duì)多導(dǎo)彈協(xié)同攻擊空中的高速大機(jī)動(dòng)目標(biāo)問(wèn)題，提出了一種帶有期望攻擊角約束的多彈協(xié)同制導(dǎo)律，建立了多彈協(xié)同制導(dǎo)模型， 并利用數(shù)值仿真驗(yàn)證了本文設(shè)計(jì)模型的有效性。

1 協(xié)同制導(dǎo)模型建立

為了簡(jiǎn)化協(xié)同制導(dǎo)模型，將導(dǎo)彈和目標(biāo)視為質(zhì)點(diǎn)?；谠摷僭O(shè)，本文針對(duì)3枚導(dǎo)彈協(xié)同攻擊高速大機(jī)動(dòng)目標(biāo)的問(wèn)題開(kāi)展研究，彈-目幾何關(guān)系如圖1所示。

圖中，Mi和T為第i枚導(dǎo)彈和目標(biāo)，ri為Mi與T的相對(duì)距離，qi為Mi的視線角，vmi和vt為Mi和T的速度，θmi和θt為Mi和T的彈道傾角，ami和at為Mi和T的法向過(guò)載，其中，i=1， 2， 3， i為導(dǎo)彈數(shù)量。協(xié)同制導(dǎo)方程如下[12-13]：

式中： r¨i為彈-目距離對(duì)時(shí)間的2階導(dǎo)數(shù);q¨i為彈-目視線角對(duì)時(shí)間的2階導(dǎo)數(shù);uri和wri分別為導(dǎo)彈和目標(biāo)的加速度在視線方向的分量;uqi和wqi分別為導(dǎo)彈和目標(biāo)的加速度在視線法向方向的分量;uri，wri，uqi，wqi變量的具體表達(dá)式可參見(jiàn)相關(guān)文獻(xiàn)。令x1i=ri，x2i=r·i，x3i=qi-qdi=qei，x4i=q·i，其中，qdi為導(dǎo)彈的期望視線角;qei為Mi的彈目實(shí)際視線角與期望視線角偏差，則由式（1）可得多導(dǎo)彈協(xié)同制導(dǎo)模型為[12-13]

根據(jù)式（2），在導(dǎo)彈視線方向和視線法向分別設(shè)計(jì)uri和uqi，使得x1i及x2i在有限時(shí)間內(nèi)達(dá)到一致，x3i及x4i在有限時(shí)間內(nèi)收斂到0， 最后，將得到的uri和uqi代入式（2）中，并通過(guò)仿真驗(yàn)證所設(shè)計(jì)制導(dǎo)律的有效性。

2 協(xié)同制導(dǎo)律設(shè)計(jì)

2.1 視線方向上的加速度指令設(shè)計(jì)

基于二階多智能體系統(tǒng)有限時(shí)間一致性原理，設(shè)計(jì)視線方向上的加速度指令，以保證所有的導(dǎo)彈同時(shí)擊中機(jī)動(dòng)目標(biāo)[12]。針對(duì)式（2），通過(guò)變量代換轉(zhuǎn)化為[13]

式中： xi=x1i，vi=x2i，ui=x1ix24i-uri。

引理1： 考慮式（2）所表示的系統(tǒng)，己知其通信拓?fù)浣Y(jié)構(gòu)圖G是無(wú)向連通圖，根據(jù)無(wú)向連通圖原理，采用如下的控制律可以保證系統(tǒng)狀態(tài)在有限時(shí)間內(nèi)趨于一致[11]：

式中： 0<αi<1（i=1， 2， 3）， 2（α1-α3）=（1+α1）α2;sig（x）α1=|x|α1sgn（x）;ψ1和ψ2是奇函數(shù);存在正數(shù)bi（i=1，2），yψi（y）>0（y≠0∈R）和ψi（y）=biy+o（y）（y≠0∈R）（在0的域內(nèi)）成立。

由引理1可知，設(shè)計(jì)式（5）的制導(dǎo)律可以使得式（2）的狀態(tài)在有限時(shí)間內(nèi)趨于一致，即實(shí)現(xiàn)視線方向的協(xié)同制導(dǎo)：

式中： wri為外部擾動(dòng)的總和，可由下式所構(gòu)成的觀測(cè)器進(jìn)行估計(jì)[12-14]：

式中： λji，μji>0（j=1，2，3，i=1，2，3）; w^ri（t）是系統(tǒng)外界總干擾wri的估計(jì)值。

2.2 視線法線方向上的加速度指令設(shè)計(jì)

為方便制導(dǎo)律設(shè)計(jì)，給出如下引理：

引理2： 對(duì)于任意給定的實(shí)數(shù)λ1>0，λ2>0，0

V·（x）+λ1V（x）+λ2Vk（x）≤0（7）

對(duì)于式（7），穩(wěn)定時(shí)間能夠用下式進(jìn)行估計(jì)：

Tk≤1λ1（1-k）lnλ1V1-k（x0）+λ2λ2（8）

即系統(tǒng)狀態(tài)有限時(shí)間收斂。 針對(duì)制導(dǎo)系統(tǒng)模型式（2），設(shè)計(jì)非奇異終端滑模面：

s=x3+βsigα（x4）+λx4（9）

式中： β，λ>0; 1<α<2;sigα（x4）=|x4|αsgn（x4）。

本文針對(duì)高速大機(jī)動(dòng)目標(biāo)， 在設(shè)計(jì)趨近律時(shí)，考慮如下含有符號(hào)函數(shù)的快速冪次趨近律：

s·=-k1s-k2sigγ（s）-εsgn（s） （10）

式中： k1，k2>0;0<γ<1;ε需要滿足的條件將在后文給出。

對(duì)式（9）求導(dǎo)，可得

s·=x3+（αβ|x4|α-1+λ）-2x2x1x4-uqx1+wqx1（11）

結(jié)合式（10）～（11），得

式中： b=αβ|x4|α-1+λ。

定理1： 對(duì)于協(xié)同制導(dǎo)模型式（2），在設(shè)計(jì)的制導(dǎo)律式（12）的作用下，可使系統(tǒng)狀態(tài)x3，x4在有限時(shí)間內(nèi)到達(dá)滑模面。在滑模面s=0上，根據(jù)快速非奇異終端滑模的基本特性，系統(tǒng)狀態(tài)x3，x4能夠在有限時(shí)間內(nèi)收斂到0，即q·（t）→0，q（t）→qd。

證明： 構(gòu)造如下的Lyapunov函數(shù)：

根據(jù)引理1可知，系統(tǒng)狀態(tài)將在有限時(shí)間內(nèi)收斂到本文設(shè)計(jì)的滑模面，且收斂時(shí)間為

在滑模面s=0上，根據(jù)滑模面的性質(zhì)，系統(tǒng)狀態(tài)x3，x4將在有限時(shí)間收斂到平衡點(diǎn)，即視線角收斂到期望值，視線角速率收斂到0，證畢。

本文在設(shè)計(jì)滑模面式（9）時(shí)，引入了λx4這一項(xiàng)。這一項(xiàng)的作用是避免系統(tǒng)狀態(tài)x4趨于0時(shí)，系統(tǒng)出現(xiàn)奇異現(xiàn)象。如果λ=0，則x4→0，b=αβ·x4α-1→0，進(jìn)而有b-1→∞，從而系統(tǒng)出現(xiàn)奇異現(xiàn)象。在設(shè)計(jì)滑模面時(shí)，通過(guò)引入λx4這一項(xiàng)，使得x4→0，b=αβx4α-1+λ→λ，b-1→λ-1，從而有效避免了奇異現(xiàn)象。

3 仿真分析

為了驗(yàn)證本文所采用的視線方向制導(dǎo)律式（5）與設(shè)計(jì)的基于快速非奇異終端滑模視線法線方向制導(dǎo)律式（12）在協(xié)同制導(dǎo)時(shí)的有效性，設(shè)計(jì)仿真場(chǎng)景： 3枚導(dǎo)彈在有限的攻擊距離內(nèi)，協(xié)同攻擊1個(gè)高速大機(jī)動(dòng)目標(biāo)。目標(biāo)的初始位置為（30 km， 29 km），速度為500 m/s，初始航向角為160°，目標(biāo)機(jī)動(dòng)過(guò)載為10gsin（t）。3枚導(dǎo)彈的初始參數(shù)如表1所示，導(dǎo)彈最大可用過(guò)載30g。

多彈協(xié)同既要保證導(dǎo)彈之間實(shí)現(xiàn)信息互通、共享其自組織網(wǎng)絡(luò)能力，還要建立多層面、異構(gòu)與同構(gòu)并存的混合體系架構(gòu)。多彈能夠根據(jù)態(tài)勢(shì)，采用動(dòng)態(tài)分配的方式實(shí)現(xiàn)內(nèi)部網(wǎng)絡(luò)成員之間的協(xié)同攻擊。本文采用有中心的組織架構(gòu)設(shè)計(jì)，彈間通訊結(jié)構(gòu)如圖2所示，M2即編號(hào)為2的導(dǎo)彈， 是3枚導(dǎo)彈的彈群中心，M2與M1和M3均能夠傳遞信息。M1和M3即編號(hào)為1和3的導(dǎo)彈，僅能夠與M2傳遞信息，M1和M3之間不進(jìn)行信息傳遞，整個(gè)彈群由指揮中心對(duì)任務(wù)進(jìn)行分配，彈間通訊權(quán)系數(shù)矩陣可表示為

導(dǎo)彈視線方向的制導(dǎo)律參數(shù)選取如下： ψ1（x）=ψ2（x）=x， α1=0.6， α2=0.75。導(dǎo)彈視線法向方向的制導(dǎo)律參數(shù)選取如下： α=5/3， β=5， λ=1.2， k1=3，k2=2。

根據(jù)上述仿真參數(shù)和本文所設(shè)計(jì)的協(xié)同制導(dǎo)律，典型目標(biāo)和攻擊態(tài)勢(shì)的仿真結(jié)果如圖3～7所示。3枚導(dǎo)彈的制導(dǎo)時(shí)間和脫靶量情況見(jiàn)表2。

圖3為慣性系下的彈-目運(yùn)動(dòng)軌跡曲線，從圖中可以看出，采用本文設(shè)計(jì)的協(xié)同制導(dǎo)律，3枚導(dǎo)彈能夠同時(shí)命中目標(biāo)，彈道無(wú)交叉，導(dǎo)彈不會(huì)相互碰撞。圖4為彈-目相對(duì)距離曲線，同樣可以看出，在3枚導(dǎo)彈初始彈-目距離不相同的情形下，

采用本文設(shè)計(jì)的制導(dǎo)律，在大約15 s左右，3枚導(dǎo)彈實(shí)現(xiàn)了彈-目距離協(xié)同，3枚導(dǎo)彈能夠同時(shí)命中目標(biāo)。圖5為導(dǎo)彈的速度曲線，可以看出，為了實(shí)現(xiàn)對(duì)高速大機(jī)動(dòng)目標(biāo)的協(xié)同攻擊，3枚導(dǎo)彈的速度需要按照視線方向過(guò)載指令進(jìn)行變化。圖6為導(dǎo)彈視線方向過(guò)載曲線，可以看出，導(dǎo)彈在視線方向的最大過(guò)載為10g， 在攻擊高速大機(jī)動(dòng)目標(biāo)時(shí)，導(dǎo)彈在視線方向的過(guò)載沒(méi)有發(fā)生飽和，且視線方向過(guò)載曲線變化平穩(wěn)。圖7

為導(dǎo)彈視線法向方向過(guò)載曲線，在視線方向的法向，導(dǎo)

彈最大過(guò)載為30g，在初始階段，導(dǎo)彈為了盡快調(diào)整彈-目相對(duì)位置關(guān)系，視線法向方向過(guò)載發(fā)生了大約5 s的過(guò)載飽和現(xiàn)象，5 s之后視線法向方向過(guò)載變化非常平穩(wěn)，在制導(dǎo)末段接近為0，保證能夠命中目標(biāo)。從表2可以看出，3枚導(dǎo)彈的脫靶量均很小。因此，根據(jù)上述仿真結(jié)

果可以看出，針對(duì)高速大機(jī)動(dòng)目標(biāo)，利用本文設(shè)計(jì)的協(xié)同制導(dǎo)律，能夠?qū)崿F(xiàn)對(duì)目標(biāo)的精確協(xié)同攻擊。

4 結(jié)論

本文針對(duì)多導(dǎo)彈協(xié)同攻擊高速大機(jī)動(dòng)目標(biāo)問(wèn)題進(jìn)行了研究，設(shè)計(jì)了能夠支持多彈協(xié)同攻擊的協(xié)同制導(dǎo)律，并進(jìn)行了理論建模與仿真分析。仿真結(jié)果表明，3枚導(dǎo)彈不但能夠命中目標(biāo)，還能夠協(xié)同并同時(shí)命中高速大機(jī)動(dòng)目標(biāo)。因此，本文所設(shè)計(jì)的協(xié)同制導(dǎo)律具有先進(jìn)性，能夠提高多導(dǎo)彈制導(dǎo)控制的精確性，提升多導(dǎo)彈協(xié)同作戰(zhàn)的效能。

參考文獻(xiàn)：

[1] 樊會(huì)濤， 閆俊. 自主化——機(jī)載導(dǎo)彈重要的發(fā)展方向[J]. 航空兵器， 2019， 26（1）： 1-10.

Fan Huitao， Yan Jun. The Important Development Direction of Airborne Missile： Autonomization[J]. Aero Weaponry， 2019， 26（1）： 1-10. （in Chinese）

[2] 郭正玉，劉琪. 從空海聯(lián)合作戰(zhàn)看美國(guó)空空導(dǎo)彈發(fā)展[J]. 航空兵器，2018（6）： 11-15.

Guo Zhengyu， Liu Qi. Evolution of Air-Sea Battle and the Development of US Air-to-Air Missile[J]. Aero Weaponry， 2018（6）： 11-15. （in Chinese）

[3] Jeon I S， Lee J I， Tahk M J. Impact-Time-Control Guidance Law for Anti-Ship Missiles[J]. IEEE Transactions on Control Systems Technology， 2006， 14（2）： 260-266.

[4] Lee J I， Jeon I S， Tahk M J. Guidance Law to Control Impact Time and Angle[C]∥ International Conference on Control and Automation， Budapest， Hungary， 2005.

[5] Harl N， Balakrishnan S N. Impact Time and Angle Guidance with Sliding Mode Control[J]. IEEE Transactions on Control Systems Technology， 2012， 20（6）：1436-1449.

[6] Zhao S Y， Zhou R. Cooperative Guidance for Multi Missile Salvo Attack[J]. Chinese Journal of Aeronautics， 2008， 21（6）： 533-539.

[7] Jeon I S ， Lee J I ， Tahk M J . Homing Guidance Law for Cooperative Attack of Multiple Missiles[J]. Journal of Guidance， Control， and Dynamics， 2010， 33（1）： 275-280.

[8] Zhao S Y， Zhou R， Wei C， et al. Design of Time-Constrained Guidance Laws via Virtual Leader Approach[J]. Chinese Journal of Aeronautics， 2010， 23（1）： 103-108.

[9] 李強(qiáng). 多導(dǎo)彈協(xié)同作戰(zhàn)目標(biāo)分配和制導(dǎo)律研究[D]. 哈爾濱： 哈爾濱工業(yè)大學(xué)， 2017.

Li Qiang. Research on Target Assignment and Guidance Law for Cooperative Engagement of Multi-Missile[D]. Harbin： Harbin Institute of Technology， 2017. （in Chinese）

[10] Zhai C， He F H， Hong Y G， et al. Coverage-Based Interception Algorithm of Multiple Interceptors Against the Target Involving Decoys[J]. Journal of Guidance， Control， and Dynamics， 2016， 39（7）：1-7.

[11] Shaferman V， Oshman Y. Stochastic Cooperative Interception Using Information Sharing Based on Engagement Staggering[J]. Journal of Guidance， Control， and Dynamics， 2016， 39（9）： 2127-2141.

[12] Wang L， Yao Y， He F H， et al. A Novel Cooperative Mid-Course Guidance Scheme for Multiple Intercepting Missiles[J]. Chinese Journal of Aeronautics， 2017， 30（3）： 1140-1153.

[13] 宋俊紅， 宋申民， 徐勝利. 帶有攻擊角約束的多導(dǎo)彈協(xié)同制導(dǎo)律[J]. 中國(guó)慣性技術(shù)學(xué)報(bào)， 2016， 24（4）： 554-560.