CN114185362A - 一种基于郊狼信息熵的无人机集群任务动态分配方法 - Google Patents

一种基于郊狼信息熵的无人机集群任务动态分配方法 Download PDF

Info

Publication number
CN114185362A
CN114185362A CN202111485575.3A CN202111485575A CN114185362A CN 114185362 A CN114185362 A CN 114185362A CN 202111485575 A CN202111485575 A CN 202111485575A CN 114185362 A CN114185362 A CN 114185362A
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
task
attack
target
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111485575.3A
Other languages
English (en)
Other versions
CN114185362B (zh
Inventor
段海滨
彭雅兰
魏晨
邓亦敏
吴江
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202111485575.3A priority Critical patent/CN114185362B/zh
Publication of CN114185362A publication Critical patent/CN114185362A/zh
Application granted granted Critical
Publication of CN114185362B publication Critical patent/CN114185362B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10—Simultaneous control of position or course in three dimensions
    • G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/104—Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开一种基于郊狼信息熵的无人机集群任务动态分配方法，步骤如下：步骤一：搭建无人机集群典型任务场景；步骤二：搭建单体无人机运动学模型；步骤三：分配所有无人机为侦察模式；步骤四：攻击编队初分配；步骤五：生成有效攻击编队组合；步骤六：根据信息熵选择打击编队；步骤七：攻击编队中的无人机控制指令生成；步骤八：输出多无人机动态任务分配结果。本发明可实现多无人机的动态任务分配过程；综合考虑了资源均衡、目标威胁与目标点动态变化的情况，能够适应复杂多变的任务场景，且在资源均衡度与任务完成效率方面均表现出较为明显的优势，增强了集群协同执行任务的时效性，能够支持高效、持久的多无人机协同任务。

Description

一种基于郊狼信息熵的无人机集群任务动态分配方法
技术领域
本发明是一种基于郊狼信息熵的无人机集群任务动态分配方法，属于无人机自主控制领域。
背景技术
多无人机动态任务分配(Multi-UAV Task Dynamic Allocation)技术是多无人机集群自主控制关键技术之一，不仅为无人机编队自主控制决策提供基础，同时也为编队控制管理提供飞行目标及任务清单。通过多无人机动态任务分配算法，集群系统内部可以实现自主的协调控制，将复杂的任务流程拆分成若干个简单的任务，通过多无人机之间多个任务的分配和协调控制，实现对复杂任务更加高效可靠地完成。
在信息化革命的推动下，无人机集群的任务目标对象从单目标到多目标、体系目标发展，相对应的执行任务方式也发生变化。针对复杂的任务目标，无人机间的协同成为集群是否能够完成复杂任务的重要条件。对于集群协同系统，如果缺乏行之有效的协同控制策略，不仅无法体现集群较之于单机的优势，甚至会出现机间的冲突，甚至碰撞的严重情况。为了确保无人机集群安全、高效地执行任务，需要根据任务需求、无人机飞行特性及任务载荷等关键因素进行合适的任务规划，确保集群中每个无人机都能充分发挥自身的作用，且实现“一加一大于二”的集群加成，提高整体的任务执行效能。集群协同多任务分配问题的目标是在执行任务无人机集群性能允许的范围内，针对不同用户提出的多任务需求，确定出各无人机的任务执行序列，包括目标集合、执行任务的类型以及不同任务的执行顺序和时间等，确保集群在多任务执行过程中的协同性。研究无人机集群高效率协同多任务分配是提升无人机集群未来飞行能力和执行任务效能的重要途径，具有较大的理论和实际意义。
真正具备完全自主能力的多无人机系统还没有出现，但随着技术进步，已经具备了部分自主能力。目前，对多无人机协同动态任务规划研究主要分为集中式任务分配和分布式任务分配两类。在集中式算法设计方面，加州理工大学学者采用混合整数线性规划方法，建立目标函数和约束条件，求解无人机集群协同攻击固定目标的时序安排问题；俄亥俄州大学学者采用匈牙利算法对集群协同目标分配问题进行研究，结合约束规划方法得到变量集合内的所有变量对应的值域，进而求解对应的组合优化问题。在分布式算法设计方面，美国空军研究室学者利用分布式博弈决策模型，实现了群对群对抗中的目标分配；法国里昂大学学者采用拍卖算法将“招标-投标-中标-确认”四个交互阶段转为竞价模式转卖，定义参与方、拍卖品、收益函数和赢家策略等要素，实现分布式集群协作指派。
综上，上述传统的大规模群体决策方法无法智能处理单一方案决策问题，针对于此，本发明提出了一种基于郊狼信息熵的无人机集群动态任务分配方法，在不确定的任务环境中，针对多无人机协同任务动态规划问题，建立了基于郊狼信息熵的无人机集群协同任务分配流程，优化了多无人机协同任务分配的时效性和任务完成效率，且同时考虑了集群内资源分配的均衡性与目标威胁属性，降低了对同一目标的重复执行次数，保证了集群执行任务的持久性，有效提高了多无人机集群的自主性。
发明内容
本发明目的在于提出一种基于郊狼信息熵的无人机集群任务动态分配方法，以提供一种高效可行的无人机集群动态任务分配方法，旨在适应动态多变的任务环境，符合无人机实际飞行约束，提高集群执行任务的效率与有效性，可以支持持久、高效的任务执行过程，为进一步提升无人机集群的自主水平奠定基础。
针对无人机集群协同搜索打击问题，本发明提出了一种基于郊狼信息熵的无人机集群任务动态分配方法，该方法的流程图如图1所示，具体实施步骤如下：
步骤一：搭建无人机集群典型任务场景
如图2所示，对于任务区域Ω，将其划分为Dx×Dy×Dz个同等大小的三维离散网格，Dx、Dy和Dz分别表示在x、y和z三个方向上的网格数量，第d个网格的三维边界坐标为(xdmin,xdmax)×(ydmin,ydmax)×(zdmin,zdmax)。
在给定的任务空间Ω内，存在我方无人机N架，该集群对K个目标执行协同搜索打击任务。考虑多种无人机类型，各类无人机具备的侦察与打击目标的能力不同，设无人机集群中共有NU类无人机，分别采用概率矩阵PR与PA分别表示各无人机对目标的侦察与打击概率，概率矩阵PR与PA为：
Figure BDA0003396378690000031
其中，p为无人机类型索引，PR(p)表示p类无人机对目标的侦察概率，PA(p)表示p类无人机对目标的打击概率。
设置K个任务目标随机分布在任务区域内，每个任务目标具有一定的攻击价值Rk(t)，其随时间变化的函数可表示为：
Figure BDA0003396378690000032
其中，k为任务目标编号Rk(t)为任务目标k在t时刻具有的攻击价值，
Figure BDA0003396378690000033
为任务目标k在初始时刻具有的攻击价值，0＜η＜1为目标攻击价值随时间的衰减因子。
步骤二：搭建单体无人机运动学模型
单体无人机运动模型选用如下所示的二阶系统模型，每架无人机都配置了速度、航向和航迹三通道的自动驾驶仪：
Figure BDA0003396378690000041
其中，
Figure BDA0003396378690000042
Figure BDA0003396378690000043
与
Figure BDA0003396378690000044
为无人机i的飞行速度、偏航角与航迹角的输入控制指令，Vi、χi与γi为第i架无人机的飞行速度、偏航角与航迹角，g为重力加速度g＝9.8m/s2，Vmin与Vmax为无人机飞行速度限制，ωmax为无人机最大航向角速度的绝对值，γmin与γmax为无人机航迹角限制，τV、τχ与τγ分别为无人机自动驾驶仪飞行速度、偏航角与航迹角的控制响应时间。
步骤三：分配所有无人机为侦察模式
任务开始阶段，所有无人机的初始任务模式都设置为侦察状态，即在任务区域Ω内寻找目标点。无人机判断是否侦察到目标的公式为：
Figure BDA0003396378690000045
其中，
Figure BDA0003396378690000046
表示第i架无人机对第k个目标的二进制探测函数，1表示探测成功，0则表示未成功。d(i)、d(k)分别表示第i架无人机和第k个目标所处的三维任务空间位置；rand为[0,1]内的随机数，p(i)表示第i架无人机的类型。若
Figure BDA0003396378690000047
则第i架无人机向任务区域内其他无人机个体广播探测成功信息，并跳转至步骤四，进入攻击编队预分配；否则，继续执行当前的侦察任务。
步骤四：攻击编队初分配
当
Figure BDA0003396378690000051
即第i架无人机成功侦察到任务目标k，并负责广播对任务目标k的探测成功信息。处于侦察状态的第j(j≠i)架无人机在接收到关于任务目标k的相关信息后，若其尚未成功探测到其他任务目标信息，则按式(5)计算其打击应激概率：
Figure BDA0003396378690000052
其中，
Figure BDA0003396378690000053
表示第j架无人机选择攻击任务目标k的应激概率；
Figure BDA0003396378690000054
表示t时刻任务目标k的打击激励函数；
Figure BDA0003396378690000055
为第j架无人机选择攻击任务目标k的打击应激阈值，其由式(6)确定
Figure BDA0003396378690000056
其中，
Figure BDA0003396378690000057
为第j架无人机选择攻击任务目标k的打击应激阈值；a1、a2、b均为模型参数；δjk表示第j架无人机对任务目标k的打击能力；PA,min为允许的最小打击概率；Δτjk表示第j架无人机到达任务目标k所处位置的预计飞行时间，由式(7)确定:
Figure BDA0003396378690000058
其中，{xk,yk,zk}分别为任务目标k的三维位置；Δdjk表示第j架无人机与任物目标k的相对距离，
Figure BDA0003396378690000059
表示向上取整，α＞1为延迟系数。第j架无人机依据打击应激概率
Figure BDA00033963786900000510
选择是否攻击任务目标k。若是，则第j架无人机向第i架无人机返回应答信息及其无人机类型p(j)；否则，继续执行侦察任务。第i架无人机在一个决策周期Ts内接收来自其他无人机的应答信息，得到应答无人机序列Qi。以此为基础构建针对任务目标k的攻击编队预分配队列Ck＝{Qi,i}，从而进入打击编队分配阶段。
打击应激机制能够使无人机集群快速反应任务目标侦察事件，及时发布和回馈具备有效打击能力的无人机列表，在一个决策周期内即可开展对时敏目标的协同打击，而无须等待其他武器平台支援，从而提高了侦察信息时效性和任务执行效率。
步骤五：生成有效攻击编队组合
首先，第i架无人机根据任务目标k的预分配攻击队列Ck，对于Ck中的第n个子集
Figure BDA0003396378690000061
定义
Figure BDA0003396378690000062
为其具备的攻击总能量，
Figure BDA0003396378690000063
为其能量需求，
Figure BDA0003396378690000064
与
Figure BDA0003396378690000065
分别由式(8)计算：
Figure BDA0003396378690000066
其中，
Figure BDA0003396378690000067
表示任务目标k的期望编队打击概率。若满足
Figure BDA0003396378690000068
则表示子集
Figure BDA0003396378690000069
中的所有无人机具有的总打击能量满足对任务目标k的打击任务需求，将其加入有效编队组合Σk。对Ck的每个子集进行轮询后判断Σk是否为空。若是，则表明当前不存在针对任务目标k的有效编队组合，第i架无人机将任务目标k暂时搁置，清除探测标志
Figure BDA00033963786900000610
并向Ck中的其他无人机发送分配无效信息，使其继续执行侦察任务；否则，对于有效编队组合Σk中的每个子集，按步骤六计算其信息熵。
步骤六：根据信息熵选择打击编队
对于Σk中的第n个子集
Figure BDA00033963786900000611
集合
Figure BDA00033963786900000612
中的每架无人机对于任务目标k都存在一个攻击评价值矢量
Figure BDA00033963786900000613
其含义是第i架无人机对任务目标k执行攻击任务的偏好度，L为影响偏好度的指标个数。按下式计算集合
Figure BDA00033963786900000614
中两两无人机之间的决策相聚度
Figure BDA00033963786900000615
Figure BDA0003396378690000071
其中，1＜a＜+∞，1＜b＜+∞，且
Figure BDA0003396378690000072
引入阈值ρ，0≤ρ≤1，设置条件：
Figure BDA0003396378690000073
以此来判断第i架无人机与第j架无人机针对任务目标k的决策是否可以相聚，当它们的偏好矢量之间的相聚度大于ρ时，则表示它们意见统一，可以相聚，否则认为它们不可以相聚。
判断子集
Figure BDA0003396378690000074
中各个成员是否相聚，由此，剔除成员间不相聚的子集，剩余的所有子集为可用的备选子集，所有备选子集组成备选编队组合Σ′k。假设针对任务目标k，共存在M个备选子集，对于第m个备选子集
Figure BDA0003396378690000075
定义其对任务目标k的攻击评价矢量即攻击偏好度为：
Figure BDA0003396378690000076
对
Figure BDA0003396378690000077
进行标准化处理，得到单位矢量，并记为
Figure BDA0003396378690000078
标准化处理方法为：
Figure BDA0003396378690000079
由此可得到每个备选子集的攻击评价标准化矢量
Figure BDA00033963786900000710
定义其熵为：
Figure BDA00033963786900000711
根据熵的含义，备选子集的熵越大，该子集越接近最优解。按照式(13)计算每个备选子集的熵，选择熵值最大的备选子集作为对任务目标k的打击编队Αk。对于每架已分配无人机j∈Αk，将其任务模式转换为攻击任务模式。
步骤七：攻击编队中的无人机控制指令生成
第i架无人机向Αk中其他无人机广播分配关于任务目标k的有效信息，向未被选择作为打击编队的无人机发送无效分配信息，接受到无效分配信息的无人机转为侦察任务模式，跳转至步骤二。打击编队Αk中的所有无人机，需要同时到达任务目标k点或一定攻击范围内，使得针对任务目标k的整个打击编队同时到达时间满足式(14)：
Figure BDA0003396378690000081
为实现这一目的，采用分散化控制方法，即将偏航角输入控制指令
Figure BDA0003396378690000082
设置为任务目标k相对于第i架无人机的视线角，将航迹角输入控制指令
Figure BDA0003396378690000083
设置为任务目标k相对于第i架无人机的高低角，飞行速度输入控制指令
Figure BDA0003396378690000084
满足一致性分散化控制策略。
步骤八：输出多无人机动态任务分配结果
实时输出无人机集群与任务目标空间分布情况，在到达仿真时间后，输出打击收益变化曲线与任务分配甘特图。
本发明提出的一种基于郊狼信息熵的无人机集群任务动态分配方法，其优点在于：本发明可实现多无人机的动态任务分配过程。该方法综合考虑了资源均衡、目标威胁与目标点动态变化的情况，建立了基于郊狼信息熵的能量分配模型，能够适应复杂多变的任务场景，且在资源均衡度与任务完成效率方面均表现出较为明显的优势，增强了集群协同执行任务的时效性，能够支持高效、持久的多无人机协同任务。
附图说明
图1基于信息熵改进的仿狼群智能多无人机动态任务分配方法流程图
图2三维网格化任务场景地图
图3郊狼信息熵任务分配仿真起始时刻场景
图4郊狼信息熵任务分配仿真300秒时刻场景
图5郊狼信息熵任务分配仿真最终时刻场景
图6任务收益变化曲线
图7任务分配甘特图
图中标号及符号说明如下：
O——坐标系原点
x——坐标系横轴
y——坐标系纵轴
z——坐标系竖轴
t——当前仿真时间
ts——单次仿真迭代时间
具体实施方式
下面通过一个具体的无人机集群动态任务分配实例来验证本发明所提方法的有效性。实验计算机配置为Inter Core i7-8750H处理器，2.20GHz主频，16G内存，软件为MATLAB2020a版本。
本方法的具体步骤如下：
步骤一：搭建无人机集群典型任务场景
任务区域Ω为20千米×20千米×600米，将其划分为10×10×5个同等大小的三维离散网格，即x、y和z三个方向上的网格数量分别为10个、10个和4个。
在任务空间Ω内，存在我方无人机30架，该集群对25个目标执行协同搜索打击任务。考虑3种无人机类型，每种类型的无人机各10架，各种无人机具备的侦察与打击目标的能力不同，设各类无人机对目标的侦察与打击概率矩阵PR与PA为：
Figure BDA0003396378690000091
设置15个任务目标随机分布在任务区域内，每个任务目标具有一定的攻击价值Rk(t)，其随时间变化的函数可表示为：
Figure BDA0003396378690000101
其中，k为任务目标索引，Rk(t)为任务目标k在t时刻具有的攻击价值，设置任务目标k在初始时刻具有的攻击价值
Figure BDA0003396378690000102
η＝0.003为目标攻击价值随时间的衰减因子。
步骤二：搭建单体无人机运动学模型
单体无人机运动模型选用如下所示的二阶系统模型，每架无人机都配置了速度、航向和航迹三通道的自动驾驶仪：
Figure BDA0003396378690000103
其中，
Figure BDA0003396378690000104
Figure BDA0003396378690000105
与
Figure BDA0003396378690000106
为无人机i的飞行速度、偏航角与航迹角的输入控制指令，Vi、χi与γi为第i架无人机的飞行速度、偏航角与航迹角，g为重力加速度g＝9.8m/s2，Vmin＝60千米/小时，Vmax＝600千米/小时，ωmax＝0.8，
Figure BDA0003396378690000107
Figure BDA0003396378690000108
τV＝τχ＝τγ＝5s。
步骤三：分配所有无人机为侦察模式
任务开始阶段，所有无人机的初始任务模式均设置为侦察状态，即在任务区域Ω内寻找目标点。无人机判断是否侦察到目标的公式为：
Figure BDA0003396378690000111
其中，
Figure BDA0003396378690000112
表示第i架无人机对第k个目标的二进制探测函数，1表示探测成功，0则表示未成功。d(i)、d(k)分别表示第i架无人机和第k个目标所处的三维任务空间位置；rand为[0,1]内的随机数，p(i)表示第i架无人机的类型。假设
Figure BDA0003396378690000113
即第i架无人机成功探测到任务目标k，并向任务区域内其他无人机个体广播探测成功信息，并跳转至步骤四，进入攻击编队预分配。
步骤四：攻击编队预分配
Figure BDA0003396378690000114
即第i架无人机成功侦察到任务目标k，并负责广播对任务目标k的探测成功信息。处于侦察状态的第j(j≠i)架无人机在接收到关于任务目标k的相关信息后，若其尚未成功探测到其他任务目标信息，则按式(19)计算其打击应激概率：
Figure BDA0003396378690000115
其中，
Figure BDA0003396378690000116
表示第j架无人机选择攻击任务目标k的应激概率；
Figure BDA0003396378690000117
表示t时刻任务目标k的打击激励函数；
Figure BDA0003396378690000118
为第j架无人机选择攻击任务目标k的打击应激阈值，其由式(20)确定
Figure BDA0003396378690000119
其中，
Figure BDA00033963786900001110
为第j架无人机选择攻击任务目标k的打击应激阈值；a1＝a2＝333，b＝0.01；δjk表示第j架无人机对任务目标k的打击能力；允许的最小打击概率PA,min＝0.7；Δτjk表示第j架无人机到达目标k所处位置的预计飞行时间，由式(21)确定:
Figure BDA0003396378690000121
其中，{xk,yk,zk}分别为任务目标k的三维位置；Δdjk表示第j架无人机与任物目标k的相对距离，
Figure BDA0003396378690000122
表示向上取整，延迟系数α＝1.2。第j架无人机依据打击应激概率
Figure BDA0003396378690000123
选择是否攻击任务目标k。若是，则第j架无人机向第i架无人机返回应答信息及其无人机类型p(j)；否则，继续执行侦察任务。第i架无人机在一个决策周期Ts＝2s内接收来自其他无人机的应答信息，得到应答无人机序列Qi。以此为基础构建针对任务目标k的攻击编队预分配队列Ck＝{Qi,i}，从而进入打击编队分配阶段。
打击应激机制能够使无人机集群快速反应任务目标侦察事件，及时发布和回馈具备有效打击能力的无人机列表，在一个决策周期内即可开展对时敏目标的协同打击，而无须等待其他武器平台支援，从而提高了侦察信息时效性和任务执行效率。
步骤五：生成有效攻击编队组合
首先，第i架无人机根据任务目标k的预分配攻击队列Ck，对于Ck中的第n个子集
Figure BDA0003396378690000124
定义
Figure BDA0003396378690000125
为其具备的攻击总能量，
Figure BDA0003396378690000126
为其能量需求，
Figure BDA0003396378690000127
与
Figure BDA0003396378690000128
分别由式(22)计算：
Figure BDA0003396378690000129
其中，
Figure BDA00033963786900001210
表示任务目标k的期望编队打击概率。若满足
Figure BDA00033963786900001211
则表示子集
Figure BDA00033963786900001212
中的所有无人机具有的总打击能量满足目标k的打击任务需求，将其加入有效编队组合Σk。对Ck的每个子集进行轮询后判断Σk是否为空。若是，则表明当前不存在针对任务目标k的有效编队组合，第i架无人机将任务目标k暂时搁置，清除探测标志
Figure BDA00033963786900001213
并向Ck中的其他无人机发送分配无效信息，使其继续执行侦察任务；否则，对于有效编队组合Σk中的每个子集，按步骤六计算其信息熵。
步骤六：根据信息熵选择打击编队
对于Σk中的第n个子集
Figure BDA0003396378690000131
集合
Figure BDA0003396378690000132
中的每架无人机对于任务目标k都存在一个攻击评价值矢量
Figure BDA0003396378690000133
其含义是第i架无人机对任务目标k执行攻击任务的偏好度，影响偏好度的指标个数L＝4。按下式计算集合
Figure BDA0003396378690000134
中两两无人机之间的决策相聚度
Figure BDA0003396378690000135
Figure BDA0003396378690000136
其中，a＝b＝2，
Figure BDA0003396378690000137
引入阈值ρ＝0.75，设置条件：
Figure BDA0003396378690000138
以此来判断第i架无人机与第j架无人机针对任务目标k的决策是否可以相聚，当它们的偏好矢量之间的相聚度大于ρ时，则表示它们意见统一，可以相聚，否则认为它们不可以相聚。
判断子集
Figure BDA0003396378690000139
中各个成员是否相聚，由此，剔除成员间不相聚的子集，剩余的所有子集为可用的备选子集，所有备选子集组成备选编队组合Σ′k。假设针对目标k，共存在M个备选子集，对于第m个备选子集
Figure BDA00033963786900001310
定义其对目标k的攻击评价矢量即攻击偏好度为：
Figure BDA00033963786900001311
对
Figure BDA00033963786900001312
进行标准化处理，得到单位矢量，并记为
Figure BDA00033963786900001313
标准化处理方法为：
Figure BDA00033963786900001314
由此可得到每个备选子集的攻击评价标准化矢量
Figure BDA0003396378690000141
定义其熵为：
Figure BDA0003396378690000142
根据熵的含义，备选子集的熵越大，该子集越接近最优解。按照式(27)计算每个备选子集的熵，选择熵值最大的备选子集作为对任务目标k的打击编队Αk。对于每架已分配无人机j∈Αk，将其任务模式转换为攻击任务模式。
步骤七：攻击编队中的无人机控制指令生成
第i架无人机向Αk中其他无人机广播分配关于任务目标k的有效信息，向未被选择作为打击编队的无人机发送无效分配信息，接受到无效分配信息的无人机转为侦察任务模式，跳转至步骤二。打击编队Αk中的所有无人机，需要同时到达任务目标点k或一定攻击范围内，使得针对任务目标k的整个打击编队同时到达时间满足式(28)：
Figure BDA0003396378690000143
为实现这一目的，采用分散化控制方法，即将偏航角输入控制指令
Figure BDA0003396378690000144
设置为任务目标k相对于第i架无人机的视线角，将航迹角输入控制指令
Figure BDA0003396378690000145
设置为任务目标k相对于第i架无人机的高低角，飞行速度输入控制指令
Figure BDA0003396378690000146
满足一致性分散化控制策略。
步骤八：输出多无人机动态任务分配结果
实时输出无人机集群与任务目标空间分布情况，图3至图5分别为起始时刻、150秒与最终时刻的仿真结果。在到达仿真时间后Tmax＝300秒，输出任务收益变化曲线如图6与任务分配甘特图如图7。

Claims (6)

1.一种基于郊狼信息熵的无人机集群任务动态分配方法，其特征在于：该方法步骤如下：
步骤一：搭建无人机集群典型任务场景
设置K个任务目标随机分布在任务区域内，每个任务目标具有一定的攻击价值Rk(t)，其随时间变化的函数可表示为：
Figure FDA0003396378680000011
其中，k为任务目标索引，Rk(t)为任务目标k在t时刻具有的攻击价值，
Figure FDA0003396378680000012
为任务目标k在初始时刻具有的攻击价值，0＜η＜1为目标攻击价值随时间的衰减因子；
步骤二：搭建单体无人机运动学模型；
步骤三：分配所有无人机为侦察模式；
若无人机判断侦察到目标，则第i架无人机向任务区域内其他无人机个体广播探测成功信息，并跳转至步骤四，进入攻击编队预分配；否则，继续执行当前的侦察任务；
步骤四：攻击编队初分配
设定打击应激机制，使无人机集群快速反应任务目标侦察事件，发布和回馈具备有效打击能力的无人机列表，在一个决策周期内开展对时敏目标的协同打击；
步骤五：生成有效攻击编队组合
步骤六：根据信息熵选择打击编队
步骤七：攻击编队中的无人机控制指令生成
第i架无人机向打击编队中其他无人机广播分配关于任务目标k的有效信息，向未被选择作为打击编队的无人机发送无效分配信息，接受到无效分配信息的无人机转为侦察任务模式，跳转至步骤二；打击编队中的所有无人机，需要同时到达任务目标k点或一定攻击范围内；
步骤八：输出多无人机动态任务分配结果。
2.根据权利要求1所述的一种基于郊狼信息熵的无人机集群任务动态分配方法，其特征在于：
当第i架无人机成功侦察到任务目标k，并负责广播对任务目标k的探测成功信息；处于侦察状态的第j(j≠i)架无人机在接收到关于任务目标k的相关信息后，若其尚未成功探测到其他任务目标信息，则按式(5)计算其打击应激概率：
Figure FDA0003396378680000021
其中，
Figure FDA0003396378680000022
表示第j架无人机选择攻击任务目标k的应激概率；
Figure FDA0003396378680000023
表示t时刻任务目标k的打击激励函数；
Figure FDA0003396378680000024
为第j架无人机选择攻击任务目标k的打击应激阈值，其由式(6)确定
Figure FDA0003396378680000025
其中，
Figure FDA0003396378680000026
为第j架无人机选择攻击任务目标k的打击应激阈值；a1、a2、b均为模型参数；δjk表示第j架无人机对任务目标k的打击能力；PA,min为允许的最小打击概率；Δτjk表示第j架无人机到达任务目标k所处位置的预计飞行时间，由式(7)确定:
Figure FDA0003396378680000027
其中，{xk,yk,zk}分别为任务目标k的三维位置；Δdjk表示第j架无人机与任物目标k的相对距离，
Figure FDA0003396378680000028
表示向上取整，α＞1为延迟系数；第j架无人机依据打击应激概率
Figure FDA0003396378680000029
选择是否攻击任务目标k；若是，则第j架无人机向第i架无人机返回应答信息及其无人机类型p(j)；否则，继续执行侦察任务；第i架无人机在一个决策周期Ts内接收来自其他无人机的应答信息，得到应答无人机序列Qi；以此为基础构建针对任务目标k的攻击编队预分配队列Ck＝{Qi,i}，从而进入打击编队分配阶段。
3.根据权利要求1所述的一种基于郊狼信息熵的无人机集群任务动态分配方法，其特征在于：所述步骤五的具体过程如下：
首先，第i架无人机根据任务目标k的预分配攻击队列Ck，对于Ck中的第n个子集
Figure FDA0003396378680000031
定义
Figure FDA0003396378680000032
为其具备的攻击总能量，
Figure FDA0003396378680000033
为其能量需求，
Figure FDA0003396378680000034
与
Figure FDA0003396378680000035
分别由式(8)计算：
Figure FDA0003396378680000036
其中，
Figure FDA0003396378680000037
表示任务目标k的期望编队打击概率；若满足
Figure FDA0003396378680000038
则表示子集
Figure FDA0003396378680000039
中的所有无人机具有的总打击能量满足对任务目标k的打击任务需求，将其加入有效编队组合Σk；对Ck的每个子集进行轮询后判断Σk是否为空；若是，则表明当前不存在针对任务目标k的有效编队组合，第i架无人机将任务目标k暂时搁置，清除探测标志
Figure FDA00033963786800000310
并向Ck中的其他无人机发送分配无效信息，使其继续执行侦察任务；否则，对于有效编队组合Σk中的每个子集，按步骤六计算其信息熵。
4.根据权利要求1所述的一种基于郊狼信息熵的无人机集群任务动态分配方法，其特征在于：所述步骤六的具体过程如下：
对于Σk中的第n个子集
Figure FDA00033963786800000311
集合
Figure FDA00033963786800000312
中的每架无人机对于任务目标k都存在一个攻击评价值矢量
Figure FDA00033963786800000313
其含义是第i架无人机对任务目标k执行攻击任务的偏好度，L为影响偏好度的指标个数；按下式计算集合
Figure FDA00033963786800000314
中两两无人机之间的决策相聚度
Figure FDA00033963786800000315
Figure FDA0003396378680000041
其中，1＜a＜+∞，1＜b＜+∞，且
Figure FDA0003396378680000042
引入阈值ρ，0≤ρ≤1，设置条件：
Figure FDA0003396378680000043
以此来判断第i架无人机与第j架无人机针对任务目标k的决策是否可以相聚，当它们的偏好矢量之间的相聚度大于ρ时，则表示它们意见统一，可以相聚，否则认为它们不可以相聚；
判断子集
Figure FDA0003396378680000044
中各个成员是否相聚，由此，剔除成员间不相聚的子集，剩余的所有子集为可用的备选子集，所有备选子集组成备选编队组合Σ′k；假设针对任务目标k，共存在M个备选子集，对于第m个备选子集
Figure FDA0003396378680000045
定义其对任务目标k的攻击评价矢量即攻击偏好度为：
Figure FDA0003396378680000046
对
Figure FDA0003396378680000047
进行标准化处理，得到单位矢量，并记为
Figure FDA0003396378680000048
标准化处理方法为：
Figure FDA0003396378680000049
由此可得到每个备选子集的攻击评价标准化矢量
Figure FDA00033963786800000410
定义其熵为：
Figure FDA00033963786800000411
根据熵的含义，备选子集的熵越大，该子集越接近最优解；按照式(13)计算每个备选子集的熵，选择熵值最大的备选子集作为对任务目标k的打击编队Αk；对于每架已分配无人机j∈Αk，将其任务模式转换为攻击任务模式。
5.根据权利要求1所述的一种基于郊狼信息熵的无人机集群任务动态分配方法，其特征在于：所述步骤七，需使得针对任务目标k的整个打击编队同时到达，其时间满足式(14)：
Figure FDA0003396378680000051
6.根据权利要求5所述的一种基于郊狼信息熵的无人机集群任务动态分配方法，其特征在于：为使得针对任务目标k的整个打击编队同时到达，采用分散化控制方法，即将偏航角输入控制指令设置为任务目标k相对于第i架无人机的视线角，将航迹角输入控制指令设置为任务目标k相对于第i架无人机的高低角，飞行速度输入控制指令满足一致性分散化控制策略。
CN202111485575.3A 2021-12-07 2021-12-07 一种基于郊狼信息熵的无人机集群任务动态分配方法 Active CN114185362B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111485575.3A CN114185362B (zh) 2021-12-07 2021-12-07 一种基于郊狼信息熵的无人机集群任务动态分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111485575.3A CN114185362B (zh) 2021-12-07 2021-12-07 一种基于郊狼信息熵的无人机集群任务动态分配方法

Publications (2)

Publication Number Publication Date
CN114185362A true CN114185362A (zh) 2022-03-15
CN114185362B CN114185362B (zh) 2023-06-06

Family

ID=80603655

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111485575.3A Active CN114185362B (zh) 2021-12-07 2021-12-07 一种基于郊狼信息熵的无人机集群任务动态分配方法

Country Status (1)

Country Link
CN (1) CN114185362B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115454146A (zh) * 2022-10-29 2022-12-09 中国电子科技集团公司第十五研究所 基于相对收益机制的无人机集群协同任务分配方法
CN115525068A (zh) * 2022-10-29 2022-12-27 中国电子科技集团公司第十五研究所 基于迭代优化的无人机集群协同任务分配方法
CN115630514A (zh) * 2022-10-29 2023-01-20 中国电子科技集团公司第十五研究所 一种无人机集群协同任务分配方法与装置
CN115840463A (zh) * 2022-11-23 2023-03-24 北京华如科技股份有限公司 一种用于无人机集群协同侦察的数据处理方法及装置
CN116185035A (zh) * 2023-02-28 2023-05-30 南开大学 基于改进仿生狼群的无人集群动态任务分配方法及系统
CN116860395A (zh) * 2023-07-13 2023-10-10 四川大学 一种分布式分层合同网火力分配方法
CN117148868A (zh) * 2023-10-31 2023-12-01 南京信息工程大学 基于群体熵度量的无人机集群控制方法、装置
CN116860395B (zh) * 2023-07-13 2024-06-11 四川大学 一种分布式分层合同网火力分配方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110069076A (zh) * 2019-04-23 2019-07-30 北京航空航天大学 一种基于猛狼围捕行为的无人机集群空战方法
CN111722643A (zh) * 2020-06-12 2020-09-29 北京航空航天大学 仿狼群合作狩猎机制的无人机集群动态任务分配方法
US20210294356A1 (en) * 2020-03-19 2021-09-23 Beihang University Dynamic recovery method and system for uavs and storage medium
CN113741508A (zh) * 2021-06-29 2021-12-03 南京航空航天大学 基于改进狼群算法的无人机任务分配方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110069076A (zh) * 2019-04-23 2019-07-30 北京航空航天大学 一种基于猛狼围捕行为的无人机集群空战方法
US20210294356A1 (en) * 2020-03-19 2021-09-23 Beihang University Dynamic recovery method and system for uavs and storage medium
CN111722643A (zh) * 2020-06-12 2020-09-29 北京航空航天大学 仿狼群合作狩猎机制的无人机集群动态任务分配方法
CN113741508A (zh) * 2021-06-29 2021-12-03 南京航空航天大学 基于改进狼群算法的无人机任务分配方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
汪汝根;李为民;刘永兰;刘金松;: "无人机集群组网任务分配方法研究", 系统仿真学报 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115454146A (zh) * 2022-10-29 2022-12-09 中国电子科技集团公司第十五研究所 基于相对收益机制的无人机集群协同任务分配方法
CN115525068A (zh) * 2022-10-29 2022-12-27 中国电子科技集团公司第十五研究所 基于迭代优化的无人机集群协同任务分配方法
CN115630514A (zh) * 2022-10-29 2023-01-20 中国电子科技集团公司第十五研究所 一种无人机集群协同任务分配方法与装置
CN115630514B (zh) * 2022-10-29 2023-07-04 中国电子科技集团公司第十五研究所 一种无人机集群协同任务分配方法与装置
CN115840463A (zh) * 2022-11-23 2023-03-24 北京华如科技股份有限公司 一种用于无人机集群协同侦察的数据处理方法及装置
CN116185035A (zh) * 2023-02-28 2023-05-30 南开大学 基于改进仿生狼群的无人集群动态任务分配方法及系统
CN116185035B (zh) * 2023-02-28 2023-09-19 南开大学 基于改进仿生狼群的无人集群动态任务分配方法及系统
CN116860395A (zh) * 2023-07-13 2023-10-10 四川大学 一种分布式分层合同网火力分配方法
CN116860395B (zh) * 2023-07-13 2024-06-11 四川大学 一种分布式分层合同网火力分配方法
CN117148868A (zh) * 2023-10-31 2023-12-01 南京信息工程大学 基于群体熵度量的无人机集群控制方法、装置
CN117148868B (zh) * 2023-10-31 2024-02-27 南京信息工程大学 基于群体熵度量的无人机集群控制方法、装置

Also Published As

Publication number Publication date
CN114185362B (zh) 2023-06-06

Similar Documents

Publication Publication Date Title
CN114185362A (zh) 一种基于郊狼信息熵的无人机集群任务动态分配方法
CN111722643B (zh) 仿狼群合作狩猎机制的无人机集群动态任务分配方法
CN108632831B (zh) 一种基于动态航迹的无人机群频谱资源分配方法
CN106990792B (zh) 混合引力搜索算法的多无人机协同时序耦合任务分配方法
CN113791634B (zh) 一种基于多智能体强化学习的多机空战决策方法
CN110412869B (zh) 一种多星协同的空间动态目标实时跟踪方法
CN108717572B (zh) 基于结构化神经网络的成像任务可调度性预测方法
CN112784445B (zh) 一种飞行控制智能体的并行分布式计算系统及方法
CN113741508B (zh) 基于改进狼群算法的无人机任务分配方法
CN111142553B (zh) 基于生物捕食能量模型的无人机集群自主任务分配方法
CN113009934A (zh) 一种基于改进粒子群算法的多无人机任务动态分配方法
CN107330588B (zh) 一种多基地异构无人机协同侦察的任务规划方法
CN107180309B (zh) 一种空天地观测资源的协同规划方法
CN112733251B (zh) 一种多无人飞行器协同航迹规划方法
CN113448703B (zh) 一种基于感知阵型的无人机蜂群动态侦察任务调度系统及方法
CN114326827B (zh) 一种无人机集群多任务动态分配方法及系统
CN108428004A (zh) 基于蚁群算法的飞行物冲突解脱路径规划方法
CN116307535A (zh) 一种基于改进差分进化算法的多星协作成像任务规划方法
CN113268087A (zh) 多约束复杂环境下基于改进蚁群算法的多无人机协同工作的航迹规划方法
CN109242290B (zh) 一种无人机群行动方案自动生成方法
CN110530373A (zh) 一种机器人路径规划方法、控制器及系统
CN114417735A (zh) 一种跨区域联合作战中的多无人机协同任务规划方法
CN113487220A (zh) 面向静态目标观测的空天异构对地观测资源协同调度方法
CN111797966B (zh) 一种基于改进羊群算法的多机协同全局目标分配方法
CN115564374A (zh) 协同多任务重分配方法、装置、设备及可读存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant