CN107103164A - 无人机执行多任务的分配方法及装置 - Google Patents
无人机执行多任务的分配方法及装置 Download PDFInfo
- Publication number
- CN107103164A CN107103164A CN201710392402.4A CN201710392402A CN107103164A CN 107103164 A CN107103164 A CN 107103164A CN 201710392402 A CN201710392402 A CN 201710392402A CN 107103164 A CN107103164 A CN 107103164A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- unmanned plane
- chromosome
- flight path
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Abstract
本发明实施例公开了一种无人机执行多任务的分配方法及装置。方法包括:获取无人机和多个目标点的位置信息,以及无人机和风场的运动参数;根据位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群;根据初始种群、无人机和风场的运动参数确定无人机飞行状态和完成欧式飞行路径的航迹段的航行时间,根据航迹段的航行时间和SUAV‑VS‑EVRP模型获取初始种群中染色体对应的航行时间;基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的欧式飞行路径作为无人机的最优飞行路径。本发明实施例将无人机航迹规划问题与无人机实际飞行环境相结合,使规划得到的最优飞行路径方案优于无人机速度恒定的方案。
Description
技术领域
本发明实施例涉及无人机技术领域,具体涉及一种无人机执行多任务的分配方法及装置。
背景技术
当前,无人机UAV(Unmanned Aerial Vehicle)在军民领域有着广泛的应用,可完成目标侦察、目标跟踪、情报收集、震后救援和地质勘探等多种类型任务。例如在多架UAV协同侦察目标时,既要最合理地为每架UAV分配其所需侦察的目标,还要为其规划最优的飞行航迹。该问题是一个受多因素约束的任务分配与航迹规划联合优化问题,也是非确定性问题。
随着UAV研究的深入,环境因素被逐渐纳入问题的研究,特别是UAV任务分配、航迹规划和飞行控制等问题中,在环境因素的影响下如何降低耗能、控制UAV的飞行状态从而使UAV消耗最少的燃料执行最多的任务、具备更好的任务执行状态和更高的安全性是当前UAV研究的主要工作。当前常用于解决UAV任务分配与任务规划问题的模型有:TSP模型,TOP模型和VRP模型,其中,TSP模型是在只有单一旅行者的条件下,使得旅行者通过所有给定的目标点之后,从而使其路径成本最小的模型;TOP模型是在存在多个成员的条件下,使得每个成员尽可能访问更多的目标点,从而使得所有成员的总收益最大的模型;VRP模型是在车辆数量固定的条件下,使得车辆访问一定数量目标点,且在此过程中每个目标点只能被访问一次,最终使得UAV航行的总距离或总时间最短的模型。
在实现本发明实施例的过程中,发明人发现现有的技术方案在实际操作中,一般是假设模型中在恒定时间内无人机的速度是恒定的。然而这个假设显然是不现实的,导致模型无法精确模拟出无人机的实际运动状态,进而无法进行最优的航迹规划。
发明内容
本发明实施例的一个目的是解决现有技术由于在进行航迹规划是设定无人机的速度是恒定的,导致模型无法精确模拟出无人机的实际运动状态,进而无法给出的最优的航迹规划。
本发明实施例提出了一种无人机执行多任务的分配方法,包括:
S1、获取无人机和多个目标点的位置信息,以及所述无人机和风场的运动参数;
S2、根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群;
S3、根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间;
S4、基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径。
可选的,根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群包括:
根据预设遗传算法的编码方式进行染色体编码生成预定规模的初始种群;所述染色体由目标点信息组成;其中所述目标点属于集合T0表示UAVs的起点,NT表示目标点数量;
所述染色体为所述目标点的随机全排列。
可选的,根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间包括:
根据其目标点被访问顺序将每个染色体对应的欧式飞行路径分为多个航迹段;
根据每个航迹段对应的起始点的坐标以及终止点的坐标,结合风场参数确定无人机飞行状态,获取所述无人机完成所述航迹段的航行时间;
根据每个航迹段对应的航行时间获取所述染色体对应的航行时间。
可选的,根据所述每个航迹段对应的起始点的坐标以及终止点的坐标,结合风场参数确定无人机飞行状态,获取所述无人机完成所述航迹段的航行时间包括:
采用以下公式计算获取无人机由目标点Tj出发飞至目标点Tk航迹段的航行时间:
其中,Tj为起始点,Tk为终止点,T表示目标点的集合,Vg为无人机在上述两目标点间的地速;
采用以下公式计算获取无人机的地速:
其中,Va表示空速大小,βa表示空速航向角,Vg表示地速的大小,βg表示地速航向角,Vw表示风速大小,βw表示风向。
采用以下公式计算获取无人机在Tj和Tk两点间的欧氏距离:
其中,X,Y分别表示对应目标点横、纵坐标。
可选的,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间包括:
根据SUAV-VS-EVRP模型获取航行时间:
其约束条件为:
其中.表示无人机由目标点Tj出发飞至目标点Tk的航行时间,是一个二元决策变量,且当UAV经Tj飞行至Tk时,则的值为1,否则的值为0,NT表示目标点的数量。
可选的,基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径,包括:
步骤1、使用所述编码方法生成初始解,并生成预定规模的初始种群并根据种群中每个染色体对应的航行时间计算其适应度;
步骤2、使用轮盘赌方法选择父代种群中的两个个体(A,B)进行交叉,交叉规则为先随机选择个体A中交叉开始位置,并选取自交叉开始位置及其以后的基因为交叉基因,然后查找个体B中与个体A交叉基因对应目标值相同的基因,将染色体A中交叉基因的目标值排列顺序放入染色体B的对应基因并将染色体B中交叉位置基因的目标值排列顺序放入染色体A的对应基因得到新的染色体C和D,并判断染色体C、D的适应度是否优于A、B,并将适应度较优的两个染色体放入种群中,然后不断迭代更新步骤1种群,得到新的子代种群;
步骤3、使用轮盘赌方法选择步骤2种群中一条染色体进行变异,即进行目标点变异;
整个染色体变异的步骤包括:首先,若染色体A的第一行顺序变异,则随机选取当前染色体A一半基因位数量的基因位,并将当前选取基因位所对应的目标点编码进行随机全排列,并按照上述随机全排列的顺序依次放入对应的基因位完成变异得到染色体B,并判断染色体B的适应度是否优于A,并将适应度较优的染色体放入种群中,不断迭代更新步骤2种群,得到新的子代种群;
步骤4、计算子代种群适应度并选取本次迭代中所有解中的最优解;
步骤5、判断当前的迭代次数是否达到预设值,若判断否,则对步骤3中的子代种群和父代种群按照一定比例组合形成新的父代种群返回步骤2;若判断为是,则结束迭代,将最终获得的最优解作为无人机的任务分配与航迹规划结果。
本发明实施例提出了一种无人机执行多任务的分配装置,包括:
获取模块,用于获取无人机和多个目标点的位置信息,以及所述无人机和风场的运动参数;
第一处理模块,用于根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群;
第二处理模块,用于根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间;
第三处理模块,用于基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径。
可选的,所述第一处理模块,用于根据预设遗传算法的编码方式进行染色体编码生成预定规模的初始种群;所述染色体由目标点信息组成;其中所述目标点属于集合T0表示UAVs的起点,NT表示目标点数量;
所述染色体为所述目标点的随机全排列。
可选的,所述第二处理模块,用于根据其目标点被访问顺序将每个染色体对应的欧式飞行路径分为多个航迹段;根据每个航迹段对应的起始点的坐标以及终止点的坐标,结合风场参数确定无人机飞行状态,获取所述无人机完成所述航迹段的航行时间;根据每个航迹段对应的航行时间获取所述染色体对应的航行时间。
可选的,所述第二处理模块,用于采用以下公式计算获取无人机由目标点Tj出发飞至目标点Tk航迹段的航行时间:
其中,Tj为起始点,Tk为终止点,T表示目标点的集合,Vg为无人机在上述两目标点间的地速。
采用以下公式计算获取无人机的地速:
其中,Va表示空速大小,βa表示空速航向角,Vg表示地速的大小,βg表示地速航向角,Vw表示风速大小,βw表示风向。
采用以下公式计算获取无人机在Tj和Tk两点间的欧氏距离:
其中,X,Y分别表示对应目标点横、纵坐标。
由上述技术方案可知,本发明实施例提出的一种无人机执行多任务的分配方法及装置首先通过对风场和无人机的运动参数进行分析,获取无人机在风场中的实际飞行状态,然后基于实际飞行状态进行飞行路径的规划,与现有技术中设定无人机速度恒定的方案相比,能根据不确定环境中风场的状态精确计算无人机在所有可能飞行路径上的航行时间,进而选择出最优的飞行路径。
附图说明
通过参考附图会更加清楚的理解本发明的特征和优点,附图是示意性的而不应理解为对本发明进行任何限制,在附图中:
图1示出了本发明一实施例提供的一种无人机执行多任务的分配方法的流程示意图;
图2示出了本发明一实施例提供的计算飞行路径的航行时间的流程示意图;
图3示出了本发明一实施例提供的遗传算法的流程示意图;
图4a-图4c示出了本发明一实施例提供遗传算法中的算子的示意图;
图5示出了本发明一实施例提供的风向示意图;
图6示出了本发明一实施例提供的速度矢量关系示意图;
图7示出了本发明一实施例提供的UAV由S飞往A点受风场影响的分析示意图;
图8示出了本发明一实施例提供的对飞行路径进行分段的示意图;
图9示出了本发明一实施例提供的一种无人机执行多任务的航迹规划的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1示出了本发明一实施例提供的一种无人机执行多任务的航迹规划的流程示意图,参见图1,该方法可由处理器实现,具体包括如下步骤:
110、获取无人机和多个目标点的位置信息,以及所述无人机和风场的运动参数;
需要说明的是,在进行任务分配和航迹规划之前,技术人员可设定或者根据实际情况测出无人机和多个目标点的位置信息,然后将其输入至处理器中。
另外,无人机的运动参数可以是技术人员根据实际飞行需要设定的,风场的运动参数可以是技术人员测量得出或者是根据实际情况设定的。
120、根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群;
130、根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间;
140、基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径。
不难理解的是,每次交叉、变异的迭代可能都有新的个体的出现,然后基于步骤130对新的染色体进行的航行时间的计算,因此,每个欧式飞行路径对应一个航行时间。
可见,本实施例首先通过对风场和无人机的运动参数进行分析,获取无人机在风场中的实际飞行状态,然后基于实际飞行状态进行飞行路径的规划,与现有技术相比,本实施例将无人机航迹规划问题与无人机实际飞行环境相结合,使规划得到的最优飞行路径方案优于无人机速度恒定的方案,进而达到能精确计算无人机在所有可能飞行路径上的航行时间,进而选择出最优的飞行路径。
下面对本发明实施例中的各步骤进行详细说明:
首先,对步骤120进行详细说明:
根据预设遗传算法的编码方式进行染色体编码生成预定规模的初始种群;所述染色体由目标点信息组成;其中所述目标点属于集合T0表示UAVs的起点,NT表示目标点数量;
所述染色体为所述目标点的随机全排列。
然后,参见图2,下面对步骤130进行详细说明:
210、根据其目标点被访问顺序将每个染色体对应的欧式飞行路径分为多个航迹段;
220、根据每个航迹段对应的起始点的坐标以及终止点的坐标,结合风场参数确定无人机飞行状态,获取所述无人机完成所述航迹段的航行时间;
230、根据每个航迹段对应的航行时间获取所述染色体对应的航行时间。
其中,步骤220包括:
采用以下公式计算获取无人机由目标点Tj出发飞至目标点Tk航迹段的航行时间:
其中,Tj为起始点,Tk为终止点,T表示目标点的集合,Vg为无人机在上述两目标点间的地速;
采用以下公式计算获取无人机的地速:
其中,Va表示空速大小,βa表示空速航向角,Vg表示地速的大小,βg表示地速航向角,Vw表示风速大小,βw表示风向。
采用以下公式计算获取无人机在Tj和Tk两点间的欧氏距离:
其中,X,Y分别表示对应目标点横、纵坐标。
步骤230具体包括:
根据SUAV-VS-EVRP模型获取航行时间:
其约束条件为:
其中.表示无人机由目标点Tj出发飞至目标点Tk的航行时间,是一个二元决策变量,且当UAV经Tj飞行至Tk时,则的值为1,否则的值为0,NT表示目标点的数量。
下面对步骤140进行详细说明:
步骤1、使用所述编码方法生成初始解,并生成预定规模的初始种群并根据种群中每个染色体对应的航行时间计算其适应度;
步骤2、使用轮盘赌方法选择父代种群中的两个个体(A,B)进行交叉,交叉规则为先随机选择个体A中交叉开始位置,并选取自交叉开始位置及其以后的基因为交叉基因,然后查找个体B中与个体A交叉基因对应目标值相同的基因,将染色体A中交叉基因的目标值排列顺序放入染色体B的对应基因并将染色体B中交叉位置基因的目标值排列顺序放入染色体A的对应基因得到新的染色体C和D,并判断染色体C、D的适应度是否优于A、B,并将适应度较优的两个染色体放入种群中,然后不断迭代更新步骤1种群,得到新的子代种群;
步骤3、使用轮盘赌方法选择步骤2种群中一条染色体进行变异,即进行目标点变异;
整个染色体变异的步骤包括:首先,若染色体A的第一行顺序变异,则随机选取当前染色体A一半基因位数量的基因位,并将当前选取基因位所对应的目标点编码进行随机全排列,并按照上述随机全排列的顺序依次放入对应的基因位完成变异得到染色体B,并判断染色体B的适应度是否优于A,并将适应度较优的染色体放入种群中,不断迭代更新步骤2种群,得到新的子代种群;
步骤4、计算子代种群适应度并选取本次迭代中所有解中的最优解;
步骤5、判断当前的迭代次数是否达到预设值,若判断否,则对步骤3中的子代种群和父代种群按照一定比例组合形成新的父代种群返回步骤2;若判断为是,则结束迭代,将最终获得的最优解作为无人机的任务分配与航迹规划结果。
下面参见图3对本发明的采用的遗传算法的原理进行详细说明:
1、开启;
2、基于技术人员的设定,生成包括指定数量染色体的种群,指定数量可具体为100个;
其中,每个染色体与图1对应实施例中的步骤130中的欧式飞行路径一一对应。
需要说明的是,染色体编码代表了问题的一种可行的解决方案。
编码举例:
如图4a所示,染色体A表示在稳定风场下UAV访问五个目标点的一种可行方案,即一号UAV从起始点S(0,0)出发,依次访问目标点3、目标点1、目标点2、目标点4和目标点5。
3、计算每个染色体的适应度;
需要说明的是,采用图1对应实施例中的步骤140的计算方法,计算无人机完成每个欧式飞行路径的航行时间,并基于航行时间计算染色体的适应度,例如:航行时间与适应度成反比关系。
不难理解的是,按照上述步骤2中的编码方式生成规定数量的种群后进行适应度的计算,本文发明中适应度的计算以目标函数为依据,其计算过程如下:
4、选择操作
根据J’通过轮盘赌的方法进行选择操作。
5、交叉操作
通过对父代染色体进行交叉,可以继承父代中比较优良的基因,获得更优的子代。针对SUAV-VS-EVRP问题本文根据当前的编码方式采用以下方法进行交叉,随机选择个体A中交叉开始位置,并选取自交叉开始位置及其以后的基因为交叉基因,然后查找个体B中与个体A交叉基因对应目标值相同的基因,将染色体A中交叉基因的目标值排列顺序放入染色体B的对应基因并将染色体B中交叉位置基因的目标值排列顺序放入染色体A的对应基因得到新的染色体C和D。
交叉操作举例:
参见图4b,有父代Parent A和Parent B,在Parent A上随机产生进行交叉的基因位为3,并根据交叉基因位3-5上目标点的值找到Parent B上对应相同目标点的交叉基因位分别为1,3和5,接下来将Parent A中交叉基因位上目标点的排列顺序2→4→5替换ParentB中交叉基因位1,3,5上的值得到OffSpring B为:2→3→4→1→5。同样的,值得到OffSpring A为:3→1→4→5→2。
6、变异操作
变异是为了防止遗传算法陷入局部最优。针对SUAV-VS-EVRP问题本文根据当前的编码方式采用以下方法进行变异,若判断得到染色体A变异,则随机选取当前染色体A一半基因位数量的基因位,并将当前选取基因位所对应的目标点编码进行随机全排列,并按照上述随机全排列的顺序依次放入对应的基因位完成变异得到染色体B。
变异操作举例:
如图4c所示,有父代Parent A,在Parent A上进行目标点变异,在进行变异前首先判断变异是否发生,在判断得到目标点变异发生时,随机选取一半总基因位数量的变异基因位,本例中选取的变异基因位是1、2和4,随后将被选取的变异基因位上的目标值进行随机全排列,并用得到的随机全排列顺序替换原顺序得到OffSpring A为:4→3→2→1→5。
7、更新操作
8、选取最优分配方案
9、判断是否终止
10、获得最优分配方案
11、结束
需要说明的是,上述步骤与图1对应实施例中的部分步骤相对应,故,相似之处此处不再赘述,具体请查看图1对应的实施例中的相关内容。
下面结合上述的遗传算法对本发明的设计原理进行详细说明:
步骤一,为避免问题过于复杂,本发明采用区域固定风场进行风场建模,即在规定区域内,其风场的风速和风向是不变的。
已知区域的风场状态可表示为:
其中,Vw表示风场中的风速,βw表示风向。
风速Vw是指风相对于地面单位时间内移动的距离,单位为m/s;风向βw是指风吹来的方向,风向的测量单位一般用方位来表示,如陆地上,一般用16个方位表示,海上多用36个方位表示,而在高空则用角度表示,即把圆周分成360度,本文规定西风(W)是0度(即360度),南风(S)是90度,东风(E)是180度,北风(N)是270度,如图5所示。
步骤二,配置UAV
用U表示四旋翼UAV,UAV在空中的配置定义为:
q=(x,y,ψ) (4)
其中,
其中,和表示的是一架UAV在笛卡尔惯性参考系中的坐标;Vg表示UAV的地速ψ是指UAV的航向角。
为使问题简化,本文提出以下关于UAV在执行任务过程中需满足的运动约束的假设:
(3)考虑UAV在固定的高度飞行;
(4)根据UAV的飞行包线,UAV在指定高度固定载荷下的飞行速度存在上下界,即Va_min和Va_max分别表示在某高度下UAV空速的最小值和最大值;
(5)UAV由出发点出发并在执行完成任务后不返回出发点。
步骤三,计算UAV的实际飞行状态
考虑风影响的UAV实际速度定义为UAV的地速大小为Vg,此时UAV的航向角为βg,UAV地速矢量将不考虑风影响的UAV理论速度定义为UAV的空速大小为Va,此时UAV的航向角为βa,UAV空速矢量UAV空速地速与风场中风速的矢量关系如图6所示。
上述速度与角度关系为:
在无风时,即UAV空速与地速相等。
下面结合图7进行实例说明:
UAV由S(0,0)飞往A(50,300),空速为8m/s,该UAV所处的环境是风速为5m/s、风向为南风(Vw=5m/s,βw=90°),根据式(7)可得到Uav在该过程中的空速和地速如表4-1所示。
表4-1四旋翼Uav在无风与南风环境下空速、地速对比表
空速 | 地速 | |
南风环境 | 28.80km/h,35.5° | 23.49km/h,80.5° |
无风环境 | 28.80km/h,80.5° | 28.80km/h,80.5° |
步骤四,目标点配置
NT个目标点的集合可表示为:
其中,集合中所有的目标点的位置和任务量均已知。在本发明中,每一个目标点上都可能有不同类型的任务需要被UAV执行,且在此过程中每架UAV只能执行一个目标点上的一个任务,即每个目标点都要被不同的UAV访问,每架UAV只能访问某个目标点一次。
步骤五,计算航行时间
在以飞行时间作为目标的UAV任务分配与航迹规划问题中,UAV的任务分配方案决定UAV访问目标点的顺序,根据UAV目标点访问顺序进行航迹规划,由航迹规划的结果计算UAV飞行时间进而由UAV飞行时间决定当前UAV任务分配与航迹规划方案是否优于已知方案。
由于Uav在两点间的航行轨迹是欧氏距离,因而Uav在两点间的航行方向固定,进而固定风场下Uav在两点间的地速是不变的,但固定风场下Uav由目标点出发在多个目标点场景中地速是变化的,其飞行时间计算如下:
其中,
表示Uav在Tj和Tk两点间的欧氏距离,Uav在两点间的地速由公式(7)求得。
从而,可根据公式(8)计算UAV的欧式路径的航行时间。
其中,表示UAV在Tj、Tk两点的航行时间;
是一个二元决策变量,且当UAV经Tj飞行至Tk时,则的值为1,否则的值为0;
J中j值取0表示UAV由起始点出发。
在求解过程中,还需满足以下几个约束条件:
上述条件保证所有的目标点都能被访问到且只能被访问一次。
上述条件保证由起始点出发UAV数量的UAV路线,并有UAV数量的UAV路径指向同一点。
上述条件在其它约束条件的基础上保证有UAV数量的路线且每一架UAV的路径是一条闭合的环,即UAV的航行轨迹是一条有序的路线,并最终回到起始点。
可见,基于上式可得到每个飞行路径的航行时间,进而从中选取出航行时间最短的飞行路径。
下面对本发明进行具体实例的详细说明:
首先,所有的仿真实验均是在4G内存、3.4GHz CPU的硬件上、在Matlab R2014a的环境中运行的。具体说明如下:
UAV模型基于UAV的数学模型,其空速为8米/秒,两架UAV均从出发点S(0,0)起飞,在完成访问任务后返回点S(0,0);风场环境是固定风场,即在一次实验过程中风速和风向都是不变的,并且为了保证UAV能够安全飞行,风速大小为5米/秒,风向取东,即180°,UAV需要访问的三个目标点坐标分别为:A(100,300)、B(200,150)、C(350,50)、D(500,150)和E(650,100)。
根据上述本发明提出的模型和算法,本文在东风风场环境和试验场景下进行实验,并得到各风场环境下无人机飞行时间最短的任务分配与航迹规划方案如表3-1所示(参见图8)。
3-1
风场环境 | 无人机任务分配与航迹规划方案 | 飞行时间 |
东风环境 | U:S→A→B→C→D→E | 235.6427 |
图9示出了本发明一实施例提供的一种无人机执行多任务的分配装置的结构示意图,参见图9,该装置包括:获取模块101、第一处理模块102、第二处理模块103以及第三处理模块104,其中:
获取模块101,用于获取无人机和多个目标点的位置信息,以及所述无人机和风场的运动参数;
第一处理模块102,用于根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群;
第二处理模块103,用于根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间;
第三处理模块104,用于基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径。
可见,本实施例首先通过对风场和无人机的运动参数进行分析,获取无人机在风场中的实际飞行状态,然后基于实际飞行状态进行飞行路径的规划,与现有技术相比,本实施例将无人机航迹规划问题与无人机实际飞行环境相结合,使规划得到的最优飞行路径方案优于无人机速度恒定的方案,进而达到能精确计算无人机在所有可能飞行路径上的航行时间,进而选择出最优的飞行路径。
下面对本装置的各功能模块进行详细说明:
第一处理模块102用于根据预设遗传算法的编码方式进行染色体编码生成预定规模的初始种群;所述染色体由目标点信息组成;其中所述目标点属于集合T0表示UAVs的起点,NT表示目标点数量;
所述染色体为所述目标点的随机全排列。
第二处理模块103,用于根据其目标点被访问顺序将每个染色体对应的欧式飞行路径分为多个航迹段;根据每个航迹段对应的起始点的坐标以及终止点的坐标,结合风场参数确定无人机飞行状态,获取所述无人机完成所述航迹段的航行时间。
采用以下公式计算获取无人机由目标点Tj出发飞至目标点Tk航迹段的航行时间:
其中,Tj为起始点,Tk为终止点,T表示目标点的集合,Vg为无人机在上述两目标点间的地速。
采用以下公式计算获取无人机的地速:
其中,Va表示空速大小,βa表示空速航向角,Vg表示地速的大小,βg表示地速航向角,Vw表示风速大小,βw表示风向。
采用以下公式计算获取无人机在Tj和Tk两点间的欧氏距离:
其中,X,Y分别表示对应目标点横、纵坐标。
第三处理模块104,用于执行以下步骤:步骤1、使用所述编码方法生成初始解,并生成预定规模的初始种群并根据种群中每个染色体对应的航行时间计算其适应度;
步骤2、使用轮盘赌方法选择父代种群中的两个个体(A,B)进行交叉,交叉规则为先随机选择个体A中交叉开始位置,并选取自交叉开始位置及其以后的基因为交叉基因,然后查找个体B中与个体A交叉基因对应目标值相同的基因,将染色体A中交叉基因的目标值排列顺序放入染色体B的对应基因并将染色体B中交叉位置基因的目标值排列顺序放入染色体A的对应基因得到新的染色体C和D,并判断染色体C、D的适应度是否优于A、B,并将适应度较优的两个染色体放入种群中,然后不断迭代更新步骤1种群,得到新的子代种群;
步骤3、使用轮盘赌方法选择步骤2种群中一条染色体进行变异,即进行目标点变异;
整个染色体变异的步骤包括:首先,若染色体A的第一行顺序变异,则随机选取当前染色体A一半基因位数量的基因位,并将当前选取基因位所对应的目标点编码进行随机全排列,并按照上述随机全排列的顺序依次放入对应的基因位完成变异得到染色体B,并判断染色体A的适应度是否优于B,并将适应度较优的两个染色体放入种群中,不断迭代更新步骤2种群,得到新的子代种群;
步骤4、计算子代种群适应度并选取本次迭代中所有解中的最优解;
步骤5、判断当前的迭代次数是否达到预设值,若判断否,则对步骤3中的子代种群和父代种群按照一定比例组合形成新的父代种群返回步骤2;若判断为是,则结束迭代,将最终获得的最优解作为无人机的任务分配与航迹规划结果。
对于装置实施方式而言,由于其与方法实施方式基本相似,所以描述的比较简单,相关之处参见方法实施方式的部分说明即可。
应当注意的是,在本发明的装置的各个部件中,根据其要实现的功能而对其中的部件进行了逻辑划分,但是,本发明不受限于此,可以根据需要对各个部件进行重新划分或者组合。
本发明的各个部件实施方式可以以硬件实现,或者以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本装置中,PC通过实现因特网对设备或者装置远程控制,精准的控制设备或者装置每个操作的步骤。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如,计算机程序和计算机程序产品)。这样实现本发明的程序可以存储在计算机可读介质上,并且程序产生的文件或文档具有可统计性,产生数据报告和cpk报告等,能对功放进行批量测试并统计。应该注意的是上述实施方式对本发明进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施方式。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中,这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种无人机执行多任务的分配方法,其特征在于,包括:
S1、获取无人机和多个目标点的位置信息,以及所述无人机和风场的运动参数;
S2、根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群;
S3、根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间;
S4、基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径。
2.根据权利要求1所述的方法,其特征在于,根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群包括:
根据预设遗传算法的编码方式进行染色体编码生成预定规模的初始种群;所述染色体由目标点信息组成;其中所述目标点属于集合T0表示UAVs的起点,NT表示目标点数量;
所述染色体为所述目标点的随机全排列。
3.根据权利要求1或2所述的方法,其特征在于,根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间包括:
根据其目标点被访问顺序将每个染色体对应的欧式飞行路径分为多个航迹段;
根据每个航迹段对应的起始点的坐标以及终止点的坐标,结合风场参数确定无人机飞行状态,获取所述无人机完成所述航迹段的航行时间;
根据每个航迹段对应的航行时间获取所述染色体对应的航行时间。
4.根据权利要求3所述方法,其特征在于,根据所述每个航迹段对应的起始点的坐标以及终止点的坐标,结合风场参数确定无人机飞行状态,获取所述无人机完成所述航迹段的航行时间包括:
采用以下公式计算获取无人机由目标点Tj出发飞至目标点Tk航迹段的航行时间:
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Length</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<msub>
<mi>V</mi>
<mi>g</mi>
</msub>
</mfrac>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
其中,Tj为起始点,Tk为终止点,T表示目标点的集合,Vg为无人机在上述两目标点间的地速;
采用以下公式计算获取无人机的地速:
<mrow>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<msub>
<mi>&beta;</mi>
<mi>g</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&beta;</mi>
<mi>g</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<msub>
<mi>V</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&beta;</mi>
<mi>a</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&beta;</mi>
<mi>w</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&beta;</mi>
<mi>a</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>sin&beta;</mi>
<mi>w</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>a</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>w</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Va表示空速大小,βa表示空速航向角,Vg表示地速的大小,βg表示地速航向角,Vw表示风速大小,βw表示风向;
采用以下公式计算获取无人机在Tj和Tk两点间的欧氏距离:
<mrow>
<msub>
<mi>Length</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>q</mi>
<mi>r</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Y</mi>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
其中,X,Y分别表示对应目标点横、纵坐标。
5.根据权利要求3所述方法,其特征在于,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间包括:
根据SUAV-VS-EVRP模型获取航行时间:
<mrow>
<mi>J</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<msub>
<mi>t</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
</mrow>
其约束条件为:
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
其中.表示无人机由目标点Tj出发飞至目标点Tk的航行时间,是一个二元决策变量,且当UAV经Tj飞行至Tk时,则的值为1,否则的值为0,NT表示目标点的数量。
6.根据权利要求4所述方法,其特征在于,基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径,包括:
步骤1、使用所述编码方法生成初始解,并生成预定规模的初始种群并根据种群中每个染色体对应的航行时间计算其适应度;
步骤2、使用轮盘赌方法选择父代种群中的两个个体(A,B)进行交叉,交叉规则为先随机选择个体A中交叉开始位置,并选取自交叉开始位置及其以后的基因为交叉基因,然后查找个体B中与个体A交叉基因对应目标值相同的基因,将染色体A中交叉基因的目标值排列顺序放入染色体B的对应基因并将染色体B中交叉位置基因的目标值排列顺序放入染色体A的对应基因得到新的染色体C和D,并判断染色体C、D的适应度是否优于A、B,并将适应度较优的两个染色体放入种群中,然后不断迭代更新步骤1种群,得到新的子代种群;
步骤3、使用轮盘赌方法选择步骤2种群中一条染色体进行变异,即进行目标点变异;
整个染色体变异的步骤包括:首先,若染色体A的第一行顺序变异,则随机选取当前染色体A一半基因位数量的基因位,并将当前选取基因位所对应的目标点编码进行随机全排列,并按照上述随机全排列的顺序依次放入对应的基因位完成变异得到染色体B,并判断染色体B的适应度是否优于A,并将适应度较优的染色体放入种群中,不断迭代更新步骤2种群,得到新的子代种群;
步骤4、计算子代种群适应度并选取本次迭代中所有解中的最优解;
步骤5、判断当前的迭代次数是否达到预设值,若判断否,则对步骤3中的子代种群和父代种群按照一定比例组合形成新的父代种群返回步骤2;若判断为是,则结束迭代,将最终获得的最优解作为无人机的任务分配与航迹规划结果。
7.一种无人机执行多任务的分配装置,其特征在于,包括:
获取模块,用于获取无人机和多个目标点的位置信息,以及所述无人机和风场的运动参数;
第一处理模块,用于根据所述无人机和所述多个目标点的位置信息和预设遗传算法,构建以欧式飞行路径为个体的初始种群;
第二处理模块,用于根据所述初始种群、无人机和风场的运动参数确定无人机飞行状态和无人机完成欧式飞行路径的航迹段的航行时间,根据所述航迹段的航行时间和SUAV-VS-EVRP模型获取初始种群中染色体对应的航行时间;
第三处理模块,用于基于遗传算法,对初始种群中染色体进行交叉、变异处理,并在达到预定迭代次数后,选取航行时间最短的染色体对应的欧式飞行路径作为所述无人机的最优飞行路径。
8.根据权利要求7所述的装置,其特征在于,所述第一处理模块,用于根据预设遗传算法的编码方式进行染色体编码生成预定规模的初始种群;所述染色体由目标点信息组成;其中所述目标点属于集合T0表示UAVs的起点,NT表示目标点数量;
所述染色体为所述目标点的随机全排列。
9.根据权利要求7或8所述的装置,其特征在于,所述第二处理模块,用于根据其目标点被访问顺序将每个染色体对应的欧式飞行路径分为多个航迹段;执行第一步骤和第二步骤;
所述第一步骤包括:采用以下公式计算获取无人机由目标点Tj出发飞至目标点Tk航迹段的航行时间:
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Length</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<msub>
<mi>V</mi>
<mi>g</mi>
</msub>
</mfrac>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
其中,Tj为起始点,Tk为终止点,T表示目标点的集合,Vg为无人机在上述两目标点间的地速;
采用以下公式计算获取无人机的地速:
<mrow>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<msub>
<mi>&beta;</mi>
<mi>g</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&beta;</mi>
<mi>g</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<msub>
<mi>V</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&beta;</mi>
<mi>a</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&beta;</mi>
<mi>w</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&beta;</mi>
<mi>a</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>sin&beta;</mi>
<mi>w</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>a</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>w</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Va表示空速大小,βa表示空速航向角,Vg表示地速的大小,βg表示地速航向角,Vw表示风速大小,βw表示风向;
采用以下公式计算获取无人机在Tj和Tk两点间的欧氏距离:
<mrow>
<msub>
<mi>Length</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>q</mi>
<mi>r</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Y</mi>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
其中,X,Y分别表示对应目标点横、纵坐标。
所述第二步骤包括:
根据SUAV-VS-EVRP模型获取航行时间:
<mrow>
<mi>J</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<msub>
<mi>t</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
</mrow>
其约束条件为:
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</munderover>
<msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>&Element;</mo>
<mi>T</mi>
</mrow>
其中.表示无人机由目标点Tj出发飞至目标点Tk的航行时间,是一个二元决策变量,且当UAV经Tj飞行至Tk时,则的值为1,否则的值为0,NT表示目标点的数量。
10.根据权利要求9所述方法,其特征在于,所述第三处理模块用于执行如下步骤:步骤1、使用所述编码方法生成初始解,并生成预定规模的初始种群并根据种群中每个染色体对应的航行时间计算其适应度;步骤2、使用轮盘赌方法选择父代种群中的两个个体(A,B)进行交叉,交叉规则为先随机选择个体A中交叉开始位置,并选取自交叉开始位置及其以后的基因为交叉基因,然后查找个体B中与个体A交叉基因对应目标值相同的基因,将染色体A中交叉基因的目标值排列顺序放入染色体B的对应基因并将染色体B中交叉位置基因的目标值排列顺序放入染色体A的对应基因得到新的染色体C和D,并判断染色体C、D的适应度是否优于A、B,并将适应度较优的两个染色体放入种群中,然后不断迭代更新步骤1种群,得到新的子代种群;步骤3、使用轮盘赌方法选择步骤2种群中一条染色体进行变异,即进行目标点变异;整个染色体变异的步骤包括:首先,若染色体A的第一行顺序变异,则随机选取当前染色体A一半基因位数量的基因位,并将当前选取基因位所对应的目标点编码进行随机全排列,并按照上述随机全排列的顺序依次放入对应的基因位完成变异得到染色体B,并判断染色体B的适应度是否优于A,并将适应度较优的染色体放入种群中,不断迭代更新步骤2种群,得到新的子代种群;步骤4、计算子代种群适应度并选取本次迭代中所有解中的最优解;步骤5、判断当前的迭代次数是否达到预设值,若判断否,则对步骤3中的子代种群和父代种群按照一定比例组合形成新的父代种群返回步骤2;若判断为是,则结束迭代,将最终获得的最优解作为无人机的任务分配与航迹规划结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710392402.4A CN107103164B (zh) | 2017-05-27 | 2017-05-27 | 无人机执行多任务的分配方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710392402.4A CN107103164B (zh) | 2017-05-27 | 2017-05-27 | 无人机执行多任务的分配方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107103164A true CN107103164A (zh) | 2017-08-29 |
CN107103164B CN107103164B (zh) | 2020-08-28 |
Family
ID=59659980
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710392402.4A Active CN107103164B (zh) | 2017-05-27 | 2017-05-27 | 无人机执行多任务的分配方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107103164B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107807665A (zh) * | 2017-11-29 | 2018-03-16 | 合肥工业大学 | 无人机编队探测任务协同分配方法及装置 |
CN108235743A (zh) * | 2017-11-24 | 2018-06-29 | 深圳市道通智能航空技术有限公司 | 任务执行方法、装置、可移动物体及计算机可读存储介质 |
CN110046857A (zh) * | 2019-04-23 | 2019-07-23 | 南京邮电大学 | 一种基于遗传算法的无人机物流系统及配送方法 |
CN111047087A (zh) * | 2019-09-18 | 2020-04-21 | 合肥工业大学 | 无人机与车辆协同下路径的智能优化方法和装置 |
CN111399533A (zh) * | 2020-02-10 | 2020-07-10 | 合肥工业大学 | 异构多无人机协同任务分配与路径优化方法 |
CN111832725A (zh) * | 2019-04-15 | 2020-10-27 | 电子科技大学 | 一种基于改进遗传算法的多机器人多任务分配方法及装置 |
CN111860984A (zh) * | 2020-07-06 | 2020-10-30 | 燕山大学 | 一种基于改进遗传算法的多异构无人机任务分配方法 |
CN112069731A (zh) * | 2020-09-02 | 2020-12-11 | 中国航空无线电电子研究所 | 一种基于遗传算法的飞行仿真机动控制优化方法 |
CN112580893A (zh) * | 2020-12-29 | 2021-03-30 | 浙江量大智能科技有限公司 | 一种基于改进遗传算法的多无人机大气监测路径规划方法 |
CN112966361A (zh) * | 2020-12-29 | 2021-06-15 | 中国计量大学 | 一种面向高架和非高架点源的多无人机监测路径优化方法 |
CN115879603A (zh) * | 2022-11-17 | 2023-03-31 | 武汉大学 | 一种面向多目标点的多无人机协同数据采集方法及装置 |
CN116011695A (zh) * | 2023-03-27 | 2023-04-25 | 湖南胜云光电科技有限公司 | 一种获取无人机目标路径的数据处理系统 |
CN116245346A (zh) * | 2023-05-12 | 2023-06-09 | 深圳大学 | 基于多种群遗传算法和局部搜索的多无人机任务分配方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105929848A (zh) * | 2016-06-28 | 2016-09-07 | 南京邮电大学 | 一种三维环境中的多无人机系统的航迹规划方法 |
CN106600147A (zh) * | 2016-12-15 | 2017-04-26 | 合肥工业大学 | 面向可分解任务的多无人机任务分配方法和装置 |
-
2017
- 2017-05-27 CN CN201710392402.4A patent/CN107103164B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105929848A (zh) * | 2016-06-28 | 2016-09-07 | 南京邮电大学 | 一种三维环境中的多无人机系统的航迹规划方法 |
CN106600147A (zh) * | 2016-12-15 | 2017-04-26 | 合肥工业大学 | 面向可分解任务的多无人机任务分配方法和装置 |
Non-Patent Citations (4)
Title |
---|
JOSE ALFREDO GUERRERO等: "UAV Path Planning for Structure Inspection in Windy Environments", 《JOURNAL OF INTELLIGENT & ROBOTIC SYSTEMS》 * |
MANGAL KOTHARI等: "UAV Path Following in Windy Urban Environments", 《JOURNAL OF INTELLIGENT & ROBOTIC SYSTEMS》 * |
屈耀红等: "基于风场信息的无人机在线航迹规划方法", 《西北工业大学学报》 * |
王国强等: "无人机编队协同任务规划仿真系统研究", 《系统仿真学报》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108235743A (zh) * | 2017-11-24 | 2018-06-29 | 深圳市道通智能航空技术有限公司 | 任务执行方法、装置、可移动物体及计算机可读存储介质 |
WO2019100326A1 (zh) * | 2017-11-24 | 2019-05-31 | 深圳市道通智能航空技术有限公司 | 任务执行方法、装置、可移动物体及计算机可读存储介质 |
US10705524B2 (en) | 2017-11-24 | 2020-07-07 | Autel Robotics Co., Ltd. | Task execution method and device, moveable object and computer readable storage medium |
CN108235743B (zh) * | 2017-11-24 | 2021-08-06 | 深圳市道通智能航空技术股份有限公司 | 任务执行方法、装置、可移动物体及计算机可读存储介质 |
CN107807665A (zh) * | 2017-11-29 | 2018-03-16 | 合肥工业大学 | 无人机编队探测任务协同分配方法及装置 |
CN111832725A (zh) * | 2019-04-15 | 2020-10-27 | 电子科技大学 | 一种基于改进遗传算法的多机器人多任务分配方法及装置 |
CN111832725B (zh) * | 2019-04-15 | 2023-05-26 | 电子科技大学 | 一种基于改进遗传算法的多机器人多任务分配方法及装置 |
CN110046857A (zh) * | 2019-04-23 | 2019-07-23 | 南京邮电大学 | 一种基于遗传算法的无人机物流系统及配送方法 |
CN110046857B (zh) * | 2019-04-23 | 2022-08-30 | 南京邮电大学 | 一种基于遗传算法的无人机物流系统及配送方法 |
CN111047087A (zh) * | 2019-09-18 | 2020-04-21 | 合肥工业大学 | 无人机与车辆协同下路径的智能优化方法和装置 |
CN111399533A (zh) * | 2020-02-10 | 2020-07-10 | 合肥工业大学 | 异构多无人机协同任务分配与路径优化方法 |
CN111399533B (zh) * | 2020-02-10 | 2022-07-26 | 合肥工业大学 | 异构多无人机协同任务分配与路径优化方法 |
CN111860984A (zh) * | 2020-07-06 | 2020-10-30 | 燕山大学 | 一种基于改进遗传算法的多异构无人机任务分配方法 |
CN112069731A (zh) * | 2020-09-02 | 2020-12-11 | 中国航空无线电电子研究所 | 一种基于遗传算法的飞行仿真机动控制优化方法 |
CN112580893A (zh) * | 2020-12-29 | 2021-03-30 | 浙江量大智能科技有限公司 | 一种基于改进遗传算法的多无人机大气监测路径规划方法 |
CN112966361A (zh) * | 2020-12-29 | 2021-06-15 | 中国计量大学 | 一种面向高架和非高架点源的多无人机监测路径优化方法 |
CN112966361B (zh) * | 2020-12-29 | 2023-07-18 | 中国计量大学 | 一种面向高架和非高架点源的多无人机监测路径优化方法 |
CN115879603A (zh) * | 2022-11-17 | 2023-03-31 | 武汉大学 | 一种面向多目标点的多无人机协同数据采集方法及装置 |
CN116011695A (zh) * | 2023-03-27 | 2023-04-25 | 湖南胜云光电科技有限公司 | 一种获取无人机目标路径的数据处理系统 |
CN116245346A (zh) * | 2023-05-12 | 2023-06-09 | 深圳大学 | 基于多种群遗传算法和局部搜索的多无人机任务分配方法 |
CN116245346B (zh) * | 2023-05-12 | 2023-07-28 | 深圳大学 | 基于多种群遗传算法和局部搜索的多无人机任务分配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107103164B (zh) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107238388B (zh) | 多无人机任务分配与航迹规划联合优化方法及装置 | |
CN107169608A (zh) | 多无人机执行多任务的分配方法及装置 | |
CN107103164A (zh) | 无人机执行多任务的分配方法及装置 | |
CN107145161A (zh) | 无人机访问多目标点的航迹规划方法及装置 | |
Lamont et al. | UAV swarm mission planning and routing using multi-objective evolutionary algorithms | |
CN106843235A (zh) | 一种面向无人自行车的人工势场路径规划法 | |
Hong et al. | Energy-efficient online path planning of multiple drones using reinforcement learning | |
CN109508035A (zh) | 基于分布式控制的多区域分级式无人机编队路径规划方法 | |
CN106843236A (zh) | 基于权重改进粒子群算法的无人自行车路径规划方法 | |
Rossmann et al. | Virtual robotic testbeds: A foundation for e-robotics in space, in industry-and in the woods | |
CN109947126A (zh) | 四旋翼无人机的控制方法、装置、设备及可读介质 | |
CN104570736B (zh) | 一种星‑臂耦合系统的动力学参数在轨辨识方法和装置 | |
CN107633105B (zh) | 一种基于改进混合蛙跳算法的四旋翼无人机参数辨识方法 | |
CN114840020A (zh) | 一种基于改进鲸鱼算法的无人机飞行轨迹规划方法 | |
CN108549210A (zh) | 基于bp神经网络pid控制的多无人机协同飞行方法 | |
CN115562357B (zh) | 一种面向无人机集群的智能路径规划方法 | |
CN111813146A (zh) | 基于bp神经网络预测航程的再入预测-校正制导方法 | |
CN109870906A (zh) | 一种基于bbo优化人工势场的高速旋翼飞行器路径规划方法 | |
CN113848974A (zh) | 一种基于深度强化学习的飞行器轨迹规划方法及系统 | |
CN108008645A (zh) | 六自由度仿真建模方法 | |
Nguyen et al. | Nonmyopic planning for long-term information gathering with an aerial glider | |
Becce et al. | Optimal path planning for autonomous spraying uas framework in precision agriculture | |
Hohmann et al. | Three-dimensional urban path planning for aerial vehicles regarding many objectives | |
CN116560384A (zh) | 一种基于深度强化学习的变体飞行器鲁棒控制方法 | |
CN112161626B (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 |