成忠濤 李聚峰 彭高祥 劉磊 王永驥

摘 要：研究了一種基于Bézier曲線的碰撞角約束的制導(dǎo)律。通過(guò)模型轉(zhuǎn)化將制導(dǎo)指令設(shè)計(jì)的問(wèn)題轉(zhuǎn)化為二次Bézier曲線形式的航跡角設(shè)計(jì)問(wèn)題。首先，利用Bézier曲線的性質(zhì)設(shè)計(jì)了導(dǎo)彈速度大小不變的制導(dǎo)律;然后進(jìn)一步對(duì)導(dǎo)彈速度時(shí)變的末端角度約束的制導(dǎo)律進(jìn)行了研究，同時(shí)對(duì)導(dǎo)彈速度時(shí)變情形下的剩余飛行時(shí)間進(jìn)行了估計(jì);最后通過(guò)不同情況下的數(shù)值仿真，驗(yàn)證了所提出的制導(dǎo)律的有效性。

關(guān)鍵詞：Bézier曲線;航跡角;碰撞角;制導(dǎo)律;剩余飛行時(shí)間

中圖分類號(hào)：TJ765.3 ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

飛行器制導(dǎo)律的設(shè)計(jì)一直是近幾十年來(lái)研究的熱點(diǎn)話題。隨著導(dǎo)彈任務(wù)的多樣化和作戰(zhàn)場(chǎng)景的復(fù)雜化，對(duì)制導(dǎo)律提出了新的要求，如視場(chǎng)約束、碰撞時(shí)間控制和碰撞角控制[1-3]。其中，特定的碰撞角具有重要意義，第一，碰撞角控制制導(dǎo)律（IACG）可以引導(dǎo)導(dǎo)彈到達(dá)目標(biāo)的最薄弱點(diǎn);第二，特定任務(wù)的特定碰撞角度可以使殺傷能力最大化;最后，IACG可以引導(dǎo)導(dǎo)彈繞過(guò)防御系統(tǒng)。

IACG的最早研究可以追溯到1973年，Kim將其應(yīng)用于再入飛行器[4]。他通過(guò)應(yīng)用最優(yōu)控制理論解決了碰撞角控制問(wèn)題。Ryoo等針對(duì)不同的導(dǎo)彈動(dòng)力學(xué)特性，提出了一種狀態(tài)反饋形式的廣義最優(yōu)制導(dǎo)律[5]。進(jìn)一步在一種新的飛行時(shí)間估計(jì)方法的基礎(chǔ)上，提出了一種新的等速導(dǎo)彈最優(yōu)碰撞角控制制導(dǎo)律[6]。文獻(xiàn)[7]提出了一種新的線性最優(yōu)IACG，與傳統(tǒng)的框架不同，線性化不是在初始視線角附近進(jìn)行的，而是在標(biāo)稱圓軌跡附近進(jìn)行的。

Lyapunov穩(wěn)定性和滑?？刂频确蔷€性控制理論，也被應(yīng)用到IACG的設(shè)計(jì)中。文獻(xiàn)[8]在Lyapunov候選函數(shù)中增加了一個(gè)碰撞角誤差項(xiàng)，使之以特定的角度擊中目標(biāo)。文獻(xiàn)[9]提出了一種基于減小航向誤差角的Lyapunov候選函數(shù)的新型制導(dǎo)律，采用一種兩級(jí)IACG來(lái)全方位命中目標(biāo)。然而，在理論上，對(duì)于基于Lyapunov的制導(dǎo)律，只有當(dāng)時(shí)間接近無(wú)限時(shí)，狀態(tài)才會(huì)收斂到零。因此，其他一些研究涉及有限時(shí)間收斂IACG。文獻(xiàn)[10]從滑?？刂评碚摮霭l(fā)，提出了一種保證視線角在有限時(shí)間內(nèi)收斂的IACG。文獻(xiàn)[11]提出了一種非奇異終端滑?？刂疲∟TSMC）的有限時(shí)間內(nèi)收斂IACG，所得到的制導(dǎo)律能夠以期望碰撞角命中目標(biāo)。文獻(xiàn)[12]得到了另一種有限時(shí)間收斂的制導(dǎo)律，這種制導(dǎo)律是由傳統(tǒng)的NTSMC和二階滑?？刂葡嘟Y(jié)合得到的。

除上述方法外，還采用幾何和多項(xiàng)式方法推導(dǎo)了IACG。文獻(xiàn)[13]首次提出了多項(xiàng)式制導(dǎo)的方法。之后，文獻(xiàn)[14]提出了一種具有碰撞角和加速度約束的增廣多項(xiàng)式制導(dǎo)律，制導(dǎo)指令以多項(xiàng)式函數(shù)形式給出，系數(shù)對(duì)應(yīng)于碰撞角和加速度的約束。為了控制撞擊時(shí)間和碰撞角，文獻(xiàn)[15]提出了作為射程的函數(shù)的制導(dǎo)指令。文獻(xiàn)[16]最新研究了涉及幾何和多項(xiàng)式方法的自適應(yīng)制導(dǎo)律。

上述對(duì)IACG的研究大多是基于等速假設(shè)。事實(shí)上，導(dǎo)彈的速度會(huì)受到推力、空氣阻力及任務(wù)的進(jìn)展而發(fā)生變化。因此，本研究將在導(dǎo)彈變速的情況下解決碰撞角控制問(wèn)題。首先，通過(guò)問(wèn)題描述將制導(dǎo)指令設(shè)計(jì)問(wèn)題轉(zhuǎn)化為航跡角的設(shè)計(jì)問(wèn)題。然后，利用二次Bézier曲線擬合航跡角。最后，通過(guò)控制點(diǎn)的不同選擇來(lái)實(shí)現(xiàn)碰撞角的控制。本研究在提出勻速導(dǎo)彈制導(dǎo)律的基礎(chǔ)上，進(jìn)一步提出了變速導(dǎo)彈的制導(dǎo)律，并對(duì)其剩余飛行時(shí)間進(jìn)行了估計(jì)。

1 問(wèn)題描述

3 數(shù)值仿真分析

本節(jié)中，將通過(guò)數(shù)值仿真來(lái)驗(yàn)證所提策略的有效性。在每種情況下，仿真步長(zhǎng)均設(shè)為0.01 s。當(dāng)相對(duì)速度符號(hào)變?yōu)檎蛳鄬?duì)作用距離小于0.01 m時(shí)，仿真終止。導(dǎo)彈的初始位置為（0，0），目標(biāo)位置為（8000 m，0）。初始航跡角為30°。導(dǎo)彈初始速度為270 m/s，最大加速度為15 g，減速率為0.005。通過(guò)三個(gè)算例的仿真驗(yàn)證所提策略的有效性。

算例1 勻速導(dǎo)彈在相同初始條件下產(chǎn)生不同碰撞角

在此情況下，假設(shè)導(dǎo)彈的速度是恒定的，導(dǎo)彈期望碰撞角范圍是-90°到0。從期望范圍中選擇四個(gè)典型的期望碰撞角。仿真結(jié)果如圖2所示。

在圖2中，實(shí)線和另外三種虛線分別代表不同的碰撞角。由圖2d可知，所提出的制導(dǎo)律可以實(shí)現(xiàn)每個(gè)期望碰撞角。圖2a為導(dǎo)彈軌跡，圖2b為橫向加速度，由圖可知符合預(yù)定的加速度約束。圖2c中，航向誤差在終端時(shí)刻減小到零，保證了目標(biāo)的脫靶量為零。仿真結(jié)果驗(yàn)證了所提制導(dǎo)律在勻速情況下是有效的。

算例2 勻速情況下不同的初始航跡角產(chǎn)生相同的碰撞角

該算例為證明在所提出的制導(dǎo)律下，不同的初始航跡角可以獲得相同的碰撞角。給出三種不同的初始航跡角，分別為30°、60°和90°。仿真結(jié)果如圖3所示。

在圖3中，實(shí)線和另外兩種虛線分別代表不同的航跡角。圖3a為導(dǎo)彈軌跡，圖3b為橫向加速度，可以看出，當(dāng)導(dǎo)彈的初始航跡角較大時(shí)，在初始時(shí)刻會(huì)產(chǎn)生較大的加速度指令。由圖3c可知，所有航向角誤差均減小到零。從圖3d可知，導(dǎo)彈的初始航跡角不同，可以獲得相同的期望碰撞角。因此，仿真結(jié)果驗(yàn)證了所提制導(dǎo)律在勻速情況下是有效的。

算例3 變速情況下的不同碰撞角問(wèn)題

在此情況下，導(dǎo)彈速度不再被假定為常數(shù)。導(dǎo)彈的期望撞擊角范圍為-90°到90°。選擇四個(gè)典型的期望碰撞角進(jìn)行仿真，結(jié)果如圖4所示。

在圖4中，實(shí)線和另外三種虛線分別代表不同的碰撞角。由圖4d可知，所提出的制導(dǎo)律可以實(shí)現(xiàn)每個(gè)期望碰撞角。圖4b為橫向加速度，由圖可知其符合預(yù)定的加速度約束。由圖4c可知，航向誤差在終端時(shí)刻減小到零，保證了目標(biāo)的脫靶量為零。因此，仿真結(jié)果驗(yàn)證了該制導(dǎo)律在變速度情況下是有效的。

4 結(jié) 論

提出了一種速度不變的碰撞角約束制導(dǎo)律，并進(jìn)一步提出了速度時(shí)變的碰撞角約束制導(dǎo)律。該制導(dǎo)律是一條控制點(diǎn)待定的Bézier曲線形式的航跡角曲線。所提出的函數(shù)同時(shí)滿足航跡角的初始條件和終止條件，并通過(guò)求解積分方程確定了未知系數(shù)，同時(shí)對(duì)提出速度時(shí)變制導(dǎo)律的剩余飛行時(shí)間進(jìn)行了估計(jì)。通過(guò)不同的數(shù)值仿真，驗(yàn)證了所提出的制導(dǎo)律的有效性。在今后的工作中，制導(dǎo)律的設(shè)計(jì)還應(yīng)考慮切向加速度的影響。

參考文獻(xiàn)

[1] 趙曜，廖選平，遲學(xué)謙，等. 終端約束末制導(dǎo)律綜述[J]. 航天控制，2017，35 （2）： 89-98.

[2] 黎克波，廖選平，梁彥剛. 基于純比例導(dǎo)引的攔截碰撞角約束制導(dǎo)策略[J]. 航空學(xué)報(bào)，2020，41（S2）：79-88.

[3] 吳放，常思江. 攻擊時(shí)間和攻擊角度控制的非奇異終端滑模制導(dǎo)律[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2021，53（6）：94-103.

[4] KIM M，GRIDER K V. Terminal guidance for impact attitude angle constrained flight trajectories[J]. IEEE Transactions on Aerospace and Electronic Systems，1973，AES9 （6）： 852-859.

[5] RYOO C K，CHO H，TAHK M J. Optimal guidance laws with terminal impact angle constraint[J]. Journal of Guidance Control & Dynamics，2005，28 （4）：724-732.

[6] RYOO C K，CHO H，TAHK M J. Timetogo weighted optimal guidance with impact angle constraints[J]. IEEE Transactions on Control Systems Technology，2006，14 （3）：483-492.

[7] TSALIK R，SHIMA T. Optimal guidance around circular trajectories for impactangle interception[J]. Journal of Guidance Control & Dynamics，2016，39 （6）：1-14.

[8] KIM M，KIM Y. Lyapunovbased pursuit guidance law with impact angle constraint[C]. IFAC Proceedings Volumes，2014，47 （3）：2509-2514.

[9] CHENG Zhongtao，LIU Lei，WANG Yongji. Lyapunovbased switchedgain impact angle control guidance[J]. Chinese Journal of Aeronautics，2018，31（4）：765-775.

[10]ZHANG Yunxi，SUN Mingwei，CHEN Zengqiang. Finitetime convergent guidance law with impact angle constraint based on slidingmode control[J]. Nonlinear Dynamics，2012，70 （1）： 619-625.

[11]HOU Zhiwei，LIU Lei，WANG Yongji，et al. Terminal impact angle constraint guidance with dual sliding surfaces and modelfree target acceleration estimator[J]. IEEE Transactions on Control Systems Technology，2016，25（1）：85-100.

[12]HE Shaoming，LIN Defu，WANG Jiang. Continuous secondorder sliding mode based impact angle guidance law[J]. Aerospace Science & Technology，2015，41（2）：199-208.

[13]MIN B M，TAHK M J，SHIM H C D，et al. Guidance law for visionbased automatic landing of UAV[J]. International Journal of Aeronautical & Space Sciences，2007，8 （1）：4653.

[14]LEE C H，KIM T H，TAHK M J，et al. Polynomial guidance laws considering terminal impact angle and acceleration constraints[J]. Aerospace & Electronic Systems IEEE Transactions on，2013，49 （1）：7492.

[15]KIM T H，LEE C H，JEON I S，et al. Augmented polynomial guidance with impact time and angle constraints[J]. IEEE Transactions on Aerospace & Electronic Systems，2013，49 （4）：2806-2817.

[16]CHENG Zhongtao，WANG Bo，LIU Lei，et al. Adaptive polynomial guidance with impact angle constraint under varying velocity[J]. IEEE Access，2019，10（7）：104210-104217.

[17]馬帥，王旭剛，王中原，等. 帶初始前置角和末端攻擊角約束的偏置比例導(dǎo)引律設(shè)計(jì)以及剩余飛行時(shí)間估計(jì)[J]. 兵工學(xué)報(bào)，2019，40（1）：68-78.

[18]查世政，王永驥，蘇茂.基于kriging濾波的高精度視線轉(zhuǎn)率估計(jì)方法[J].計(jì)算技術(shù)與自動(dòng)化，2014，33 （4）：20-25.