考虑落角约束和气动不确定性的固定时间制导控制一体化设计方法

胡诚伟, 韦悠, 汪学平

弹箭与制导学报 ›› 2023, Vol. 43 ›› Issue (4) : 98-104.

PDF(1893 KB)
文章检索
PDF(1893 KB)
弹箭与制导学报 ›› 2023, Vol. 43 ›› Issue (4) : 98-104. DOI: 10.15892/j.cnki.djzdxb.2023.04.015

考虑落角约束和气动不确定性的固定时间制导控制一体化设计方法

作者信息 +

Fixed-time Integrated Guidance and Control for Impact Angle Constrained Interception with Multiple Uncertainties

Author information +
文章历史 +

摘要

文中研究落角约束和气动不确定性条件下的机动目标制导控制一体化设计问题。首先,将目标的加速度和气动不确定性考虑为有界的外部扰动。其次,将典型的反步法与动态曲面控制技术相结合,提出了一种固定时间滑模控制方法来解决制导控制一体化设计问题,并通过引入固定时间微分器有效避免反步法中存在的“微分膨胀”问题。然后,借助于固定时间稳定性理论对闭环系统进行了稳定性分析,表明所提方法能够使所有误差信号在任意初始条件下实现半全局固定时间一致最终有界稳定。最后,通过仿真验证了所提方法的有效性。

Abstract

This paper considers the integrated guidance and control (IGC) problem for maneuvering target interception with impact angle constraint and uncertain aerodynamics. The target’s acceleration and uncertain aerodynamics are considered as bounded external disturbances. By combining the typical backstepping idea with dynamic surface control (DSC) technique, a fixed-time sliding mode control approach is presented to address the IGC law design problem and the inherent problem of “explosion of complexity” in typical backstepping design is efficiently avoided by introducing a fixed-time differentiator. Furthermore, fixed-time stability theory is presented to carry out the stability analysis showing that with the proposed control scheme all signals in the sliding mode system can achieve semi-globally fixed-timely uniformly ultimately bounded stabilization independent of initial conditions. Finally, the superiority of the proposed IGC law is demonstrated by simulation results.

关键词

制导控制一体化设计 / 落角约束 / 机动目标拦截 / 固定时间稳定性 / 气动不确定性

Key words

integrated guidance and control / impact angle constraint / maneuvering target interception / fixed-time stability / aerodynamic uncertainties

引用本文

导出引用
胡诚伟, 韦悠, 汪学平. 考虑落角约束和气动不确定性的固定时间制导控制一体化设计方法[J]. 弹箭与制导学报, 2023, 43(4): 98-104 https://doi.org/10.15892/j.cnki.djzdxb.2023.04.015
HU Chengwei, WEI You, WANG Xueping. Fixed-time Integrated Guidance and Control for Impact Angle Constrained Interception with Multiple Uncertainties[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2023, 43(4): 98-104 https://doi.org/10.15892/j.cnki.djzdxb.2023.04.015
中图分类号: TJ765.3   

0 引言

现代战争中更加能动的战场战术和精度要求的提高使得经典控制系统设计方法在很多应用中不能满足要求。通常来说,制导系统和控制系统可以通过解耦分为两个独立的子系统,从而形成典型的两回路控制方案[1-2]。基于这种设计思路,已经有文献提出了包括比例导引[3-4]、鲁棒制导律[5-6]和滑模制导律[7-8]等多种有效的制导控制方法。上述方法虽然已成功应用于目标拦截,但采用分离设计的制导和控制系统,其拦截高机动目标时性能可能会严重下降。
为了提高导弹的杀伤效能,制导控制一体化设计已成为战术导弹制导与控制系统设计的新趋势。一般来说,现有制导控制一体化设计方法可以分为基于自动驾驶补偿的制导律设计方法和完全制导控制一体化设计方法[9]。对于第一类方案,通常将自动驾驶仪假设为一阶或二阶环节,然后在制导律设计中加入补偿来处理自动驾驶仪带来的响应滞后[10-11]。尽管这些算法能够实现目标拦截,但其性能严重依赖于对自动驾驶仪近似估计的准确性。然而,对于大空域机动的导弹,其空气动力学特性本身就存在着较大的不确定性,这使得自动驾驶仪很难被精确地近似为简单的一阶或二阶环节。而完全制导控制一体化设计方法能够避免对自动驾驶仪的近似估计,直接针对完全非线性导弹动力学模型设计制导控制方法,其中典型的方法包括反馈线性化方法[12]、状态依赖Riccati方法[13]θ-D方法[14]、模型预测静态规划方法[15]和有限时域最优控制[16]。然而,完全制导控制一体化设计方法计算较为复杂,且大都无法保证闭环系统的鲁棒性[17]
由于导弹制导和控制系统具有严格的反馈形式[18],因此反步设计思想在制导控制一体化设计中也被广泛应用。然而,在反步设计过程中对虚拟控制量进行重复微分会导致“微分膨胀”问题[19]。为了解决这个问题,动态面控制技术被用于制导控制一体化设计中。文献[20]针对慢机动目标,基于自适应动态面控制给出了一种制导控制一体化设计方法。文献[21]通过引入了模糊系统来估计和补偿不确定性,提出了一种基于自适应模糊动态面的制导控制方法。文献[22]提出了一种二阶滑模微分器来估计虚拟控制量的导数,实现了考虑气动不确定性和落角约束的制导控制一体化设计。然而,上述基于反步控制和动态面思想给出的方法只能实现半全局渐近最终有界稳定,而非有限时间稳定[23]。因此研究有限时间收敛的制导控制一体化设计方法具有重要的实际意义。作为有限时间稳定性理论的扩展,固定时间稳定性引起收敛时间上界不依赖于初始条件而被广泛的关注。然而,固定时间制导方法仅适用于二阶系统,无法直接应用到导弹的制导控制一体化设计问题。
基于上述分析,文中研究考虑落角约束和气动不确定性的固定时间制导控制一体化设计问题。基于反步法提出了固定时间滑模制导控制一体化设计方法,并借助于固定时间微分器解决了传统反步设计中存在的“微分膨胀”问题,最终实现了闭环系统的半全局固定时间一致最终有界稳定。主要价值在于:1)提出了一种固定时间制导控制一体化设计方法,并实现了不受初始条件影响的半全局固定时间一致最终有界稳定;2)在制导控制一体化设计中考虑了目标机动不确定性和导弹气动不确定性,使得所提出方法更加适用于实际应用。

1 问题描述和预备知识

考虑铅垂面内落角约束的制导控制一体化设计问题。图1为导弹与目标相对运动示意图,其中VMVT为导弹和目标速度且满足VT/VM>1;aMaT为导弹和目标垂直于速度方向的加速度;γMγT为导弹和目标弹道倾角;rλ分别表示导弹与目标之间的相对距离和视线角。为简单起见,导弹和目标被视为匀速运动的点质量。
图1 导弹与目标相对运动

Fig.1 Relative motion between missile and target

Full size|PPT slide

导弹和目标间相对运动可以表示为:
r˙=VTcos(γTλ)VMcos(γMλ)
(1)
rλ˙=VTsin(γTλ)VMsin(γMλ)
(2)
γ·M=aMVM
(3)
γ·T=aTVT
(4)
λd=γTarctan2sin γdcos γdVT/VM
(5)
文中所研究的带落角约束的制导控制一体化设计问题可以表述为设计控制输入δz使得导弹以期望视线角λd与目标交会。
为了给出主要结果,考虑以下系统:
x˙(t)=f(x(t)),x(0)=x0,xRn
(6)
假设V(x):Rn→R是一个正定函数。如果存在正实数k1k2一个任意小的正实数ε,正奇数p1,q1,p2,q2满足p1>q1,p2<q2和一个闭集X0,使得 V·≤-k1 Vp 1/q 1-k2 Vp 2/q 2+ε,x0X0,那么式(6)是半全局固定时间一致最终有界的,并且可以在固定时间内稳定到原点附近任意小的邻域,其上界为:
Tf1k1q1p1-q1+1k2(2p 2/q 2-1)q2p2-q2
(7)

2 集成制导与控制律设计及稳定性分析

2.1 制导控制一体化设计

为了给出制导控制一体化设计方法,需借助于以下假设:
存在两个正实数 r¯, r_,使得在末制导阶段弹目相对距离r满足 r_<r< r¯,其中下界 r_可以根据目标物理尺寸得出。
基于上述假设,下面将根据典型的反步设计过程逐步给出制导控制方法。
步骤1:令s1=λ-λd,s1对时间求导可得:
s˙1=Vλλ˙d=Vλ+Δ1
(8)
其中,Δ1=-aT/VT被视为外部干扰,其上界d1>0,即|Δ1|≤d1
为式(8)选择以下虚拟控制律:
Vλ*=(k1|s1|pq+c1|s1|qp+ρ1)sgn(s1)
(9)
其中:k1,c1,p,q是正实数;ρ1≥d1
在传统的反步控制设计中,虚拟控制量 Vλ*会被反复微分,导致“微分膨胀”问题。为了解决这个问题,可以借助具有时间常数τ1>0的固定时间微分器Vλ,d Vλ*进行估计:
τ1V·λ,d=-|Vλ,d-Vλ*|pqsgn(Vλ,d-Vλ*)-|Vλ,d-Vλ*|qpsgn(Vλ,d-Vλ*)Vλ,d(0)=Vλ*0
(10)
步骤2:令s2=Vλ-Vλ,d,s2对时间求导得出:
s˙2=aTcos(γTλ)raMcos(γMλ)r2r˙λ˙rV˙λ,d=2r˙λ˙rV˙λ,d+gcos(γMλ)cosλMrQScy,nαcos(γMλ)mrα+Δ2
(11)
其中Δ2=aTcos(γT-λ)/r-QS cy,uααcosM-γ)/mr被具有上边界d2>0的外部扰动,即|Δ2|≤d2。对于子系统式(11),考虑以下虚拟控制律:
α=mrQScy,nαcos(γMλ)(2r˙γ˙rV˙γ,d+s1+gcos(γMλ)cosγMr+k2|s2|pqsgn(s2)+k2|s2|pqsgn(s2)+c2|s2|qpsgn(s2)+ρ2sgn(s2))
(12)
其中:k2和c2为正数;ρ2≥d2。引入新状态变量αd其动力学由式(13)给出:
τ2α·d=-|αd-α*|pqsgn(αd-α*)-|αd-α*|qpsgn(αd-α*)αd(0)=α*0
(13)
其中τ2>0为时间常数。
步骤3:令s3=α-αd,并对s3求导得出:
s·3=ωz-QScy,nααmVM+gcos γMVM-α·d+Δ3
(14)
其中Δ3=-QS cy,uαα/mVM被认为是具有上边界d3>0的外部扰动,即|Δ3|≤d3。对于子系统式(13),考虑以下虚拟控制律:
ωz=k3|s3|pqsgn(s3)c3|s2|qpsgn(s3)+α˙dρ3sgn(s3)+QScy,nααmVMgcosγMVM+QScy,nαcos(γMλ)mrs2
(15)
其中k3和c3为正数,且ρ3≤d3。引入新状态变量ωd,其动力学满足:
τ3ω·d=-|ωd-ω*|pqsgn(ωd-ω*)-|ωd-ω*|qpsgn(ωd-ω*)ωd(0)=ω*0
(16)
其中τ3>0为时间常数。
步骤4:令s4zz,d,并对s4求导得:
s·4=QSlmz,nααJz+QSl2mz,nωzωzVMJz+QSlmz,nδzJzδz- ω·z,d+Δ4
(17)
其中 Δ4=QSlmz,uααJz+QSl2mz,uωzωzVMJz+QSlmz,uδzδzJz视为具有上边界d4>0的外部干扰,即|Δ4|≤d4
对于子系统式(16),考虑控制律:
δz=JzQSLmz,nδδ(k4|s4|pqsgn(s4)c4|s4|qpsgn(s4)+ω˙z, ds3ρ4sgn(s4)QSlmz,nααJzQSl2mz,nωzωzJzVM)
(18)
其中:k4和c4是正数;ρ4≥d4

2.2 稳定性分析

给出稳定性分析,以证明提出方法具有半全局固定时间稳定性。首先定义以下误差变量:
e1=Vλ,d-Vλ*e2=αd-α*e3=ωz,d-ωz*
(19)
对于使用文中所提控制律的闭环系统,存在足够大的正常数ki,ci(i=1,2,3,4)、足够小的正常数τi(i=1,2,3)和正奇数p,q,满足p>q且(p+q)/2是正奇数,使得闭环系统是半全局一致最终有界的,同时在任意初始条件下所有误差均可以在固定时间Ts内稳定在原点任意小的邻域中。
考虑以下李雅普诺夫函数:
V1=12i=14si2+12i=13ei2
(20)
选择设计参数使得:
k3> γ1-pq+1pq+1,k2> δγ2-pq+1pq+1,k3> γ3-pq+1pq+1,k4>0,
1τ1>max γ1qp+1qp+1,η1pq+1pq+1, 1τ2>max δγ2qp+1qp+1,η2pq+1pq+1,
1τ3>max γ3qp+1qp+1,η3pq+1pq+1和ci>0 (i=1,2,3,4)。
σ1=min{k1γ1(pq+1)pq+1,k2δγ2(pq+1)pq+1,k3γ3(pq+1)pq+1,k4,1τiηipq+1pq+1},i=1,2,3,4
(21)
σ2=minci,1τ1-γ1qp+1qp+1,1τ2-δγ2qp+1qp+1,1τ3-γ3qp+1qp+1
(22)
σ3=i=13ηi-qp+1qp+1|Θi|qp+1
(23)
计算所定义李雅普诺夫函数的导数。
假设Ω>0为闭环系统的初始条件的边界,即V1(0)≤Ω。因此总可选择合适的参数,使得以下不等式成立:
-2 pq+12·7 1-pq2σ1 V1 pq+12-2 qp+12σ2 V1 qp+12+σ30
(24)
因此,{V1≤Ω}是一个不变集,并且对于所有t>0,V1(t)≤Ω都成立,这表明整个闭环系统是最终有界的。闭环系统可以固定时间内稳定在原点的一个小邻域中,收敛时间上界为:
Ts=12 pq+12·7 1-pq2σ1·2qp-q+12 qp+12σ22 p+q2p-1·2pp-q
(25)

3 仿真结果

3.1 参数设置

通过数值仿真验证文中提出的制导控制一体化设计方法的有效性。假设目标加速度为aT=10·sin(0.1πt)。导弹的气动参数为: QScy,nαmVM=0.3487; QSLmz,nαJz=-17.801; QSL2mz,nωzJzVM=-0.2741; QSLmz,nδzJz=-56.267。
假设空气动力系数相对于标称值有30%的不确定性。不失一般性,设计参数可选择为p=9,q=6,k1=1,k2,=2,k3=5,k4=10,ci=0.08(i=1,2,3,4),τi=0.01(i=1,2,3)。导弹和目标的初始条件见表1
表1 初始条件

Table 1 Initial condition

Parameters Value
Missile initial position/m 0
0
Missile initial velocity/(m/s) 500
Target initial velocity/(m/s) 250
Initial ballistic inclination
angle of missile/(°)
45
Initial ballistic inclination
angle of target/(°)
120
Initial relative distance/m 10 000
Initial LOS angle/(°) 30
Initial angle of attack of missile/(°) 0
Initial pitch angular velocity of
missile/((°)/s)
0

3.2 不同落角约束仿真

给出不同落角约束下机动目标拦截的仿真结果。期望落角γd分别选择为0°、30°、60°和90°。仿真结果如图2所示。
图2 不同落角约束下的仿真结果

Fig.2 Simulation result of different impact angle constraint

Full size|PPT slide

末端脱靶量和落角误差如表2所示。从表2图2可以看出,文中所提出的制导控制一体化设计方法针对不同落角约束下的目标拦截问题,都能够取得较小的脱靶量和落角误差。具体而言,在所有拦截场景中落角误差和脱靶量分别小于0.004°和0.4 m。从图2还可以看出,在拦截过程中舵偏角在允许的最大值范围内且较为平滑。
表2 脱靶量与落角误差

Table 2 The miss distance and impact angle error

Expected impact
angle/(°)
Miss
distance/m
Impact angle
error/(°)
0 0.107 45 0.001 60
30 0.302 95 0.003 55
60 0.266 29 0.000 84
90 0.269 40 0.002 29

3.3 对比仿真

为了突出文中所提出的制导控制一体化设计方法的优越性,给出了与非奇异终端滑模制导律文献[7]和固定时间非奇异终端滑模制导律文献[8]的对比仿真。加速度约束为|aM|≤200 m/s2,期望落角为γd=45°仿真结果如图3所示。
图3 对比仿真

Fig.3 Comparative simulation

Full size|PPT slide

从图中可以看出与文献[7-8]中提出的方法相比,文中所设计的制导控制一体化方法能够以相对更小的加速度和更短的攻击时间获得相对较高的命中精度。

4 结论

文中针对考虑气动不确定性和落角约束的机动目标拦截问题,基于固定时间稳定性理论提出了一种制导控制一体化设计方法。利用典型的反步设计思想构建了滑模控制器,实现半全局固定时间一致最终有界稳定;同时通过引入固定时间微分器避免了典型反步设计存在的“微分膨胀”问题;并给出了李雅普诺夫函数框架下的稳定性分析。仿真结果与对比实验充分证明了文中所提出的制导控制一体化设计方法的有效性。

参考文献

[1]
YAN H, TAN S P, HE Y Z. A small-gain method for integrated guidance and control in terminal phase of reentry[J]. Acta Astronautica, 2017, 132: 282-292.
[2]
WANG X H, WANG J Z. Partial integrated guidance and control with impact angle constraints[J]. Journal of Guidance, Control, and Dynamics, 2015: 38(5): 925-936.
[3]
YANG C D, YANG C C. Analytical solution of three-dimensional realistic true proportional navigation[J]. Journal of Guidance, Control, and Dynamics, 2012, 19(3): 569-577.
[4]
OHIMEYER E. Root-mean-square miss distance of proportional navigation missile against sinusoidal target[J]. Journal of Guidance, Control, and Dynamics, 2015, 19(3): 563-568.
[5]
ZHOU D, MU C, XU W. Adaptive sliding-mode guidance for a homing missile[J]. Journal of Guidance, Control, and Dynamics, 1999, 22(4): 589-594.
[6]
MOON J K, KIM Y D. Design of missile guidance law via variable structure control[J]. Journal of Guidance, Control, and Dynamics, 2001, 24(4): 659-664.
[7]
KUMAR S, RAO S, GHOSE D. Nonsingular terminal sliding mode guidance with impact angle constraints[J]. Journal of Guidance, Control, and Dynamics, 2014, 37(4): 610-621.
[8]
ZHANG Y, TANG S J, GUO J. Adaptive terminal angle constraint interception against maneuvering targets with fast fixed-time convergence[J]. International Journal of Robust and Nonlinear Control, 2018, 2: 2996-3014.
[9]
GURFIL P, JODORKOVSKYM M, GUELMAN M. Finite time stability approach to proportional navigation systems analysis[J]. Journal of Guidance, Control, and Dynamics, 1998, 21(6): 853-861.
[10]
CHWA D, CHOI J. Adaptive nonlinear guidance law considering control loop dynamics[J]. IEEE Transactions on Aerospace and Electronic Systems, 2003, 39(4): 1134-1143.
[11]
SUN S, ZHOU D, HOU W T. A guidance law with finite time convergence accounting for autopilot lag[J], Aerospace Science and Technology, 2013, 25(1): 132-137.
[12]
MENON P, OHLMEYER E. Integrated design of agile missile guidance and autopilot systems[J]. Control Engineering Practice, 2001, 9(10): 1095-1106.
[13]
VADDO S, MENON P, OHLMEYER E. Numerical state-dependent Raccati equation approach for missile integrated guidance control[J]. Journal of Guidance, Control, and Dynamics, 2009, 32(2): 699-703.
[14]
XIN M, BALAKRISHNAN S N, OHLMEYER E J. Integrated guidance and control missiles with method[J]. IEEE Transactions on Control Systems Technology, 2006, 14(6): 981-992.
[15]
RADHAKANT P, CHARU C, PRIYA D G. Partial integrated guidance and control of interceptors for high-speed ballistic targets[J]. Journal of Guidance, Control, and Dynamics, 2014, 37(1): 149-163.
[16]
KIM J H, WHANG I H, KIM B M. Finite horizon integrated guidance and control for terminal homing in vertical plane[J], Journal of Guidance, Control, and Dynamics, 2016, 39(5): 1104-1112.
[17]
HE S M, LIN D F, WANG J. Continuous second-order sliding mode based impact angle guidance law[J]. Aerospace Science and Technology, 2015, 41: 199-208.
[18]
SHIN H S, HWANG T, TSOPURDOS A, et al. Integrated intercept missile guidance and control with terminal angle constraint[C]// ICAS. Proceedings of the 26th International Congress of the Aeronautical Sciences. Brisbane: ICAS, 2008: 1875-1884.
[19]
SHEN Y C, YU J Q, LUO G C, et al. Observer-based adaptive sliding mode backstepping output-feedback DSC for spin-stabilized canard-controlled projectiles[J]. Chinese Journal Aeronautics, 2017, 30(3): 1115-1126.
[20]
HOU M Z, LIANG X L, DUAN G R. Adaptive block dynamic surface control for integrated missile guidance and autopilot[J]. Chinese Journal of Aeronautics, 2013, 26(3): 741-750.
[21]
RAN M P, WANG Q, HOU D L, et al. Backstepping design of missile guidance and control based on adaptive fuzzy sliding mode control[J]. Chinese Journal of Aeronautics, 2014, 27(3): 634-642.
[22]
HE S M, SONG T, LIN D F. Impact angle constrained integrated guidance and control for maneuvering target interception[J]. Journal of Guidance, Control, and Dynamics, 2017, 40(10): 2652-2660.
[23]
NI J K, LIU L, LIU C X, et al. Fast fixed-time nonsingular terminal sliding mode control and its application to chaos suppression in power system[J]. IEEE Transactions on Circuits and Systems, 2017, 64(2): 151-155.
[24]
NI J K, LIU L, LIU C X, et al. Fixed-time dynamic surface high-order sliding mode control for chaotic oscillation in power system[J]. Nonlinear Dynamics, 2016, 86(1): 401-420.

PDF(1893 KB)

27

Accesses

0

Citation

Detail

段落导航
相关文章

/