CN115388719B - 一种基于激光制导武器的飞行器激光照射航路规划方法 - Google Patents

Abstract

本申请提供一种基于激光制导武器的飞行器激光照射航路规划方法，包括如下步骤：S1、确定输入信息和输出信息；S2、确定飞行方案；S3、根据飞行器做圆周运动的方向以及飞行器朝向目标飞行的方向确定搜索空间；S4、对飞行器和激光发射器的性能约束进行建模，将飞行器和激光发射器的性能约束转换成搜索区域各维度的约束空间模型；S5、确定评价函数；S6、基于约束空间模型和评价函数，采用四个粒子群分别在四个搜索区域内进行搜索，得到四个最优飞行路径；S7、基于评价函数，从四个最优飞行路径中筛选出最终的飞行器激光照射路径。本申请解决了飞行器激光照射航路规划中搜索空间跨度大，单个粒子群无法跨区域搜索的问题。

Description

一种基于激光制导武器的飞行器激光照射航路规划方法

技术领域

本申请涉及飞行器航路规划技术领域，尤其涉及一种基于激光制导武器的飞行器激光照射航路规划方法。

背景技术

激光制导是一种先进的制导技术，在实际应用中激光制导武器通常具有较高的命中率。激光制导的基本原理就是用激光器发射激光束照射目标，而装在弹体的激光接收装置接收目标反射的激光信号，计算出弹体偏离反射激光束的程度，不断调整飞行轨迹，使战斗部沿着反射激光前进，最终命中目标。激光制导方式可分为半主动寻的式、主动寻的式和波束式三种。目前激光制导武器中大多采用半主动激光制导方式，即导引头与激光照射装置分开配置于两地，前者随弹飞行，后者置于地面、车船或飞机上。

目前关于激光制导的攻击区计算、光斑检测和复杂环境下的抗干扰等方面的研究较多，但是关于飞行器时如何使用激光制导武器的相关研究非常少。当飞行器上装载激光照射器，使用半主动的激光制导方式时，需要预先规划飞行器的飞行路径，以确保在激光制导武器从导引头开机到命中目标的这一时间段内，飞行器能一直提供有效的激光信号。

现有技术中，用于飞行器航路规划的方法主要有人工势场法、模拟退火算法、遗传算法、粒子群算法等。人工势场法是一种反馈控制策略，具有一定的鲁棒性，但是容易产生局部极小值。模拟退火算法过程简单通用，但是收敛速度慢，算法性能对参数变化敏感。与其他算法相比，粒子群算法具有实现简单、收敛速度快、设置参数少等优点。

然而，在飞行器进行激光照射的场景下，由于激光制导武器本身的特性，飞行器具有多个可选择的飞行区域和飞行方式，因此各个搜索区域跨度较大，单个粒子群算法无法完成搜索，即存在搜索不充分的问题。为此，本申请提出一种基于激光制导武器的飞行器激光照射航路规划方法。

发明内容

本申请的目的是针对以上问题，提供一种基于激光制导武器的飞行器激光照射航路规划方法。

本申请提供一种基于激光制导武器的飞行器激光照射航路规划方法，所述方法包括如下步骤：

S1、确定输入信息和输出信息；所述输入信息包括飞行器投放位置、飞行器姿态、飞行器投放时刻、目标位置、激光照射开启时间、激光照射结束时间、激光发射器照射能力模型、飞行器飞行姿态约束；所述输出信息为飞行路径的关键航路点；

S2、确定飞行方案；所述飞行方案为飞行器按照Dubins机动路径切入到绕目标做圆周运动；

S3、根据飞行器做圆周运动的方向以及飞行器朝向目标飞行的方向确定搜索空间；所述搜索空间包括四个搜索区域；

S4、对飞行器和激光发射器的性能约束进行建模，将飞行器和激光发射器的性能约束转换成所述搜索区域各维度的约束空间模型；

S5、确定评价函数；

S6、基于约束空间模型和所述评价函数，采用四个粒子群分别在四个所述搜索区域内进行搜索，得到四个最优飞行路径；

S7、基于所述评价函数，从四个所述最优飞行路径中筛选出最终的飞行器激光照射路径。

根据本申请某些实施例提供的技术方案，步骤S2中，确定的所述飞行方案具体包括：

从武器投放点P₀按照Dubins机动路径切入到绕目标做圆周运动的P₁点；

从切入点P₁点绕目标做半径为r_xz的圆周运动飞到照射起始点P₂；

从照射起始点P₂飞到照射结束点P₃。

根据本申请某些实施例提供的技术方案，步骤S3中，各所述搜索区域表示为：

其中，v₀为P₀点到P₁点的规划速度；v₁为P₁点到P₂点的规划速度；r_xz为水平面转弯半径；为P₂点相对圆心的方位角；/>为P₁点相对圆心的方位角。

根据本申请某些实施例提供的技术方案，步骤S4中，对飞行器和激光发射器的性能约束进行建模，将飞行器和激光发射器的性能约束转换成所述搜索区域各维度的约束空间模型，具体包括：

确定激光发射器的激光照射时间约束；

确定激光发射器的激光照射角度约束，具体包括确定和/>的约束范围；

确定飞行器的飞行性能约束，具体包括确定v₀、v₁以及r_xz的约束范围。

根据本申请某些实施例提供的技术方案，步骤S5中，确定的评价函数公式为：

其中，r_xz为水平面转弯半径；Δt为期望开启照射时间与规划的路径中实际开启照射时间的差；w₁为权重；f₁(r_xz)＝|r_xz-r_xzmax*0.8|；f₂(Δt)＝|Δt|；r_xzmax为r_xz的约束范围的最大值。

根据本申请某些实施例提供的技术方案，步骤S6中，各粒子群在对应的所述搜索区域内进行搜索的方法具体包括：

通过在所述搜索区域中均匀取点的方式对粒子进行初始化，得到各粒子的初始位置以及初始速度；

迭代更新各粒子在所述搜索区域中的位置和速度，并根据更新后的位置确定新的搜索区域，再根据新的搜索区域对位置和速度进行修正；

根据飞行器在投放点P₀和切入点P₁点的方位调用Dubins算法，寻找飞行器的Dubins飞行路径，确定各粒子位置对应的飞行航路；

计算规划航路中的实际开启照射时间，并计算各粒子的当前评价值，更新各粒子的历史最优位置及最优值，并择优更新粒子群的最优位置及最优值；

当搜索过程满足迭代终止条件时退出搜索步骤。

与现有技术相比，本申请的有益效果：本申请提供的基于激光制导武器的飞行器激光照射航路规划方法，使用4个粒子群对飞行器可能飞行的所有区域分别进行搜索，最后综合所有粒子群的输出结果确定飞行器的激光照射路径，解决了飞行器激光照射航路规划中搜索空间跨度大，单个粒子群无法跨区域搜索的问题。该方法使用飞行器通过Dubins路径接近目标，绕其圆周运动进行照射的飞行方案，减少了照射过程中飞行器姿态和位置的影响，简化了激光发射器的照射能力模型，然后考虑飞行器和激光照射器的性能约束，通过坐标系转换、运动学方程等方式，确定了对算法搜索空间进行建模。该方法在粒子群初始化时，对各个变量的搜索空间均匀取点，根据每个维度的取点个数确定粒子数量，使粒子群能充分搜索全部空间，从而减少多次运行该方法输出的评价函数值的方差，保证了每次输出的路径的稳定性。

附图说明

图1为本申请实施例提供的基于激光制导武器的飞行器激光照射航路规划方法的流程图；

图2为四个粒子群的搜索区域示意图；

图3为单个粒子群算法的使用流程图。

具体实施方式

为了使本领域技术人员更好地理解本申请的技术方案，下面结合附图对本申请进行详细描述，本部分的描述仅是示范性和解释性，不应对本申请的保护范围有任何的限制作用。

应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。

本实施例提供一种基于激光制导武器的飞行器激光照射航路规划方法，所述方法的流程图如图1所示，所述方法包括如下步骤：

S1、确定输入信息和输出信息；所述输入信息包括飞行器投放位置、飞行器姿态、飞行器投放时刻、目标位置、激光照射开启时间、激光照射结束时间、激光发射器照射能力模型、飞行器飞行姿态约束；所述输出信息为飞行路径的关键航路点。

其中，激光发射器照射能力模型由激光发射器本身性能决定，飞行器飞行姿态约束是由其自身性能决定，包括飞行倾斜角度限制范围等；飞行路径的关键航路点包括组成飞行器飞行路径的若干离散点，其可以包括从什么时候开始拐弯、从什么时候退出拐弯、什么时候可以直线飞行等。

S2、确定飞行方案；所述飞行方案为飞行器按照Dubins机动路径切入到绕目标做圆周运动。

确定的所述飞行方案具体包括：

飞行器从武器投放点P₀按照Dubins机动路径切入到绕目标做圆周运动的P₁点；

飞行器从切入点P₁点绕目标做半径为r_xz的圆周运动，飞到照射起始点P₂；

开启激光照射器，飞行器从照射起始点P₂飞到照射结束点P₃。

武器挂在飞行器上，当飞行器飞行至武器投放点P₀时，武器投放，飞行器开始从切入点P₁绕目标做圆周运动，当其飞行至P₂点时，开启激光照射，当其飞行至P₃点时，停止激光照射，其中，P₁点、P₂点以及P₃点均在圆周上。

S3、根据飞行器做圆周运动的方向以及飞行器朝向目标飞行的方向确定搜索空间；所述搜索空间包括四个搜索区域；

各所述搜索区域表示为：

其中，v₀为P₀点到P₁点的规划速度；v₁为P₁点到P₂点的规划速度，也为P₂点到P₃点的规划速度，两段路径的规划速度相同；r_xz为水平面转弯半径；为P₂点相对圆心的方位角；为P₁点相对圆心的方位角。

由于照射角度的限制，取值范围为两个相互独立的连续区间，故对其进行拆分处理，在每个区间内，飞行器又可以分为逆时针飞行和顺时针飞行，因此分别使用四个粒子群在各自连续区间内进行搜索；即搜索空间由四个搜索区域组成，四个搜索区域的确定是根据飞行器朝向目标飞行时的左右方向、以及飞行器绕圆的顺时针和逆时针方向两两组合得到的，四个粒子群的搜索区域的示意图如图2所示。

S4、对飞行器和激光发射器的性能约束进行建模，将飞行器和激光发射器的性能约束转换成所述搜索区域各维度的约束空间模型。

步骤S4中，对飞行器和激光发射器的性能约束进行建模，将飞行器和激光发射器的性能约束转换成所述搜索区域各维度的约束空间模型，具体包括：

S41、确定激光发射器的激光照射时间约束；

定义激光照射开启、结束时间分别为t_s和t_e，则激光发射器的激光照射时间约束为t_e-t_s。

S42、确定激光发射器的激光照射角度约束，即确定和/>的约束范围；

定义激光发射器指向目标时相对飞行器机体系的俯仰角为α，目标相对飞行器机体系的离轴角为β；激光对目标的入射角度与导引头相对目标的方位角为 满足

其中，步骤S1中确定的激光发射器照射能力模型，包含了激光照射航段角度约束和/>这里/>有两种取值，一种为：/>另一种为：/>通过将的两种取值分别带入上述不等式后，便可以得到两个不等式，/>是在投放武器时就确定的，从而分别可以确定/>和/>的约束范围，即分别获得/>的约束空间模型。

S43、确定飞行器的飞行性能约束，具体包括确定v₀、v₁以及r_xz的约束范围。

飞行器本身的飞行性能限制包括飞行器滚转限制γ_max、飞行速度限制v_min和v_max、盘旋角速度限制ω_max。此外，激光发射器的限制包括照射距离d和俯仰限制α_max。对规划路径做如下假设：航段P₀P₁的Dubins路径为匀速运动，速度大小为v₀；航段P₁P₂及P₂P₃为匀速运动，速度大小为v₁，且俯仰角近似θ≈0。

因此，可以确定速度v₀、v₁的限制范围，即二者的取值分别在v_min和v_max之间；然后根据飞行器的位置姿态、飞行器和激光发射器的性能限制，确定水平面转弯半径r_xz的约束空间模型，r_xz的范围公式为：

其中，y₀为目标高度，y_T为飞行器航路高度，g为重力加速度。

计算所有单个粒子的上述搜索空间后，规定各维度的速度最大值为各自位置区间长度的0.2倍。

S5、确定评价函数。

步骤S5中，确定的评价函数公式为：

其中，r_xz为水平面转弯半径；Δt为期望开启照射时间与规划的路径中实际开启照射时间的差；w₁为权重，w₁的取值范围为0～1；f₁(r_xz)＝|r_xz-r_xzmax*0.8|；r_xzmax为r_xz的约束范围的最大值；/>为Dubins路径的长度。

为保证飞行器在输入的开启照射时间和结束照射时间点准时到达并开始照射，并且在照射的过程中尽可能远离目标以保证飞行器的安全，设计了上述评价函数，评价函数的结果作为评价飞行方案好坏的指标，评价函数的结果即评价值越小，说明飞行方案越好。

S6、基于约束空间模型和所述评价函数，采用四个粒子群分别在四个所述搜索区域内进行搜索，得到四个最优飞行路径。

粒子群搜索方法采用现有技术中常规的方法，在采用粒子群算法的时候，需要将步骤S1中确定的飞行器投放位置、飞行器姿态、飞行器投放时刻、目标位置、激光照射开启时间、激光照射结束时间等输入进去。各粒子群在对应的所述搜索区域内进行搜索的方法具体包括：

S61、通过在所述搜索区域中均匀取点的方式对粒子进行初始化，得到各粒子的初始位置以及初始速度；

S62、迭代更新各粒子在所述搜索区域中的位置和速度，并根据更新后的位置确定新的搜索区域，再根据新的搜索区域对位置和速度进行修正；

S63、根据飞行器在投放点P₀和切入点P₁点的方位调用Dubins算法，寻找飞行器的Dubins飞行路径，确定各粒子位置对应的飞行航路；

S64、计算规划航路中的实际开启照射时间，并计算各粒子的当前评价值，更新各粒子的历史最优位置及最优值，并择优更新粒子群的最优位置及最优值；

S65、当搜索过程满足迭代终止条件时退出搜索步骤。

具体地，采用在各个搜索区域中均匀取点的方式对粒子进行初始化，从而使粒子遍布整个搜索空间，提高搜索效果。假设在每个搜索维度x₀,x₁,x₂,x₃,x₄的空间内分别均匀取k₀,k₁,k₂,k₃,k₄个点，即在维度x₀的空间内取k₀个点，在维度x₁的空间内取k₁个点，在维度x₂的空间内取k₂个点，在维度x₃的空间内取k₃个点，在维度x₄的空间内取k₄个点，k₀、k₁、k₂、k₃、k₄这5个数的取值均大于1，表示为k_j＞1，(j＝0、1、2、3、4)，则一个搜索区域内的总的粒子个数N＝Πk_j；然后依次在各维度的搜索空间内均匀取点，确定粒子的初始位置，再随机生成粒子各维度的初始速度；然后根据飞行器在P₀和P₁点的方位调用Dubins算法，寻找飞行器的Dubins飞行路径，再根据评价函数计算该粒子的当前评价值，并记为粒子历史最优值，选择各粒子中的最优值作为种群历史最优值。

对于第i(i＝0、1、2......N)个粒子的第j(j＝0、1、2、3、4)个变量的位置的更新公式为/>由于变量x₂,x₃,x₄的约束区间分别与x₁,x₂,x₃有关，因此，从x₂开始，每次更新x_j的位置后，需要根据x_j-1的位置重新确定x_j的约束区间，再对x_j进行修正，然后再对下一个维度x_j+1的位置进行更新。

对于第i个粒子的第j个变量的速度的更新公式为其中，w_v为惯性权重，其取值范围为0～1，其数值大小随迭代次数增加而线性递减，/>为第i个粒子的第j个变量的局部最优，即第i个粒子历史最优值对应的第j个变量取值，gbest_j为第j个变量的全局最优，即是所有粒子中的评价函数最优值对应的第j个变量取值，c₁、c₂为学习因子，其通常取值分别为2，rand(·)为0～1之间的随机数。由于粒子速度的空间约束随粒子变量位置的空间约束变化而变化，因此根据新的粒子速度的约束空间/>对/>进行修正。

迭代的过程中若满足以下任意一个条件，则终止迭代过程：1)迭代次数大于等于最大迭代次数；2)在连续K次迭代中，评价函数的历史最优值变化量均小于Δ。其中，最大迭代次数人为设定，一般设定为4；Δ也为认为设定，其取值比如可以为10^-3。

综上，单个粒子群算法使用流程图如图3所示，其依次包括：确定粒子的空间约束、均匀初始化粒子的位置和速度、保存粒子群和粒子的历史最优值、更新粒子位置、更新粒子空间约束、修正粒子位置并更新粒子速度、更新粒子群和粒子的历史最优值、判断是否满足结束条件、若满足则输出粒子群历史最优对应的路径，若不满足则返回至更新粒子位置的步骤，并重复上述步骤直至满足结束条件为止。

S7、基于所述评价函数，从四个所述最优飞行路径中筛选出最终的飞行器激光照射路径。

通过之前的步骤，可以确定飞行器从P₀点飞到P₁点的飞行路径、P₁点到P₂点、P₃点的路径也知道了，因此可以根据路径的长度计算出飞行时间，同时可以获得P₁点、P₂点、P₃点所在的圆周半径r_xz，将飞行时间和半径代入到步骤S5中的评价函数中，即可以计算出评价函数的值。

因为采用了四个粒子群，因此可以计算出四个评价函数的值，选取最小的评价函数的值所对应的粒子群的规划结果，将其作为最终的飞行器激光照射路径。

本申请实施例提供的基于激光制导武器的飞行器激光照射航路规划方法，采用由Dubins机动路径和绕目标做圆周运动组成的飞行照射方案，建立了一个关于时间和距离的评价函数，然后对激光照射器和飞行器的性能约束进行建模，确定各变量的搜索范围，然后在粒子群初始化时对各个变量的搜索空间均匀取点，使用粒子群算法进行搜索，从而生成飞行器的激光照射关键航路点。该方法提供了一种飞行器激光照射的飞行方案，充分考虑了实际使用过程中的各种限制条件，具有较强的实用性，能够提供较为稳定的输出结果。此外，该方法智能化地提供了飞行路径，既考虑了完成任务的准时性，又考虑了飞行器的安全性。

本申请在飞行器使用激光制导武器的场景下，为飞行器提供多个关键航路点，从而保证飞行器在投放该武器后，在激光导引头开机前，准时到达目标的照射区域，使用激光照射器对目标进行持续照射，直到武器命中目标；本申请解决了现有技术中单个粒子群算法搜索不充分的问题，同时也解决了激光照射航路规划建模困难的问题。

本申请提供的基于激光制导武器的飞行器激光照射航路规划方法，具有以下优点：第一，该方法提供了一种飞行器进行激光照射的飞行方案；第二，该方法充分考虑了激光照射过程中的各种限制条件，具有较强的实用性。第三，该方法能对搜索空间进行充分的搜索，提高了输出结果的稳定性；第四，该方法提供的飞行路径即考虑了完成任务的准时性，又使飞行器尽可能和目标保持一定距离，考虑了飞行器的安全性。

本文中应用了具体个例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的方法及其核心思想。以上所述仅是本申请的优选实施方式，应当指出，由于文字表达的有限性，而客观上存在无限的具体结构，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进、润饰或变化，也可以将上述技术特征以适当的方式进行组合；这些改进润饰、变化或组合，或未经改进将发明的构思和技术方案直接应用于其他场合的，均应视为本申请的保护范围。

Claims (2)

1.一种基于激光制导武器的飞行器激光照射航路规划方法，其特征在于，所述方法包括如下步骤：

S1、确定输入信息和输出信息；所述输入信息包括飞行器投放位置、飞行器姿态、飞行器投放时刻、目标位置、激光照射开启时间、激光照射结束时间、激光发射器照射能力模型、飞行器飞行姿态约束；所述输出信息为飞行路径的关键航路点；

S2、确定飞行方案；所述飞行方案为飞行器按照Dubins机动路径切入到绕目标做圆周运动；

S3、根据飞行器做圆周运动的方向以及飞行器朝向目标飞行的方向确定搜索空间；所述搜索空间包括四个搜索区域；

S4、对飞行器和激光发射器的性能约束进行建模，将飞行器和激光发射器的性能约束转换成所述搜索区域各维度的约束空间模型；

S5、确定评价函数；

S6、基于约束空间模型和所述评价函数，采用四个粒子群分别在四个所述搜索区域内进行搜索，得到四个最优飞行路径；

S7、基于所述评价函数，从四个所述最优飞行路径中筛选出最终的飞行器激光照射路径；

步骤S2中，确定的所述飞行方案具体包括：

从武器投放点P₀按照Dubins机动路径切入到绕目标做圆周运动的P₁点；

从切入点P₁点绕目标做半径为r_xz的圆周运动飞到照射起始点P₂；

从照射起始点P₂飞到照射结束点P_3；

步骤S3中，各所述搜索区域表示为：

其中，v₀为P₀点到P₁点的规划速度；v₁为P₁点到P₂点的规划速度；r_xz为水平面转弯半径；为P₂点相对圆心的方位角；/>为P₁点相对圆心的方位角；

步骤S4中，对飞行器和激光发射器的性能约束进行建模，将飞行器和激光发射器的性能约束转换成所述搜索区域各维度的约束空间模型，具体包括：

确定激光发射器的激光照射时间约束；

确定激光发射器的激光照射角度约束，具体包括确定和/>的约束范围；

确定飞行器的飞行性能约束，具体包括确定v₀、v₁以及r_xz的约束范围；

步骤S5中，确定的评价函数公式为：

其中，r_xz为水平面转弯半径；Δt为期望开启照射时间与规划的路径中实际开启照射时间的差；w₁为权重；f₁(r_xz)＝|r_xz-r_xzmax*0.8|；f₂(Δt)＝|Δt|；r_xzmax为r_xz的约束范围的最大值。

2.根据权利要求1所述的基于激光制导武器的飞行器激光照射航路规划方法，其特征在于，步骤S6中，各粒子群在对应的所述搜索区域内进行搜索的方法具体包括：

通过在所述搜索区域中均匀取点的方式对粒子进行初始化，得到各粒子的初始位置以及初始速度；

迭代更新各粒子在所述搜索区域中的位置和速度，并根据更新后的位置确定新的搜索区域，再根据新的搜索区域对位置和速度进行修正；

根据飞行器在投放点P₀和切入点P₁点的方位调用Dubins算法，寻找飞行器的Dubins飞行路径，确定各粒子位置对应的飞行航路；

计算规划航路中的实际开启照射时间，并计算各粒子的当前评价值，更新各粒子的历史最优位置及最优值，并择优更新粒子群的最优位置及最优值；

当搜索过程满足迭代终止条件时退出搜索步骤。