CN103413186B - 一种基于混合优化算法的多飞行器协同目标分配方法 - Google Patents
一种基于混合优化算法的多飞行器协同目标分配方法 Download PDFInfo
- Publication number
- CN103413186B CN103413186B CN201310361554.XA CN201310361554A CN103413186B CN 103413186 B CN103413186 B CN 103413186B CN 201310361554 A CN201310361554 A CN 201310361554A CN 103413186 B CN103413186 B CN 103413186B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- flight
- fuel consumption
- node
- 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.)
- Expired - Fee Related
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005457 optimization Methods 0.000 title claims abstract description 25
- 239000000446 fuel Substances 0.000 claims abstract description 69
- 239000002245 particle Substances 0.000 claims abstract description 63
- 238000013178 mathematical model Methods 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000002068 genetic effect Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000002153 concerted effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Landscapes
- Traffic Control Systems (AREA)
Abstract
本发明提供了一种基于混合优化算法的多飞行器协同目标分配方法,包括以下步骤:(11)采集我方飞行器的飞行燃料消耗量和飞行时间长度;(12)建立离散粒子群算法中的适应度函数;(13)采用离散粒子群算法对我方飞行器进行目标分配。本发明基于混合优化算法的多飞行器协同目标分配方法,可以应对瞬息万变的战场态势,快速的确定出目标的分配方案,保证了该算法的实时性;同时,该分配方案能够保证与实际吻合,在随后的路径规划中得到印证,保证了该算法的适用性。
Description
技术领域
本发明涉及计算机辅助空战以及电子模拟空战领域,特别是一种基于混合优化算法的多飞行器协同目标分配方法。
背景技术
目标分配算法是多飞行器协同作战的关键技术之一,其主要是依据战场态势,合理对敌方目标进行分配,有效提高战场的整体作战效能。目标分配是一个约束众多而复杂的组合优化问题,其解空间随着武器总数和任务总数的增加而呈几何级数增长,使其成为一个多参数、多约束的NP问题。
传统的求解方法是将目标分配问题简化为数学规划模型,然后采用穷举法、动态规划、分支定界等方法进行求解。然而问题的求解受维数的限制,当问题的维数增加时,运算复杂度将成指数增加,很难得到最优解。后来由于智能算法的发展,提出了大量用智能方法求解的思想,主要有遗传算法、遗传模拟退火算法、蚁群算法、粒子群算法及它们的混合算法。
潘寒尽等针对多机协同空战中目标分配的实时性和复杂性问题,提出了一种基于分布式计算的协同空战目标分配算法。该算法参考空中作战常见的三级指挥体制,构建了相应的分层任务分解结构和多智能体目标分配模型,将传统的集中控制、集中优化问题分布到多个作战控制单元进行计算,从而提高优化的性能。(潘寒尽,王威,邱学军等.基于分布式计算的多机协同空战目标分配算法.电光与控制,2013,20(1):33-35)。
李红亮等针对反舰导弹预先目标分配无法应对战场环境动态变化的难题,提出了基于领弹、数据链、卫星定位技术的在线目标分配的设想。综合考虑影响反舰作战效能的各种因素,构建了目标分配优化模型和作战效能主要评价指标的技术模型。为了满足反舰导弹协同攻击在线目标分配的实时性要求,提出了一种改进的粒子群算法。(李红亮,孔丽,宋贵宝.基于改进粒子群算法的反舰导弹在线目标分配.战术导弹技术,2013,(1):64-70)。
专利申请“一种空战多目标分配的启发式量子遗传方法”(专利申请号:201210147056.0)将协同多目标攻击空战决策问题的威胁经验公式进行变形转换,对我方武器分配方案进行量子位编码,扩大了可行解的表示范围;根据每个武器的多于攻击分配值,提出并设计出优先攻击分配值向量,使染色体按照分配值向量启发式的修正量子染色体,加快收敛速度。
总结国内外公开的基于混合优化算法的多飞行器协同目标分配方法相关文献和研究成果,存在以下两个方面的不足:
1)由于没有对分配策略本身做出改进,单纯以任务代价为优化目标得到的分配方案适应环境变化的能力较差。同时,现代战场上,需要应对瞬息万变的战场态势,态势的改变可能导致原有分配方案效率降低甚至不可行。
虽然目标分配问题的求解速度与精度得到很大提高,但是优化后的结果可能并不和实际吻合,或者说优化后的结果可能在随后的路径规划中得到否认。
发明内容
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种基于混合优化算法的多飞行器协同目标分配方法。
为了解决上述技术问题,本发明公开了一种基于混合优化算法的多飞行器协同目标分配方法,包括以下步骤:
(11)采集我方飞行器的飞行燃料消耗量和飞行时间长度;
(12)建立离散粒子群算法中的适应度函数;
(13)采用离散粒子群算法对我方飞行器进行目标分配。
本发明中,所述采集我方飞行器的飞行燃料消耗量和飞行时间长度(11)步骤具体为:
21)采用动态规划法进行轨迹优化:以我方飞行器燃料消耗量最小为性能指标,采用动态规划法对我方第ii架飞行器飞行至敌方第jj架飞行器的飞行轨迹进行优化,ii=1,2,…,m,jj=1,2,…,n,m为我方飞行器的数量,n为敌方飞行器的数量;
22)采集我方飞行器的飞行燃料消耗量和飞行时间长度:飞行轨迹优化结束后,计算我方第ii架飞行器飞行至敌方第jj架飞行器所需要的飞行燃料消耗量与飞行时间长度。
本发明中所述混合优化算法包括离散粒子群算法与动态规划法的混合。
本发明中,所述采用动态规划法进行轨迹优化(21)步骤中,每架我方飞行器轨迹优化方法包括以下步骤:
31)建立网格坐标系:在笛卡尔直角坐标系o-xyz中,按照我方飞行器的出发点和目标点建立网格坐标系,o为我方飞行器的出发点,X为横向,Y为高度方向,Z为纵向;
32)确定我方飞行器的移动空间走廊(X,Y,Z):
式中:Xmin、Xmax分别为我方飞行器在笛卡尔直角坐标系X轴上的最小、最大距离,Ymin、Ymax分别为我方飞行器在笛卡尔直角坐标系Y轴上的最小、最大距离,Zmin、Zmax分别为我方飞行器在笛卡尔直角坐标系Z轴上的最小、最大距离;
33)对我方飞行器的移动空间走廊进行网格划分:在我方飞行器的移动空间走廊内,进行网格划分,将我方飞行器轨迹用网格上的点P(ik,jk,kk)表示,ik,jk,kk分别代表我方飞行器在网格坐标系中X,Y,Z方向的坐标,ik,jk,kk分别为自然数;
34)建立我方飞行器的数学模型:
式中:P为我方飞行器的推力、为我方飞行器的阻力、为我方飞行器的升力、cx、cy分别为我方飞行器的阻力系数,升力系数、ρ为我方飞行器的所处高度的空气密度、V为我方飞行器的速度、S为我方飞行器的特征面积、θ为我方飞行器的轨迹倾角、ψv为我方飞行器的轨迹偏角、h为我方飞行器的高度、x为我方飞行器的横向距离、z为我方飞行器的侧向距离、α为我方飞行器的迎角、β为我方飞行器的侧滑角、m为我方飞行器的机体质量、t为我方飞行器的飞行时间;
35)确定性能指标J:我方飞行器燃料消耗量mf最小,即:
J=minmf;
36)计算我方飞行器在网格中各个节点段飞行的性能指标:基于我方飞行器的数学模型,分别计算我方飞行器在网格内各个节点段的性能指标,其中各个节点段包括:从原点到ik=1层每个节点,从ik=1层每个节点到ik=2层每个节点,直到从ik=ikmax-1层每个节点到ik=ikmax层每个节点,从ik=ikmax层每个节点到目标点;
37)优化我方飞行器轨迹:根据36)计算出的网格内各个节点段的性能指标,以我方飞行器燃料消耗量mf最小为性能指标,我方飞行器进行轨迹优化。
本发明中,步骤(37)具体包括:
41)计算我方飞行器从ik=ikmax-1层中每个节点到目标点的最优飞行轨迹及燃料消耗量,针对ik=ikmax-1层中每个节点,比较我方飞行器从每个节点到目标点的所有飞行轨迹的燃料消耗量,得出并存储燃料消耗量最小的节点到目标点的最优飞行轨迹与燃料消耗量;
42)计算我方飞行器从ik=ikmax-L层中每个节点到目标点的最优飞行轨迹及燃料消耗量:针对ik=ikmax-L层中每个节点,我方飞行器从该节点按照步骤41)得出的在ik=ikmax-(L-1)层中每个节点到目标点的最优飞行轨迹飞行,对上述所有飞行轨迹的燃料消耗量进行比较,得出并存储燃料消耗量最小的节点到目标点的最优轨迹及燃料消耗量,L=2,3,…,ikmax-1;
43)判断ik=1层是否完成,如果是,转到步骤44);否则,返回步骤42);
44)计算我方飞行器从起点到目标点的最优飞行轨迹及燃料消耗量:我方飞行器从起点按照步骤42)得出的在ik=1层中每个节点到目标点的最优飞行轨迹飞行,对上述所有飞行轨迹的燃料消耗量进行比较,得出并存储燃料消耗量最小的节点到目标点的最优轨迹及燃料消耗量,该飞行轨迹及燃料消耗量即为我方飞行器从起点到目标点的最优飞行轨迹及性能指标。
本发明中,所述建立离散粒子群算法中的适应度函数(12)步骤中,所述的适应度函数具体为:
51)计算我方飞行器的燃料消耗程度C1:我方飞行器的燃料消耗程度是指我方各飞行器在安全航路飞行,命中目标的前提下,通过最小化我方飞行器攻击敌方飞行器的燃料消耗量,引导目标分配向着减少我方飞行器燃料消耗程度的方向进行:
式中,miijj为第ii架我方飞行器以步骤44)的飞行轨迹攻击敌方第jj架飞行器所用的燃料消耗量,m_max为攻击敌方第jj架飞行器中我方飞行器所用的最大的燃料消耗量;
52)计算我方飞行器飞行时间长度比值C2:我方飞行器飞行时间长度比值是指我方各飞行器在安全航路飞行,命中目标的前提下,通过最小化我方飞行器攻击敌方飞行器的飞行时间引导目标分配向着减少我方飞行器飞行时间的方向进行:
式中,tiijj为第ii架我方飞行器以步骤44)的飞行轨迹攻击敌方第jj架飞行器所用的飞行时间,t_max为攻击敌方第jj架飞行器中我方飞行器所用的最长的飞行时间;
53)计算目标的价值收益C3:目标价值收益最大指标通过对我方飞行器执行任务时所获取的目标价值的评估,来引导目标分配的优化和决策向着使作战效能最大化的方向进行,该指标使我方飞行器趋向于攻击高价值目标,目标的价值收益包括目标的价值v、我方飞行器对杀伤概率PS,假设第ii架我方飞行器攻击敌方第jj架飞行器时计算公式为:
54)计算我方飞行器使用均衡度C4:假设m架我方飞行器攻击n个敌方飞行器,则计算公式为:
其中:Nii表示第ii架我方飞行器攻击敌方飞行器的架次,表示我方飞行器攻击敌方飞行器架次的平均数,我方飞行器使用均衡度越小,表示我方各型号飞行器使用越均衡;
55)建立离散粒子群算法中的适应度函数:根据我方飞行器的燃料消耗程度最小化指标,我方飞行器飞行时间的长度最小化指标,目标价值收益最大化指标,最小化我方飞行器使用均衡度,多飞行器协同目标分配的适应度函数C为:
C=ω1C1+ω2C2-ω3C3+ω4C4,
其中:ω1,ω2,ω3,ω4为权系数,且ω1=ω2=ω3=ω4=0.25。
本发明中,所述采用离散粒子群算法对我方飞行器进行目标分配(13)步骤具体为:
61)设置种群大小、迭代次数;
62)随机初始化种群中各粒子的位置Xi和速度Vi:Xi=[xi,1,xi,2,…,xi,d],Vi=[vi,1,vi,2,…,vi,d];其中:i为粒子的序号,d为搜索空间的维度,维度的数量为敌方飞行器的架数,d=1,2,…,n;
63)计算种群中每个粒子的适应度:评价种群中所有粒子,将当前各粒子的位置和目标值存储于各粒子的Pi=[pi,1,pi,2,…,pi,d]中,将所有Pi=[pi,1,pi,2,…,pi,d]目标值最优的个体位置和目标值存储于Pg=[pg,1,pg,2,…,pg,d]中;其中:Pi=[pi,1,pi,2,…,pi,d]为第i个粒子的最优值,Pg=[pg,1,pg,2,…,pg,d]为群体最优值;
64)按下式更新各个粒子的速度vi,j和位置xi,j:
其中:t_t表示粒子的迭代次数,ω=0.729为惯性权因子,b1=2.05和b2=2.05为正的加速常数,r1和r2为在0和1之间均匀分布的随机数,j表示维度;
对于m架我方飞行器攻击敌方n个目标,在第t_t次迭代时,第i个粒子的第j维的位置xi,j和速度vi,j取值范围为:
65)计算种群中每个粒子的适应度,评价种群中所有粒子;
66)更新最优值Pi:比较种群中每个粒子当前目标值与其最优值Pi的目标值,若当前目标值更优,则用粒子的当前位置和目标值更新最优值Pi;
67)更新最优值Pg:比较当前所有最优值Pi和最优值Pg的目标值,更新最优值Pg;
68)若迭代次数达到最大值,则输出最优值Pg及其目标值并结束,否则返回步骤64)。
本发明基于由离散粒子群算法与动态规划法组合的混合优化算法的多飞行器协同目标分配方法,可以应对瞬息万变的战场态势,快速的确定出目标的分配方案,保证了该算法的实时性;同时,该分配方案能够保证与实际吻合,在随后的路径规划中得到印证,保证了该算法的适用性。
本发明与现有技术相比:具有以下显著的特点:
1)为了保证分配算法的实时性,采用离散粒子群算法用于多飞行器协同目标分配方案的制定;
2)为了保证分配算法的适用性,采用动态规划法粗略优化计算每架飞行器到每个目标的路径。
下面结合附图对本发明作进一步详细描述。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1是基于混合优化算法的多飞行器协同目标分配方法的流程图。
图2是采集我方飞行器的飞行燃料消耗量和飞行时间长度的流程图。
图3是采用动态规划法,对我方每架飞行器分别进行轨迹优化的总体框架流程图。
图4是采用动态规划法,优化我方飞行器轨迹的流程图。
图5是我方飞行器的燃料消耗程度,飞行时间的长度,目标价值收益。
图6是采用离散粒子群算法对我方飞行器进行目标分配的流程图。
具体实施方式
结合图1,本发明基于混合优化算法的多飞行器协同目标分配方法,步骤如下:
1)采集我方飞行器的飞行燃料消耗量和飞行时间长度,具体步骤如图2所示。
21)采用动态规划法进行轨迹优化:以我方飞行器燃料消耗量最小为性能指标,采用动态规划法对我方第ii架飞行器飞行至敌方第jj架飞行器的飞行轨迹进行优化,ii=1,2,…,m,jj=1,2,…,n,m为我方飞行器的数量,n为敌方飞行器的数量;例如:我方飞行器的数量:m=6,敌方飞行器的数量:n=12。具体步骤如图3所示。
31)建立网格坐标系:在笛卡尔直角坐标系o-xyz中,按照我方飞行器的出发点和目标点建立网格坐标系,o为我方飞行器的出发点,X为横向,Y为高度方向,Z为纵向。
32)确定我方飞行器的移动空间走廊:
式中:Xmin、Xmax分别为我方飞行器在笛卡尔直角坐标系X轴上的最小、最大距离,Ymin、Ymax分别为我方飞行器在笛卡尔直角坐标系Y轴上的最小、最大距离,Zmin、Zmax分别为我方飞行器在笛卡尔直角坐标系Z轴上的最小、最大距离。
网格的大小(ΔX,ΔY,ΔZ)是由飞行器的机动能力确定的,因为发动机的推力大小是受限制的,飞行空间走廊越来越大,因为随着发动机工作时间的增加,飞行器质量减小,使得飞行器机动能力增强。
例如:
第一架我方飞行器的移动空间走廊为:
0≤X≤10km,5km≤Y≤15km,-5km≤Z≤5km;
第二架我方飞行器的移动空间走廊为:
0≤X≤10km,5km≤Y≤15km,-5km≤Z≤5km;
第三架我方飞行器的移动空间走廊为:
0≤X≤10km,5km≤Y≤15km,-5km≤Z≤5km;
第四架我方飞行器的移动空间走廊为:
0≤X≤10km,5km≤Y≤15km,-5km≤Z≤5km。
33)对我方飞行器的移动空间走廊进行网格划分:在我方飞行器的移动空间走廊内,进行网格划分。将我方飞行器轨迹用网格上的点P(ik,jk,kk)表示,ik,jk,kk分别代表我方飞行器在网格坐标系中X,Y,Z方向的坐标,ik,jk,kk分别为自然数。
网格划分越小,所得到的轨迹精度越高,但付出的代价是时间的花费。所以,应根据实际情况,在网格划分和时间的代价上折衷处理。
例如:
第一架我方飞行器在移动空间走廊网格划分为:1000×1000×1000;
第二架我方飞行器在移动空间走廊网格划分为:1000×1000×1000;
第三架我方飞行器在移动空间走廊网格划分为:1000×1000×1000;
第四架我方飞行器在移动空间走廊网格划分为:1000×1000×1000。
34)建立我方飞行器的数学模型:
式中:P为我方飞行器的推力、为我方飞行器的阻力、为我方飞行器的升力、cx、cy分别为我方飞行器的阻力系数,升力系数、ρ为我方飞行器的所处高度的空气密度、V为我方飞行器的速度、S为我方飞行器的特征面积、θ为我方飞行器的轨迹倾角、ψv为我方飞行器的轨迹偏角、h为我方飞行器的高度、x为我方飞行器的横向距离、z为我方飞行器的侧向距离、α为我方飞行器的迎角、β为我方飞行器的侧滑角、m为我方飞行器的机体质量、t为我方飞行器的飞行时间。
35)确定性能指标:我方飞行器燃料消耗量mf最小,即:
J=minmf。
36)计算我方飞行器在网格中各个节点段飞行的性能指标:基于我方飞行器的数学模型,分别计算我方飞行器在网格内各个节点段的性能指标,其中各个节点段包括:从原点到ik=1层每个节点,从ik=1层每个节点到ik=2层每个节点,…,从ik=ikmax-1层每个节点到ik=ikmax层每个节点,从ik=ikmax层每个节点到目标点;
37)采用动态规划法,优化我方飞行器轨迹:根据36)计算出的网格内各个节点段的性能指标,以我方飞行器燃料消耗量mf最小为性能指标,采用动态规划法对我方飞行器进行轨迹优化。具体步骤如图4所示:
41)计算我方飞行器从ik=ikmax-1层中每个节点到目标点的最优飞行轨迹及性能指标:针对ik=ikmax-1层中每个节点,比较我方飞行器从该节点到目标点的所有飞行轨迹的性能指标,得出并存储该节点到目标点的最优飞行轨迹及性能指标;
42)计算我方飞行器从ik=ikmax-L层中每个节点到目标点的最优飞行轨迹及性能指标:针对ik=ikmax-L层中每个节点,我方飞行器从该节点按照41)得出的在ik=ikmax-(L-1)层中每个节点到目标点的最优飞行轨迹飞行,对上述所有飞行轨迹的性能指标进行比较,得出并存储该节点到目标点的最优轨迹及性能指标,L=2,3,…,ikmax-1;
43)判断ik=1层是否完成。如果是,转到44);否则,转到42)。
44)计算我方飞行器从起点到目标点的最优飞行轨迹及性能指标:我方飞行器从起点按照42)得出的在ik=1层中每个节点到目标点的最优飞行轨迹飞行,对上述所有飞行轨迹的性能指标进行比较,得出并存储该节点到目标点的最优轨迹及性能指标,该飞行轨迹及性能指标即为我方飞行器从起点到目标点的最优飞行轨迹及性能指标。
22)采集我方飞行器的飞行燃料消耗量和飞行时间长度:飞行轨迹优化结束后,对我方第ii架飞行器飞行至敌方第jj架飞行器所需要的飞行燃料消耗量与飞行时间长度两数据进行采集。该步骤采集的数据将作为步骤51)我方飞行器的燃料消耗程度和52)我方飞行器飞行时间长度比值的输入数据。
2)建立离散粒子群算法中的适应度函数,适应度函数具体为:
多飞行器协同目标分配是以整体作战效能最优为目的的,我方飞行器的燃料消耗程度、我方飞行器飞行时间的长度比值、目标的价值收益、我方飞行器的使用均衡度是评价目标分配好坏的主要指标。因此,衡量多飞行器协同目标分配方案优劣的适应度函数主要包括以下四方面:
51)我方飞行器的燃料消耗程度。
我方飞行器的燃料消耗程度是指我方各飞行器在安全航路飞行,命中目标的前提下,通过最小化我方飞行器攻击敌方飞行器的燃料消耗量引导目标分配向着减少我方飞行器燃料消耗程度的方向进行。该指标是通过动态规划法计算飞行器航迹规划得到的。设第ii架我方飞行器攻击敌方第jj架飞行器所用的燃料消耗量为miijj,则:
式中,m_max为攻击敌方第jj架飞行器中我方飞行器所用的最大的燃料消耗量。miijj和m_max的数据由步骤22)采集我方飞行器的飞行燃料消耗量和飞行时间长度输入。
52)我方飞行器飞行时间长度比值
我方飞行器飞行时间长度比值是指我方各飞行器在安全航路飞行,命中目标的前提下,通过最小化我方飞行器攻击敌方飞行器的飞行时间引导目标分配向着减少我方飞行器飞行时间的方向进行。该指标是通过动态规划法计算我方飞行器航迹规划得到的。设第ii架我方飞行器攻击敌方第jj架飞行器所用的飞行时间为:tiijj,则:
式中,t_max为攻击敌方第jj架飞行器中我方飞行器所用的最长的飞行时间。tiijj和t_max的数据由步骤22)采集我方飞行器的飞行燃料消耗量和飞行时间长度输入。
53)目标的价值收益
目标价值收益最大指标通过对我方飞行器执行任务时所获取的目标价值的评估,来引导目标分配的优化和决策向着使作战效能最大化的方向进行。该指标使我方飞行器趋向于攻击高价值目标。综合考虑目标的价值v、杀伤概率PS,则第ii架我方飞行器攻击敌方第jj架飞行器时,目标的价值收益为:
式中:目标的价值v的取值范围:0<v<1,杀伤概率PS的取值范围:0<PS<1。
54)我方飞行器使用均衡度
我方飞行器使用均衡度C4表示我方各飞行器利用率的均衡程度,定义为:
假设m架我方飞行器攻击n个敌方飞行器,则:
其中:Nii表示第ii架我方飞行器攻击敌方飞行器的架次,表示我方飞行器攻击敌方飞行器架次的平均数。我方飞行器使用均衡度越小,表示我方各型号飞行器使用越均衡。
例如:我方飞行器的燃料消耗程度、我方飞行器飞行时间的长度比值、目标的价值收益如图5所示。
55)建立离散粒子群算法中的适应度函数:根据我方飞行器的燃料消耗程度最小化指标,我方飞行器飞行时间的长度最小化指标,目标价值收益最大化指标,最小化我方飞行器使用均衡度,多飞行器协同目标分配的适应度函数为:
C=ω1C1+ω2C2-ω3C3+ω4C4
其中:ω1,ω2,ω3,ω4为权系数,且ω1=ω2=ω3=ω4=0.25,反映了每个指标的重要程度。
3)采用离散粒子群算法对我方飞行器进行目标分配,具体步骤如图6:
多飞行器协同目标分配的关键在于确定任务目标由哪架飞行器来执行。因此,本发明采用自然数编码方式来表达。
61)设置种群大小、迭代次数;
例如:种群大小为50个、迭代次数为100次
62)随机初始化种群中各粒子的位置和速度:Xi=[xi,1,xi,2,…,xi,d],Vi=[vi,1,vi,2,…,vi,d];其中:i为粒子的序号,d为搜索空间的维度,维度的数量为敌方飞行器的架数,d=1,2,…,n。
63)计算种群中每个粒子的适应度:评价种群中所有粒子,将当前各粒子的位置和目标值存储于各粒子的Pi=[pi,1,pi,2,…,pi,d]中,将所有Pi=[pi,1,pi,2,…,pi,d]中目标值最优的个体位置和目标值存储于Pg=[pg,1,pg,2,…,pg,d]中;其中:
Pi=[pi,1,pi,2,…,pi,d]为第i个粒子的最优值,Pg=[pg,1,pg,2,…,pg,d]为群体最优值。
64)按下式更新各个粒子的速度和位置:
其中:t_t表示粒子的迭代次数,ω=0.729为惯性权因子,c1=2.05和c2=2.05为正的加速常数,r1和r2为在0和1之间均匀分布的随机数;
对于m架我方飞行器攻击敌方n个目标,在第t_t次迭代时,第i个粒子的第j维的位置xi,j和速度vi,j取值范围为:
65)计算种群中每个粒子的适应度,评价种群中所有粒子;
66)更新Pi:比较种群中每个粒子当前目标值与其Pi的目标值,若当前目标值更优,则用粒子的当前位置和目标值更新其Pi;
67)更新Pg:比较当前所有Pi和Pg的目标值,更新Pg;
68)若迭代次数达到最大值,则输出Pg及其目标值并停止算法,否则转向步骤64)。
例如:最优分配结果为:510134345220。
本发明提供了一种基于混合优化算法的多飞行器协同目标分配方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (1)
1.一种基于混合优化算法的多飞行器协同目标分配方法,其特征在于,包括以下步骤:
(11)采集我方飞行器的飞行燃料消耗量和飞行时间长度;
(12)建立离散粒子群算法中的适应度函数;
(13)采用离散粒子群算法对我方飞行器进行目标分配;
所述采集我方飞行器的飞行燃料消耗量和飞行时间长度(11)步骤具体为:
21)采用动态规划法进行轨迹优化:以我方飞行器燃料消耗量最小为性能指标,采用动态规划法对我方第ii架飞行器飞行至敌方第jj架飞行器的飞行轨迹进行优化,
ii=1,2,…,m,jj=1,2,…,n,m为我方飞行器的数量,n为敌方飞行器的数量;
22)采集我方飞行器的飞行燃料消耗量和飞行时间长度:飞行轨迹优化结束后,计算我方第ii架飞行器飞行至敌方第jj架飞行器所需要的飞行燃料消耗量与飞行时间长度;
所述采用动态规划法进行轨迹优化步骤21)中,每架我方飞行器轨迹优化方法包括以下步骤:
31)建立网格坐标系:在笛卡尔直角坐标系o-xyz中,按照我方飞行器的出发点和目标点建立网格坐标系,o为我方飞行器的出发点,X为横向,Y为高度方向,Z为纵向;
32)确定我方飞行器的移动空间走廊(X,Y,Z):
式中:Xmin、Xmax分别为我方飞行器在笛卡尔直角坐标系X轴上的最小、最大距离,Ymin、Ymax分别为我方飞行器在笛卡尔直角坐标系Y轴上的最小、最大距离,Zmin、Zmax分别为我方飞行器在笛卡尔直角坐标系Z轴上的最小、最大距离;
33)对我方飞行器的移动空间走廊进行网格划分:在我方飞行器的移动空间走廊内,进行网格划分,将我方飞行器轨迹用网格上的点P(ik,jk,kk)表示,ik,jk,kk分别代表我方飞行器在网格坐标系中X,Y,Z方向的坐标,ik,jk,kk分别为自然数;
34)建立我方飞行器的数学模型:
式中:P为我方飞行器的推力、为我方飞行器的阻力、为我方飞行器的升力、cx、cy分别为我方飞行器的阻力系数,升力系数、ρ为我方飞行器的所处高度的空气密度、V为我方飞行器的速度、S为我方飞行器的特征面积、θ为我方飞行器的轨迹倾角、ψv为我方飞行器的轨迹偏角、h为我方飞行器的高度、x为我方飞行器的横向距离、z为我方飞行器的侧向距离、α为我方飞行器的迎角、β为我方飞行器的侧滑角、m为我方飞行器的机体质量、t为我方飞行器的飞行时间;
35)确定性能指标J:我方飞行器燃料消耗量mf最小,即:
J=minmf;
36)计算我方飞行器在网格中各个节点段飞行的性能指标:基于我方飞行器的数学模型,分别计算我方飞行器在网格内各个节点段的性能指标,其中各个节点段包括:从原点到ik=1层每个节点,从ik=1层每个节点到ik=2层每个节点,直到从ik=ikmax-1层每个节点到ik=ikmax层每个节点,从ik=ikmax层每个节点到目标点;
37)优化我方飞行器轨迹:根据36)计算出的网格内各个节点段的性能指标,以我方飞行器燃料消耗量mf最小为性能指标,对我方飞行器进行轨迹优化;
步骤37)具体包括:
41)计算我方飞行器从ik=ikmax-1层中每个节点到目标点的最优飞行轨迹及燃料消耗量,针对ik=ikmax-1层中每个节点,比较我方飞行器从每个节点到目标点的所有飞行轨迹的燃料消耗量,得出并存储燃料消耗量最小的节点到目标点的最优飞行轨迹与燃料消耗量;
42)计算我方飞行器从ik=ikmax-L层中每个节点到目标点的最优飞行轨迹及燃料消耗量:针对ik=ikmax-L层中每个节点,我方飞行器从该节点按照步骤41)得出的在ik=ikmax-(L-1)层中每个节点到目标点的最优飞行轨迹飞行,对上述所有飞行轨迹的燃料消耗量进行比较,得出并存储燃料消耗量最小的节点到目标点的最优轨迹及燃料消耗量,L=2,3,…,ikmax-1;
43)判断ik=1层是否完成,如果是,转到步骤44);否则,返回步骤42);
44)计算我方飞行器从起点到目标点的最优飞行轨迹及燃料消耗量:我方飞行器从起点按照步骤42)得出的在ik=1层中每个节点到目标点的最优飞行轨迹飞行,对上述所有飞行轨迹的燃料消耗量进行比较,得出并存储燃料消耗量最小的节点到目标点的最优轨迹及燃料消耗量,该飞行轨迹及燃料消耗量即为我方飞行器从起点到目标点的最优飞行轨迹及性能指标;
所述建立离散粒子群算法中的适应度函数(12)步骤中,所述的适应度函数具体为:
51)计算我方飞行器的燃料消耗程度C1:我方飞行器的燃料消耗程度是指我方各飞行器在安全航路飞行,命中目标的前提下,通过最小化我方飞行器攻击敌方飞行器的燃料消耗量,引导目标分配向着减少我方飞行器燃料消耗程度的方向进行:
式中,miijj为第ii架我方飞行器以步骤44)的飞行轨迹攻击敌方第jj架飞行器所用的燃料消耗量,m_max为攻击敌方第jj架飞行器中我方飞行器所用的最大的燃料消耗量;
52)计算我方飞行器飞行时间长度比值C2:我方飞行器飞行时间长度比值是指我方各飞行器在安全航路飞行,命中目标的前提下,通过最小化我方飞行器攻击敌方飞行器的飞行时间引导目标分配向着减少我方飞行器飞行时间的方向进行:
式中,tiijj为第ii架我方飞行器以步骤44)的飞行轨迹攻击敌方第jj架飞行器所用的飞行时间,t_max为攻击敌方第jj架飞行器中我方飞行器所用的最长的飞行时间;
53)计算目标的价值收益C3:目标价值收益最大指标通过对我方飞行器执行任务时所获取的目标价值的评估,来引导目标分配的优化和决策向着使作战效能最大化的方向进行,该指标使我方飞行器趋向于攻击高价值目标,目标的价值收益包括目标的价值v、我方飞行器对杀伤概率PS,假设第ii架我方飞行器攻击敌方第jj架飞行器时计算公式为:
54)计算我方飞行器使用均衡度C4:假设m架我方飞行器攻击n个敌方飞行器,则计算公式为:
其中:Nii表示第ii架我方飞行器攻击敌方飞行器的架次,表示我方飞行器攻击敌方飞行器架次的平均数,我方飞行器使用均衡度越小,表示我方各型号飞行器使用越均衡;
55)建立离散粒子群算法中的适应度函数:根据我方飞行器的燃料消耗程度最小化指标,我方飞行器飞行时间的长度最小化指标,目标价值收益最大化指标,最小化我方飞行器使用均衡度,多飞行器协同目标分配的适应度函数C为:
C=ω1C1+ω2C2-ω3C3+ω4C4,
其中:ω1,ω2,ω3,ω4为权系数,且ω1=ω2=ω3=ω4=0.25;
所述采用离散粒子群算法对我方飞行器进行目标分配(13)步骤具体为:
61)设置种群大小、迭代次数;
62)随机初始化种群中各粒子的位置Xi和速度Vi:Xi=[xi,1,xi,2,…,xi,d],
Vi=[vi,1,vi,2,…,vi,d];其中:i为粒子的序号,d为搜索空间的维度,维度的数量为敌方飞行器的架数,d=1,2,…,n;
63)计算种群中每个粒子的适应度:评价种群中所有粒子,将当前各粒子的位置和目标值存储于各粒子的Pi=[pi,1,pi,2,…,pi,d]中,将所有Pi=[pi,1,pi,2,…,pi,d]目标值最优的个体位置和目标值存储于Pg=[pg,1,pg,2,…,pg,d]中;其中:Pi=[pi,1,pi,2,…,pi,d]为第i个粒子的最优值,Pg=[pg,1,pg,2,…,pg,d]为群体最优值;
64)按下式更新各个粒子的速度vi,j和位置xi,j:
其中:t_t表示粒子的迭代次数,ω=0.729为惯性权因子,b1=2.05和b2=2.05为正的加速常数,r1和r2为在0和1之间均匀分布的随机数,j表示维度;
对于m架我方飞行器攻击敌方n个目标,在第t_t次迭代时,第i个粒子的第j维的位置xi,j和速度vi,j取值范围为:
65)计算种群中每个粒子的适应度,评价种群中所有粒子;
66)更新最优值Pi:比较种群中每个粒子当前目标值与其最优值Pi的目标值,若当前目标值更优,则用粒子的当前位置和目标值更新最优值Pi;
67)更新最优值Pg:比较当前所有最优值Pi和最优值Pg的目标值,更新最优值Pg;
68)若迭代次数达到最大值,则输出最优值Pg及其目标值并结束,否则返回步骤64)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310361554.XA CN103413186B (zh) | 2013-08-19 | 2013-08-19 | 一种基于混合优化算法的多飞行器协同目标分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310361554.XA CN103413186B (zh) | 2013-08-19 | 2013-08-19 | 一种基于混合优化算法的多飞行器协同目标分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103413186A CN103413186A (zh) | 2013-11-27 |
CN103413186B true CN103413186B (zh) | 2016-01-20 |
Family
ID=49606193
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310361554.XA Expired - Fee Related CN103413186B (zh) | 2013-08-19 | 2013-08-19 | 一种基于混合优化算法的多飞行器协同目标分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103413186B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104090819B (zh) * | 2014-07-22 | 2017-07-18 | 中国科学院空间科学与应用研究中心 | 一种空间天文卫星巡天扫描多目标任务规划方法 |
CN105739303B (zh) * | 2015-12-29 | 2019-11-22 | 沈阳航空航天大学 | 基于滚动时域法的多无人机协同打击任务分配方法 |
CN105469139B (zh) * | 2016-01-11 | 2017-12-01 | 昆明理工大学 | 一种基于嵌入式控制的无人机空中实时协同制导方法 |
CN105975651B (zh) * | 2016-04-12 | 2018-09-07 | 哈尔滨工业大学 | 基于遗传粒子群多学科设计优化算法的导弹参数设计方法 |
CN105929848A (zh) * | 2016-06-28 | 2016-09-07 | 南京邮电大学 | 一种三维环境中的多无人机系统的航迹规划方法 |
CN106681358A (zh) * | 2017-03-03 | 2017-05-17 | 合肥工业大学 | 集中式无人机编队任务分配方法和装置 |
CN107037827B (zh) * | 2017-04-14 | 2020-02-21 | 合肥工业大学 | 无人机航空作业任务分配与航迹规划联合优化方法及装置 |
CN109640258B (zh) * | 2019-01-29 | 2021-02-19 | 广东交通职业技术学院 | 一种无线传感网络三维定位方法及装置 |
CN109976164B (zh) * | 2019-04-25 | 2022-11-29 | 南开大学 | 一种多旋翼无人机能量优化视觉覆盖轨迹规划方法 |
CN110163502B (zh) * | 2019-05-21 | 2021-03-05 | 北京壹氢科技有限公司 | 一种多弹协同多阶段目标分配方法 |
CN113031650B (zh) * | 2021-03-04 | 2022-06-10 | 南京航空航天大学 | 一种不确定环境下的无人机集群协同目标分配设计方法 |
CN113805569B (zh) * | 2021-09-23 | 2024-03-26 | 北京理工大学 | 基于多智能体技术的对抗系统、方法、终端及存储介质 |
CN114240170B (zh) * | 2021-12-17 | 2024-05-24 | 中国航空研究院 | 一种面向多域作战的武器目标分配方法 |
-
2013
- 2013-08-19 CN CN201310361554.XA patent/CN103413186B/zh not_active Expired - Fee Related
Non-Patent Citations (5)
Title |
---|
二级入轨空天飞机上升轨迹优化;黄国强 等;《宇航学报》;20100330;第31卷(第3期);641-647 * |
多UCAV协同任务分配模型及粒子群算法求解;杜继永 等;《控制与决策》;20121115;第27卷(第11期);1751-1755 * |
多UCAV协同目标分配算法研究;叶文 等;《系统工程与电子技术》;20100115;第32卷(第1期);105、107 * |
小推力深空探测轨道全局优化设计;黄国强 等;《航空学报》;20100725;第31卷(第7期);1366-1368、1371 * |
某型电视制导空地导弹仿真研究;陈蕾 等;《系统仿真技术》;20120728;第8卷(第3期);192-196 * |
Also Published As
Publication number | Publication date |
---|---|
CN103413186A (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103413186B (zh) | 一种基于混合优化算法的多飞行器协同目标分配方法 | |
CN112733421B (zh) | 一种针对有人无人机协同对地作战的任务规划方法 | |
CN106020230B (zh) | 一种能耗约束下的多无人机任务分配方法 | |
CN102506863B (zh) | 一种基于万有引力搜索的无人机航路规划方法 | |
CN103744290B (zh) | 一种多无人机编队分层目标分配方法 | |
CN110766254A (zh) | 一种基于改进遗传算法的多无人机协同任务分配方法 | |
CN103324993B (zh) | 一种基于多飞行器协同作战的轨迹优化方法 | |
CN105739303B (zh) | 基于滚动时域法的多无人机协同打击任务分配方法 | |
CN102436604B (zh) | 一种基于多目标进化方法的多弹协同航路计算方法 | |
CN103471592A (zh) | 一种基于蜂群协同觅食算法的多无人机航迹规划方法 | |
CN113392521B (zh) | 面向海空联合作战任务的资源编组模型构建方法及系统 | |
CN101908097B (zh) | 一种空战决策的粒子群优化方法 | |
CN108153328B (zh) | 一种基于分段贝塞尔曲线的多导弹协同航迹规划方法 | |
CN107832885A (zh) | 一种基于自适应迁移策略bbo算法的舰艇编队火力分配方法 | |
CN106529674A (zh) | 多无人机协同多目标分配方法 | |
CN105302153A (zh) | 异构多无人机协同察打任务的规划方法 | |
CN108958292A (zh) | 一种基于rrt*算法的飞行器突防轨迹规划方法 | |
Zeng et al. | Modeling for UAV resource scheduling under mission synchronization | |
CN106599537A (zh) | 一种基于多目标克隆进化算法的大规模武器‑目标分配方法 | |
CN110186328A (zh) | 一种多目标拦截的导弹分配方法及系统 | |
CN106681358A (zh) | 集中式无人机编队任务分配方法和装置 | |
CN111773722B (zh) | 一种模拟环境中的战斗机规避机动策略集生成方法 | |
CN111157002A (zh) | 基于多智能体进化算法的飞行器3d路径规划方法 | |
CN114815891A (zh) | 一种基于per-idqn的多无人机围捕战术方法 | |
Xiong et al. | Multi-uav 3d path planning in simultaneous attack |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 |