导弹动态可拦截区基于目标机动策略的计算.pdf

总第347期1引言随着现代战争模式逐渐向着高技术、信息化、远程控制方向发展，空中武器的使用愈加频繁，各国对武器的精确性和实时控制性有着越来越高的要求1。在面对敌方来袭目标飞行器时，空空导弹的拦截是防御系统的重要组成部分2，能有效提高自卫飞机或友军目标的生存能力，在导弹群协同一体化拦截对方多飞行器目标时，可以依照每个拦截弹的动态可拦截区实时在线地分配拦截目标，对集群的目标分配有重大意义3。因此，对空空导弹的拦截特性动态可拦截区的研究有着非常重要的意义4。计算可拦截区时，综合考虑到各种因素。第一，导弹自身的参数，例如导弹的结构参数、发动机特性、气动特性、导弹机动能力5。第二，导弹和目标之间关系，例如导弹与目标的相对位置、目标特性、目标的速度矢量、目标的机动能力6。第三，外收稿日期：2022年11月8日，修回日期：2022年12月17日作者简介：冯开，男，硕士，研究方向：飞行器动力学与控制。南英，男，博士，教授，研究方向：飞行力学与控制、飞行器设计与飞行仿真。王昕，男，硕士，研究方向：飞行力学与控制。晁泽睿，男，硕士，研究方向：飞行器任务规划。导弹动态可拦截区基于目标机动策略的计算冯开南英王昕晁泽睿（南京航空航天大学自然控制实验室南京211100）摘要为能在飞行器的攻防对抗中快速精确地拦截动态目标，需要对导弹的射后动态可拦截区进行求解。论文解算导弹可拦截区边界点时，在平移数值算法的基础上，提出一种导弹可拦截区的实时快速高精度拟合方法，将每次平移的数值基于脱靶量来确定，提高了仿真计算速度。考虑到实际战场上的随机干扰，针对目标突防机动的多种典型的飞行轨迹，都能进行动态、实时、有效的拦截，并可以实时动态地更新记录可拦截区的边界与导弹各项数据。仿真数据结果表明，该方法得到的导弹射后可拦截区与实际计算结果误差较小，具有很大的工程应用价值。关键词机动突防；导弹动态可拦截区；攻防对抗；飞行器拦截中图分类号E927DOI：10.3969/j.issn.1672-9730.2023.05.009Analysis of Missile Dynamic Interceptible Zone Based on theCalculation of Target Maneuver StrategyFENG KaiNAN YingWANG XinCHAO Zerui（Natural Control Laboratory，Nanjing University of Aeronautics and Astronautics，Nanjing211100）AbstractIn order to intercept dynamic targets quickly and accurately in the attack and defense confrontation of aircraft，it isnecessary to solve the dynamic interceptable zone of missile after firing.In this paper，a real-time，fast and high-precision fittingmethod of missile interceptable zone is proposed based on the translation numerical algorithm when solving the boundary points ofmissile interceptable zone.The value of each translation is determined based on the miss distance，which improves the simulationspeed.Considering the random interference in the actual battlefield，various typical flight trajectories of target penetration maneuvercan be intercepted dynamically，real-time and effectively，and the boundary of the interceptable area and various missile data canbe dynamically updated and recorded in real time.The simulation results show that the error between the interceptable zone obtainedby this method and the actual calculation results is small，which has great engineering application value.Key Wordsmobile penetration，missile dynamic interceptor zone，offensive and defensive confrontation，aircraft interceptionClass NumberE927舰 船 电 子 工 程Ship Electronic Engineering总第 347 期2023 年第 5 期Vol.43 No.540舰 船 电 子 工 程2023 年第 5 期部随机因素的干扰，例如，随机风场干扰、各系统中的随机误差和控制误差7。空空导弹的动态可拦截区是导弹在拦截过程中实时可拦截的范围2，与传统的射后可拦截区相比，具有实时性、高效性、准确性、适用性等多种优势8。传统算法计算可拦截区包括二分法或黄金比例法等9，运算量较大，周期较长，运算速度受到限制，很难实时满足战场瞬息万变的局势10。本文采用一种基于脱靶量的平移数值算法，大大减少了计算量，提高了计算速度11，可以做到快速并高精度地优化可拦截区的边界12。并且在拦截各种目标时，针对不同机动方式的目标，分别计算出其可拦截区的内外边界进行对比13。可以看出该方法计算可拦截区，具有实时、在线、快速等特点的同时14，对不同类型目标都可适用，具有很强的普适性15。2问题描述空空导弹的动态可拦截区是指空空拦截导弹发射后，沿自身弹道飞行过程中，在发射tM时间后，在诸多外部限制条件下，能够有概率命中目标的区域，即：Rmax=Rmax(hvvnTN1N2N3.)|t=tMRmin=Rmin(hvvnTN1N2N3.)|t=tMLmax=Lmax(hvvnTN1N2N3.)|t=tMDmax=Dmax(hvvnTN1N2N3.)|t=tM（1）其中，()RmaxRminLmaxDmax|t=tM分别表示空空导弹射后tM时刻可攻击区的外边界、内边界、左侧边界、右侧边界；()hvv|t=tM分别表示空空导弹沿着自身弹道飞行tM时间后的所在高度，飞行速度，轨迹偏角，轨迹倾角；(nT)|t=tM表示机动目标沿 着 自 身 弹 道 飞 行tM时 间 后 的 机 动 过 载；N1N2N3.为其他限制条件，如目标飞行器的机动类型、目标飞行器的最大机动过载、导弹最大可用过载、导弹推进系统、风场作用力的影响及系统中其他的干扰因素等。机动目标有其自身的突防对策，即机动过载nT随时间的变化策略。本文主要研究机动目标的突防对策对拦截该目标的空空导弹的动态可拦截区各边界()RmaxRminLmaxDmax的影响。3数学模型3.1拦截导弹的运动学与动力学模型对于导弹3自由度的运动方程，主要描述的是导弹的质心运动状态，主要考虑气动力、推力以及重力等。在讨论导弹绕质心的运动特性或对称姿态运动时，一些次要因素是可以忽略的，如地球自转、扁率、地球半径的变化等。给出导弹一般运动方程组如下：dxmdt=Vmcosmcosvm（2）dymdt=Vmsinm（3）dzmdt=-Vmcosmsinvm（4）式中：()xmymzm分别表示导弹位置的经度、高度、纬度；()Vmmvm分别表示导弹飞行时的速度、弹道倾角、弹道偏角。3.2制导律模型拦截导弹的制导律采用比例导引法，在速度坐标系中，利用式（5）（7）可以计算出拦截导弹在导弹速度坐标系的指令加速度：axv=n1|r qxv+kmaTxv（5）ayv=n2|r qyv+kmaTyv（6）azv=n3|r qzv+kmaTzv（7）式中，n1，n2，n3，km是比例系数；r是拦截弹与目标相对距离对时间的偏导数，即拦截弹与目标之间的相对速度；qxv，qyv，qzv是导弹导引头的转动速度投影在导弹速度坐标轴上的分量；aTxv，aTyv，aTzv是目标的飞行的加速度投影在导弹速度坐标轴上的分量；qxv，qyv，qzv是导弹导引头的转动速度投影在导弹速度坐标轴上的分量，拦截导弹在导弹速度坐标系的指令加速度可以通过矩阵转化为地面坐标系的指令加速度。3.3目标的数学模型VT=aTxT=aTyVT VT=aTzVTcosT（8）xT=VTcosTcosVT（9）hT=VTsinT（10）zT=-VTcosTsinVT（11）式中，下标T代表目标，aTx，aTy，aTz是目标的飞行的加速度投影在笛卡尔坐标轴上的分量。通过设定其随时间变化的函数，即可产生不同类型的机动突防飞行轨迹。4导弹动态可拦截区计算方法4.1确定可拦截区边界点根据空空拦截导弹和不同机动方式的目标的41总第347期数学模型与飞行状态，改进传统的平移数值算法，将每次平移的数值基于脱靶量来确定，对导弹可拦截区的可拦截内外边界进行初步解算，建立导弹可拦截区边界点的数据库。将拦截导弹与目标之间的距离初始值设置为明显大于导弹最大拦截距离或明显小于导弹最小拦截距离，保证可以让导弹脱靶。将空空导弹相对目标的接近速度、空空导弹的工作时间、过载、视线俯仰角、视线方位角等指标取最大或最小值，作为最大或最小可拦截区的约束条件。根据脱靶量判向内或向外平移目标飞行器，平移的距离依据脱靶量的大小设置，作为下次迭代计算的初始值位置，如此循环计算求出边界，计算流程图如图1所示。图1可拦截区计算流程4.2多种目标的拦截轨迹及攻击区4.2.1动态可拦截区的实时变化与传统的射前可拦截区相比，动态可拦截区能做到实时计算，实时变化，针对目标的不同运动状态，不同时刻，都能实时计算出可拦截区。选取迎向攻击直线运动目标为例，计算出拦截弹发射时刻与发射一段时间后的动态可拦截区，仿真结果如图2。其中红色与黑色线条分别表示拦截导弹与目标的运动轨迹，在其轨迹上取蓝色、橙色、黑色、品红色四点，分别为发射时刻、发射后5s、发射后20s、发射后40s的导弹与目标所在位置，并有颜色相同的可拦截区与之对应。仿真计算表明，可拦截区的边界在飞行拦截过程中不断发生变化，所以进行实时计算与更新是有必要的。20高度/km151050100-10-20Z轴/kmX轴/km1020303040500发射时刻可拦截区边界发射后5s可拦截区边界发射后20s可拦截区边界发射后40s可拦截区边界蓝色橙色黑色品红色品红色蓝色橙色黑色图2迎击目标直线运动动态可拦截区与轨迹4.2.2内外边界与拦截轨迹考虑到目标机动方式的多样性，选取几种典型的目标机动方式对其进行拦截，通过C+程序进行计算，得到动态可拦截区，取其发射时刻的内外边界画出图像，并选取其中一个拦截目标的算例，得到将拦截弹与目标的轨迹图像，利用Matlab画出图例如图38。20高度/km15105400-20-40Z轴/kmX轴/km10205030406007080目标轨迹导弹轨迹导弹发射点可拦截区内边界可拦截区外边界蓝色线条红色线条蓝色圆点品红色圆点黑色圆点红色线条黑色圆点品红色圆点蓝色圆点蓝色线条图3迎击目标直线运动可拦截区与轨迹10高度/km15105200-10-20Z轴/kmX轴/km510251520300354540红色线条黑色圆点品红色圆点蓝色圆点蓝色线条目标轨迹导弹轨迹导弹发射点可拦截区内边界可拦截区外边界蓝色线条红色线条蓝色圆点品红色圆点黑色圆点图4尾击目标直线运动可拦截区与轨迹如图3图8所示，蓝色圆点表示拦截导弹发射点，设置为（0，10km，0），蓝色线条表示目标飞行轨迹，红色线条表示拦截导弹飞行轨迹，品红色圆点表示计算得到的发射时刻可拦截区内边界，黑色圆点表示发射时刻可拦截区的外边界，选取的拦截轨冯开等：导弹动态可拦截区基于目标机动策略的计算42舰 船 电 子 工 程2023 年第 5 期迹算例是在外边界计算时获取的。20高度/km15105400-20-40Z轴/kmX轴/km10402030500608070目标轨迹导弹轨迹导弹发射点可拦截区内边界可拦截区外边界蓝色线条红色线条蓝色圆点品红色圆点黑色圆点红色线条黑色圆点品红色圆点蓝色圆点蓝色线条图5迎击目标U型运动可拦截区与轨迹5高度/km141210864100-5-10Z轴/kmX轴/km51025152030035目标轨迹导弹轨迹导弹发射点可拦截区内边界可拦截区外边界蓝色线条红色线条蓝色圆点品红色圆点黑色圆点红色线条黑色圆点品红色圆点蓝色圆点蓝色线条图6尾击目标U型运动可拦截区与轨迹10高度/km15105200-10-20Z轴/kmX轴/km51025152030035目标轨迹导弹轨迹导弹发射点可拦截区内边界可拦截区外边界蓝色线条红色线条蓝色圆点品红色圆点黑色圆点红色线条黑色圆点品红色圆点蓝色圆点蓝色线条图7目标sin型机动可拦截区与轨迹20高度/km15105400-20-40Z轴/kmX轴/km10402030500目标轨迹导弹轨迹导弹发射点可拦截区内边界可拦截区外边界蓝色线条红色线条蓝色圆点品红色圆点黑色圆点红色线条黑色圆点品红色圆点蓝色圆点蓝色线条图8目标螺旋机动可拦截区与轨迹4.2.3拦截不同机动目标时导弹数据对比如图9图12所示，分别表示导弹在射后至燃料耗尽时间内速度、加速度、轨迹倾角、轨迹偏角随时间的变化过程，由图可以看出，在追击不同机动目标时，导弹速度、加速度、轨迹倾角、轨迹偏角会根据目标的不同运动状态自适应地调节，以便更精确地追踪目标。图9导弹速度随时间变化曲线图10导弹加速度随时间变化曲线图11导弹轨迹倾角随时间变化曲线图12导弹轨迹偏角随时间变化曲线43总第347期4.2.4不同机动方式内外边界对比由图1315明显看出，在其他条件相同时，可拦截区的边界大小与目标的机动方式有关，当目标速度一定时，可拦截区由内而外，由小到大的排序依次为尾击直线运动、尾击U型机动、sin型机动、螺旋机动、迎击U型机动、迎击直线运动。这个现象对分析可拦截区的大小边界有着引导性作用。图13目标不同机动可拦截区X-Y平面对比图14目标不同机动可拦截区X-Z平面对比20高度/km1614121086400-20-40Z轴/kmX轴/km10402030500607080发射点迎击直线尾击直线迎击U型尾击U型sin型机场螺旋机动图15目标不同机动可拦截区三维对比5结语本文为了更高精度、更快速地计算导弹动态可拦截区，基于数值平移法的原理进行优化，提出了根据脱靶量快速迭代计算可拦截区的边界值，大大减少了计算量，计算结果与实际仿真解算的可拦截区边界点误差很小，且能够进行实时、快速、在线计算。并针对六种不同机动轨迹的目标进行了有效拦截，并对拦截各种机动目标的导弹轨迹以及数据进行分析比较。在实际作战过程中，可以首先判断突防飞行器的机动策略，做出相应的拦截策略，实时判断可拦截区与拦截时间。在作战过程中的导弹攻防对抗、群体导弹拦截一体化协同等方面，有很强的工程实际意义。参 考 文 献1周须峰，唐硕，高小五.空间目标可拦截区与可遭遇区的确定 J.航天控制，2006（01）：33-382胡朝晖，吕跃.反空空导弹弹道及拦截区解算研究 J.弹道学报，2020，32（04）：1-6.3葛鲁亲，孙旺，南英.对空导弹射后动态可攻击区的快速高精度拟合方法研究 J.航空兵器，2019，26（05）：41-47.4吴胜亮，南英.空空导弹射后动态可攻击区计算 J.弹箭与制导学报，2013，33（05）：49-54.5吴胜亮，南英，黄国强，等.随机风场与升阻比误差对空空导弹可攻击区的敏感度研究 J.弹道学报，2012，24（02）：25-30.6黄国强，南英，陆宇平.空空导弹擦肩发射全向三维可攻击区计算 J.弹道学报，2011，23（02）：37-41.7惠耀洛，南英，陈哨东，等.空空导弹动态攻击区的高精度快速算法研究 J.弹道学报，2015，27（02）：39-45.8李强，夏群利，崔莹莹，等.空地制导炸弹可攻击区的快速解算方法 J.兵工学报，2012，33（04）：390-394.9王凯，李望西，轩永波.基于BP神经网络的空空导弹攻击区解算及其仿真 J.弹箭与制导学报，2010，30（01）：75-77.10周灿辉，周德云，聂志强，等.中程空空导弹可攻击区解算的新方法J.火力与指挥控制，2012，37（07）：54-56，60.11王海涛，佟惠军，王洋.基于改进的BP神经网络空空导弹攻击区解算方法 J.电子设计工程，2014，22（03）：28-30，33.12孙旺.空空导弹对突防目标群的协同拦截研究 D.南京：南京航空航天大学，2020.13安彬.拦截临近空间飞行器制导与仿真研究 D.南京：南京航空航天大学，2017.14苗涛，杨毅，南英.导弹动态可攻击区实时在线算法研究 J.飞行力学，2018，36（02）：39-43.15钱杏芳，林瑞雄，赵亚男.导弹飞行力学 M.北京：北京理工大学出版社，2000.冯开等：导弹动态可拦截区基于目标机动策略的计算44