1 巡飞弹编组任务分配建模
1.1 问题描述
表1 符号说明Tab.1 Symbol Description |
| 类别 | 属性 |
|---|---|
| 巡飞弹 | 巡飞弹总数n |
| 巡飞弹位置Poin =(xi,yi) | |
| 巡飞速度Speedm | |
| 巡飞弹续航Durationi | |
| 巡飞弹i执行任务CK消耗的资源Consumeik | |
| 单巡飞弹对目标j最大侦察时间Timej | |
| 巡飞弹i健康状态Statusi | |
| 目标 | 目标数量m |
| 目标位置Poin =(xj,yj) | |
| 任务总个数C=2*m | |
| 目标价值Valuej |
1.2 主要任务分配指标
1.2.1 任务价值收益
1.2.2 任务损耗代价
1.3 目标函数
1.4 约束条件
2 自适应模拟退火遗传算法设计
2.1 编码
2.2 初始化
表2 初始化操作流程Tab.2 Initialization Operation Process |
| 初始化流程 |
|---|
| 步骤1:在目标集合中选择目标,每个目标有两次选中机会,分别表示侦察任务与打击任务:如果为侦察任务,跳转步骤2;如果为打击任务,跳转步骤3; |
| 步骤2:在侦察巡飞弹集合中选择执行侦察任务的巡飞弹,优先选择执行任务少的巡飞弹,跳转步骤4; |
| 步骤3:在打击巡飞弹集合中选择执行打击任务的巡飞弹,并在侦察巡飞弹集合和打击巡飞弹集合中删除该巡飞弹,跳转步骤5; |
| 步骤4:设置巡飞弹侦察阶段的时长,如果目标被侦察总时长与侦察所需时长相等,跳转步骤1;否则,跳转步骤2; |
| 步骤5:如果目标被打击总枚数与毁伤所需枚数相等,跳转步骤1;否则,跳转步骤3。 |
2.3 基于罚函数的适应度函数
2.4 自适应交叉、变异操作
表3 交叉操作流程Tab.3 Cross operation process |
| 交叉算子 |
|---|
| 步骤1:计算出最大个体适应度fitmax、种群平均适应度fitavg与代交叉的两个染色体中较大的适应度fitp |
| 步骤2:将待交叉染色体中较大的适应度与种群平均适应度比较,根据结果不同设置不同交叉概率,如下: PC= |
| 步骤3:结束。 |
表4 变异算子流程Tab.4 Mutation Operator Process |
| 变异算子 |
|---|
| 步骤1:算出最大个体适应度fitmax、种群平均适应度fitavg与代变异染色体的适应度fitp; |
| 步骤2:将待变异染色体的适应度与种群平均适应度比较,根据结果不同设置不同变异概率,如下: Pm= |
| 步骤3:结束。 |
2.5 模拟退火算法迭代
表5 模拟退火算法步骤Tab.5 Steps of Simulated Annealing Algorithm |
| 模拟退火算法 |
|---|
| 步骤1:在当前温度下得到有遗传算法优化的新解; |
| 步骤2:计算新解与上一温度的解之间的目标函数差异; |
| 步骤3:如果目标函数值改善,则接受新解作为当前解,执行步骤5;如果目标函数值较差,执行步骤4; |
| 步骤4:计算温度T下的概率P |
| P=eΔfit/KT-1 |
| 生成一个随机数a,如果a<P,那么接受这个解,否则放弃这个解; |
| 步骤5:降低温度,执行遗传算法下一次迭代。 |
3 案例分析
3.1 想定概述
表6 目标属性值Tab.6 Target Attribute Values |
| 目标编号 | 目标位置 | 侦察所需时长 | 目标毁伤概率 | 目标价值 |
|---|---|---|---|---|
| 1 | (10,15) | 360 | 0.82 | 10 |
| 2 | (13,14) | 420 | 0.9 | 9 |
| 3 | (5,6) | 400 | 0.71 | 7 |
| 4 | (5,10) | 300 | 0.9 | 10 |
| 5 | (10,10) | 500 | 0.82 | 8 |
| 6 | (12,7) | 450 | 0.71 | 9 |
| 7 | (11,15) | 450 | 0.82 | 8 |
3.2 结果分析
表7 巡飞弹编队任务分配结果Tab.7 Task allocation results of patrol missile formation |
| 侦察阶段 | 打击阶段 | |
|---|---|---|
| 巡飞弹1 | 侦察目标1(176秒) | 目标1 |
| 巡飞弹2 | 侦察目标2(121秒)→侦察目标7(26秒) | |
| 巡飞弹3 | 侦察目标1(119秒)→侦察目标6(40秒) | 目标4 |
| 巡飞弹4 | 侦察目标5(42秒)→侦察目标1(32秒)→侦察目标7(10秒) | |
| 巡飞弹5 | 侦察目标5(47秒)→侦察目标2(45秒)→侦察目标7(42秒) | 目标7 |
| 巡飞弹6 | 侦察目标5(125秒) | 目标5 |
| 巡飞弹7 | 侦察目标7(165秒) | |
| 巡飞弹8 | 侦察目标4(71秒)→侦察目标2(11秒)→侦察目标6(64秒) | 目标3 |
| 巡飞弹9 | 侦察目标5(37秒)→侦察目标6(101秒) | 目标2 |
| 巡飞弹10 | 侦察目标6(26秒)→侦察目标7(84秒) | |
| 巡飞弹11 | 侦察目标3(64秒)→侦察目标5(28秒)→侦察目标1(19秒)→侦察目标2(30秒) | 目标2 |
| 巡飞弹12 | 侦察目标5(146秒) | 目标5 |
| 巡飞弹13 | 侦察目标3(86秒)→侦察目标4(62秒) | 目标4 |
| 巡飞弹14 | 侦察目标3(107秒)→侦察目标4(54秒) | 目标4 |
| 巡飞弹15 | 侦察目标3(53秒)→侦察目标2(105秒) | 目标2 |
| 巡飞弹16 | 侦察目标4(63秒)→侦察目标6(86秒)→侦察目标7(24秒) | 目标7 |
| 巡飞弹17 | 侦察目标5(32秒)→侦察目标6(44秒)→侦察目标7(49秒) | 目标7 |
| 巡飞弹18 | 侦察目标3(41秒)→侦察目标4(48秒)→侦察目标6(85秒) | 目标6 |
| 巡飞弹19 | 侦察目标3(62秒)→侦察目标5(11秒)→侦察目标2(35秒)→侦察目标7(45秒) | |
| 巡飞弹20 | 侦察目标5(57秒)→侦察目标2(51秒)→侦察目标1(20秒) | 目标1 |
中国指挥与控制学会会刊 