随着无人机(Unmanned Aerial Vehicle，UAV)技术和人工智能的快速发展，利用大量具有自主作战能力且成本低廉的UAV组成UAV集群突破对手防御体系，对目标实施饱和打击以及对入侵机群进行空中拦截是UAV集群作战的重要手段^[1-2]。美军已经开展了多项关于UAV集群研究及试验验证，并将UAV“蜂群”战术进一步向实战推进^[3]。UAV集群作战系统在高对抗的战场环境中，如何充分发挥整体协调优势，高效合理进行攻击火力资源规划，实现多个作战单元协同攻击，同时尽可能减少作战单元消耗，使得作战效能最佳，是UAV集群遂行作战任务需要研究解决的重要问题。火力分配问题，即武器-目标分配(Weapon-Target Assignment, WTA)问题，历来是作战指挥辅助决策研究中的核心内容之一，其解空间随着武器数目和目标总数的增加而呈指数级递增，是多参数、多约束的离散非确定性多项式完全问题^[4]。求解火力资源分配问题的算法，分为传统算法和智能优化算法。传统算法如文献[5]采用线性规划解决了舰空火力分配问题，其他算法还包括分支定界法、动态规划法等；智能优化算法如文献[6-7]分别采用改进的遗传算法，文献[8-9]采用蚁群算法，文献[10]提出了自适应灰狼优化算法；混合优化策略如文献[11-12] 以粒子群算法为基础分别混合模拟退火、引力搜索算法，文献[13]以自适应蚁群优化(ACO)算法为基础结合了遗传算法和粒子群算法特点设计了融合算法，文献[14]采用了动态差分改进的蝙蝠算法等。这些算法都能够获得满意解，但不同程度存在以下缺陷：易早熟、进化速度慢或者算法设计实现困难。在实际允许的时间内求解其最优解是不现实的，只能根据假设的作战原则求其满意解。

狼群算法(Wolf Pack Algorithm，WPA)作为群体智能优化算法^[15]，已成功在高维复杂函数优化、背包问题^[16-17]、无人机航迹规划^[18]等优化问题领域得到了很好的应用，但文献[19-20]指出了WPA还存在寻优精度较低、易陷局部最优、效率不高的弊端，并进行了有效改进。WPA在搜索能力方面还存在上升空间，本文引入了蚁群算法^[21-22]中信息素启发规则改进WPA，针对游走行为中人工狼的更新规则，提出使用信息素引导搜索过程，通过不断更新信息素形成狼群对历次捕猎的“记忆”，并应用于狼群更新机制，以改善算法的全局搜索能力。本文将这种新的信息素启发狼群算法(Pheromone Heuristic Wolf Pack Algorithm，PHWPA)应用于UAV集进攻中的火力分配问题，构建了问题的数学模型，进而给出了算法实现过程，并通过仿真实验验证了该方法的可行性及有效性。

1.   构建火力分配数学模型

UAV集群在进攻作战中突出系统整体优势，可根据目标情况进行调整，快速适应任务要求，进而协同作战，以较少的消耗达成作战目的。假设UAV集群作战攻击多个目标的具体场景如下：

集群共有v型UAV，第z型UAV数量为m_z(z=1，2，…，v)，且，每架UAV视为一个作战单元，仅携带一种火力载荷；待攻击目标总数为n个，第g个目标的价值W_g(g=1，2，…，n)；p_zg为第z型UAV对第g个目标攻击的单次杀伤概率；决策变量x_zg为用于攻击第g个目标的第z型UAV的数量，x_zg≥0，第z型UAV对第g个目标的杀伤概率P_zg为

所有v型UAV对第g个目标的综合杀伤概率P_g为

火力分配是以UAV集群作战系统的整体作战效能最优为目标的，即在保证任务完成的情况下，消耗较少的UAV作战单元，使得攻击后的目标价值收益F最大。本文采用了带有杀伤概率门限的火力分配数学模型^[11]。

式中：P_g为P_g的均值。

模型约束条件包括：

1) 任务完成约束。每个目标至少分配1架UAV对其进行攻击。

式中：A_g为第g个目标分配的UAV总数。

2) 有效杀伤约束。该约束条件要求对每个目标的综合杀伤概率大于预设的杀伤概率门限。

式中：P_dg为第g个目标的预设杀伤门限，可根据具体情况指定。

3) 攻击消耗约束。用于攻击的任何一型UAV不能超过其数量限制。

整个模型的特征如下：

1) 保证每个目标均能被有效杀伤。通过衡量目标综合杀伤概率是否超过预设杀伤概率门限，判定目标是否被有效杀伤，若低于预设杀伤概率门限，则认为对目标的分配为无效分配。

2) 目标价值W_g可以保证高价值目标被优先分配。

3) A_g的大小比p_zg对P_g值的影响更大，A_g值越小，综合杀伤概率的平均值就越大，所以模型可以保证使用较少的UAV火力单元。

4) P_g保证目标g在分配UAV火力单元目相同的情况下，即A_g相同，选择P_g大的UAV组合，使综合杀伤概率尽可能大。

上述模型带有非线性约束，采用罚函数法将其转化为一个无约束优化问题来求解。该问题属于非线性整数规划问题，这里采用PHWPA进行解决。

2.   算法实现

2.1   狼群算法

狼群组织严密，分工明确，通过各自履行责任进而完成共同协作，保证狼群的生存和发展。WPA模拟狼群分工协作式捕猎行为、猎物分配规则，通过狼群个体对猎物气味、环境信息的探知、人工狼群相互间信息的共享和交互以及人工狼基于自身职责的个体行为决策最终实现了狼群捕猎的全过程，相互关系如图 1所示^[15]。

图  1  狼群的捕猎模型

Figure  1.  Hunting model of wolf pack

下载: 全尺寸图片 幻灯片

根据狼群特征，WPA采用基于人工狼主体的自下而上的设计方法和基于职责分工的协作搜索路径结构，将人工狼区分为头狼、探狼、猛狼，整个捕猎过程抽象为游走、召唤、围攻3种智能行为以及“胜者为王”的头狼产生规则和“强者生存”的狼群更新机制，相关含义参见文献[15]。

2.2   求解火力分配问题的WPA

围绕火力分配模型，首先进行了火力分配问题的算法描述，并借鉴文献[17]中运动算子设计方法，重新设计了游走、召唤2个算子改进WPA，以便于问题的求解；进而给出了头狼产生规则、游走行为、召唤行为、围攻行为、狼群更新机制等智能行为的详细描述以及新算子在具体智能行为中的应用。

在解决WTA问题中采用了基于整数的编码方式，分配方案X=(x₁，x₂，…，x_j，…，x_m)，变量x_j为0~n之间的整数，x_j=t表示将第j架UAV分配给第t个目标；x_j=0表示第j架UAV没有分配给任一目标^[6]。用人工狼的位置代表一种候选分配方案，设人工狼位置矢量维度为m(UAV总数量)，N为人工狼总数，在N×m的欧式空间中人工狼i的位置X=(x_i1, x_i2, …，x_ij, …, x_im)，x_ij为第i(i=1，2，…，N)匹人工狼在第j(j=1，2，…，m)维变量空间中所处的位置。人工狼感知到的猎物气味浓度Y=f(X)，即目标函数值；人工狼p与人工狼q之间距离d_pq为两者位置编码的Manhattan距离：

定义1  游走算子Ω(X_i，M_a，r)，人工狼i的位置为X_i=(x_i1, x_i2, …，x_ij, …, x_im)，M_a= {1，2，…，m}为编码位集合，可理解为人工狼的可活动范围，r为进行改变的编码位的数目，可理解为人工狼的游走步长。游走算子表示在M_a中随机选择r个编码位形成集合R，将x_ij(j∈R)改变为第j维变量空间中选取的随机数。

定义2  召唤算子Ψ(X_i，M_b，r)，人工狼i的位置为X_i=(x_i1，x_i2，…，x_ij，…，x_im)，M_b为人工狼位置X_i和头狼位置X_d不相同编码位的集合且不为空集，集合M_b为

式中：j=1，2，…，m；k的初值为1；null表示空值；x_dj表示头狼位置第j维取值。召唤算子Ψ为在M_b中随机选择r个编码位形成集合R，并将x_ij值按式(9)进行改变，可理解为人工狼i接收到了头狼传递的部分信息。

WPA的规则和智能行为^[15-17]如下：

1) 头狼产生规则。算法中具有最优目标函数值的人工狼为头狼，迭代过程中根据目标函数进行头狼更替，头狼不执行游走、召唤、围攻行为。

2) 游走行为。选取最优人工狼为头狼，除头狼外所有的人工狼视为探狼，探狼以游走步长step_a向h个方向试探性走一步，即对探狼位置X_i执行h次游走算子X_ip^new= Ω(X_i，M_a，step_a)(p=1，2，…，h)。而后将所得的h个新位置{X_ip^new}进行贪婪决策，即选择具有最优目标函数值且优于原有位置的为该人工狼新位置。设探狼i和头狼位置对应的目标函数值分别为Y_i和Y_lead，游走次数为T，最大游走次数为T_max，重复以上过程直至Y_i＞Y_lead或T≥T_max。

3) 召唤行为。猛狼(除头狼外所有的人工狼)经头狼召唤迅速以较大奔袭步长step_b向头狼所在位置X_d靠拢，即对猛狼位置X_i执行1次召唤算子X_i^new=Ψ(X_i，M_b，step_b)。若M_b为空集时执行游走算子Ψ(X_i，M_a，1)。若Y_i＞Y_lead，则Y_lead=Y_i，猛狼i替代头狼进行召唤；若Y_i < Y_lead，则猛狼i继续奔袭直到猛狼i与头狼的距离小于判定距离d_near。

4) 围攻行为。将头狼所在位置X_d视为猎物的位置，参与围攻的人工狼i的位置X_i依式(10)进行位置变换得到新位置：比较人工狼实施围攻行为前后在新旧位置所感知到的猎物气味浓度并进行贪婪决策。

上述智能行为所涉及的游走步长step_a、奔袭步长step_b、围攻步长step_c皆为整数，表示人工狼搜索的精细程度。

5) 狼群更新机制。按照狼群更新机制进行群体更新，即淘汰N_new匹人工狼，N_new为[N/(2U)，N/U]之间的随机整数，U为更新比例因子。算法模拟自然界狼群繁衍方式，新人工狼作为头狼子女继承头狼的优良基因，即头狼的部分编码位，新的人工狼位置X_new由式(11)计算得到：

式中：X_d为头狼所在位置；M_a={1，2，…，m}；编码位改变的数目L由式(12)计算得到：

式中：*为向下取整；k为算法迭代次数；k_max为算法的最大迭代次数；N为人工狼总数。

2.3   蚁群算法中的信息素启发规则

蚁群算法模拟自然界中蚁群觅食机制，采用了分布式正反馈并行计算机制^[21-22]。WTA问题中蚂蚁各自构建其目标分配方案，蚂蚁从第1架UAV开始在所有目标中选择一个目标分配给该UAV；接着对第2架UAV，蚂蚁s在当前运行分配的目标集合中选择一个目标分配给该UAV；依此顺序分配，指导完成全部目标分配^[8]。算法中信息素启发的具体规则如下：

1) 状态转移规则。蚂蚁s依据伪随机规则选择目标g分配给第u架UAV。

式中：τ_ut(k)为k时第u架UAV与目标t之间的信息素，k即迭代次数；η_ut为与问题相关的启发信息，在火力分配问题中，η_ut根据数学模型的最优准则给定，本文设为第u架UAV对目标t单次杀伤概率与目标t价值的积；α和β分别为τ_ut和η_ut的相对重要性；q为0~1之间均匀分布的随机数；q₀(0≤q₀≤1)为一个指定常数；G为依据如下随机比例规则从allowed_s中选择一个目标。

式中: allowed_s为蚂蚁s当前可分配的目标集； η_ug为第u架UAV对目标g单次杀伤概率与目标g价值的积；τ_ug(k)为k时第u架UAV与目标g之间的信息素。

2) 信息素更新规则。当一次迭代过程结束后，根据当前最优目标函数值F_max及其对应解X_best，按式(15)全局更新规则对信息素进行调整：

式中: Δτ_ug为第u架UAV选择目标g相应信息素增量；ρ为信息素挥发因子，0 < ρ < 1。

式中：Q为信息素强度，是一个常数；F_max(k)为第k次迭代目标函数最优值；X_best(k)为对应的最优解。

2.4   算法改进策略

为强化对解空间的学习，避免陷入局部最优，提高算法寻优效率，在求解火力分配的WPA基础上提出了PHWPA，针对WPA的游走行为和更新机制做出如下改进：

1) 在人工狼搜索过程中，增加信息素引导，可以理解为狼群搜捕猎物的过程形成了一定记忆。在游走行为中，对人工狼位置执行游走算子时，相应编码位由选取随机数改为按照式(13)、式(14)选取新位置。

2) 在更新狼群过程中使用信息素，这一过程可以理解为：一方面新的人工狼继承了头狼的优良基因，另一方面采用信息素引导使新的人工狼传承了整个狼群在捕猎中形成的“智慧”，符合自然界种群繁衍进化的特点。在狼群更新机制中，新人工狼按照式(12)产生时，同样其编码位按照式(13)、式(14)计算其改变值。

信息素更新的时机选择在每次迭代时狼群完成围猎之后、更新狼群之前，利用头狼位置信息及其目标函数值按式(15)对信息素进行全局更新，不断更新的信息素代表着狼群在整个捕猎过程中形成的捕猎“智慧”，有利于该智能优化算法对解空间进行更好的学习。

2.5   PHWPA的实现步骤

步骤1  初始化。人工狼总数N，初始化每匹人工狼的位置X_i，最大迭代次数k_max，更新比例因子U，最大游走次数T_max，各步长step_a、step_b、step_c，判定距离d_near，信息素的启发因子α，期望启发因子β，信息素强度Q，信息素挥发因子ρ以及常数q₀，信息素初值τ_initial，最大值τ_max，最小值τ_min。

步骤2  计算各人工狼的目标函数值。选取头狼，其余人工狼执行游走行为，利用式(13)、式(14)更新，直到某匹人工狼感知的猎物气味浓度(目标函数值)大于头狼感知的气味浓度，或达到最大游走限制次数，转入步骤3。

步骤3  除头狼外的所有人工狼执行召唤行为，向发出召唤的头狼进行奔袭。若奔袭过程中，人工狼感知的气味浓度大于头狼所感知的气味浓度，则进行头狼更替，之后人工狼向新的头狼位置奔袭，直到人工狼与头狼的距离d≤d_near，转入步骤4。

步骤4  除头狼外的所有人工狼执行围攻行为。将头狼所在位置视为猎物的位置，参与围攻的人工狼位置根据式(10)进行变换并进行贪婪决策。

步骤5  全局信息素更新。根据当前头狼位置信息及其所感知的猎物气味浓度，根据式(15)、式(16)对全局信息素进行更新。

步骤6  狼群更新。按照狼群更新机制，更新N_new匹人工狼，新人工狼变化的编码位相应数值按照式(13)、式(14)计算得到。

步骤7  判断是否结束。判断是否达到优化精度要求或最大迭代次数k_max，若达到则输出头狼的位置，即所求问题的最优解，否则转步骤2。

综上，可得出基于PHWPA的火力分配流程，如图 2所示。

图  2  基于PHWPA的火力分配流程

Figure  2.  WTA flowchart based on PHWPA

下载: 全尺寸图片 幻灯片

3.   算法验证及仿真

为验证PHWPA求解火力分配问题的可行性和有效性，进行了仿真实验，算例参考文献[11]。

设UAV集群共有4种型号，每型4架，需要协同攻击10个目标。UAV编号与目标标号如表 1所示。

表  1  无人机编号与目标编号

Table  1.  UAV number and target number

UAV型号	UAV编号	目标编号
Ⅰ型	V1~V4	T1~T10
Ⅱ型	V5~V8	T1~T10
Ⅲ型	V9~V12	T1~T10
Ⅳ型	V13~V16	T1~T10

下载: 导出CSV 

| 显示表格

UAV集群中作战单元对目标的杀伤概率矩阵为P。

各目标的杀伤概率门限P_dj均设为0.9。各目标的价值矩阵W=[0.6，0.7，0.3，0.5，0.6，0.35，0.65，0.55，0.4，0.75]。

采用本文的PHWPA、WPA、ACO算法、最大最小蚁群(MMAS)算法、模拟退火离散粒子群算法(SA-DPSO)^[11]以及仅在更新狼群时使用信息素的PHWPA1、仅在探狼搜索过程使用信息素引导的PHWPA2分别对该火力分配问题进行优化求解。最大迭代次数均为100，涉及WPA的参数：狼群规模N_w=32，更新比例因子U_w=6，游走次数T_max=10，游走步长step_a=2，奔袭步长step_b=4，围攻步长step_c=1，判定距离d_near=2step_b；涉及ACO算法的参数：蚁群规模N_a=32，α=2，β=5，Q=20，ρ=0.4，q₀=0.8，信息素初始值τ_initial=10，信息素最大值τ_max=10，信息素最小值τ_min=2。

图 3表示采用PHWPA得出的最优方案，该火力分配方案共使用4型13架UAV攻击10个目标，具体UAV与目标对应分配情况如表 2所示。表 3中各目标杀伤概率均达到杀伤概率门限指标要求，方案表明在保证任务完成的情况下减少了UAV作战单元的消耗，有利于保持UAV集群整体优势。

图  3  无人机-目标最优分配方案

Figure  3.  Optimal UAV-target assignment

下载: 全尺寸图片 幻灯片

表  2  最优攻击分配方案

Table  2.  Optimal UAV-target assignment

UAV型号	UAV编号	对应分配目标编号
Ⅰ型	V1  V2  V3  V4	T8  T5  T6  T3
Ⅱ型	V5  V6  V7  V8	-  T4  -  T8
Ⅲ型	V9  V10  V11  V12	-  T9  T2  T10
Ⅳ型	V13  V14  V15  V16	T7  T2  T10  T1

下载: 导出CSV 

| 显示表格

表  3  各目标杀伤概率

Table  3.  Kill probability of each target

目标编号	杀伤概率
T1	0.9
T2	0.985
T3	0.91
T4	0.91
T5	0.9
T6	0.9
T7	0.979
T8	0.91
T9	0.92
T10	0.979

下载: 导出CSV 

| 显示表格

由图 4、图 5可看出，因为初始狼群、蚁群是随机产生的，所以每次迭代会略有不同，为检验PHWPA的收敛性，并证明其有效性。对各算法分别进行了100多次仿真实验，算法收敛趋势基本一致，PHWPA与各算法目标函数值变化情况如图 5所示，可以看出，PHWPA能够迅速收敛至全局最优解，并且每次仿真均可以收敛至最优解；WPA收敛至全局最优解时间较长，收敛速度相对较慢，且部分解是局部最优解；基本ACO算法收敛速度慢，且常陷入局部最优解；MMAS算法收敛趋势与WPA相近，且部分解为局部最优解；SA-DPSO算法能够稳定收敛，但常陷入局部最优解；PHWPA1、PHWPA2收敛趋势与PHWPA较为接近，但实验中这2种算法也存在收敛速度慢、陷入局部最优解的情况，不如PHWPA收敛情况稳定。因此，本文提出的PHWPA能够有效提高对全局最优解的寻优效率。

图  4  各算法最优值迭代过程

Figure  4.  Iterative process of optimal value of each algorithm

下载: 全尺寸图片 幻灯片

图  5  PHWPA与其他算法迭代过程比较

Figure  5.  Comparison of iterative process between PHWPA and other algorithms

下载: 全尺寸图片 幻灯片

4.   结论

UAV集群作战是UAV自主作战发展的必然趋势，战场环境下UAV集群协同攻击多目标的火力分配问题研究具有重要意义。

1) 针对UAV集群进攻作战火力分配问题进行了研究，主要考虑达成有效杀伤目的、同时节约UAV火力资源且满足任务约束等要求，建立了合理的火力分配数学模型。

2) 在使用WPA进行模型求解中，重新设计了游走、召唤算子，并在智能行为中描述新算子的具体应用。

3) 针对WPA在解决火力分配问题时易陷入局部最优解的问题，提出了PHWPA，引入ACO算法中信息素启发机制，强化了人工狼群对搜索空间的认知，对狼群算法搜索环节进行启发式引导，赋予了人工狼群“记忆”，在人工狼群更新环节引入信息素启发，继承了狼群的捕猎“智慧”，提升了狼群算法的搜索效率。

仿真结果表明，该优化算法是快速且有效的，在寻优精度和稳定上都有较好的表现，为研究UAV集群作战系统火力分配问题提供了一种新的解决方案。