基于改进凸优化算法的多机编队突防航迹规划.pdf

第 卷第期 年月系统工程与电子技术 文章编号：（）网址：收稿日期：；修回日期：；网络优先出版日期：。网络优先出版地址：基金项目：国家自然科学基金（）资助课题通讯作者引用格式：刘玉杰，李樾，韩维，等基于改进凸优化算法的多机编队突防航迹规划系统工程与电子技术，（）：犚犲 犳 犲 狉 犲 狀 犮 犲犳 狅 狉犿犪 狋：，（）：基于改进凸优化算法的多机编队突防航迹规划刘玉杰，李樾，韩维，崔凯凯（海军航空大学航空基础学院，山东 烟台 ；海军招收飞行学员工作办公室，北京 ；中国人民解放军 部队，上海 ；中国人民解放军 部队，北京 ）摘要：为更好地发挥多机编队在低空突防作战中的优势，对已有的凸优化算法进行改进，提出一种多机编队低空突防航迹规划方法。首先，根据低空突防任务特点进行问题建模，包含多机编队航迹规划模型、障碍物模型以及任务分配评估模型。其次，提出基于匈牙利算法的集群任务分组概念，设计更符合战场需求和运动学规律的预规划航迹。之后，考虑新类型的外部障碍物和不同任务小组间的组间避障，通过对约束进行合理的近似与凸优化，实现多机安全飞行。最后，通过对比仿真，验证所提改进方法的可行性和有效性。结果表明，改进后的航迹规划算法能实现对低空突防任务的合理分配，求解效率和成功率均有所提高，对新加入的障碍有良好的避障效果。关键词：多机编队；低空突防；任务分配；航迹规划；凸优化；避障中图分类号：文献标志码：犇犗犐：犜 狉 犪 犼 犲 犮 狋 狅 狉 狔狆 犾 犪 狀 狀 犻 狀 犵犳 狅 狉狆 犲 狀 犲 狋 狉 犪 狋 犻 狅 狀狅 犳犿狌 犾 狋 犻 犪 犻 狉 犮 狉 犪 犳 狋犳 狅 狉犿犪 狋 犻 狅 狀犫 犪 狊 犲 犱狅 狀犻 犿狆 狉 狅 狏 犲 犱犮 狅 狀 狏 犲 狓狅 狆 狋 犻 犿 犻 狕 犪 狋 犻 狅 狀犪 犾 犵 狅 狉 犻 狋 犺犿 ，（犃狏 犻 犪 狋 犻 狅 狀犉狅 狌 狀犱犪 狋 犻 狅 狀犆狅 犾 犾 犲 犵 犲，犖犪 狏 犪 犾犃犲 狉 狅 狀 犪 狌 狋 犻 犮 犪 犾犝狀 犻 狏 犲 狉 狊 犻 狋 狔，犢犪 狀 狋 犪 犻 ，犆犺 犻 狀 犪；犖犪 狏 狔犚犲 犮 狉 狌 犻 狋 犻 狀犵犗犳犳 犻 犮 犲犳 狅 狉犛 狋 狌犱 犲 狀 狋犘 犻 犾 狅 狋 狊，犅犲 犻 犼 犻 狀犵 ，犆犺 犻 狀 犪；犝狀 犻 狋 狅 犳狋 犺 犲犘犔犃，犛犺 犪 狀犵犺 犪 犻 ，犆犺 犻 狀 犪；犝狀 犻 狋 狅 犳狋 犺 犲犘犔犃，犅犲 犻 犼 犻 狀犵 ，犆犺 犻 狀 犪）犃犫 狊 狋 狉 犪 犮 狋：，犓犲 狔狑狅 狉 犱 狊：；系统工程与电子技术第 卷引言现代化战争正在向全方位、立体化、大纵深的多兵种合成作战方向发展，防空技术的快速进步使得对空战的要求不断提升。应运而生的低空突防技术是多机编队作战以及网络化作战的重要技术，有效的突防部署可使多机编队更有效地躲避威胁障碍，提高任务完成效率，保证飞行安全。面对复杂多变的低空战场环境，只依赖人工规划航迹或飞行员实时操纵已无法满足现代战场多变的要求，需要综合考虑自然环境、电磁防御、炮火防御等多种威胁以及飞机本身的性能指标，并按照科学有效的航迹规划算法，规划出合理的飞行轨迹，从而提高飞机的作战效能。因而，多机编队低空突防的重要基础是确定航迹规划方法，规划方法的优劣直接决定了低空突防任务完成的效果。航迹规划根据建模方式的不同，通常分成两类，一类是基于转弯角、航段长度而建立代数运动方程，包括势场法、图论法、随机树法 等，此类建模方式适合大尺度的全局规划问题。另一类则是以过载或加速度等控制量为基础，建立微分运动方程模型，此类建模方式更贴近工程实际，与控制理论也更为契合。本文采取后一类方式进行建模，该类模型的求解方式又可分为直接法和间接法。前者的原理是将最优控制问题转化为有限维参数优化问题，用非线性规划方法进行求解，其缺点是求解精度不高；间接法包括伪谱法、混合整数规划算法 以及启发式算法 等，其原理是将航迹规划问题视为最优控制问题进行求解，存在不易收敛的不足。当前，作为直接法的典型代表之一，凸优化方法已成为求解航迹规划问题的一种高效稳定算法，这与凸优化理论的迅速发展有密切关系。但是，凸优化建模方法需要基于线性动力学假设，并受限于建模对象的不同特点，因而普适性一般。对此，序列凸优化方法采取逐次凸化的方式，将原模型中的非凸问题分解为若干凸优化子问题，从而达到求解一般非凸约束航迹规划问题的目的，并已在多机编队避障、卫星自主交会、行星着陆 等航迹规划案例中得到校验。特别地，文献 针对异构多机系统中各飞机种类的任务特点，设计了不同的优化目标函数，并通过凸优化近似各类约束条件，为多机编队协同航迹规划问题提供了一种新的求解思路，特别是证明了凸优化算法在解决高维协同航迹规划问题上比伪谱法等经典算法有更大的优势。研究发现，将凸优化算法直接应用到多机编队低空突防航迹规划过程中，还应对以下几个方面进行改进和完善：初始规划航迹的改进。在序列凸优化迭代求解的过程中，初始解即首次输入的初始航迹采用单纯的初始点、目标点连线方式确定，这不符合飞行的运动学规律，虽不影响求解，但导致迭代收敛缓慢，同时，并未结合任务需求进行合理的初始任务分配，“协同作战”概念不足，含“战”量低。障碍模型及其近似凸化过程的完善。在已有文献中，对障碍物过于简化，忽略了其所代表的实际意义，仅考虑了所谓的圆形障碍物，未考虑其他构型的障碍，更没有讨论其他构型障碍的近似与凸优化过程。多机执行任务时的飞行安全问题需进一步考虑。由于低空突防的飞机数量更多，子任务不尽相同，应考虑建立“任务小组”的概念，在完成任务分配的基础上，重点研究航迹规划中不同任务小组之间的避障问题，以进一步保证多机飞行安全。基于上述分析，本文以固定翼飞机多机编队为研究对象，以低空突防为任务输入，对已有的凸优化算法进行改进，提出一种基于改进凸优化算法的多机编队低空突防航迹规划方法。本文的创新点包括以下点：首先，提出了基于匈牙利算法的集群任务分组概念，设计了更符合运动学规律、更贴近战场实际需求的预规划航迹；其次，在凸优化过程中加入了更贴合战场环境的多边形障碍，用以模拟禁飞区，并解决了航迹规划离散点之间可能与多边形障碍发生碰撞的问题；最后，针对不同任务小组之间可能发生的碰撞问题，引入“虚拟多边形”的概念以进行凸优化处理，以进一步确保各任务小组中飞机的飞行安全。问题建模 多机编队航迹规划模型的建立多机编队航迹规划模型是基于飞行器运动模型建立的，设定多机编队在执行任务过程中无侧滑飞行，且发动机推力与飞行速度方向一致，从而可建立单机三维空间点质量运动模型，如下所示：狓犻犞犻 犻 犻狔犻犞犻 犻 犻犺犻犞犻 犻犞犻犵（狀狓，犻 犻）犻犵狀狔，犻犞犻 犻犻犵（狀狕，犻 犻）犞烅烄烆犻（）式中：（狓犻，狔犻，犺犻）表示飞机犻在地面坐标系下的三维坐标；（犻，犻）表示飞机犻的航迹倾角和航向角；犞犻表示飞机犻的地速；（狀狓，犻，狀狔，犻，狀狕，犻）分别表示飞机犻的水平轴向过载、法向过载的水平分量和法向过载的垂直分量；狊犻和狌犻分别表示飞机犻的状态量及控制量，其中，狊犻（狓犻，狔犻，犺犻，犞犻，犻，犻），狌犻（狀狓，犻，狀狔，犻，狀狕，犻）。由于本文所研究的多机为同一型号，因而编队中各机的运动模型一致，用下标犻进行区分。在完成运动模型建立的基础上，多机航迹规划问题还需要考虑若干约束条件，具体包含飞机初末状态量约束、状态控制量的幅值约束、控制量的变化率约束、机间防碰第期刘玉杰等：基于改进凸优化算法的多机编队突防航迹规划 撞约束和与外部障碍物防碰撞约束等，上述约束的数学表达式分别如下所示：狊犻（狋）狊犻，狊犻（狋犳）狊犻，烅烄烆犳（）狊犻，狊犻（狋）狊犻，狌犻，狌犻（狋）狌犻，烅烄烆，狋狋狋犳（）狌犻，狌犻（狋）狌犻，（）犉狊犻（狋）狊犼（狋）犚狋，犻犼；狋狋狋犳（）犈狊犻（狋）狆 ，犿狉犿，犿，犕；狋狋狋犳（）式中：狋和狋犳分别表示航迹规划的初始时刻和终端时刻；矩阵犉犐，犚狋是机间安全距离；犈犐，狆 ，犿和狉犿分别是外部障碍的中心坐标和对应的安全距离；犕为静态威胁个数。约束条件确定后，需设定航迹规划模的目标函数，结合低空突防任务快速、精确打击的特点要求，设定以最短时间到达目标点为目标函数，具体如下所示：狊犻（狋），狌犻（狋）狋犳狋狋 式（）式（）（）障碍物模型的建立低空突防过程中，在空域内往往会遇到雷达、高炮阵及禁飞区等障碍，在对上述障碍进行建模时应既考虑其共性，又考虑其特性，同时也还要进行合理的简化和提炼，降低计算复杂度。值得一提的是，虽然本文所研究的低空突防航迹规划问题是基于三维空间的，但借鉴之前凸优化求解的理论成果，将障碍物视为高度不限的三维柱体，因此在障碍物建模时只需考虑二维截面。其中，敌方雷达和高炮阵地可以近似为圆形障碍，以雷达为例，在突防过程中，雷达的覆盖面积可用以雷达所处位置为圆心、探测范围为半径的圆形表示。图给出了飞机与雷达位置关系的示意图，当满足犇犚狉犱狊时，可认为飞机实现了避障，图中犚狉为雷达探测半径，犱狊为安全裕度。图飞机与雷达位置关系图 低空突防过程中所遇到的禁飞区多为不规则形状，将不规则图形简化为可用数学表达式表达的规则图形，当满足下式时，飞机实现了对禁飞区的避障：狓犻狓狅犪犚犱（）狊犽狔犻狔狅犫犚犱（）狊犽（）式中：（狓狅，狔狅）为禁飞区中心坐标；（犪，犫）为禁飞区的宽度与长度；犽表示禁飞区所近似多边形的形状参数。根据式（），当犽时，禁飞区表示菱形；当犽时，禁飞区表示椭圆形；当犽时，禁飞区表示矩形。任务分配评估模型的建立任务分配的合理性直接关系到我方突防效率及最终空战结果，合理的任务分配可使多机编队在最小的损耗代价下获得对敌方最大的空战打击优势。当今，比较主流的任务分配方法包括启发式方法、最优化方法和基于市场交易思想的拍卖算法等。其中，启发式算法求解的效率较高，包括遗传算法、鱼群算法等。最优化方法包含匈牙利算法、动态规划方法等。这类算法应用更为广泛，但处理某些特殊的数据时易出现不收敛的情况。拍卖算法 假设任务为拍卖品，任务分配对象根据出价策略与收益函数竞拍或拍卖某一任务，通过任务分配对象竞价的方式实现任务分配。低空突防任务中包含近距轰炸、远程打击、电子战压制、空中加油等多个子任务。匈牙利算法是目前公认的一种有效解决任务分配的方法，其适合求解“一对一”的匹配问题，在解决小规模任务分配问题时，具有求解迅速、可靠性高的特点。本文所研究的对象是不具备多目标攻击能力的三代机，因而在低空突防任务中，不能采取一架飞机对多个任务的分配模式；另一方面，由于我方出动的战机数量通常大于任务数量，因而可以根据任务数量对我方战机进行对应数量的分组，形成若干“任务小组”，将“多对一”匹配问题转化为“一对一”的匹配问题，从而可采用匈牙利算法进行计算和求解。由于匈牙利算法的原理较为经典，受篇幅限制，在此不再详细展开介绍，该算法的关键是要输入一个科学有效的评估矩阵。为此，本节将重点放在对低空突防过程中子任务的评估建模上，以更好地提高任务分配的合理性，从而进一步保证多机航迹规划中初始航迹输入的有效性。本文在对子任务评估建模时，提出了创新的评价标准，重点考虑了两个方面，即子任务的威胁和子任务的固有价值，前者又由敌我相对态势指标和敌我相对战斗力指标共同决定。由于引入了“任务小组”的概念，因而在评估建模时，需假设小组内各飞机正以同速、同航向巡航飞行，其机间距离相对子任务与任务小组距离可忽略不计，即在计算相对态势时，将任务小组视为质点，取小组中长机状态量进行计算。其中，长机是任务小组中的某一架飞机，通常由领导者驾驶，负责整个任务小组的指挥、下达命令等，是预先设定好的。同时，由于执行子任务的位置通常距离子任务实际位置仍有一段距离，如远程打击子任务并不需要实际达到目标点上空，因而还需提前设定实际的任务执行点，同样只取长机状态量即可。系统工程与电子技术第 卷首先，对敌我相对态势指标犜进行建模，犜由相对距离态势参数犜狉、相对高度态势参数犜犺及相对角度态势参数犜犪这个参数决定，其对应的数学表达式如下所示：犜狉，狉犚犪；狉犚犫 狉犚犫犚犪犚犫，犚犫狉犚犪，犚犫狉犚犪，（犚犪，犚犫）狉犚犮，烅烄烆其他（）犜犺，犺犺犺犺犺，犺犺犺 ，犺犺烅烄烆（）犜犪 狇犫狇犪（）式中：狉为任务小组长机距实际任务执行点的距离；犚犫为子任务目标的防御攻击距离；犚犪为任务小组中所携导弹的最大射程；犚犮为任务小组雷达最大跟踪距离；犺为任务小组长机相对实际任务执行点的高度差；犺和犺 均为高度经验值；狇犪表示任务小组中长机的速度矢量与长机的质心和实际任务执行点的质心之间连线的夹角；狇犫表示实际任务执行点的速度矢量与长机的质心和实际任务执行点的质心之间连线的夹角。上述部分物理量的关系如图所示。图敌我相对态势指标中各物理量示意图 个参数确定后，可得敌我相对态势指标犜：犜犜犪犜狉犜犺（）再来分析敌我相对战斗力指标犆，通俗来讲，犆是指我方任务小组和敌方战斗力的差。对战斗力的评估，可以参考文献 所提的评估公式，其受个要素的影响，分别是机动能力、打击能力、侦查能力、操纵性能、生存指标、飞行距离以及电子对抗能力，则本文中犆的表达式为犆犆犫犆犪 犅犫 犃犫（）犃（）犫犫犫犫犫 犅犪 犃犪（）犃（）犪犪犪犪犪（）式中：犆犪、犆犫分别为我方任务小组和子任务目标的战斗力指标；犅为机动能力参数；犃为打击能力参数；犃为侦查能力参数；表示操纵性能参数；表示生存指标参数；表示飞行距离参数；表示电子对抗能力参数。本文中，对式（）进行进一步拓展，引入飞行员综合素质参数，其值设定如下：，优秀飞行员，良好飞行员，烅烄烆合格飞行员（）从而有：犆犆犫犫犆犪犪（）通常，可根据经验数据库提前确定犆犪、犆犫，在后续的仿真中，在此也不再详细展开的具体计算过程。最后，子任务的固有价值犌可由我方指挥官利用打分原理确定，通常犜、犆、犌都需要进行归一化处理，得到犜、犆、犌。其中，归一化方法采取较为简单的处理方式，即分别除以各自指标的最大值即可。至此，可得到子任务小组犻相对任务犼的子任务的评估值犠犻 犼。犠犻 犼犌犻 犼（犜犻 犼犆犻 犼）（）式中，、为比例系数。在此基础上得到评估矩阵犠，可表示为犠犠 犠 犠犼犠 犠 犠犼犠犻犠犻犠犻熿燀燄燅犼（）凸优化算法的改进首先对式（）式（）进行离散化，将最优控制问题转化为参数优化问题；之后基于基准轨迹对约束进行近似和凸化，其中外部障碍约束仅考虑了圆形障碍；最后利用序列凸优化迭代的思想对模型进行求解。本文将重点在预规划航迹、外部障碍处理和任务小组避障约束等方面对凸优化算法进行改进。预规划航迹的改进对预规划航迹的改进主要体现以下两个方面：一方面是利用第 节构建的子任务评估模型，结合匈牙利算法，对各任务小组进行任务预分配，得到大尺度规模下任务小组的预规划航迹。由于航迹规划过程中包含末状态量约束，因而各任务小组还应根据所分配的子任务特点，确定末状态编队队形，合理的队形选择能够有效提高子任务的完成率。以三机以上的任务小组为例，当末状态编队构型为横向平行队形时，适和完成大面积侦查和搜索任务；为一字纵向队形时，适用于完成近距轰炸任务，或者伴随掩护任务；为楔形队形时，适和完成巡航或空投任务；为菱形护卫队形时，适用于保护某架友机，该飞机可能出现油量不足或损伤的突发情况。第期刘玉杰等：基于改进凸优化算法的多机编队突防航迹规划 另一方面，在既往研究所述的序列凸优化迭代求解过程中，首次输入的基准航迹采用简单的初始点、目标点连线方式进行确定，而初始离散时刻的航迹点则根据离散时刻数进行平均分配，这显然不符合飞行的运动学规律，虽不影响求解，但易导致迭代收敛缓慢，甚至无解的情况发生。为此，在任务分配完成的基础上，本文对组内各机的预规划航迹进行改进，利用初始序列凸优化迭代的方式求取预规划航迹。在这一过程中，重点保证的是飞机的运动特性和求解效率，因而只需考虑飞机的初末状态量约束、状态控制量的幅值约束以及控制量的变化率约束，而机间防碰撞约束和与外部障碍物防碰撞约束则暂时忽略，初始序列凸优化迭代的优化模型如下：狊犻狊犻，狊犻犓狊犻，烅烄烆犳，犻，犖（）狊 狊犻犽狊，犽，犓；犻，犖（）犝 狌犻犽犝，犽，犓；犻，犖（）犝 狌犻犽犝，犽，犓；犻，犖（）狊犻（犽），狌犻（犽），狋，犻，犖；犽，犓狋 犃 犣犫 （犃 犣犫，）狊狇犻犽狊狇犻犽狇，犽，犓（）式中：犓为离散时刻数，是指将初始时刻到终端时刻这段时间离散成的时间段数量；犖为任务小组中飞机的数量；狇为迭代次数，为收敛域，狇每增加，将变成原来的。通过这种策略使得循环逐渐趋向收敛；等式约束犃 犣犫 包括式（），不等式约束犃 犣犫 包括式（）式（），和 为上述两项的比例系数；犣包括离散时间步长狋和集合狕犻，其中：狕犻（狊犻，狊犻 狊犻犓，狌犻，狌犻 狌犻犓），犻，犖（）至此，可得到小尺度规模下任务小组的预规划航迹，从而完成了对预规划航迹的改进。外部障碍的增加及凸优化既往研究仅考虑了所谓的圆形障碍物并忽略了其所代表的实际意义。本文在第 节建立了基于低空突防任务的环境障碍物模型，包括雷达、高炮阵及禁飞区等，除可近似为圆形障碍物的雷达、高炮阵外，还有可近似为多边形障碍的禁飞区。本节对禁飞区的避障进行研究，以式（）中犽为例进行讨论，此时禁飞区表示矩形。设我方通过侦察获取到一处可近似为矩形的禁飞区，用 表示，其个顶点从最西顶点逆时针标记为（狓犻，狔犻），其中犻，式（）中的约束条件能保证飞机不与禁飞区相撞，但该约束显然是非凸的，无法直接用于序列凸优化求解。为此，借助基准航迹，即上一次迭代所求的航迹对该约束进行凸化近似。首先，取某一时刻基准点狆犽与 形心的连线，若该连线与 的某一条边有交点，则置预警信号 ，并定义该边为基准边，若该连线与 的任何一条边都没有交点，则置预警信号 为一恒正值，同时计算狆犽到形心的距离，记为犱。其次，根据基准边两顶点的坐标，可求得该边所在直线的标准表达式，该直线划分平面为两部分，显然两部分均为凸的，其中基准点所在的半平面为该次循环狆犽的可行域。最后，将该约束作为序列凸优化模型中的一项不等式约束进行惩罚，用犆 表示，其表达式如式下 （犪狆狓狆犽犫狆狔狆犽犮狆）（犪狆狓狆犽犫狆狔狆犽犮狆）（，）（犪狆犪狆犫狆犫狆 槡）犆 犱 烅烄烆（）式中：（狓狆犽，狔狆犽）为狆犽的坐标；（狓狆犽，狔狆犽）为狆犽的坐标；（犪狆，犫狆，犮狆）为基准边所在直线的标准表达式参数。特别是犆 表达式中 犱 的设计有利于更好地引导收敛。（犪狆狓狆犽犫狆狔狆犽犮狆）（犪狆狓狆犽犫狆狔狆犽犮狆）；（，）（犪狆犪狆犫狆犫狆 槡）犆 犱 烅烄烆（）式中：（犪狆，犫狆，犮狆）为狆犽与 形心所决定的基准边所在直线标准表达式系数；犱 为狆犽到形心的距离，同理有 。在保证离散点处的避障后，还应保证航迹离散点之间的航迹不与 发生碰撞，须在式（）的基础上确保式（）中的犆 表达式成立，可保证离散点之间的飞行安全。式（）和式（）联立得到狆犽的可行域如图所示的深色区域，显然两凸集的交集还是凸的，这符合序列凸优化求解的原则。图多边形障碍的凸化可行域 系统工程与电子技术第 卷 任务小组间的避障约束及凸化在多机编队低空突防任务中，划分若干任务小组，若继续按照既往研究 所述的小组内飞机间避障约束进行优化，存在以下问题。一是飞机数量多导致的求解效率降低。由于已有的机间避障约束是令每一架飞机躲避其他所有飞机，飞机数量越多，这一约束的数量和复杂程度将越大，不利于发挥序列凸优化算法求解迅速的优势。二是各任务小组的任务不同，执行任务时间相对灵活，原机间约束的适用性降低。之前的理论成果多以多机同时到达目标点为基础条件展开，当面对更为复杂的低空突防任务时，由于子任务的评价值不同，其被执行的时序也不尽相同，因而多机同时到达目标点的理想设定将难以保证。为此，提出“虚拟多边形”的概念来解决任务小组间的避障问题。以两个任务小组为例，设在收到指挥官任务指令后，优先级更高的任务小组率先按照规划航迹进行飞行，所规划的到达子任务目标点的时间为犜狆，小组内飞机间避障按照传统的方式进行避障，其约束如下所示：（狊犻犽狊犼犽）犈犈（狊犻犽狊犼犽）犚犈（狊犻犽狊犼犽）狋（犞犻，犞犼，）犚犚烅烄烆狋，犻，犖；犼，犖；犻犼；犽，犓；（）式中：犞犻，表示飞机犻的最大飞行速度；犞犼，表示飞机犼的最大飞行速度；犚表示平面内的安全距离。在任务小组飞行至犜时刻时，任务小组开始执行任务，此时对任务小组进行航迹规划时需考虑任务小组的航迹。截取任务小组在犜犜狆时段的规划航迹，标记出这段航迹在地面坐标系下的横纵坐标最值，并确定个最值点，进而依据最值点确定一个矩形，称该矩形为组间避障的虚拟多边形，该过程的示意图如图所示。图组间避障示意图 在对任务小组进行航迹规划时，应将该虚拟多边形按照第 节所述的外部障碍进行躲避。虽然该近似过程具有一定的保守性，但这是对经典凸优化算法思想的一个延续，该步操作有利于实现多机编队简单、快速、安全的组间避障的效果，并可拓展到多任务小组协同避障的案例中。仿真校验本节对所提算法的有效性进行校验。设定仿真背景如下：经预警侦查，获取到敌方总部位于我东北方向约 处，同时侦察到重要的电磁干扰区和战略性建筑各一处，此外还有若干雷达区和禁飞区。将打击电磁干扰第期刘玉杰等：基于改进凸优化算法的多机编队突防航迹规划 区、敌方总部、战略性建筑设为个子任务，分别记为子任务、子任务和子任务，将我方 架战机分为个任务小组，分别记为任务小组、任务小组和任务小组，每组架飞机，敌我双方的初始状态如表和表所示。其中，表中的狊（狋）第一行为组内长机的参数。由于执行的是低空突防任务，因而任务执行点的目标队形高度均不超过 。优化目标均设定为多机到达时间最短，并选取一字纵向队形和楔形队形等适合攻击的队形为主要队形。障碍物模型有关信息如表所示，优化参数设置如表所示。仿真环境为 ，计算机配置为 酷睿 处理器，主频 ，内存，位 操作系统。表各任务小组的状态量犜 犪 犫 犾 犲犛 狋 犪 狋 狌 狊狊 犮 犪 犾 犲狅 犳犲 犪 犮 犺狋 犪 狊 犽犵 狉 狅 狌 狆物理量任务小组任务小组任务小组狊（狋）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）初始队形菱形横向平行队形菱形狊 （，）（，）（，）狊（，）（，）（，）狌 （，）（，）（，）狌（，）（，）（，）狌 （，）（，）（，）狌（，）（，）（，）犚犪 犚犮 表各子任务信息犜 犪 犫 犾 犲犐 狀 犳 狅 狉犿犪 狋 犻 狅 狀狅 犳犲 犪 犮 犺狊 狌 犫 狋 犪 狊 犽物理量任务小组任务小组任务小组真实位置（，）（，）（，）任务执行点（，）（，）（，）狊（狋犳）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）任务执行点队形楔形楔形一字纵向队形犚犫 犌 表障碍物模型信息犜 犪 犫 犾 犲犐 狀 犳 狅 狉犿犪 狋 犻 狅 狀狅 犳狅 犫 狊 狋 犪 犮 犾 犲犿狅 犱 犲 犾物理量雷达雷达禁飞区禁飞区中心坐标（，）（，）犚狉犱狊 边界坐标（，）（，）（，）（，）（，）（，）（，）（，）表优化参数设置犜 犪 犫 犾 犲犗狆 狋 犻 犿 犻 狕 犪 狋 犻 狅 狀狆 犪 狉 犪犿犲 狋 犲 狉狊 犲 狋 狋 犻 狀 犵阶段物理量数值预规划阶段犺 犿犪狓 续表犆狅 狀 狋 犻 狀 狌 犲 犱犜 犪 犫 犾 犲阶段物理量数值航迹规划阶段犓 犿犪狓 犳 犾 犪犵 （，）系统工程与电子技术第 卷 预规划航迹的校验根据表中长机的狊（狋）以及犚犪、犚犮和表中任务执行点的状态量及犚犫，结合表设定的预规划阶段参数，利用式（）式（）可解算出各子任务相对各任务小组的态势指标犜。受篇幅限制，不再列出犜狉、犜犺及犜犪具体值，而是直接在表中列出了求解后的犜值，同时在该表格中列出了根据文献 估算的犆值以及根据经验打分获取的犌值。根据表中数据，通过式（）和表中的参数，可计算得到评估矩阵犠，具体如下：犠 熿燀燄燅 （）将犠带入匈牙利算法进行求解犠狉，得出结果矩阵：犠狉 熿燀燄燅 （）对应的最优值为 。从而可以得到初始任务分配结果，即任务小组执行子任务，任务小组执行子任务，任务小组执行子任务。表任务小组对子任务的评价值犜 犪 犫 犾 犲犈狏 犪 犾 狌 犪 狋 犻 狅 狀狏 犪 犾 狌 犲狅 犳狋 犺 犲狊 狌 犫 狋 犪 狊 犽犫 狔狋 犺 犲狋 犪 狊 犽狋 犲 犪犿分组子任务犜犆犌子任务犜犆犌子任务犜犆犌任务小组 任务小组 任务小组 为进一步验证小尺度范围内预规划航迹的改进效果，暂时不考虑本文所提的组间避障约束和多边形外部障碍约束，即以文献 的算法结构为基础，采用已完成任务分配的个任务小组为实验对象，分别利用初始序列凸优化迭代的方式，采用简单的初始点、目标点连线方式确定初始航迹。其中，前者为实验组，后者为对照组，分别进行 次完整的航迹规划对照试验以验证改进效果。航迹规划的已知条件和约束均已在表表中列出，实验结果如表所示。表预规划航迹对比表犜 犪 犫 犾 犲犆狅犿狆 犪 狉 犻 狊 狅 狀狅 犳狆 狉 犲 狆 犾 犪 狀 狀 犲 犱狋 狉 犪 犼 犲 犮 狋 狅 狉 狔分组平均求解时间实验组对照组求解成功率实验组对照组任务小组 任务小组 任务小组 从表中可以看出，实验组的平均求解时间在任务小组和任务小组两组规划过程中用时更短，最大平均差异为 ；任务小组中两组用时接近，出现这种情况的原因是采用两种方式对任务小组进行的预规划航迹均近似接近直线，复杂度较低，差异化较小，具体如图所示，因而后续输入到完整约束的航迹规划模型中差异不大。而另外两组航迹因相对复杂，预规划航迹差异较大而使得完整的求解时间有差异，由此可以得出改进后的预规划航迹使得算法的求解效率有一定的提升。第期刘玉杰等：基于改进凸优化算法的多机编队突防航迹规划 图任务小组预规划航迹 另一方面，实验组的求解成功率也相比对照组更优，且同样有着航迹越复杂，优势越明显的特点，在任务小组的 次航迹规划求解中，实验组求解成功率高出对照组，证明了本文所提改进策略的优越性。航迹避障效果的校验在第 节任务分配的基础上，验证航迹避障效果。各任务小组的出发时序由后方的指挥中心根据战场的实际作战需求和战场态势决定，先出发的小组优先级更高，后出发的任务小组需对其进行避障。由于电磁压制任务（子任务）是保障其他两项子任务的基础，因而先出动任务小组。任务小组的航迹规划结果如图所示，图中不同颜色的线条代表组内各飞机的航迹。结合表，从图中可看出组内各飞机均可到达实际任务执行点，并形成指定队形，所求解出的狋 。同时，从图中可以直观地看到，各飞机均有效地躲避了雷达和禁飞区，验证了本文所提多边形障碍约束的有效性。图任务小组执行子任务航迹变化图 当任务小组按规划航迹飞行至 时，出动任务小组执行打击战略建筑的任务，以实现对敌方进一步的压制。由于该时刻任务小组还未抵达实际任务执行点，因而对任务小组进行航迹规划时，需考虑组间避障。通过计算可知，此时任务小组飞至第 个离散时刻点，截取第 离散时间段的航迹，按照第 节所述方法，确定虚拟多边形，其个顶点坐标分别为（，），（，），（，），（，），将其加入任务小组的航迹规划模型进行求解，可得出相应航迹。为校验组间约束的效果，去除组间约束进行对比实验，同时保留其他约束和参数设定不变，再次进行规划，两次规划的结果如图所示。系统工程与电子技术第 卷图任务小组航迹对比 从图（）和图（）两图中，可以看出添加组间约束后，任务小组为避免与任务小组可能发生的碰撞，组内飞机均采取向禁飞区西侧绕飞的方式抵达实际任务执行点，而未添加组间约束时，任务小组内架飞机，各有架分别从东、西两侧对禁飞区进行绕飞。结合图（）和图（）两图可看出，两种方式均可到达指定任务执行点，并且组内最小机间距均大于设定的安全距离，但结合图（）可知，在任务小组出动前期，任务小组正从西向东飞行，未添加组间避障的规划航迹此时从禁飞区东侧由南向北进行穿越。经计算，在同一时刻两小组最小机间距仅为 ，小于允许的最小安全距离，而添加机间避障约束的最小机间距显然始终远大于安全距离，可以证明本文所提的组间约束的可行性和有效性。当任务小组出发后的第 ，出动任务小组打击敌方总部。由于任务小组按照添加组间约束进行规划所求出的狋 ，而任务小组的狋 ，通过简单计算可知，此时任务小组已抵达实际任务执行点，任务小组仍未抵达对应的任务执行点，因而需要考虑任务小组与任务小组的组间避障。按照上一轮仿真的计算方式，截取任务小组第 离散时间段内的航迹，形成虚拟多边形，其个顶点坐标分别为（，），（，），（，），（，），则任务小组的规划航迹如图所示。由图可以看出，组内飞机按照各约束条件，安全抵达了任务执行点。图任务小组执行子任务航迹变化 第期刘玉杰等：基于改进凸优化算法的多机编队突防航迹规划 综上所述，采用本文所提算法使飞机能有效躲避多边形障碍，同时可以通过躲避虚拟多边形的方式实现有效的组间避障。此外，在执行子任务的过程中，各机严格满足始末状态量约束、状态控制量的幅值约束、控制量的变化率约束、机间防碰撞约束，可有效躲避雷达等圆形障碍物。考虑到序列凸优化算法在求解该航迹规划模型时较其他传统算法在效率上也具有明显优势，可充分证明本文算法能有效满足多机低空突防任务航迹规划需求。结束语本文对已有的凸优化算法进行改进，研究了多机编队低空突防航迹规划问题，得到如下结论：（）提出了基于匈牙利算法的集群任务分组概念，并设计了更符合运动学规律、更贴近战场实际需求的预规划航迹，仿真证明改进后的预规划航迹有助于提升模型求解效率和成功率；（）在凸优化过程中引入更为贴合战场环境的多边形障碍，以提高模拟禁飞区的真实性，并有效解决航迹规划离散点之间可能与多边形障碍发生碰撞的潜在风险；（）针对不同任务小组之间存在的碰撞风险，引入“虚拟多边形”概念进行凸化处理，可进一步提升各任务小组间的飞行安全性。后续将针对如何进一步完善任务分配模型、降低凸优化算法的保守性等问题进行更为深入的探索研究。对于突防任务而言，还应加强对动态的、无法事先确定的威胁的研究，并将本文所提算法作进一步延伸拓展。此外，为有效应对当前人工智能高速发展趋势，还应尽早针对有人无人协同作战航迹规划和无人机机群作战航迹规划特点将该算法作适应性改进研究。参考文献闵昌万飞行器航迹规划与轨迹控制研究西安：西北工业大学，：，：黄鹤，吴琨，王会峰，等基于改进飞蛾扑火算法的无人机低空突防路径规划中国惯性技术学报，（）：，（）：张志伟，蒋道刚，袁坤低空突防航线规划研究飞行力学，（）：，（）：，：，（）：丁明跃，郑昌文，周成平，等无人飞行器航迹规划北京：电子工业出版社，：，：，：韩尧，李少华基于改进人工势场法的无人机航迹规划系统工程与电子技术，（）：，（）：，：，：，：，：，：，：，（）：，（）：，（）：王祝，刘莉，龙腾，等基于罚函数序列凸规划的多无人机轨迹规划航空学报，（）：，（）：严冰，张进，罗亚中面向编队卫星的空间系绳在轨服务系统工程与电子技术，（）：，（）：系统工程与电子技术第 卷，（）：李牧东，赵辉，黄汉桥，等基于模型的三维实时平滑航迹规划系统工程与电子技术，（）：，（）：，：，李樾，韩维，陈清阳，等凸优化算法在有人无人机协同系统航迹规划中的应用宇航学报，（）：，（）：，（）：，：，邵将，徐扬，罗德林无人机多机协同对抗决策研究信息与控制，（）：，（）：，（）：杨尚君，孙永维，庞宇基于改进鱼群算法的多无人机任务分配研究计算机仿真，（）：，（）：，（）：，：张进，郭浩，陈统基于可适应匈牙利算法的武器 目标分配问题兵工学报，（）：，（）：高阳阳，余敏建，杨杰多机空战指挥引导不平衡目标分配方法研究航空工程进展，（）：，（）：朱宝鎏，朱荣昌，熊笑非作战飞机效能评估北京：航空工业出版社，：，：，：作者简介刘玉杰（），男，高级工程师，博士，主要研究方向为航空飞行、航空保障、航空心理、空域管理、航空人才选拔培养。李樾（），男，工程师，博士，主要研究方向为舰载航空保障、舰载机动力学。韩维（），男，教授，博士，主要研究方向为舰载航空保障、舰载机动力学。崔凯凯（），男，工程师，博士，主要研究方向为舰载航空保障、舰载机动力学。