专利摘要

专利摘要

本发明公开了一种基于遗传算法的舰船对岸打击火力分配方法，包括：以一艘舰船上的一种武器作为火力单元，将火力单元对岸上目标的火力分配数目作为编码位，采用随机分布生成法对每个火力单元进行编码，得到初始种群；利用适应度函数计算初始种群的种群适应度，根据种群适应度结合轮盘赌对初始种群进行选择，得到选择后种群，对选择后种群进行染色体交叉和变异，得到新的种群；利用新的种群更新初始种群后重复上一步，重复多次后，得到舰船对岸打击火力分配方案。本发明求解结果很稳定，耗时短、收敛速度快、初始种群生成容易。

权利要求

1.一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，包括：

(1)以一艘舰船上的一种武器作为火力单元，将火力单元对岸上目标的火力分配数目作为编码位，采用随机分布生成法对每个火力单元进行编码，得到初始种群；

(2)利用适应度函数计算初始种群的种群适应度，根据种群适应度结合轮盘赌对初始种群进行选择，得到选择后种群，对选择后种群进行染色体交叉和变异，得到新的种群；

(3)利用新的种群更新初始种群后重复步骤(2)，重复多次后，得到舰船对岸打击火力分配方案。

2.如权利要求1所述的一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，所述步骤(1)包括：

(1-1)获取舰船和岸上目标的相关信息，所述相关信息包括：参战每一艘舰船的经纬度、舰船载有武器的种类、每种武器的数目、每一个岸上目标的战略价值和武器对岸上目标的命中毁伤概率；

(1-2)以一艘舰船上的一种武器作为火力单元，利用相关信息得到火力单元的地理位置、武器数量、武器射程和火力单元对岸上目标的毁伤概率矩阵；

(1-3)将火力单元对岸上目标的火力分配数目作为编码位，采用随机分布生成法对每个火力单元进行编码，得到初始种群，根据每一个岸上目标的战略价值和火力单元对岸上目标的毁伤概率矩阵，以岸上目标毁伤最大为目的，建立适应度函数。

3.如权利要求2所述的一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，所述步骤(1-2)包括：

以一艘舰船上的一种武器作为火力单元，利用相关信息得到火力单元的地理位置、武器数量、武器射程；

所述地理位置为火力单元所在的经纬度，其中，东经为正，西经为负，经度的取值范围为(-180°，180°)，北纬为正，南纬为负，纬度的取值范围为(0°，180°)；

所述火力单元的数量为N，岸上目标的数量为M，火力单元i与岸上目标j的距离的Lij：

C＝sin(lati)*sin(latj)*cos(loni-lonj)+cos(lati)*cos(latj)

Lij＝R*Arccos(C)*Pi/180

其中，lati表示第i个火力单元的经度，loni表示第i个火力单元的纬度，火力单元i的武器射程为ri，对于岸上目标，武器射程如果大于等于火力单元与岸上目标之间的距离，即Lij≤ri，则火力单元对岸上目标的毁伤概率为pij，否则火力单元的毁伤概率为0，所有火力单元对岸上目标的毁伤概率pij构成火力单元对岸上目标的毁伤概率矩阵P，i＝1，2，......，N；j＝1，2，......，M。

4.如权利要求2或3所述的一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，所述编码为对每个火力单元对M个岸上目标中的(M-1)个岸上目标进行编码。

5.如权利要求4所述的一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，所述编码的具体实现方式为：

对于火力单元i，火力单元的数量为N，岸上目标的数量为M，火力单元i装载的导弹数量为mj，火力单元i有(M-1)个编码位，当xij≥0时，以随机分布生成法生成随机排列{a1，a2，......，aM}，其中as，s≤M为基因位置的随机排列，按照随机排列进行火力单元对应基因的初始化，其中，rand(k)为产生(0-k)之间的整数，若as＝M，即xiM＝rand(mi-sum)，xiM不出现在编码位中，对每一个火力单元进行上述编码操作，生成具有随机分布的初始种群。

6.如权利要求5所述的一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，所述适应度函数为：

S(x)＝F(x)*G(x)

其中，S(x)为适应度函数，F(x)为初始适应度函数，G(x)为指示函数，pij为火力单元i对岸上目标j的毁伤概率，N为火力单元的数量，M为岸上目标的数量，约束条件为wj为岸上目标j的战略价值，wM为岸上目标M的战略价值，piM为火力单元i对岸上目标M的毁伤概率，mi表示火力单元i装载的导弹数量。

7.如权利要求1或2所述的一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，所述步骤(2)包括：

利用适应度函数计算初始种群的种群适应度为Pu为种群数量，累加种群适应度得到左间隔点总适应度为随机产生Pu个(0-ga)大小的随机数，将相邻的左间隔点作为一个区间，将随机数置于对应区间，得到选择后种群；

设置交换概率为Pc，将选择后的种群中的个体随机成对作为父代，对父代以Pc概率进行交换，产生随机数r1∈[0，1]，如果r1＜Pc则发生交换，从父代的染色体中选取随机一点为交叉点，交换交叉点之后的基因，获得子代；如果父代未发生交换，则由父代直接产生子代；设置变异概率为Pm，对子代中的每一个体，以Pm概率进行变异，产生随机数r2∈[0，1]，如果r2＜Pm则发生变异，对该个体的染色体中的任意一个基因位进行变异，得到新的种群。

8.如权利要求1或2所述的一种基于遗传算法的舰船对岸打击火力分配方法，其特征在于，所述步骤(3)包括：

利用新的种群更新初始种群后重复步骤(2)，重复多次，在每一代迭代过程中，记录该次迭代时的最大适应度值Sbest、迭代次数bestq和最优编码方案Xbest，Sbest＝S(Xbest)，对于当前迭代次数cur，若当前迭代时的适应度值Scur＝S(Xcur)，Scur＞Sbest，更新当前迭代下的最大适应度值Sbest、迭代次数bestq和最优编码方案Xbest，即Xbest＝Xcur，Sbest＝S(Xcur)，bestq＝cur，若cur-bestq＞预设迭代次数，那么就终止遗传算法，得到目标编码方案，通过目标编码方案得到舰船对岸打击火力分配方案。

说明书

技术领域

本发明属于调度分配领域，更具体地，涉及一种基于遗传算法的舰船对岸打击火力分配方法。

背景技术

火力分配问题为当今作战系统工程的研究热点问题，由于该种问题为NP完全问题，没有多项式复杂度算法完成最优解的求取。舰船对岸打击火力分配问题，就是对于舰队编队，利用舰载的多种武器对岸上多个目标进行打击，获取较优作战效果的武器分配方案。获取较优解，可以用到各种现代数值求解算法例如传统遗传算法、蛙跳算法、模拟退火算法等。

在传统遗传算法求取火力分配方案，是一种可行的方案，定义适应度函数后可以获得较好的武器分配方案。但是传统遗传算法求解结果很不稳定，编码位过多会造成求解时间长、收敛速度慢、初始种群生成困难等问题。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提供了一种基于遗传算法的舰船对岸打击火力分配方法，由此解决传统遗传算法求解结果很不稳定，编码位过多会造成求解时间长、收敛速度慢、初始种群生成困难的技术问题。

为实现上述目的，本发明提供了一种基于遗传算法的舰船对岸打击火力分配方法，包括：

(1)以一艘舰船上的一种武器作为火力单元，将火力单元对岸上目标的火力分配数目作为编码位，采用随机分布生成法对每个火力单元进行编码，得到初始种群；

(2)利用适应度函数计算初始种群的种群适应度，根据种群适应度结合轮盘赌对初始种群进行选择，得到选择后种群，对选择后种群进行染色体交叉和变异，得到新的种群；

(3)利用新的种群更新初始种群后重复步骤(2)，重复多次后，得到舰船对岸打击火力分配方案。

进一步地，步骤(1)包括：

(1-1)获取舰船和岸上目标的相关信息，所述相关信息包括：参战每一艘舰船的经纬度、舰船载有武器的种类、每种武器的数目、每一个岸上目标的战略价值和武器对岸上目标的命中毁伤概率；

(1-2)以一艘舰船上的一种武器作为火力单元，利用相关信息得到火力单元的地理位置、武器数量、武器射程和火力单元对岸上目标的毁伤概率矩阵；

(1-3)将火力单元对岸上目标的火力分配数目作为编码位，采用随机分布生成法对每个火力单元进行编码，得到初始种群，根据每一个岸上目标的战略价值和火力单元对岸上目标的毁伤概率矩阵，以岸上目标毁伤最大为目的，建立适应度函数。

进一步地，步骤(1-2)包括：

以一艘舰船上的一种武器作为火力单元，利用相关信息得到火力单元的地理位置、武器数量、武器射程；

所述地理位置为火力单元所在的经纬度，其中，东经为正，西经为负，经度的取值范围为(-180°，180°)，北纬为正，南纬为负，纬度的取值范围为(0°，180°)；

所述火力单元的数量为N，岸上目标的数量为M，火力单元i与岸上目标j的距离的Lij：

C＝sin(lati)*sin(latj)*cos(loni-lonj)+cos(lati)*cos(latj)

Lij＝R*Arccos(C)*Pi/180

其中，lati表示第i个火力单元的经度，loni表示第i个火力单元的纬度，火力单元i的武器射程为ri，对于岸上目标，武器射程如果大于等于火力单元与岸上目标之间的距离，即Lij≤ri，则火力单元对岸上目标的毁伤概率为pij，否则火力单元的毁伤概率为0，所有火力单元对岸上目标的毁伤概率pij构成火力单元对岸上目标的毁伤概率矩阵P，i＝1，2，......，N；j＝1，2，......，M。

进一步地，编码为对每个火力单元对M个岸上目标中的(M-1)个岸上目标进行编码。

进一步地，编码的具体实现方式为：

对于火力单元i，火力单元的数量为N，岸上目标的数量为M，火力单元i装载的导弹数量为mi，火力单元i有(M-1)个编码位，当 xij≥0时，以随机分布生成法生成随机排列{a1，a2，......，aM}，其中as，s≤M为基因位置的随机排列，按照随机排列进行火力单元对应基因的初始化， 其中， rand(k)为产生(0-k)之间的整数，若as＝M，即xiM＝rand(mi-sum)，XiM不出现在编码位中，对每一个火力单元进行上述编码操作，生成具有随机分布的初始种群。

进一步地，适应度函数为：

S(x)＝F(x)*G(x)

其中，S(x)为适应度函数，F(x)为初始适应度函数，G(x)为指示函数，pij为火力单元i对岸上目标j的毁伤概率，N为火力单元的数量，M为岸上目标的数量，约束条件为 i＝1，2，......，N，xij≥0，j＝1，2，......，M-1，wj为岸上目标j的战略价值，wM为岸上目标M的战略价值，piM为火力单元i对岸上目标M的毁伤概率，mi表示火力单元i装载的导弹数量。

进一步地，步骤(2)包括：

利用适应度函数计算初始种群的种群适应度为 Pu为种群数量，累加种群适应度得到左间隔点 总适应度为 随机产生Pu个(0-ga)大小的随机数，将相邻的左间隔点作为一个区间，将随机数置于对应区间，得到选择后种群；

设置交换概率为Pc，将选择后的种群中的个体随机成对作为父代，对父代以Pc概率进行交换，产生随机数r1∈[0，1]，如果r1＜Pc则发生交换，从父代的染色体中选取随机一点为交叉点，交换交叉点之后的基因，获得子代；如果父代未发生交换，那么由父代直接产生子代；设置变异概率为Pm，对子代中的每一个体，以Pm概率进行变异，产生随机数r2∈[0，1]，如果r2＜Pm则发生变异，对该个体染色体中的任意一个基因位进行变异，得到新的种群。

进一步地，步骤(3)包括：

利用新的种群更新初始种群后重复步骤(2)，重复多次，在每一代迭代过程中，记录该次迭代时的最大适应度值Sbest、迭代次数bestq和最优编码方案Xbest，Sbest＝S(Xbest)，对于当前迭代次数cur，若当前迭代时的适应度值Scur＝S(Xcur)，Scur＞Sbest，更新当前迭代下的最大适应度值Sbest、迭代次数bestq和最优编码方案Xbest，即Xbest＝Xcur，Sbest＝S(Xcur)，bestq＝cur，若cur-bestq＞预设迭代次数，那么就终止遗传算法，得到目标编码方案，通过目标编码方案得到舰船对岸打击火力分配方案。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：

(1)本发明采用随机分布生成法对每个火力单元进行编码，得到初始种群；这样生成的初始种群很容易，同时本发明使用的改进的遗传算法为利用适应度函数计算初始种群的种群适应度，根据种群适应度结合轮盘赌对初始种群进行选择，得到选择后种群，对选择后种群进行染色体交叉和变异，得到新的种群；使得对此迭代后的求解结果很稳定。

(2)本发明采用减少一位编码的方式进行编码，使得迭代过程耗时短、收敛速度快，本发明的适应度函数充分考虑了减少一位编码的操作，同时结合约束函数，使得最后得到的舰船对岸打击火力分配方案准确、稳定性强。

附图说明

图1是本发明实施例提供的一种基于遗传算法的舰船对岸打击火力分配方法的流程图；

图2是本发明实施例提供的遗传算法的流程图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

如图1所示，一种基于遗传算法的舰船对岸打击火力分配方法，包括：

(1)以一艘舰船上的一种武器作为火力单元，将火力单元对岸上目标的火力分配数目作为编码位，采用随机分布生成法对每个火力单元进行编码，得到初始种群；

(2)利用适应度函数计算初始种群的种群适应度，根据种群适应度结合轮盘赌对初始种群进行选择，得到选择后种群，对选择后种群进行染色体交叉和变异，得到新的种群；

(3)利用新的种群更新初始种群后重复步骤(2)，重复多次后，得到舰船对岸打击火力分配方案。

如图2所示，步骤(1)包括：

(1-1)获取舰船和岸上目标的相关信息，所述相关信息包括：参战每一艘舰船的经纬度、舰船载有武器的种类、每种武器的数目、每一个岸上目标的战略价值和武器对岸上目标的命中毁伤概率；

(1-2)以一艘舰船上的一种武器作为火力单元，利用相关信息得到火力单元的地理位置、武器数量、武器射程；

所述地理位置为火力单元所在的经纬度，其中，东经为正，西经为负，经度的取值范围为(一180°，180°)，北纬为正，南纬为负，纬度的取值范围为(0°，180°)；

所述火力单元的数量为N，岸上目标的数量为M，火力单元i与岸上目标j的距离的Lij：

C＝sin(lati)*sin(latj)*cos(loni-lonj)+cos(lati)*cos(latj)

Lij＝R*Arccos(C)*Pi/180

其中，lati表示第i个火力单元的经度，loni表示第i个火力单元的纬度，火力单元i的武器射程为ri，对于岸上目标，武器射程如果大于等于火力单元与岸上目标之间的距离，即Lij≤ri，则火力单元对岸上目标的毁伤概率为pij，否则火力单元的毁伤概率为0，所有火力单元对岸上目标的毁伤概率pij构成火力单元对岸上目标的毁伤概率矩阵P，i＝1，2，......，N；j＝1，2，......，M。

(1-3)将火力单元对岸上目标的火力分配数目作为编码位，采用随机分布生成法对每个火力单元进行编码，得到初始种群，对于火力单元i，火力单元的数量为N，岸上目标的数量为M，火力单元i装载的导弹数量为mi，火力单元i有(M-1)个编码位，当 xij≥0时，以随机分布生成法生成随机排列{a1，a2，......，aM}，其中as，s≤M为基因位置的随机排列，按照随机排列进行火力单元对应基因的初始化， 其中， rand(k)为产生(0-k)之间的整数，若as＝M，即xiM＝rand(mi-sum)，xiM不出现在编码位中，对每一个火力单元进行上述编码操作，生成具有随机分布的初始种群，设置种群容量为Pu＝100。

根据每一个岸上目标的战略价值和火力单元对岸上目标的毁伤概率矩阵，以岸上目标毁伤最大为目的，建立适应度函数。

适应度函数为：

S(x)＝F(x)*G(x)

其中，S(x)为适应度函数，F(x)为初始适应度函数，G(x)为指示函数，pij为火力单元i对岸上目标j的毁伤概率，N为火力单元的数量，M为岸上目标的数量，约束条件为 i＝1，2，......，N，xij≥0，j＝1，2，......，M-1，wj为岸上目标j的战略价值，wM为岸上目标M的战略价值，piM为火力单元i对岸上目标M的毁伤概率，mi表示火力单元i装载的导弹数量。

步骤(2)包括：

利用适应度函数计算初始种群的种群适应度为 Pu为种群数量，累加种群适应度得到左间隔点 总适应度为 随机产生Pu个(0-ga)大小的随机数，将相邻的左间隔点作为一个区间，将随机数置于对应区间，得到选择后种群；具有高适应度的解具有更高的几率被选择到，用于后面的交叉和变异过程中，以生成子代。种群容量保持Pu＝100。

设置交换概率为Pc＝0.95，将选择后的种群中的个体随机成对作为父代，对父代以Pc概率进行交换，产生随机数r1∈[0，1]，如果r1＜Pc则发生交换，从父代的染色体中选取随机一点为交叉点，交换交叉点之后的基因，获得子代；如果父代未发生交换，那么由父代直接产生子代；设置变异概率为Pm，对子代中的每一个体，以Pm＝0.05概率进行变异，产生随机数r2∈[0，1]，如果r2＜Pm则发生变异，对该个体染色体中的任意一个基因位进行变异，得到新的种群。

具体地：

交叉为一点交叉。

变异：采用一点变异策略，对于要变异的个体xu，随机选取两点1≤i1≤N，1≤j1≤M-1，则i1对应的基因 那么更新 这里rand(y)为取0-y的整数。

步骤(3)包括：

利用新的种群更新初始种群后重复步骤(2)，重复多次，在每一代迭代过程中，记录该次迭代时的最大适应度值Sbest、迭代次数bestq和最优编码方案Xbest，Sbest＝S(Xbest)，对于当前迭代次数cur，若当前迭代时的适应度值Scur＝S(Xcur)，Scur＞Sbest，更新当前迭代下的最大适应度值Sbest、迭代次数bestq和最优编码方案Xbest，即Xbest＝Xcur，Sbest＝S(Xcur)，bestq＝cur，若cur-bestq＞100，那么就终止遗传算法，得到目标编码方案，通过目标编码方案得到舰船对岸打击火力分配方案。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

一种基于遗传算法的舰船对岸打击火力分配方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。