CN107037827B - Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device - Google Patents

Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device Download PDF

Info

Publication number
CN107037827B
CN107037827B CN201710245797.5A CN201710245797A CN107037827B CN 107037827 B CN107037827 B CN 107037827B CN 201710245797 A CN201710245797 A CN 201710245797A CN 107037827 B CN107037827 B CN 107037827B
Authority
CN
China
Prior art keywords
unmanned aerial
farmland
aerial vehicle
time
chromosomes
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.)
Active
Application number
CN201710245797.5A
Other languages
Chinese (zh)
Other versions
CN107037827A (en
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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201710245797.5A priority Critical patent/CN107037827B/en
Publication of CN107037827A publication Critical patent/CN107037827A/en
Application granted granted Critical
Publication of CN107037827B publication Critical patent/CN107037827B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Genetics & Genomics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Artificial Intelligence (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Physiology (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention relates to a method and a device for joint optimization of task allocation and track planning of aerial work of an unmanned aerial vehicle. The method provided by the invention can automatically obtain the task and the flight path planning of the unmanned aerial vehicle, so that the unmanned aerial vehicle can automatically execute the operation task according to the task and the flight path planning, and the influence of manual operation is avoided. In addition, the unmanned aerial vehicle based on this result execution job task also can be in the operation of limited duration of flight to the polylith farmland and obtain the biggest gross profit and can improve the efficiency of operation effectively.

Description

Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device
Technical Field
The embodiment of the invention relates to the technical field of unmanned aerial vehicles, in particular to a method and a device for joint optimization of unmanned aerial vehicle aerial work task allocation and flight path planning.
Background
With the continuous deepening of the agricultural mechanization degree, the unmanned aerial vehicle rapidly becomes an important mode in the agricultural operation process with the advantages of high operation efficiency, low labor intensity, low comprehensive cost and the like, and has wide application in agricultural aviation operations such as precision seeding, vegetation detection, pesticide spraying and the like. For example, the germination condition and the weed degree of herbaceous plants can be detected by an unmanned aerial vehicle, or the planthoppers can be controlled by spraying pesticides on rice fields by using the unmanned aerial vehicle, and the like.
Current drones can be broadly classified into two broad categories, multi-rotor (e.g., quad-rotor, six-rotor, or eight-rotor drones, etc.) and fixed-wing. Wherein fixed wing unmanned aerial vehicle is comparatively wide application in agricultural operation with advantages such as long, the area of cruising is big, flying speed is fast, height.
However, in the process of implementing the present invention, the inventor finds that, because the current operation of the fixed-wing unmanned aerial vehicle is mainly performed by artificial remote control, the effect of the actual operation is greatly influenced by the operation level of the operator, and the flight path planned by the artificial and instant looking mode deviates from the theoretical flight path seriously, which leads to the high operation leakage rate and repetition rate of the unmanned aerial vehicle. And when multiple operators operate multiple fixed wing drones, there is also a lack of coordination before each other.
In addition, in the process of carrying out pesticide spraying aviation operation by using the fixed-wing unmanned aerial vehicle, a plurality of factors influencing operation execution are provided, such as the time for executing work, the flight time of the unmanned aerial vehicle and the like. For example, when the working time of a day is divided into two sections, how to spray pesticides on a plurality of farmlands by a plurality of fixed-wing unmanned aerial vehicles under the condition of limited flight time and how to carry out the most reasonable flight path planning makes the total income of the farmlands after spraying maximum (that is, the total sum of pesticide effects of the farmlands completing the spraying task maximum) become a problem to be solved urgently.
Disclosure of Invention
An embodiment of the invention provides a method and a device for joint optimization of mission allocation and flight path planning of unmanned aerial vehicle aerial work, which are used for solving the problems that in the prior art, the flight of the unmanned aerial vehicle is greatly influenced by manual operation, and when the working time of one day is divided into two sections, how to obtain the maximum total profit for operation of a plurality of farmlands by utilizing a plurality of fixed-wing unmanned aerial vehicles within the limited flight time length is high.
In a first aspect, an embodiment of the present invention provides a joint optimization method for mission allocation and flight path planning of unmanned aerial vehicle, where when multiple fixed-wing unmanned aerial vehicles perform a mission on multiple rectangular farmlands, the method includes:
acquiring two time windows for executing the task, information of the fixed-wing unmanned aerial vehicle and farmland information of the candidate farmlands;
coding the time window information, the unmanned aerial vehicle information and the farmland information, and randomly generating a plurality of chromosomes; screening out chromosomes meeting preset constraint conditions defined by a preset MUAV-MTW-DTOP model from the plurality of chromosomes generated randomly, and constructing an initial solution set of the MUAV-MTW-DTOP model; the MUAV-MTW-DTOP model is an objective function which enables the multi-frame fixed-wing unmanned aerial vehicle flying in a preset flying mode to obtain the maximum total benefit in the task; the preset constraints comprise flight time constraints and two time window constraints of each fixed-wing unmanned aerial vehicle;
and continuously crossing and mutating the initial solution set based on a preset genetic algorithm until the initial solution set meets the constraint of the iteration times, ending the crossing and the mutation, selecting the optimal solution which enables the model to obtain the maximum total income from the mutated solution set, and taking the optimal solution as the task allocation and track planning result of the operation.
In a second aspect, a further embodiment of the present invention provides a joint optimization apparatus for mission allocation and flight path planning for unmanned aerial vehicles, which includes, when multiple fixed-wing unmanned aerial vehicles perform a mission on multiple rectangular fields:
the information acquisition unit is used for acquiring two time windows for executing the task, information of the fixed-wing unmanned aerial vehicle and farmland information of the candidate farmlands;
the initial solution acquisition unit is used for encoding the time window information, the information of the unmanned aerial vehicle and the farmland information and randomly generating a plurality of chromosomes; screening out chromosomes meeting preset constraint conditions defined by a preset MUAV-MTW-DTOP model from the plurality of chromosomes generated randomly, and constructing an initial solution set of the MUAV-MTW-DTOP model; the MUAV-MTW-DTOP model is an objective function which enables the multi-frame fixed-wing unmanned aerial vehicle flying in a preset flying mode to obtain the maximum total benefit in the task; the preset constraints comprise flight time constraints and two time window constraints of each fixed-wing unmanned aerial vehicle;
and the optimal solution calculation unit is used for continuously crossing and mutating the initial solution set based on a preset genetic algorithm until the crossing and the mutation are finished after the constraint of the iteration times is met, selecting an optimal solution which enables the model to obtain the maximum total income from the mutated solution set, and taking the optimal solution as a task allocation and flight path planning result of the operation.
The invention provides a joint optimization method for task allocation and flight path planning of unmanned aerial vehicle aerial work, aiming at the condition that a plurality of fixed-wing unmanned aerial vehicles execute a task on a plurality of rectangular farmlands, the method comprises the steps of firstly obtaining two time windows for executing the task, information of the fixed-wing unmanned aerial vehicles and farmland information of a plurality of candidate farmlands, then obtaining an optimal solution which can enable the model to obtain the maximum total income according to the information based on a preset MUAV-MTW-DTOP model and a genetic algorithm, and taking the optimal solution as a task allocation and flight path planning result of the operation. Compared with the existing manual remote control mode, the method provided by the invention can automatically obtain the task and flight path planning of the unmanned aerial vehicle in the operation according to the preset model and algorithm, so that the unmanned aerial vehicle can automatically execute the operation task according to the task and flight path planning, and the influence of manual operation is avoided. In addition, the method provided by the invention takes the optimal solution of the preset farmland maximization profit model as a flight path planning result, so that the unmanned aerial vehicle executing the operation task based on the result can operate a plurality of farmlands within the limited flight time and obtain the maximum total profit, thereby effectively improving the operation efficiency and enabling the operation form of the unmanned aerial vehicle to be applied to wider farmland operation tasks.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a flowchart of an embodiment of a joint optimization method for task allocation and flight path planning in unmanned aerial vehicle aerial work provided by the invention;
FIGS. 2(a) -2(c) are schematic diagrams of a farmland spraying mode provided by an embodiment of the invention;
fig. 3 is a schematic view of a flight path between farmlands of an unmanned aerial vehicle under a dynamic constraint condition, provided by an embodiment of the invention;
FIG. 4 is a flow chart of a genetic algorithm provided by an embodiment of the present invention;
FIG. 5 is an example of chromosome coding provided by an embodiment of the invention;
FIG. 6 is an example of chromosome crossing provided by an embodiment of the invention;
FIG. 7 is an example of chromosomal variation provided by an embodiment of the invention;
FIG. 8 is a schematic diagram of the distribution of 20 farmland provided by an embodiment of the invention;
FIG. 9 is a schematic diagram of the convergence of the optimal solution provided by the embodiment of the present invention;
fig. 10 is a schematic structural diagram of an embodiment of an apparatus for joint optimization of unmanned aerial vehicle aerial work task allocation and flight path planning, provided by the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In a first aspect, an embodiment of the present invention provides a method for joint optimization of mission allocation and flight path planning for unmanned aerial vehicle aerial work, where when multiple fixed-wing unmanned aerial vehicles perform a mission on multiple rectangular farmlands, as shown in fig. 1, the method includes:
s101, obtaining two time windows for executing the task, information of the fixed-wing unmanned aerial vehicle and farmland information of the candidate farmlands;
s102, coding the time window information, the unmanned aerial vehicle information and the farmland information, and randomly generating a plurality of chromosomes; screening out chromosomes meeting preset constraint conditions defined by a preset MUAV-MTW-DTOP model from the plurality of chromosomes generated randomly, and constructing an initial solution set of the MUAV-MTW-DTOP model; the MUAV-MTW-DTOP model is an objective function which enables the multi-frame fixed-wing unmanned aerial vehicle flying in a preset flying mode to obtain the maximum total benefit in the task; the preset constraints comprise flight time constraints and two time window constraints of each fixed-wing unmanned aerial vehicle;
s103, continuously crossing and mutating the initial solution set based on a preset genetic algorithm until the initial solution set meets the constraint of iteration times, ending crossing and mutating, selecting an optimal solution which enables the model to obtain the maximum total income from the mutated solution set, and taking the optimal solution as a task allocation and flight path planning result of the operation.
The invention provides a joint optimization method for task allocation and flight path planning of unmanned aerial vehicle aerial work, aiming at the condition that a plurality of fixed-wing unmanned aerial vehicles execute a task on a plurality of rectangular farmlands, the method comprises the steps of firstly obtaining two time windows for executing the task, information of the fixed-wing unmanned aerial vehicles and farmland information of a plurality of candidate farmlands, then obtaining an optimal solution which can enable the model to obtain the maximum total income according to the information based on a preset MUAV-MTW-DTOP model and a genetic algorithm, and taking the optimal solution as a task allocation and flight path planning result of the operation. Compared with the existing manual remote control mode, the method provided by the invention can automatically obtain the task and flight path planning of the unmanned aerial vehicle in the operation according to the preset model and algorithm, so that the unmanned aerial vehicle can automatically execute the operation task according to the task and flight path planning, and the influence of manual operation is avoided. In addition, the method provided by the invention takes the optimal solution of the preset farmland maximization profit model as a flight path planning result, so that the unmanned aerial vehicle executing the operation task based on the result can operate a plurality of farmlands within the limited flight time and obtain the maximum total profit, thereby effectively improving the operation efficiency and enabling the operation form of the unmanned aerial vehicle to be applied to wider farmland operation tasks.
In practical implementation, it can be understood that the objective function and the constraint condition included in the MUAV-MTW-DTOP model in the above method are important bases for obtaining the optimal planning result of the present invention, and may be set in various ways, and one of the alternative setting ways is described in detail below.
(1) Unmanned plane
Denoted by U is N for executing a task to be sprayedUSet of isomorphic drones
Figure BDA0001270604300000061
Each unmanned aerial vehicle can only carry one pesticide; during flight, all drones have the same minimum turning radius RUAnd a flying speed V, and all carry a spraying radius RDThe spray head of (1).
Considering the characteristic that the unmanned aerial vehicle executes the pesticide spraying task, the embodiment of the invention makes the following assumptions:
(a) the unmanned aerial vehicles have the capability of automatically avoiding obstacles, and can adopt an autonomous avoidance control strategy under the condition of collision, so that the generated path deviation is very small and negligible relative to the total flight path length;
(b) the unmanned aerial vehicles all fly at the same cruising speed and cruising height, so that the influence of the factors on the spraying effect is not considered;
(c) influence of external environment on the flight track of the unmanned aerial vehicle is not considered in the flight process of the unmanned aerial vehicle;
(d) the unmanned aerial vehicle can carry pesticides required by executing tasks in the flight process, but the fuel is limited;
(2) farmland
Is provided with
Figure BDA0001270604300000062
The starting point and the end point of the unmanned aerial vehicle are respectively, and the starting point and the end point are the same in the embodiment of the invention;
Figure BDA0001270604300000063
for N to be sprayed with pesticidesARectangular farmland, and farmland AiIs side length of Lix,LiyThe vertex coordinate is (A)i1,Ai2,Ai3,Ai4) Area is SiA rectangle of (2); the set of the starting point, the terminal point and the farmland area of the unmanned aerial vehicle is
Figure BDA0001270604300000064
When unmanned plane UuFor farmland AiWhen the pesticide is sprayed In the coverage mode, the entry point of the unmanned aerial vehicle flying into the farmland is IniuThe departure point of the flying away from the farmland is OutiuAnd assume that the drone must completely spray the entire field before leaving. Meanwhile, the pesticide can be sprayed on each farmland only once at most.
(3) Flight path of preset flight mode
In the process of executing agricultural aviation operation tasks by the unmanned aerial vehicle, the agricultural operation tasks are required to be completed by spraying the agricultural material in a covering mode in farmlands, and the unmanned aerial vehicle also needs to fly among different farmlands to realize switching among the tasks, so that two types of flight paths are generated, namely flight paths in farmlands and in farmlands.
Flight path inside the field: in farmland AiInside, unmanned aerial vehicle flies under the dynamics constraint and will have turning radius to use the parallel scanning strategy to carry out the overlay pesticide and spray. In the process, the unmanned plane UuFrom farmland AiIn (2) ofiuPoint entering, the path after entering the farmland is parallel to a certain edge of the farmland, and then from OutiuPoint away, at this moment, unmanned plane UuThe spent time for spraying the farmland is tiwu
The parallel scanning strategy is an optimal strategy for the unmanned aerial vehicle to carry out coverage type pesticide spraying on a rectangular farmland area, is also the most conveniently realized strategy in control, and is widely applied to the execution process of an area coverage task. This strategy has two common implementations, namely flight parallel to the short sides of the rectangle and flight parallel to the long sides of the matrix. For example, there are two ways to overlay scan the rectangular area in fig. 2(a) using a parallel scanning strategy, with the flight trajectory of the drone parallel to the short side of the rectangle (as shown in fig. 2 (b)) and the flight trajectory of the drone parallel to the long side of the rectangle (as shown in fig. 2 (c)). Meanwhile, the position of the unmanned aerial vehicle entering the rectangular area and the actual turning radius of the unmanned aerial vehicle on the rectangular boundary can influence the length of the flight path of the unmanned aerial vehicle in the farmland at the same time. Although the unmanned aerial vehicle can enter from any point of the boundary of the farmland, the entry point when the unmanned aerial vehicle covers the shortest flight path of the rectangular area is that the distance from the vertex of the farmland is RDPoint (2) of (c). Therefore, the embodiment of the invention discretizes the entry point of the unmanned aerial vehicle into the rectangular area, and selects the distance R between the 4 vertexes of the rectangle on the long side and the short side of the rectangle D8 points of (1) { R ]D1,RD2...RD8As an entry point for the drone into the rectangular area. For a rectangular area, the shortest path for a drone to fly inside the area must enter the rectangular area from one of the 8 entry points, and the exit point Out of the rectangular area for the drone to exit on this shortest pathiuI.e. uniquely determined.
Flight path between farmlands:
the unmanned aerial vehicle finishes the execution of the farmland AiAfter the pesticide spraying task, if the farmland A needs to be sprayedjSpraying pesticide, then the unmanned aerial vehicle must also be in two farmlands Ai,AjFly along the Dubins path. The starting point of the path is unmanned plane UuIn farmland AiOut ofiuThe termination point is the unmanned plane UuIn farmland AjIn (b) at the entry point InjuThe time taken by the Dubins path is denoted as tijwu
According to the generation principle of the Dubins path, the shortest Dubins path between two points can be generated by combining an arc-segment path and a straight-segment path, and there are six cases where D ═ { LSL, RSR, RSL, LSR, RLR, LRL }. Wherein, L represents a section of arc that unmanned aerial vehicle turned to the left, and R represents a section of arc that unmanned aerial vehicle turned to the right, and S represents that unmanned aerial vehicle flies with the straight line mode.
For example, fig. 3 depicts a flight path between fields for a drone under dynamic constraints. According to the generation principle of the Dubins path, the unmanned plane follows A of the rectangular areaiThe point starts to turn right and fly along an arc path (R), then fly along a straight path (S), and finally turn left and fly along an arc path (L) to reach A of a rectangular areajThe path of the point. Because a parallel scanning strategy is adopted in the farmland, the angle of the unmanned aerial vehicle entering or leaving a rectangular farmland area is necessarily the boundary vertical to the farmland, so that in the optimization process of the shortest Dubins path, the starting angle of the starting point and the entering angle of the ending point of the unmanned aerial vehicle are determined, and the changing factor is that the U of the unmanned aerial vehicle is a UuIn farmland AiOut ofiuAnd field AjIn (b) at the entry point Inju
In order to complete the task allocation problem, the embodiment of the invention selects all the profits of all the unmanned aerial vehicles when completing the task as the maximized benefit function.
The objective function of the MUAV-MTW-DTOP model is as follows:
Figure BDA0001270604300000081
the constraint conditions of the MUAV-MTW-DTOP model are as follows:
Figure BDA0001270604300000083
Figure BDA0001270604300000084
Figure BDA0001270604300000085
Figure BDA0001270604300000091
Figure BDA0001270604300000092
Figure BDA0001270604300000093
wherein N isUFor unmanned plane UuThe number of the cells; n is a radical ofAFor farmland AiThe number of the cells; w is the w-th time window, and is 1 or 2; n th0,NA+1The block farmland represents the starting point and the end point of the unmanned aerial vehicle; SQiIndicating farmland AiArea (unit hectare); piwIndicating that field A is completed in the w-th time windowiThe benefits of the task, i.e., the profits; siwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time to start spraying the pesticide; t is tijwuRefers to unmanned plane UuFor farmland A in the w time windowi,AjThe time of flying according to a preset flying mode; sjwuRefers to unmanned plane UuFor farmland A in the w time windowjThe time to start spraying the pesticide; m is a preset value; o isiwu,CiwuAre respectively unmanned aerial vehicle UuFor farmland A in the w time windowiA start time and an end time at which a task can be executed; t is tiwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time for executing the task according to a preset flight mode; ewuRepresenting the maximum flight duration of the drone. If xiwuIndicate unmanned plane U as 1uTo farmland A in the w time windowiComplete the drug spraying task, otherwise the unmanned aerial vehicle does notIs directed to farmland AiExecuting the task; if yijwuIndicate unmanned plane U as 1uThrough farmland Ai,AjOtherwise, the unmanned aerial vehicle does not pass through the farmland Ai,Aj
And the target function formula (1) is used for maximizing the total income of the farmland after the pesticide spraying task is completed. Constraint (2) is to ensure that all routes, i.e. all drones, start at a0End point is
Figure BDA0001270604300000094
In this disclosure, the coordinates of the starting point and the ending point are the same, and the number of routes is the sum of the routes of the unmanned aerial vehicle in two time windows. Constraint (3) allows each edge to have connectivity. The constraint equation (4) is to explain the service time required for accessing the farmland. Constraint (5) ensures that each field is visited at most once. Constraint (6) entails that the time to perform the task must be within two time windows and that the drone may not be within a time window from the starting point to the field and from the field back to the starting point. Constraint (7) translates the limit on fuel required by the drone into a limit on the flight time of each drone. Constraint (8) is the definition of variables such as target, path, etc.
It is understood that, after the MUAV-MTW-DTOP model is obtained, the method provided by the embodiment of the present invention can solve the optimal solution of the MUAV-MTW-DTOP model according to a preset genetic algorithm. The predetermined genetic algorithm for finding the optimal solution can be implemented by various methods, and one of the alternative ways is described in detail below.
The preset genetic algorithm is used for solving the problem of joint optimization of the unmanned aerial vehicle aviation operation task allocation and flight path planning, and the feasible solution is a chromosome. The population is composed of a plurality of chromosomes, the scale of the population is defined according to problems, the population is updated through the crossing and variation of the chromosomes to form a new population, wherein the crossing means that two parent chromosomes form two new child chromosomes according to the crossing probability, the variation means that one chromosome forms one new chromosome according to the variation probability, and the loop iteration is continuously carried out, so that the optimal child, namely the optimal solution, is selected. The specific flow is shown in fig. 4. The coding, interleaving, mutation, fitness function rules involved in fig. 4 are as follows:
(1) encoding
The MDD coding method comprises coding of time information of a target area farmland, an unmanned aerial vehicle, a target area entry point and an unmanned aerial vehicle accessing a first farmland in each route. And the first line of the coded chromosome is a target area, the second line is an unmanned aerial vehicle for executing tasks, the third line is an entry point when the unmanned aerial vehicle executes the tasks on the target area, and the fourth line is the time for the unmanned aerial vehicle to visit the first farmland in each route. Wherein the target region belongs to the set {1, 2.. NAN, the drone belongs to the set {0,1UAnd f, enabling the entry point of the farmland in the target area to belong to a set {1, 2.. 8}, and enabling the time for accessing the first farmland in each route to belong to any time window.
For example, in FIG. 5, there are four paths for the chromosome in two time windows, namely, the UAV U1And unmanned plane U2The represented route. At [9:00:00,12:00:00]In the time window, unmanned plane U1First at time 9.926 from RD6Point-entering farmland A2From R to RD7Point-entering farmland A6Finally returning to the starting point A0(ii) a Unmanned plane U2First at time 9.071 from RD2Point-entering farmland A7From R to RD1Point-entering farmland A1Finally returning to the starting point A0(ii) a At [14:00:00,17:00:00]In the time window, unmanned plane U1From R at time 15.57D3Point-entering farmland A5And then returns to the starting point a0(ii) a Unmanned plane U2From R at time 15.53D4Point-entering farmland A4And then returns to the starting point a0(ii) a And object A3Has not been accessed.
(2) Fitness function and selection
The fitness of a chromosome is the total profit for all visited target regions, which is related to the profit variation function of the target region, i.e. to the visited region, the area of the region. Fitness can be obtained by equation (1) and Fitness ═ Profit. The selection of chromosomes is made by roulette, wherein the likelihood of each chromosome being selected is the same.
(3) Crossing
Firstly, randomly determining the crossing position, crossing the genes in the two chromosomes of the parent, and because the fourth row of the genes represents the starting time of the unmanned aerial vehicle accessing the first farmland in the same route, in the crossed offspring chromosomes, replacing the fourth rows of other genes related to the same route with new time.
For example, in fig. 6, parents a and B are two chromosomes to be crossed selected randomly, the crossing position 7 is randomly determined in parent a, and then the seventh column of parent a is replaced with the 4 th column of parent B, for parent a, the unmanned plane U is used in the second time window1The time of the fifth column and the fourth row is replaced by 16.75, so that the time of the fifth column and the fourth row is also replaced by 16.75 to keep the time consistency of the first farmland accessed in the same route in the same time, thereby obtaining the offspring A; similarly, in the parent B, 9.926 is required to replace the time in the fourth row of the third column to obtain the child B.
(4) Variation of
The mutation in the genetic algorithm is to prevent the genetic algorithm from falling into local optima, so that the genetic algorithm has the possibility of gene mutation, the mutation can be one gene or a plurality of genes, and the chromosome mutation in the embodiment of the invention mainly has the following conditions: the sequence variation of the farmland, the variation of the unmanned aerial vehicle, the entry point variation of the target area, and the time variation of the unmanned aerial vehicle accessing the first farmland in each route.
For example, as shown in fig. 7, when the chromosome a is mutated by the drone, the seventh row of the chromosome a changes from 1 to 2, and as a result, the field 6 is originally transformed by the drone U1Access change to drone U2And (6) accessing. Simultaneously in order to guarantee unmanned aerial vehicle U1The start time of the first farmland visit in the first time window is consistent, and the seventh column and the fourth row of the chromosome A are updated to 9.071 on the basis of the variation
To show the superiority of the method provided by the embodiment of the present invention, how to solve the MUAV-MTW-DTOP model by using the genetic algorithm according to the above function setting, and thus obtain the final task allocation result, will be described in detail below with reference to several specific examples.
Specifically, the genetic algorithm is used for solving the MUAV-MTW-DTOP model in the environment of MATLAB 2013, and experiments are carried out.
Suppose that two unmanned aerial vehicles spray pesticides on 20 farmland areas in two time windows [9,12], [14,17], and use the genetic algorithm to obtain a distribution scheme, wherein the crossover probability of the genetic algorithm is taken to be 0.9, the variation probability is 0.5, the population scale is 600, and the iteration number is 100. The specific parameters involved in the experimental procedure are described below:
(1) unmanned plane
In the experiment of the embodiment of the invention, the specific configuration of the unmanned aerial vehicle is shown in table 1, the speed of the unmanned aerial vehicle is 4m/s, the maximum spraying radius is 5m, and the maximum duration is 2600 s. And the course angles of the unmanned aerial vehicle at the starting point and the ending point are both 0.
Table 1 basic parameter configuration table for unmanned aerial vehicle
Unmanned aerial vehicle parameters A0\AN+1 V RD Ewu
Unmanned aerial vehicle information (0,0) 4m/s 5m 2600s
(2) Farmland area and income thereof
This field area had 20 fields to be sprayed, as shown in detail in fig. 8. Specific farmland coordinates and earnings are shown in table 2.
TABLE 2 Farmland coordinate information and yield
Figure BDA0001270604300000121
Figure BDA0001270604300000131
The optimal solution obtained by using the genetic algorithm for the above scenario has the benefit of 15.601, and has converged at the 38 th generation, and the convergence rate is higher, as shown in fig. 9. One of the optimal allocation schemes is shown in table 3, where in all fields, fields 3,9,12,20 are not being visited and other fields have been completed with spraying within the time window.
TABLE 3 optimal solution Table
Figure BDA0001270604300000141
In a second aspect, an embodiment of the present invention further provides an apparatus for joint optimization of mission allocation and flight path planning for unmanned aerial vehicle aviation operation, wherein when a plurality of fixed-wing unmanned aerial vehicles perform a mission on a plurality of rectangular fields, the apparatus comprises:
an information obtaining unit 201, configured to obtain two time windows for executing the task, information of the fixed-wing drone, and farmland information of the multiple candidate farmlands;
an initial solution obtaining unit 202, configured to encode the time window information, the information of the unmanned aerial vehicle, and the farmland information, and randomly generate a plurality of chromosomes; screening out chromosomes meeting preset constraint conditions defined by a preset MUAV-MTW-DTOP model from the plurality of chromosomes generated randomly, and constructing an initial solution set of the MUAV-MTW-DTOP model; the MUAV-MTW-DTOP model is an objective function which enables the multi-frame fixed-wing unmanned aerial vehicle flying in a preset flying mode to obtain the maximum total benefit in the task; the preset constraints comprise flight time constraints and two time window constraints of each fixed-wing unmanned aerial vehicle;
and the optimal solution calculation unit 203 is configured to perform continuous crossing and variation on the initial solution set based on a preset genetic algorithm until the crossing and variation are finished after the constraint of the iteration times is satisfied, select an optimal solution which enables the model to obtain the maximum total yield from the varied solution set, and use the optimal solution as a task allocation and track planning result of the operation.
In specific implementation, the objective function of the MUAV-MTW-DTOP model is as follows:
Figure BDA0001270604300000151
the constraint conditions of the MUAV-MTW-DTOP model are as follows:
Figure BDA0001270604300000152
Figure BDA0001270604300000153
Figure BDA0001270604300000154
Figure BDA0001270604300000155
Figure BDA0001270604300000156
Figure BDA0001270604300000157
Figure BDA0001270604300000158
wherein N isUFor unmanned plane UuThe number of the cells; n is a radical ofAFor farmland AiThe number of the cells; w is the w-th time window, and is 1 or 2; n th0,NA+1The block farmland represents the starting point and the end point of the unmanned aerial vehicle; SQiIndicating farmland AiThe area of (d); piwIndicating that field A is completed in the w-th time windowiThe benefit of the task; siwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time to start spraying the pesticide; t is tijwuRefers to unmanned plane UuFor farmland A in the w time windowi,AjThe time of flying according to a preset flying mode; sjwuRefers to unmanned plane UuFor farmland A in the w time windowjThe time to start spraying the pesticide; m is a preset value; o isiwu,CiwuAre respectively unmanned aerial vehicle UuFor farmland A in the w time windowiA start time and an end time at which a task can be executed; t is tiwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time for executing the task according to a preset flight mode; ewuRepresenting the maximum flight time of the unmanned aerial vehicle; if xiwu1, then indicate unmanned plane UuTo farmland A in the w time windowiThe medicine spraying task is completed, otherwise the unmanned aerial vehicle does not spray the farmland AiExecuting the task; if yijwu1, then indicate unmanned plane UuThrough farmland Ai,AjOtherwise, the unmanned aerial vehicle does not pass through the farmland Ai,Aj
In specific implementation, the preset flight mode is as follows:
the method comprises the following steps of flying in a farmland according to a parallel scanning flying mode: entering the farmland from a first entry point on a first side in a direction perpendicular to the first side of the farmland, wherein the distance between the first entry point and the nearest vertex of the farmland is the scanning radius of the unmanned aerial vehicle, and the first side is any side of the region to be detected; when the unmanned aerial vehicle needs to turn, the unmanned aerial vehicle turns with the minimum turning radius larger than or equal to that of the unmanned aerial vehicle;
flying between farmlands according to a Dubins path flying mode, wherein the Dubins path flying mode is that the Dubins path flying mode flies according to a path formed by combining an arc path and a straight path, and the path formed by combining the arc path and the straight path comprises { LSL, RSR, RSL, LSR, RLR and LRL }; wherein, L represents the arc section route that unmanned aerial vehicle turned to the left, and R represents the arc section route that unmanned aerial vehicle turned to the right, and S represents that unmanned aerial vehicle flies with the straight line mode.
In a specific implementation, the initial solution obtaining unit is further configured to:
encoding the solution of the MUAV-MTW-DTOP model into a chromosome composed according to a preset structure by adopting an MMD encoding method to form an initial solution; the chromosome comprises information of a target area farmland, an unmanned aerial vehicle, a target area entry point and time for the unmanned aerial vehicle to access the first farmland in each route; wherein the target area farmland belongs to a set {1,2AN, the drone belonging to the set {0,1UAn entry point of the target region farmland belongs to a set {1, 2.. 8}, and the time for accessing the first farmland in each route belongs to a time window;
the first action of the chromosome is the random full arrangement of the target area farmlands, the second action is the random combination of unmanned aerial vehicles for executing tasks, the third action is the random combination of entry points when the unmanned aerial vehicles execute the tasks on the target area, the fourth action is the random combination of the time for the unmanned aerial vehicles to visit the first farmlands in each route, and the starting time for the same unmanned aerial vehicle to visit the first farmlands in each route is the same.
In specific implementation, the optimal solution calculating unit is further configured to perform the following steps:
step 1, generating an initial solution by using the MMD coding method, judging whether the initial solution meets the constraint condition of the MUAV-MTW-DTOP model, if not, continuously generating the initial solution meeting the constraint condition of the MUAV-MTW-DTOP model, generating a parent population with a preset scale and calculating the fitness of the parent population;
step 2, selecting two individuals (A, B) in a parent population to cross by using a roulette method, wherein the crossing rule is that the crossing position in the individual A is randomly selected, then the gene in the individual B, which is the same as the first line of the crossing position of the individual A, is searched, the gene at the crossing position in the chromosomes A and B is replaced, all genes, which are newly exchanged in the chromosomes A and B and are the same as the second line of the genes and the fourth line of the genes is in the same time window, are searched, and the fourth line of the searched gene is replaced by the fourth line of the newly exchanged genes in the chromosomes A and B to obtain chromosomes C and D; judging whether chromosomes C and D meet the constraint conditions of the MUAV-MTW-DTOP model, if so, replacing the chromosomes A and B in the population by using the chromosomes C and D, otherwise, ending, and continuously iterating and updating the population in the step 1 to obtain a new offspring population;
and 3, selecting one chromosome in the population in the step 2 for mutation by using a roulette method, wherein the method for mutating the chromosome is at least one of the following mutation methods, and the method comprises the following steps: carrying out random full arrangement on the first row sequence of the chromosomes; and carrying out variation on the unmanned aerial vehicle information of the second line of the chromosome, wherein the variation range is {0, 1.. NU}; and mutating an entry point of a target region farmland of a third row of chromosomes, wherein the entry point belongs to a set {1, 2.. 8 }; and a fourth temporal variation of the chromosome belonging to two time windows;
the whole chromosome variation steps include: firstly, if the first row sequence of the chromosome is mutated, a farmland sequence full array is randomly generated; selecting whether the second row is mutated or not and the position of the mutation, and if so, randomly generating a mutated unmanned aerial vehicle to replace the original position; randomly selecting whether the third row has variation and a variation position, and if the variation exists, randomly generating a variation entry point to replace the original position; randomly selecting whether the fourth line is mutated or not and a mutation position, randomly generating mutation time if the fourth line is mutated, replacing the original position, finding a gene which is the same as the second line of the mutation position and in the same time window as the fourth line, mutating the fourth line into new mutation time, and finally obtaining a new mutation chromosome;
judging whether the varied chromosomes meet the constraint conditions of the MUAV-MTW-DTOP model or not, if so, replacing the chromosomes in the population, otherwise, finishing the replacement and continuously iterating and updating the population in the step 2 to obtain a new offspring population;
step 4, calculating the population fitness of the offspring in the step 3 and selecting the optimal solution of all solutions in the iteration;
step 5, judging whether the current iteration times reach a preset value or not, if not, combining the child population and the parent population in the step 3 according to a certain proportion to form a new parent population, and returning to the step 2; if so, ending the iteration, and taking the finally obtained optimal solution as a task allocation and track planning result of the unmanned aerial vehicle.
Since the device for the joint optimization of the assignment of the unmanned aerial vehicle aerial work task and the flight path planning, which is introduced in this embodiment, is a device capable of executing the method for the joint optimization of the assignment of the unmanned aerial vehicle aerial work task and the flight path planning, which is introduced in this embodiment of the present invention, based on the method for the joint optimization of the assignment of the unmanned aerial vehicle aerial work task and the flight path planning, those skilled in the art can understand the specific implementation manner and various variations of the device for the joint optimization of the assignment of the unmanned aerial vehicle aerial work task and the flight path planning, and therefore, how to implement the method for the joint optimization of the assignment of the unmanned aerial vehicle aerial work task and the flight path planning, which is introduced in this embodiment, is not described in detail herein. As long as the person skilled in the art implements the device adopted by the method for joint optimization of unmanned aerial vehicle aerial work task allocation and flight path planning in the embodiment of the present invention, the device belongs to the protection scope of the present application.
In the description provided herein, numerous specific details are set forth. It is understood, however, that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Similarly, it should be appreciated that in the foregoing description of exemplary embodiments of the invention, various features of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. However, the disclosed method should not be interpreted as reflecting an intention that: that the invention as claimed requires more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive aspects lie in less than all features of a single foregoing disclosed embodiment. Thus, the claims following the detailed description are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate embodiment of this invention.
Those skilled in the art will appreciate that the modules in the device in an embodiment may be adaptively changed and disposed in one or more devices different from the embodiment. The modules or units or components of the embodiments may be combined into one module or unit or component, and furthermore they may be divided into a plurality of sub-modules or sub-units or sub-components. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or elements of any method or apparatus so disclosed, may be combined in any combination, except combinations where at least some of such features and/or processes or elements are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings) may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.
Furthermore, those skilled in the art will appreciate that while some embodiments herein include some features included in other embodiments, rather than other features, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the following claims, any of the claimed embodiments may be used in any combination.
Some component embodiments of the invention may be implemented in hardware, or in software modules running on one or more processors, or in a combination thereof. Those skilled in the art will appreciate that a microprocessor or Digital Signal Processor (DSP) may be used in practice to implement some or all of the functionality of some or all of the components of a gateway, proxy server, system according to embodiments of the present invention. The present invention may also be embodied as apparatus or device programs (e.g., computer programs and computer program products) for performing a portion or all of the methods described herein. Such programs implementing the present invention may be stored on computer-readable media or may be in the form of one or more signals. Such a signal may be downloaded from an internet website or provided on a carrier signal or in any other form.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The usage of the words first, second and third, etcetera do not indicate any ordering. These words may be interpreted as names.

Claims (8)

1. A joint optimization method for task allocation and flight path planning of unmanned aerial vehicle aerial work is characterized in that when a plurality of fixed-wing unmanned aerial vehicles execute a task on a plurality of rectangular farmlands, the method comprises the following steps:
acquiring two time windows for executing the task, information of the fixed-wing unmanned aerial vehicle and farmland information of a plurality of candidate farmlands;
coding the time window information, the unmanned aerial vehicle information and the farmland information, and randomly generating a plurality of chromosomes; screening out chromosomes meeting preset constraint conditions defined by a preset MUAV-MTW-DTOP model from the plurality of chromosomes generated randomly, and constructing an initial solution set of the MUAV-MTW-DTOP model; the MUAV-MTW-DTOP model is an objective function which enables the multi-frame fixed-wing unmanned aerial vehicle flying in a preset flying mode to obtain the maximum total benefit in the task; the preset constraints comprise flight time constraints and two time window constraints of each fixed-wing unmanned aerial vehicle;
continuously crossing and mutating the initial solution set based on a preset genetic algorithm until the initial solution set meets the constraint of iteration times, ending the crossing and mutation, selecting an optimal solution which enables the model to obtain the maximum total income from the mutated solution set, and taking the optimal solution as a task allocation and track planning result of the operation;
wherein the objective function of the MUAV-MTW-DTOP model is as follows:
the constraint conditions of the MUAV-MTW-DTOP model are as follows:
Figure FDA0002302487290000012
Figure FDA0002302487290000013
Figure FDA0002302487290000014
Figure FDA0002302487290000021
Figure FDA0002302487290000022
Figure FDA0002302487290000023
Figure FDA0002302487290000024
wherein N isUFor unmanned plane UuThe number of the cells; n is a radical ofAFor farmland AiThe number of the cells; w is the w-th time window, and is 1 or 2; n th0,NA+1The block farmland represents the starting point and the end point of the unmanned aerial vehicle; SQiIndicating farmland AiThe area of (d); piwIndicating that field A is completed in the w-th time windowiThe benefit of the task; siwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time to start spraying the pesticide; t is tijwuRefers to unmanned plane UuFor farmland A in the w time windowi,AjThe time of flying according to a preset flying mode; sjwuRefers to unmanned plane UuFor farmland A in the w time windowjThe time to start spraying the pesticide; m is a preset value; o isiwu,CiwuAre respectively unmanned aerial vehicle UuFor farmland A in the w time windowiA start time and an end time for executing the task; t is tiwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time for executing the task according to a preset flight mode; ewuRepresenting the maximum flight time of the unmanned aerial vehicle; if xiwu1, then indicate unmanned plane UuTo farmland A in the w time windowiThe medicine spraying task is completed, otherwise the unmanned aerial vehicle does not spray the farmland AiExecuting the task; if yijwu1, then indicate unmanned plane UuThrough farmland Ai,AjOtherwise the nobodyThe machine has not passed through the farmland Ai,Aj
2. The method according to claim 1, characterized in that the predetermined flight pattern is:
the method comprises the following steps of flying in a farmland according to a parallel scanning flying mode: entering the farmland from a first entry point on a first side in a direction perpendicular to the first side of the farmland, wherein the distance between the first entry point and the nearest vertex of the farmland is the scanning radius of the unmanned aerial vehicle, and the first side is any side of the region to be detected; when the unmanned aerial vehicle needs to turn, the unmanned aerial vehicle turns with the minimum turning radius larger than or equal to that of the unmanned aerial vehicle;
flying between farmlands according to a Dubins path flying mode, wherein the Dubins path flying mode is that the Dubins path flying mode flies according to a path formed by combining an arc path and a straight path, and the path formed by combining the arc path and the straight path comprises { LSL, RSR, RSL, LSR, RLR and LRL }; wherein, L represents the arc section route that unmanned aerial vehicle turned to the left, and R represents the arc section route that unmanned aerial vehicle turned to the right, and S represents that unmanned aerial vehicle flies with the straight line mode.
3. The method of claim 1, wherein the step of selecting chromosomes satisfying the predetermined constraints defined by the predetermined MUAV-MTW-DTOP model from among the randomly generated chromosomes, and the step of constructing the initial solution set of the MUAV-MTW-DTOP model comprises:
encoding the solution of the MUAV-MTW-DTOP model into a chromosome composed according to a preset structure by adopting an MMD encoding method to form an initial solution; the chromosome comprises information of a target area farmland, an unmanned aerial vehicle, a target area entry point and time for the unmanned aerial vehicle to access the first farmland in each route; wherein the target area farmland belongs to a set {1,2AN, the drone belonging to the set {0,1UAn entry point of the target region farmland belongs to a set {1, 2.. 8}, and the time for accessing the first farmland in each route belongs to a time window;
the first action of the chromosome is the random full arrangement of the target area farmlands, the second action is the random combination of unmanned aerial vehicles for executing tasks, the third action is the random combination of entry points when the unmanned aerial vehicles execute the tasks on the target area, the fourth action is the random combination of the time for the unmanned aerial vehicles to visit the first farmlands in each route, and the starting time for the same unmanned aerial vehicle to visit the first farmlands in each route is the same.
4. The method according to claim 3, wherein the initial solution set is continuously crossed and varied based on a preset genetic algorithm until the crossing and variation are finished after the constraint of the iteration times is satisfied, an optimal solution which enables the model to obtain the maximum total income is selected from the varied solution set, and the optimal solution is used as a task allocation and flight path planning result of the operation, and the method comprises the following steps:
step 1, generating an initial solution by using the MMD coding method, judging whether the initial solution meets the constraint condition of the MUAV-MTW-DTOP model, if not, continuously generating the initial solution meeting the constraint condition of the MUAV-MTW-DTOP model, generating a parent population with a preset scale and calculating the fitness of the parent population;
step 2, selecting an individual A and an individual B in a parent population to be crossed by using a roulette method, wherein the crossing rule is that the crossing position in the individual A is randomly selected, then the gene in the individual B, which is the same as the first line of the crossing position of the individual A, is searched, the gene at the crossing position in the chromosomes A and B is replaced, all genes, which are newly exchanged in the chromosomes A and B and are the same as the second line of the genes and the fourth line of the genes is in the same time window, are searched, and the fourth line of the searched genes is replaced by the fourth line of the newly exchanged genes in the chromosomes A and B to obtain chromosomes C and D; judging whether chromosomes C and D meet the constraint conditions of the MUAV-MTW-DTOP model, if so, replacing the chromosomes A and B in the population by using the chromosomes C and D, otherwise, ending, and continuously iterating and updating the population in the step 1 to obtain a new offspring population;
and 3, selecting one chromosome in the population in the step 2 for mutation by using a roulette method, wherein the mutation mode of the chromosome is at least one of the following mutation modesThe method comprises the following steps: carrying out random full arrangement on the first row sequence of the chromosomes; and carrying out variation on the unmanned aerial vehicle information of the second line of the chromosome, wherein the variation range is {0, 1.. NU}; and mutating an entry point of a target region farmland of a third row of chromosomes, wherein the entry point belongs to a set {1, 2.. 8 }; and a fourth temporal variation of the chromosome belonging to two time windows;
the whole chromosome variation steps include: firstly, if the first row sequence of the chromosome is mutated, a farmland sequence full array is randomly generated; selecting whether the second row is mutated or not and the position of the mutation, and if so, randomly generating a mutated unmanned aerial vehicle to replace the original position; randomly selecting whether the third row has variation and a variation position, and if the variation exists, randomly generating a variation entry point to replace the original position; randomly selecting whether the fourth line is mutated or not and a mutation position, randomly generating mutation time if the fourth line is mutated, replacing the original position, finding a gene which is the same as the second line of the mutation position and in the same time window as the fourth line, mutating the fourth line into new mutation time, and finally obtaining a new mutation chromosome;
judging whether the varied chromosomes meet the constraint conditions of the MUAV-MTW-DTOP model or not, if so, replacing the chromosomes in the population, otherwise, finishing the replacement and continuously iterating and updating the population in the step 2 to obtain a new offspring population;
step 4, calculating the population fitness of the offspring in the step 3 and selecting the optimal solution of all solutions in the iteration;
step 5, judging whether the current iteration times reach a preset value or not, if not, combining the child population and the parent population in the step 3 according to a certain proportion to form a new parent population, and returning to the step 2; if so, ending the iteration, and taking the finally obtained optimal solution as a task allocation and track planning result of the unmanned aerial vehicle.
5. The utility model provides an unmanned aerial vehicle air work task allocation and flight path planning unite optimization device which characterized in that, when many fixed wing unmanned aerial vehicles carry out a task to polylith rectangle farmland, the device includes:
the information acquisition unit is used for acquiring two time windows for executing the task, information of the fixed-wing unmanned aerial vehicle and farmland information of a plurality of candidate farmlands;
the initial solution acquisition unit is used for encoding the time window information, the information of the unmanned aerial vehicle and the farmland information and randomly generating a plurality of chromosomes; screening out chromosomes meeting preset constraint conditions defined by a preset MUAV-MTW-DTOP model from the plurality of chromosomes generated randomly, and constructing an initial solution set of the MUAV-MTW-DTOP model; the MUAV-MTW-DTOP model is an objective function which enables the multi-frame fixed-wing unmanned aerial vehicle flying in a preset flying mode to obtain the maximum total benefit in the task; the preset constraints comprise flight time constraints and two time window constraints of each fixed-wing unmanned aerial vehicle;
the optimal solution calculation unit is used for continuously crossing and mutating the initial solution set based on a preset genetic algorithm until the initial solution set meets the constraint of iteration times, ending the crossing and the mutation, selecting an optimal solution which enables the model to obtain the maximum total income from the mutated solution set, and taking the optimal solution as a task allocation and flight path planning result of the operation;
wherein the objective function of the MUAV-MTW-DTOP model is as follows:
Figure FDA0002302487290000051
the constraint conditions of the MUAV-MTW-DTOP model are as follows:
Figure FDA0002302487290000061
Figure FDA0002302487290000063
Figure FDA0002302487290000064
Figure FDA0002302487290000065
Figure FDA0002302487290000067
wherein N isUFor unmanned plane UuThe number of the cells; n is a radical ofAFor farmland AiThe number of the cells; w is the w-th time window, and is 1 or 2; n th0,NA+1The block farmland represents the starting point and the end point of the unmanned aerial vehicle; SQiIndicating farmland AiThe area of (d); piwIndicating that field A is completed in the w-th time windowiThe benefit of the task; siwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time to start spraying the pesticide; t is tijwuRefers to unmanned plane UuFor farmland A in the w time windowi,AjThe time of flying according to a preset flying mode; sjwuRefers to unmanned plane UuFor farmland A in the w time windowjThe time to start spraying the pesticide; m is a preset value; o isiwu,CiwuAre respectively unmanned aerial vehicle UuFor farmland A in the w time windowiA start time and an end time for executing the task; t is tiwuRefers to unmanned plane UuFor farmland A in the w time windowiThe time for executing the task according to a preset flight mode; ewuRepresenting the maximum flight time of the unmanned aerial vehicle; if xiwu1, then indicate unmanned plane UuTo farmland A in the w time windowiThe medicine spraying task is completed, otherwise the unmanned aerial vehicle does not spray the farmland AiExecuting the task; if yijwu1, then indicate unmanned plane UuThrough farmland Ai,AjOtherwise, theUnmanned aerial vehicle does not pass through farmland Ai,Aj
6. The device according to claim 5, characterized in that the predetermined flight pattern is:
the method comprises the following steps of flying in a farmland according to a parallel scanning flying mode: entering the farmland from a first entry point on a first side in a direction perpendicular to the first side of the farmland, wherein the distance between the first entry point and the nearest vertex of the farmland is the scanning radius of the unmanned aerial vehicle, and the first side is any side of the region to be detected; when the unmanned aerial vehicle needs to turn, the unmanned aerial vehicle turns with the minimum turning radius larger than or equal to that of the unmanned aerial vehicle;
flying between farmlands according to a Dubins path flying mode, wherein the Dubins path flying mode is that the Dubins path flying mode flies according to a path formed by combining an arc path and a straight path, and the path formed by combining the arc path and the straight path comprises { LSL, RSR, RSL, LSR, RLR and LRL }; wherein, L represents the arc section route that unmanned aerial vehicle turned to the left, and R represents the arc section route that unmanned aerial vehicle turned to the right, and S represents that unmanned aerial vehicle flies with the straight line mode.
7. The apparatus of claim 5, wherein the initial solution obtaining unit is further configured to:
encoding the solution of the MUAV-MTW-DTOP model into a chromosome composed according to a preset structure by adopting an MMD encoding method to form an initial solution; the chromosome comprises information of a target area farmland, an unmanned aerial vehicle, a target area entry point and time for the unmanned aerial vehicle to access the first farmland in each route; wherein the target area farmland belongs to a set {1,2AN, the drone belonging to the set {0,1UAn entry point of the target region farmland belongs to a set {1, 2.. 8}, and the time for accessing the first farmland in each route belongs to a time window;
the first action of the chromosome is the random full arrangement of the target area farmlands, the second action is the random combination of unmanned aerial vehicles for executing tasks, the third action is the random combination of entry points when the unmanned aerial vehicles execute the tasks on the target area, the fourth action is the random combination of the time for the unmanned aerial vehicles to visit the first farmlands in each route, and the starting time for the same unmanned aerial vehicle to visit the first farmlands in each route is the same.
8. The apparatus of claim 7, wherein the optimal solution computing unit is further configured to perform the following steps:
step 1, generating an initial solution by using the MMD coding method, judging whether the initial solution meets the constraint condition of the MUAV-MTW-DTOP model, if not, continuously generating the initial solution meeting the constraint condition of the MUAV-MTW-DTOP model, generating a parent population with a preset scale and calculating the fitness of the parent population;
step 2, selecting an individual A and an individual B in a parent population to be crossed by using a roulette method, wherein the crossing rule is that the crossing position in the individual A is randomly selected, then the gene in the individual B, which is the same as the first line of the crossing position of the individual A, is searched, the gene at the crossing position in the chromosomes A and B is replaced, all genes, which are newly exchanged in the chromosomes A and B and are the same as the second line of the genes and the fourth line of the genes is in the same time window, are searched, and the fourth line of the searched genes is replaced by the fourth line of the newly exchanged genes in the chromosomes A and B to obtain chromosomes C and D; judging whether chromosomes C and D meet the constraint conditions of the MUAV-MTW-DTOP model, if so, replacing the chromosomes A and B in the population by using the chromosomes C and D, otherwise, ending, and continuously iterating and updating the population in the step 1 to obtain a new offspring population;
and 3, selecting one chromosome in the population in the step 2 for mutation by using a roulette method, wherein the method for mutating the chromosome is at least one of the following mutation methods, and the method comprises the following steps: carrying out random full arrangement on the first row sequence of the chromosomes; and carrying out variation on the unmanned aerial vehicle information of the second line of the chromosome, wherein the variation range is {0, 1.. NU}; and mutating an entry point of a target region farmland of a third row of chromosomes, wherein the entry point belongs to a set {1, 2.. 8 }; and a fourth temporal variation of the chromosome belonging to two time windows;
the whole chromosome variation steps include: firstly, if the first row sequence of the chromosome is mutated, a farmland sequence full array is randomly generated; selecting whether the second row is mutated or not and the position of the mutation, and if so, randomly generating a mutated unmanned aerial vehicle to replace the original position; randomly selecting whether the third row has variation and a variation position, and if the variation exists, randomly generating a variation entry point to replace the original position; randomly selecting whether the fourth line is mutated or not and a mutation position, randomly generating mutation time if the fourth line is mutated, replacing the original position, finding a gene which is the same as the second line of the mutation position and in the same time window as the fourth line, mutating the fourth line into new mutation time, and finally obtaining a new mutation chromosome;
judging whether the varied chromosomes meet the constraint conditions of the MUAV-MTW-DTOP model or not, if so, replacing the chromosomes in the population, otherwise, finishing the replacement and continuously iterating and updating the population in the step 2 to obtain a new offspring population;
step 4, calculating the population fitness of the offspring in the step 3 and selecting the optimal solution of all solutions in the iteration;
step 5, judging whether the current iteration times reach a preset value or not, if not, combining the child population and the parent population in the step 3 according to a certain proportion to form a new parent population, and returning to the step 2; if so, ending the iteration, and taking the finally obtained optimal solution as a task allocation and track planning result of the unmanned aerial vehicle.
CN201710245797.5A 2017-04-14 2017-04-14 Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device Active CN107037827B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710245797.5A CN107037827B (en) 2017-04-14 2017-04-14 Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710245797.5A CN107037827B (en) 2017-04-14 2017-04-14 Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device

Publications (2)

Publication Number Publication Date
CN107037827A CN107037827A (en) 2017-08-11
CN107037827B true CN107037827B (en) 2020-02-21

Family

ID=59536349

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710245797.5A Active CN107037827B (en) 2017-04-14 2017-04-14 Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device

Country Status (1)

Country Link
CN (1) CN107037827B (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107320222B (en) * 2017-08-14 2019-01-15 侯策 A kind of intelligent artifucial limb replacement method and its system based on unmanned plane
CN109460053A (en) * 2017-09-06 2019-03-12 广东泰高新技术发展有限公司 Aerial survey of unmanned aerial vehicle turning flight course planning method
CN107886201B (en) * 2017-11-29 2021-01-12 合肥有空科技信息技术有限公司 Multi-objective optimization method and device for multi-unmanned aerial vehicle task allocation
CN108594645B (en) * 2018-03-08 2021-02-19 中国人民解放军国防科技大学 Planning method and system for single-station multi-unmanned aerial vehicle distribution and flight route
WO2019183856A1 (en) * 2018-03-28 2019-10-03 深圳市大疆软件科技有限公司 Control method and apparatus for operation of plant protection unmanned aerial vehicle
CN108445911A (en) * 2018-05-30 2018-08-24 河北清华发展研究院 A kind of sliceable unmanned aerial vehicle group control method
CN108919832A (en) 2018-07-23 2018-11-30 京东方科技集团股份有限公司 Unmanned machine operation flight course planning method, unmanned plane application method and device
CN109197278B (en) * 2018-10-18 2021-01-19 广州极飞科技有限公司 Method and device for determining operation strategy and method for determining drug spraying strategy
CN111310992B (en) * 2020-02-10 2022-03-22 合肥工业大学 Multi-unmanned aerial vehicle path optimization method for rapid evaluation after earthquake disaster
CN111461464B (en) * 2020-05-06 2021-06-08 农业农村部南京农业机械化研究所 Plant protection unmanned aerial vehicle cluster operation task allocation method and device
CN112819242B (en) * 2021-02-22 2023-06-09 西北工业大学 Civil transportation airplane flight test task allocation optimization method
CN113126650A (en) * 2021-03-03 2021-07-16 华南农业大学 Automatic weeding operation method for unmanned aerial vehicle
CN114021914B (en) * 2021-10-22 2024-07-09 北京市农林科学院信息技术研究中心 Unmanned aerial vehicle cluster flight protection scheduling method and device
CN114355985B (en) * 2022-03-18 2022-07-08 北京卓翼智能科技有限公司 Path planning method and device for unmanned aerial vehicle cluster, controller and storage medium

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9651950B2 (en) * 2012-07-18 2017-05-16 The Boeing Company Mission re-planning for coordinated multivehicle task allocation
CN103413186B (en) * 2013-08-19 2016-01-20 中国电子科技集团公司第二十八研究所 A kind of multi-aircraft based on hybrid optimization algorithm works in coordination with target assignment method
CN103699135B (en) * 2014-01-03 2016-04-06 中南大学 The flight path automatic planning in depopulated helicopter pesticide spraying farmland operation region
CN104406593A (en) * 2014-12-03 2015-03-11 广西民族大学 Method for determining optimal route of airway of unmanned aerial vehicle
KR101666132B1 (en) * 2014-12-26 2016-10-18 서울대학교산학협력단 Market-based decentralized task assignment method for UAV and recording medium recorded program to perform the method
CN104808660B (en) * 2015-03-04 2017-10-03 中南大学 Bumps mixing complex polygon farmland unmanned plane spraying operation path planning method
CN105512769A (en) * 2015-12-16 2016-04-20 上海交通大学 Unmanned aerial vehicle route planning system and unmanned aerial vehicle route planning method based on genetic programming
CN106406346B (en) * 2016-11-01 2019-04-16 北京理工大学 A kind of multiple no-manned plane collaboration rapid Cover search path planning method

Also Published As

Publication number Publication date
CN107037827A (en) 2017-08-11

Similar Documents

Publication Publication Date Title
CN107037827B (en) Unmanned aerial vehicle aerial work task allocation and flight path planning combined optimization method and device
CN108362285B (en) Plant protection unmanned aerial vehicle accurate coverage track planning method and system
CN104503464B (en) Computer-based convex polygon field unmanned aerial vehicle spraying operation route planning method
CN106969778B (en) Path planning method for cooperative pesticide application of multiple unmanned aerial vehicles
CN103699135B (en) The flight path automatic planning in depopulated helicopter pesticide spraying farmland operation region
CN107390709B (en) A kind of plant protection drone multi-machine collaborative operational method and system
CN107289950B (en) Plant protection drone operation flight course planning method and plant protection drone
Thayer et al. Multirobot routing algorithms for robots operating in vineyards
CN106598070B (en) More obstacles, barrier-avoiding method and unmanned plane under small obstacle during a kind of agricultural plant protection unmanned plane sprays
CN106842963B (en) Multiple no-manned plane detection mission is distributed and trajectory planning combined optimization method and device
CN107807665B (en) Unmanned aerial vehicle formation detection task cooperative allocation method and device
CN109459026B (en) Multi-moving-body collaborative full-coverage path planning method
CN108229719A (en) Unmanned plane formation task distributes the Multipurpose Optimal Method and device with trajectory planning
Balampanis et al. Area partition for coastal regions with multiple UAS
CN105841702A (en) Method for planning routes of multi-unmanned aerial vehicles based on particle swarm optimization algorithm
Luo et al. Optimization of Pesticide Spraying Tasks via Multi‐UAVs Using Genetic Algorithm
CN107748499B (en) Optimization method and device for multi-zone detection task planning of fixed-wing unmanned aerial vehicle
Sujit et al. Self assessment-based decision making for multiagent cooperative search
Yu et al. Balanced multi-region coverage path planning for unmanned aerial vehicles
CN116880561A (en) Optimization method and system based on improved particle swarm unmanned aerial vehicle path planning safety enhancement
Wang et al. Collaborative path planning and task allocation for multiple agricultural machines
CN107037826B (en) Unmanned plane detection mission distribution method and device
Stachniss et al. Coordinated multi-robot exploration
CN113703488A (en) Multi-operation plant protection unmanned aerial vehicle path planning method based on improved ant colony algorithm
Chen et al. Cooperative task assignment for multi-uav attack mobile targets

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