一种带终端视线角约束的反舰导弹三维协同制导律.pdf

收稿日期：2022-04-16修回日期：2022-06-10作者简介：张书森（1998），男，河北吴桥人，硕士研究生。研究方向：导航制导与控制。摘要：针对导弹的协同攻击问题，设计了一种带终端视线角约束的三维协同制导律，在视线坐标系下建立导弹-目标相对运动方程，基于一致性理论设计沿视线方向上的加速度使得多枚导弹同时攻击目标，基于有限时间理论设计沿视线法向上的加速度，使得各枚导弹以期望的视线角命中目标，并利用李雅普诺夫理论对所设计的制导律给出了稳定性证明，针对 4 枚导弹同时攻击海平面上机动的舰船这一场景进行了仿真，仿真结果表明所设计的制导律能够满足期望指标要求，具有较高的打击精度，并且所设计的制导律在导弹通讯拓扑结构发生变化时同样有效。关键词：协同制导律；终端视线角约束；有限时间理论；一致性理论；终端滑模中图分类号：TJ765.3文献标识码：ADOI：10.3969/j.issn.1002-0640.2023.06.013引用格式：张书森，孟秀云，丁晓.一种带终端视线角约束的反舰导弹三维协同制导律 J.火力与指挥控制，2023，48（6）：84-90.一种带终端视线角约束的反舰导弹三维协同制导律张书森，孟秀云，丁晓（北京理工大学宇航学院，北京100081）A Three-dimensional Cooperative Guidance Law for Anti-ship Missilewith Terminal Line of Sight Angle ConstraintZHANG Shusen，MENG Xiuyun，DING Xiao（School of Aerospace Engineering Beijing Institute of Technology，Beijing 100081，China）Abstract：Aiming at the cooperative attack problem of multiple missiles，a three-dimensionalcooperative guidance law with terminal line of sight angle constraint is designed.The relative motionequation of missile and target is established in the line of sight coordinate system.Based on theconsistency theory，the acceleration along the line of sight direction is designed to make multiplemissiles attack the target at the same time.Based on the finite time theory，the acceleration up thenormal line of sight is designed to make each missile hit the target at the expected line of sight angle，and the stability of the designed guidance law is proved with Lyapunov theory.Then，the simulation iscarried out for the scene that four missiles attack the maneuvering ships on the sea level at the sametime.The simulation results showed that the guidance law designed in this paper could meet theexpected index requirements with high striking accuracy.The designed guidance law is also effectivewhen the missile communication topology structure changes.Key words：cooperative guidance law；terminal line of sight angle constraint；finite time theory；consistency theory；terminal sliding modeCitation format：ZHANG S S，MENG X Y，DING X.A three-dimensional cooperative guidance lawfor anti-ship missile with terminal line of sight angle constraint J.Fire Control&Command Control，2023，48（6）：84-90.文章编号：1002-0640（2023）06-0084-07Vol.48，No.6Jun，2023火 力 与 指 挥 控 制Fire Control&Command Control第 48 卷第 6 期2023 年 6 月84（总第 48-）0引言现代战争态势不断发展，反导系统愈加完善，仅靠单枚导弹难以满足作战需要，多枚导弹协同攻击能够有效弥补单枚导弹突防能力不足、抗干扰能力差的缺陷，因此，多枚导弹协同攻击更符合现代化信息作战的需要1。最早的多导弹协同攻击策略是通过给各枚导弹设置相同的攻击时间实现的，导弹个体之间互不干扰、独立作战，归根结底研究的还是单枚导弹的制导问题2，严格意义上讲不属于真正的协同，真正的协同无需在导弹发射前预置飞行时间，而是通过导弹在飞行过程中不断的信息交互进而达到攻击时间一致3，为了实现对目标的全方位打击，在多导弹时间协同的基础上，还要求导弹以期望的角度攻击目标，将导弹-目标相对运动方程建立在视线坐标系下是一种较为常用的方法4-6，沿视线方向设计制导律可以协调导弹飞行的时间，沿视线法向设计制导律可以使导弹以特定的视线角攻击目标。现有的制导律大多基于二维平面内导弹-目标相对运动模型设计的，针对三维空间设计制导律的文献较少，本文针对多枚导弹同时攻击海面上机动的舰船这一场景，设计一种有终端视线角约束的三维协同制导律，在视线方向上，区别于文献 4-6采用构造关于时间的滑模面的方式，本文采用多智能体一致性理论设计时间协同制导律，制导算法相对简洁；在视线法向上，本文设计一种带终端视线角约束的非奇异自适应滑模制导律，精度高且收敛速度快，最后选取特定场景对所设计的制导律进行仿真验证。1导弹目标三维相对运动模型本文将导弹-目标相对运动方程建立在视线坐标系下，如图 1 所示，MXLYLZL为视线坐标系，OXIYIZI为惯性坐标系，q兹、q渍分别为视线倾角和视线偏角。相对速度在视线坐标系下的表达如下：（1）式中，VR为弹目接近速度，棕 为视线坐标系相对大地坐标系的旋转角速度，R 为弹目相对距离。在视线坐标系下有：（2）根据视线坐标系与大地坐标系的转换关系有：（3）将式（2）、式（3）代入式（1）可得：（4）根据牛顿第二定律有：（5）将式（3）、式（4）代入式（5）可得三维空间下导弹-目标相对运动动力学方程组形式如下：（6）其中，aM=aMxaMyaMzT为导弹在视线坐标系下三轴方向的加速度，aT=aTxaTyaTzT为目标在视线坐标系下三轴方向的加速度，本文的目的在于通过设计沿视线方向上的加速度 aMx，使得多枚导弹同时击中目标；通过设计沿视线法向上的加速度 aMy、aMz，一方面使得视线角速率在导弹接近目标的过程中快速趋近于零，另一方面控制导弹在末端以期望的视线角攻击目标。2三维协同制导律设计2.1视线方向上制导律设计定理 17：如果存在一个含有 n 个节点的系统可以用 G=（V，E，A）描述，并且节点的状态可以用一阶动态环节表达，那么对于任意给定的初值 孜i（0），坌i沂V，当 t寅肄 时，有 孜i=孜j，坌i，j沂V，则称系统节点之间的状态能够达到渐进一致。（7）式中，孜i沂Rm用于描述第 i 个节点的状态；ui沂Rm用图 1三维空间弹目相对运动几何关系图Fig.1Geometric relation diagram of relative motion betweenmissile and target in three-dimensional space张书森等：一种带终端视线角约束的反舰导弹三维协同制导律851039（总第 48-）火 力 与 指 挥 控 制2023 年第 6 期于描述第 i 个节点的输入。一种常见的一致性算法表达形式如下：（8）假设节点间的通信拓扑结构保持固定，采用式（8）的控制算法，使得系统状态能够达到渐进一致的必要条件是 G 中存在有向生成树，系统最终的一致状态由状态初值、通信拓扑结构和加权矩阵共同决定8。特别地，如果 G 是无向图，那么系统状态将会达到平均一致，一致状态为：为方便起见，令导弹同时启动发射的时刻为时间原点，记为 t=0 时刻，在此后的飞行过程中，记 Ri为 t 时刻第 i 枚导弹与目标的相对距离，vi为 t 时刻第 i 枚导弹与目标的接近速度，Ri、vi均是关于时间 t的函数。对 Ri关于时间 t 求导有：（9）式中，Ri的值可以根据导弹和目标的位置坐标求出，vi的值可以分别将导弹和目标的速度投影到视线坐标系下得到。根据求得的 Ri和 vi，即可估计出第 i 枚导弹击中目标的时刻为：（10）式中，i=Ri/vi为第 i 枚导弹在当前时刻下的期望攻击时间（expected time-of-attack，ETA）。根据式（10），显然有：（11）由式（11）可以看出，多枚导弹同时击中目标即击中时刻一致，等价于在某一时刻 t 之后各导弹的期望攻击时间一致。将式（10）两边同时对时间 t 求导，并将式（9）代入可得：（12）将式（9）两边同时对时间 t 求导：（13）将式（13）与式（6）第 1 式代入式（12）可得：（14）令，则式（14）化简为：（15）令 ui=琢i（vic-vi）i/vi，将对 aMx的控制转化为对 ui的控制，显然有（16）式（16）满足定理 1，灼i为状态变量，因此，上述所设计的算法能够使各枚导弹的击中时刻达到渐进一致，即实现多导弹协同攻击。由于 灼i与i等价，选取i作为状态变量，参照式（8）、式（15）可以得到一种简洁的分布式控制算法：（17）根据式（17），通过设计导弹之间的通讯拓扑结构和加权邻接矩阵，即可计算出各枚导弹的速度指令，从而得到所需要的视线方向加速度指令，进而实现多枚导弹期望攻击时间一致。2.2视线法向上制导律设计在进行视线方向上的制导律设计之前，首先简要介绍有限时间理论，方便后续制导律的设计与证明。定理 2 考虑如下非线性系统9：（18）假设在系统原点的邻域 U沂Rn内存在一个连续函数 V（x，t），并且函数 V（x，t）具有一阶连续导数，如果再存在一个大于零的实数 琢，使得 V（x，t）在U 上是正定的，并且使得V觶（x，t）+琢V姿（x，t）在 U 上是半负定的，其中，姿（0，1），那么称系统（18）的原点可以达到有限时间稳定。证明：因为已知V觶（x，t）+琢V姿（x，t）0，那么有（19）因为 V（x，t）是正定的，将式（19）左右两边同时对时间 t 积分，可以得到（20）且有终端时刻 tr取决于初始值 x（0）=x0，那么有V（x，0）=V（x0），上界满足（21）其中，x0是系统原点附近某一开域的任意一点，证毕。在导弹飞向目标的过程中，导弹和目标相对距离的变化率小于零，因此，只要控制视线角的变化率为零，就能保证导弹击中目标10，同时考虑到有861040（总第 48-）期望终端视线角的约束，选取视线角速率和视线角误差为状态变量，法向加速度为控制输入，具体定义如下：其中，q兹f、q渍f为期望的终端视线倾角和视线偏角。根据三维空间内导弹-目标相对式（6）可以得到状态方程表达如下：（22）对于上述系统构造滑模面：（23）式中，s1、s2的表达形式设计如下：（24）其中，茁1、茁2=const0，0琢1、琢2=const1，s1表示纵向运动滑模面，s2表示侧向运动滑模面。在趋近律方面构造一种自适应滑模趋近律，以保证系统有良好的动态特性：（25）式中，着=着1着2T。它的物理意义是在制导初期弹目相对距离 R 较大，趋近滑模面的速率要小，而随着导弹不断地逼近目标，趋近滑模面的速率不断加快，确保导弹在击中目标前系统能够进入滑动模态。将式（24）两边同时对时间 t 求导并将式（22）、式（25）代入可以得到制导律形式如下：（26）值得注意的是：控制率 u1、u2中含有 x1和 x3的指数项，当 x1和 x3趋近于 0 时会出现奇异问题，因此，采用非线性方法对其进行限幅处理，即：（27）至此，有终端视线角约束的三维制导律设计完毕，利用式（26）中的制导方法，能够保证三维制导系统（22）在有限时间内收敛到滑模面 s=0，并且在滑模面 s=0 上系统的状态也能在有限时间内收敛到零。证明：选取李雅普诺夫函数 V1=12sTs，显然 V1正定。对 V1求导，并将制导律表达式（26）代入，整理得到：（28）式中，着=min（着1，着2）由定理 2 可知，制导系统（22）能够在有限时间内收敛到滑模面 s=0 上。而在滑模面 s1=0 上，有：（29）选取新的李雅普诺夫函数 V2=x12，显然 V2正定。对 V2求导，有：（30）当 琢1（0，1）时，有（琢1+1）2（0，1），根据定理2 可知，系统状态在终端滑模面 s1=0 能够在有限时间内收敛到零，同理可证系统状态在终端滑模面s2=0 也能够在有限时间内收敛到零。3仿真结果与分析为了说明本文所设计的三维空间内有终端视线角约束的协同制导律的有效性，接下来针对 4 枚导弹同时攻击一个在海平面上作机动航行的舰船这一场景，对制导律表达式（17）、式（26）进行仿真验证，4 枚导弹之间的通讯拓扑结构由图 2 给出，4枚导弹的初始位置、初始速度以及所期望的末端攻击视线角如表 1 所示，每枚导弹能承受的最大过载为 50，g 取 9.8 m/s2，仿真步长取 10 ms，同时为了抑制滑模态出现抖振现象，用以下的饱和函数来替代符号函数 sgn（s），并取 啄=0.01。张书森等：一种带终端视线角约束的反舰导弹三维协同制导律871041（总第 48-）火 力 与 指 挥 控 制2023 年第 6 期图 3弹目三维运动轨迹Fig.33-D trajectory of missile and target movement图 4弹目相对距离Fig.4Relative distancebetween missile and target目标舰船的初始位置为（5 000 m，0 m，5 000 m），航行速度 VT=30 m/s 恒定，初始航迹偏角为 10 ，并以 aT=5sin（0.01t）m/s2的加速度作侧向机动，控制参数选取如下：琢1=琢2=0.9，茁1=0.4，茁2=0.2，着1=着2=0.2。图 2导弹间的通讯拓扑Fig.2Communication topology between missiles表 1导弹 Mi初始条件和期望的视线角（1，2，3，4）Table 1Initial conditions and expected line of sight angle ofmissile Mi（i=1，2，3，4）4 枚导弹的三维运动轨迹、弹目相对距离、预估剩余飞行时间、沿视线方向的过载 nx、纵向平面内沿视线法向的过载 ny、侧向平面内沿视线法向的过载nz、视线倾角、视线偏角、视线倾角角速率和视线偏角角速率随时间的变化曲线如图 3 图 12 所示。由图 3 和图 4 可以看出，即使起始点位置不同，但在本文设计协同制导律的控制作用下，4 枚导弹依然能够从不同的方向沿不同路径同时击中在海平面上机动的目标。由图 5 可以看出，在分布式控制算法（17）的作用下，4 枚导弹在制导初期通过信息交互协调剩余攻击时间，大约在 5 s 时各导弹的剩余攻击时间能够达成一致。图 6 反映了 4 枚导弹沿视线方向的过载随时间的变化情况，过载在前5 s 变化迅速，之后逐渐趋于稳定，这是因为在制导前期需要过载指令的变化从而控制各导弹的飞行速度，进而使各导弹的剩余攻击时间达成一致，实现对目标的协同打击，这与图 5 所示的现象吻合。图 7 和图 8 反映了 4 枚导弹沿视线法向的过载随时间的变化情况，为了使视线角速率快速收敛至零，导弹过载在前期指令较大且变化较为剧烈，之后随着视线角速率逐渐趋近于零导弹过载也随之变为平缓并收敛至零，过载曲线光滑且抖振抑制效果明显。图 9 图 12 表明，在本文所设计的协同制导律的作用下，弹目视线倾角和视线偏角能在有限时间内收敛到期望值，同时，弹目视线倾角角速率和视线偏角角速率能在有限时间内收敛到零。以上的分析结果表明，本文设计的有终端视线角约束的协同制导律，能够保证 4 枚导弹以期望的终端视线角同时击中目标，验证了所设计的协同制导律的有效性。以上仿真是假定导弹间的通讯拓扑结构固定不变的，下面考虑拓扑结构切换的情况，文献 11表明对于一致性算法式（8），即使系统存在随时间切变的拓扑结构，但只要保证每一相邻有限时间隔内所有切换子图的并集存在有向生成树，系统的状i初始位置/m初始速度/（m/s）期望视线倾角/deg期望视线偏角/deg1（0，8 000，0）300-30-402（4 000，7 000，0）320-60-453（0，5 000，8 000）260-90104（2 000，10 000，0）280-50-90图 5剩余时间Fig.5Time left881042（总第 48-）态还是能够在有限时间内收敛到一致，只不过子图切换的频率会影响系统的收敛速度。在以上仿真条件的基础上，假设第 3 枚导弹由于通讯故障在 2 s 时刻与系统失去链接，在 10 s 时刻恢复链接，即导弹之间的通讯拓扑按照图 3 所示的拓扑结构切换。4 枚导弹的弹目相对距离和预估剩余时间的变化曲线如下页图 14 和图 15 所示。由图 14 可以看出，2 s 时第 3 枚导弹与系统失图 13通讯拓扑切换Fig.13Communication topology switching图 6视线方向过载指令 nxFig.6Overload command along LOS图 7视线法向过载指令 nyFig.7Normal overload command nyof LOS图 8视线法向过载指令 nzFig.8Normal overload command nzof LOS图 9视线倾角Fig.9Line of sight inclinationangle图 10视线偏角Fig.10Line of sight deflectionangle图 11视线倾角角速度Fig.11Line of sight inclinationpalstance图 12视线偏角角速度Fig.12Line of sight deflection palstance张书森等：一种带终端视线角约束的反舰导弹三维协同制导律891043（总第 48-）火 力 与 指 挥 控 制2023 年第 6 期图 14剩余时间Fig.14Time left图 15弹目相对距离Fig.15Relative distance between missile and target去链接，导致其无法与其他 3 枚导弹达成时间一致，10 s 后恢复链接，新的拓扑图存在有向生成树，使得4 枚导弹重新达成时间一致，图 15 表明即使系统在时间段内存在无向的通讯拓扑结构图，但在整个有限时间段内系统所有子图的并图存在有向生成树，所以依然能够保证 4 枚导弹同时击中目标。4结论本文针对导弹的协同攻击问题，在视线坐标系下建立了导弹-目标相对运动模型，在此基础上基于一致性理论设计了视线方向上的制导律，基于有限时间理设计了视线法向上的非奇异滑模制导律；对所设计的制导律给出了稳定性证明；通过数学仿真实验检验制导律的性能品质。仿真结果表明，所设计制导律能够使得多枚导弹以期望的视线角同时击中目标，具有一定的工程意义。本文考虑了目标机动方式已知的情况，后续可通过设计扩张状态观测器估计目标机动，从而设计考虑目标机动未知情况下的协同制导律。参考文献：1 姚禹正，余文斌，杨立军，等.多导弹协同制导技术综述J.飞航导弹，2021，41（6）：59-67.YAO Y Z，YU W B，Y ANG L J，et al.Overview of multimissile collaborative guidance technology J.Aviation Mis-sile，2021，41（6）：59-67.2 赵建博，杨树兴.多导弹协同制导研究综述 J.航空学报，2017，38（1）：13.ZHAO J B，YANG S X.Review of multi-missile cooperativeguidance J.Acta Aeronautica et Astronautica Sinica，2017，38（1）：13.3 孙新勇，吴杰，侯至丞，等.输入饱和约束下多智能体的有限时间一致 J.控制理论与应用，2020，37（11）：2391-2397.SUN X Y，WU J，HOU Z C，et al.Finite-time consensus ofmulti-agent with input saturation constraintsJ.ControlTheory&Applications，2020，37（11）：2391-2397.4 YANG B，JING W，GAO C.Three-dimensional cooperativeguidance law for multiple missiles with impact angle con-straint J.Journal of Systems Engineering and Electronics，2020，31（6）：1286-1296.5SONG J，SONG S，XU S.Three-dimensional cooperativeguidance law for multiple missiles with finite-time conver-gence J.Aerospace Science and Technology，2017，67：193-205.6 宋俊红，宋申民，徐胜利.一种拦截机动目标的多导弹协同制导律 J.宇航学报，2016，37（12）：1306-1314.SONG J H，SONG S M，XU S L.A Cooperative guidance lawfor multiple missiles to intercept maneuvering targetsJ.Journal of Astronautics，2016，37（12）：1306-1314.7 毛宁，范军芳，李斌.基于快速一致性理论的多导弹协同制导 J.兵器装备工程学报，2021，42（5）：227-234.MAO N，FAN J F，LI B.Multi-missile cooperative guidancebased on fast consistency theory J.Journal of OrdnanceEquipment Engineering，2021，42（5）：227-234.8 董晓飞，任章，池庆玺，等.有向拓扑条件下针对机动目标的分布式协同制导律设计J.航空学报，2020，41（S1）：116-124.DONG X F，REN Z，CHI Q X，et al.Distributed cooperativeguidance for maneuvering targets with directed communica-tion topologies J.Acta Aeronautica et Astronautica Sinica，2020，41（S1）：116-124.9PUKDEBOON C，SIRICHARUANUN P.Nonsingular ter-minal sliding mode based finite-time control for spacecraftattitude tracking J.International Journal of Control，Au-tomation&Systems，2014，12（3）：530-540.10 李磊，郭建国，周军，等.一种新的有限时间收敛的末制导律设计 J.火力与指挥控制，2015，40（12）：145-148.LI L，GUO J G，ZHOU J，et al.A new terminal guidancelaw design with finite time convergence J.Fire Control&Command Control，2015，40（12）：145-148.11WEI R，BEARD R W.Consensus seeking in multiagentsystems under dynamically changing interaction topologiesJ.IEEE Transactions on Automatic Control，2005，50（5）：655-661.901044