CN111461464A - Plant protection unmanned aerial vehicle cluster operation task allocation method and device - Google Patents
Plant protection unmanned aerial vehicle cluster operation task allocation method and device Download PDFInfo
- Publication number
- CN111461464A CN111461464A CN202010370842.1A CN202010370842A CN111461464A CN 111461464 A CN111461464 A CN 111461464A CN 202010370842 A CN202010370842 A CN 202010370842A CN 111461464 A CN111461464 A CN 111461464A
- Authority
- CN
- China
- Prior art keywords
- matrix
- plant protection
- sub
- unmanned aerial
- field
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Development Economics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Genetics & Genomics (AREA)
- Educational Administration (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Physiology (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Feedback Control In General (AREA)
- Traffic Control Systems (AREA)
Abstract
The invention discloses a plant protection unmanned aerial vehicle cluster operation task allocation method which comprises the steps of dividing an operation field into N sub-fields, allocating all the sub-fields to all plant protection unmanned aerial vehicles participating in operation according to a preset allocation principle, searching an optimal operation allocation matrix by adopting an MOSE L A algorithm, determining an operation time matrix of the sub-fields according to the optimal operation allocation matrix and an aircraft operation time matrix, defining an operation sequence matrix of the sub-fields by combining a preset assembly rule, and optimizing the operation sequence matrix of the sub-fields by adopting a genetic algorithm.
Description
Technical Field
The invention relates to the technical field of cluster operation of plant protection unmanned aerial vehicles, in particular to a method and a device for distributing cluster operation tasks of plant protection unmanned aerial vehicles.
Background
The plant protection unmanned aerial vehicle has the remarkable characteristics of high speed, high efficiency, wide adaptability and the like, is flexible in operation, does not need a take-off and landing runway, can adapt to complex terrains such as hills, mountainous areas, sloping fields, paddy fields and the like, is rapidly developed in China in recent years, and the annual accumulated operation area exceeds 5 hundred million mu times.
However, the plant protection unmanned aerial vehicle is limited by the carrying capacity of the machine (less than 20kg), the single-frame operation area generally does not exceed 30 mu, and in the actual operation process, the machine needs to waste 40% of the operation time for frequently taking off and landing, replacing batteries and filling pesticides. In order to improve the operation efficiency, part of enterprises and scientific research teams develop the operation technical research of one control of multiple machines, the current research work mainly focuses on networking communication and system hardware development, the most representative of the research work is that a P-30XP type plant protection unmanned aerial vehicle is released by Guangzhou polar flight science and technology company in 2019, and a T-20 type plant protection unmanned aerial vehicle released by Dajiang innovation company in 2019 supports that one person controls more than two machines to operate simultaneously. But the operation path planning mode is single, the field blocks are simply divided according to the operation capacity of the machines, and each machine adopts an independent air route to perform sequencing operation. Only the full-coverage path planning is considered, the cooperation and the coordination among machines are rarely considered, and the situations of lifting conflict and unreasonable task allocation are easily generated in the actual operation.
Disclosure of Invention
The invention aims to provide a method and a device for distributing operation tasks of a plant protection unmanned aerial vehicle cluster, which are used for integrating two constraint conditions of time and energy consumption on the basis of ensuring the full coverage of an airway, establishing a double-layer decision method for distributing the shortest flight distance and the optimal time, reducing time conflict and improving operation efficiency.
In order to achieve the above purpose, with reference to fig. 1, the present invention provides a plant protection unmanned aerial vehicle cluster operation task allocation method, where the allocation method includes the following steps:
s1, dividing the operation field into N sub-fields according to the input field parameters and the plant protection unmanned aerial vehicle parameters, wherein the operation requirements of each sub-field are the same, and only the sub-field is required to be operatedCorresponding to one unmanned aircraft, defining the set of N sub-fields as Hj,j=1,2,3,…,N;
S2, distributing all the sub-fields to all the plant protection unmanned aerial vehicles participating in the operation according to a preset distribution principle, wherein the preset distribution principle is as follows: a single plant protection unmanned aerial vehicle is allocated with at least one sub-field block, each sub-field block can be allocated only once, and all the sub-field blocks are allocated;
s3, setting the number of unmanned planes, the number of field blocks K, the number of sub-field blocks N and the operation time matrix (T)K×NAnd the distance matrix of the in/out field is (L)1)K×N、(L2)K×NBased on the operation, an MOSE L A algorithm is adopted to find an optimal operation distribution matrix (X)K×N;
S4, distributing matrix (X) according to the optimal operationK×NAnd the aircraft operation time matrix (T)K×NDetermining an operation time matrix (T') of the sub-field blockN×1(ii) a Defining a sequence of operations matrix (SE) for the sub-fields in combination with predetermined assembly rulesN×1The preset assembly rule is that after each unmanned aircraft lands, assembly tasks including battery replacement and medicine chest replacement are carried out, but only one unmanned aircraft executes the assembly tasks at the same time;
s5, adopting genetic algorithm to work sequence matrix (SE) of the sub-farmland blocksN×1And (6) optimizing.
As a preferred example, in step S1, the process of dividing the working field into N sub-fields according to the input field parameters and plant protection unmanned aerial vehicle parameters includes the following steps:
s11, putting the fields with the same operation requirements into the same set, wherein the number of the fields is M;
s12, further field processing is carried out on all the fields in the M sets according to the one-time lifting/landing operation area of the plant protection unmanned aerial vehicle, and the field processing comprises the operations of adjacent small field merging and large field dividing;
s13, extracting all the sub-field blocks in the M sets, setting the number of the sub-field blocks to be N, and defining the N sub-field block sets to be Hj,j=1,2,3,…,N。
As a preferred example, in step S1 and step S2, the allocating all the sub-fields to all the plant protection unmanned aerial vehicles participating in the operation according to the preset allocation principle includes:
defining the operation distribution matrix of plant unmanned plane as (X)K×NWherein a parameter x is assignedijThe assignment is defined as:
where j is 1,2,3, …, N.
As a preferred example, in step S1 and step S3, the MOSE L A algorithm is used to find the optimal job distribution matrix (X)K×NComprises the following steps:
s31, inputting population scale S, sub-population number S, individual frog num of the population, total population iteration number D and population iteration number I;
s32, initializing the frog group P ═ X1,X2,…,XS};
S33, calculating adaptive values, arranging the adaptive values in a descending order, and dividing the group;
s34, intra-population evolution: each sub-population makes a first jump according to the following formula:
wherein: xbAnd XwRespectively, the optimal individual and the worst individual of the sub-group;
s35, merging of first population: and (4) performing descending arrangement according to the adaptive value, and combining the following formula, and performing second jump on the whole group:
wherein: xgAnd XwRespectively, the optimal individual and the worst individual of the whole population;
s36, second population merging: and (4) performing descending arrangement according to the adaptive value, and performing third jump on the weakest individuals of the population by combining the following formula:
(Xw)′=rand(X),X∈QM×N
wherein: rand (X) represents randomly generating a decision matrix in the feasible solution space;
s37, judging whether the current cycle times are less than the total group iteration times D, if so, turning to the step S33, otherwise, outputting the optimal individual X1;
Wherein, the related parameters of population update are defined as follows:
let matrix X1∈QM×N、X2∈QM×NWherein Q isM×NFor a feasible solution space, a matrix is definedRepresents X2And X1The matrix of the differences between the two,indicating an ectopic or operational point between matrices;
the meaning of z ═ cut (z, k) is: intercepting the front k columns of the one-dimensional array z and storing the front k columns into a new one-dimensional array z';
(X1) ' -F (B, Φ) means: updating an approximation matrix B for the matrix X according to the difference matrix B, counting element sums of all columns of B, if the element sums are larger than 0, storing the column number into a one-dimensional array z, and defining an operator phi as follows:
φ=~(B,X1,cut(z,ri(n)))
wherein: n is the length of the array z, ri (n) denotes the interval [1, n ]]Any integer within; "-" indicates that X is assigned to the column number ri (n) designated by z1Element correction in the column is consistent with the arrangement of the B element;
the meaning of (A) is: for matrix X1According to the operatorPerforming mutation operation, operatorIs defined as:
in the formula: switch (r)1,r2) Representing r in a switching matrix1Row and r2Line, swicth (c)1,c2) Representing c in the switching matrix1Column sum c2And (4) columns.
As a preferred example, in step S1, in step S33, the process of dividing the population includes the following steps:
s331, defining an adaptive value function FV of an individual t1As follows below, the following description will be given,
s332, taking the adaptive value of the individual as a reference for dividing the group according to FV1 tIn descending order, the individuals in the ith subgroup are:
{FV1 k|k∈[(i-1)×num+1,i×num]}
wherein i is 1,2,3, …, s.
As a preferred example, in step S1 and step S34, the process of intra-population evolution includes the following steps:
s341, jumping to the optimal individual of the sub-population according to the following formula:
s342, introducing a mutation operator for detecting whether the plant protection unmanned aerial vehicle which is not distributed to the field exists, and performing mutation detection and processing by adopting the mutation operator:
(Xw)″=F((Xw)′,switch([r0,rm],c))
in the formula: r is0Representation matrix (X)w) ' line element and line number 0; r ismRepresentation matrix (X)w) ' line element sum is the largest line number; c represents rmAny column number with a median of 1; switch ([ r ]0,rm]And c) is represented in matrix (X)w) In column c of `, r is exchanged0And rmWhen r is an element of0When the value of (A) is not unique, a jump formula corresponding to the first jump is adopted for carrying out multiple variations;
s343, updating the optimal and worst individuals;
and S344, repeating the steps S341 to S343 until the cycle number reaches I times.
As a preferred example, in step S1, in step S35, the process of merging the first population includes the following steps:
s351, jumping to the population-optimized individual using the following formula:
s352, introducing a mutation operator for detecting whether the plant protection unmanned aerial vehicle which is not distributed to the field exists, and performing mutation detection and processing by adopting the mutation operator:
(Xw)″=F((Xw)′,switch([r0,rm],c))
in the formula: r is0Representation matrix (X)w) ' line element and line number 0; r ismRepresentation matrix (X)w) ' line element sum is the largest line number; c represents rmAny column number with a median of 1; switch ([ r ]0,rm]And c) is represented in matrix (X)w) In column c of `, r is exchanged0And rmWhen r is an element of0When the value of (A) is not unique, a jump formula corresponding to the first jump is adopted for carrying out multiple variations;
and S353, updating the optimal and worst individuals.
As a preferable example of these, in step S4, the operation sequence matrix (SE) defining the sub-fieldsN×1Comprises the following steps:
s41, according to the airplane operation time matrix, the time matrix is (T)K×NDetermining a time matrix (T') of the operation of the fieldN×1,tj' is:
wherein j is 1,2,3, …, N;
s42, defining the operation sequence matrix of the field block as (SE)N×1,sejIs [1, N ]]A non-repeating set of integers within the interval;
s42, according to the operation Sequence (SE)N×1The time matrix after defining the sequence isIt is composed ofComprises the following steps:
s43, based on (SE)N×1Anddefining the operation time of the plant protection unmanned plane-field block as (TA)K×N’,taij′Indicates the time of j ' th operation of the ith aircraft, j ' is 1,2,3, …, N ' is max { sum [ x (i,:)]};(TA)K×N’Each column of the non-assigned items is completed by 0;
s44, considering the overlapping time of assembling batteries and medicine boxes simultaneously on a plurality of plant protection unmanned planes, defining the landing time matrix of the unmanned planes as (T L)K×N’Takeoff time matrix is (TR)K×N’,tlij′And trij′Are respectively provided withComprises the following steps:
wherein, tsij′The time required for the aircraft to fit to the next takeoff after landing, here set to a fixed length, tsij′=c。
As a preferred example, in step S5, the operation sequence matrix (SE) of the sub-field block is implemented by using genetic algorithmN×1The process of performing the optimization includes the steps of:
s51, evaluating the adaptive value of the sub-field operation:
order of operation (SE)N×1As chromosome, in determining (SE)N×1Under the condition of (3), determining the final landing time of different plant protection unmanned aerial vehicles to be tlKN’Defining an fitness value selection function FV of a chromosome2Comprises the following steps:
wherein g (SEt) max [ tl [ ]KN′(SEt,c)](ii) a t represents the position of the chromosome in the population; wt is a time adaptive value parameter, the unit is min, and the value range is [300,600,];
s52, genetic operator definition:
a. chromosome selection operations
Selecting initial population by roulette method, if fitness of a certain individual t is ftThen the probability that the individual is selected is denoted as Pt:
Wherein N represents the chromosome scale;
b. chromosome crossing operations
Carrying out chromosome crossing operation by adopting a partial matching crossing method: randomly generating two cross points and determining a matching section; generating two sub-individuals according to the mapping relation between the matching segments and the non-matching segments in the parent;
c. operation of chromosomal variation
Randomly selecting N positions, and interchanging genes corresponding to the N positions, wherein N is more than 1;
s53, updating chromosome based on GA algorithm to obtain optimal operation sequence matrix (SE) of sub-fieldN×1The method comprises the following substeps:
s531, inputting population size N, maximum iteration times max (Gen) of evolution, and cross rate pcThe rate of variation pm;
S532, generating an initial population { SE }PThe total population size is P;
s533, calculating the adaptive value of the population individuals;
s534, determining the following conditions for stopping genetic calculation: (1) the iteration number reaches max (Gen), (2) the requirement of an adaptive value is met, (3) the optimal chromosome is continuously kept for 10 generations; if any one of the conditions is satisfied, the iteration process is terminated, and the optimal individual { SE } is outputmax(f)If not, go to step S535;
s535, selecting and calculating other chromosomes except the optimal chromosome;
s536, performing cross operation on other chromosomes except the optimal chromosome;
s537, performing mutation operation on other chromosomes except the optimal chromosome, and returning to step S533.
Based on the distribution method, the invention also provides a plant protection unmanned aerial vehicle cluster operation task distribution device, which comprises an ARM embedded system, and a GPS positioning system, a temperature and humidity sensor, a data transmission module, a power supply module, an L CD touch screen and a cache module which are respectively connected with the ARM embedded system;
the temperature and humidity sensor is used for detecting the real-time temperature and the real-time humidity of the operation area and feeding back the detection result to the ARM embedded system;
the L CD touch screen is used for inputting field parameters and plant protection unmanned aerial vehicle parameters to the ARM embedded system, and the field parameters comprise boundary coordinates of each field;
the ARM embedded system combines input field parameters, plant protection unmanned aerial vehicle parameters and environment information fed back by a temperature and humidity sensor, adopts the distribution method, jointly calculates all plant protection unmanned aerial vehicle cooperative task distribution models and corresponding scheduling schemes of plant protection unmanned aerial vehicle cluster operation, sends the scheduling schemes to the flight control system of each plant protection unmanned aerial vehicle through a data transmission module, and stores the input parameters and the corresponding scheduling schemes to a cache module;
the power supply module is used for providing electric energy required by normal work of the ARM embedded system, the GPS positioning system, the temperature and humidity sensor, the data transmission module, the power supply module, the L CD touch screen and the buffer module;
the GPS is used for acquiring the position information of each plant protection unmanned aerial vehicle and displaying the acquired position information to a user through an L CD touch screen.
Compared with the prior art, the technical scheme of the invention has the following remarkable beneficial effects:
(1) on the basis of ensuring the full coverage of the air route, two constraint conditions of time and energy consumption are integrated, a double-layer decision method for distributing the shortest flight distance and the optimal time is established, time conflict is reduced, and the operation efficiency is improved.
(2) The whole operation process of all plant protection unmanned aerial vehicles is automatically planned, the manpower loss is reduced, the maintenance and the monitoring of all unmanned aerial vehicles can be completed by only one worker, and the true one-control-multiple-machine operation is realized.
It should be understood that all combinations of the foregoing concepts and additional concepts described in greater detail below can be considered as part of the inventive subject matter of this disclosure unless such concepts are mutually inconsistent. In addition, all combinations of claimed subject matter are considered a part of the presently disclosed subject matter.
The foregoing and other aspects, embodiments and features of the present teachings can be more fully understood from the following description taken in conjunction with the accompanying drawings. Additional aspects of the present invention, such as features and/or advantages of exemplary embodiments, will be apparent from the description which follows, or may be learned by practice of specific embodiments in accordance with the teachings of the present invention.
Drawings
The drawings are not intended to be drawn to scale. In the drawings, each identical or nearly identical component that is illustrated in various figures may be represented by a like numeral. For purposes of clarity, not every component may be labeled in every drawing. Embodiments of various aspects of the present invention will now be described, by way of example, with reference to the accompanying drawings, in which:
fig. 1 is a flowchart of a plant protection unmanned aerial vehicle cluster work task allocation method of the present invention.
Fig. 2 is a flow chart of a distribution method according to a first embodiment of the present invention.
FIG. 3 is a flow diagram illustrating the division of a field of the present invention.
FIG. 4 is a MOSF L A based search for an optimal job assignment matrix (X)K×NAnd (4) a flow chart.
FIG. 5 is a schematic diagram of chromosome crossover preparation.
FIG. 6 is a schematic diagram of chromosome crossing completion.
FIG. 7 is a schematic diagram of a chromosomal mutation process.
FIG. 8 is an operation order matrix (SE) for a sub-field block using a genetic algorithmN×1And carrying out an optimization flow chart.
Fig. 9 is a schematic structural diagram of a plant protection unmanned aerial vehicle cluster operation task distribution device.
Fig. 10 is a schematic view of a target work field of a plant protection unmanned aerial vehicle according to a third embodiment.
Fig. 11 is a schematic diagram of a target field block preprocessing method in the third embodiment.
FIG. 12 is a comparison of the results of the third embodiment.
Detailed Description
In order to better understand the technical content of the present invention, specific embodiments are described below with reference to the accompanying drawings.
Detailed description of the preferred embodiment
With reference to fig. 1, the present invention provides a plant protection unmanned aerial vehicle cluster work task allocation method, which includes the following steps:
s1, dividing the operation field into N sub-fields according to the input field parameters and the plant protection unmanned aerial vehicle parameters, wherein the operation requirements of each sub-field are the same, and only corresponding to one unmanned aerial vehicle, the set of N sub-fields is defined as Hj,j=1,2,3,…,N。
S2, distributing all the sub-fields to all the plant protection unmanned aerial vehicles participating in the operation according to a preset distribution principle, wherein the preset distribution principle is as follows: a single plant protection drone is assigned at least one sub-field, each sub-field can only be assigned once and all sub-fields are assigned.
S3, setting the number of unmanned planes, the number of field blocks K, the number of sub-field blocks N and the operation time matrix (T)K×NAnd the distance matrix of the in/out field is (L)1)K×N、(L2)K×NBased on the operation, an MOSE L A algorithm is adopted to find an optimal operation distribution matrix (X)K×N。
S4, distributing matrix (X) according to the optimal operationK×NAnd the aircraft operation time matrix (T)K×NDetermining an operation time matrix (T') of the sub-field blockN×1(ii) a Defining a sequence of operations matrix (SE) for the sub-fields in combination with predetermined assembly rulesN×1And after each unmanned aircraft lands, the preset assembly rule is that assembly tasks including battery replacement and medicine box replacement are carried out, but only one unmanned aircraft executes the assembly tasks at the same moment.
And S5, optimizing a work sequence matrix (SE) N × 1 of the sub-field block by adopting a genetic algorithm.
Referring to fig. 1 and 2, assuming known conditions of ① known field positions, the number of plant protection unmanned planes and take-off and landing positions, ② known flight and operation performance of the planes and unmanned plane operation parameters of each field, ③ only one labor force (one control machine) when the unmanned planes land for medicine box and battery assembly, the process of unmanned plane distribution and field operation sequencing of the operation fields under the coordination of a plurality of plant protection unmanned planes comprises the following steps:
first, repartitioning of the field pieces
Assuming that n blocks of the common target operation field exist, the dividing steps are as follows: (I) putting the field blocks with the same operation requirements (pesticide type, pesticide solution concentration, airplane operation parameters and the like) into the same set, wherein the total number of the field blocks is M; (II) according to the one-time lifting/landing operation area of the plant protection unmanned aerial vehicle, further dividing all the fields in the M sets into fields, merging adjacent small fields and dividing large fields; (III) extracting all sub-field blocks in the M sets, counting as N, redefining the N field block sets into Hj(j ═ 1,2,3, …, N). The dividing step is shown in fig. 3 below.
Secondly, distributing unmanned aerial vehicle operation field blocks based on MOSF L A
A single plant protection drone may be assigned to multiple fields (< N), but each field may only be assigned once.
Defining the operation distribution matrix of plant unmanned plane as (X)K×NWherein a parameter x is assignedijAn assignment is defined as
Where j is 1,2,3, …, N. The unmanned aircraft allocation parameters follow the constraint conditions: each field can only be operated once and all fields are assigned to plant protection drones.
Third, job allocation optimization based on MOSF L A
Given the number of airplanes, the field size K, N, and the operation time matrix (T)K×NThe distance matrix of the in/out field is (L)1)K×N、(L2)K×NThe MOSF L A algorithm is adopted to find the optimal operation distribution matrix as (X)K×N。
Through the co-evolution of multi-clan groups, the mixed frog leaping algorithm can realize the update of the group. To optimize the distribution parameter X in equation (1), the definition of the population update method is as follows.
Definition 2z ═ cut (z, k) is defined as: and intercepting the first k columns of the one-dimensional array z and storing the first k columns into a new one-dimensional array z'.
Definition 3 (X)1) ' - < F (B, φ), based on the difference matrix B (i.e.) The matrix X is updated to approximate the matrix B. And counting the element sum of each column B, and if the element sum is greater than 0, storing the column number into a one-dimensional array z. The definition of the operator φ is:
φ=~(B,X1,cut(z,ri(n))). (2)
wherein: n is the length of the array z, ri (n) denotes the interval [1, n ]]Any integer; "-" indicates: according to the column number ri (n) designated by z' (i.e. cut (z, ri (n))), X is added1Element correction in the column is consistent with B element arrangement
in the formula: switch (r)1,r2) Representing r in a switching matrix1Row and r2Line, swicth (c)1,c2) Representing c in the switching matrix1Column sum c2And (4) columns.
Frog leaping rule of MOSF L A
The rules related to the frog leaping designed in the invention are as follows
The number 1 of the jumps is the number one jump,
wherein: xbAnd XwRespectively, the best individual and the worst individual of the subgroup.
The 2 nd hop, the first hop,
wherein: xgAnd XwRespectively the best and worst individuals of the entire population.
The 3 rd time of the jump is the jump,
(Xw)′=rand(X),X∈QM×N. (6)
wherein: rand (x) represents the random generation of a decision matrix in the feasible solution space.
Therefore, the optimization capability of the worst individual in the sub-population and the total population to the optimal individual in the population and the global updating process of the total population are embodied.
The new individuals generated may have a problem with the aircraft not being assigned to the field (i.e., the aircraft is not assigned to the field)). The invention introduces mutation operators as follows to carry out mutation processing on the partial components:
(Xw)″=F((Xw)′,switch([r0,rm],c)). (7)
in the formula: r is0Representation matrix (X)w) ' line element and line number 0; r ismRepresentation matrix (X)w) ' line element sum is the largest line number; c represents rmAny column number with a median of 1; switch ([ r ]0,rm]And c) is represented in matrix (X)w) In column c of `, r is exchanged0And rmWhen r is an element of0The value of (3) is not exclusive and may be varied several times according to the expression (4).
Group division of individuals
Assuming that the total number of the frog-jump population is S, the number of the sub-populations is S, and the number of the sub-populations is num (S is S × num), an adaptive value function FV of the individual t is defined based on the formula (6)1As follows below, the following description will be given,
the present invention uses the fitness value of an individual as a reference for population division. According to FV1 tIn descending order, the individuals in the ith subgroup are respectively,
{FV1 k|k∈[(i-1)×num+1,i×num]}. (9)
wherein i is 1,2,3, …, s.
According to the algorithm design, the specific steps of the modified MOSF L a algorithm for the multi-machine cooperative allocation problem proposed by the present invention are shown in fig. 4.
Sixth, field operation ordering model based on GA
Based on the MOSF L A algorithm, we have found a superior job assignment matrix as (X)K×N. According to the aircraft operation time matrix as (T)K×NWe can determine the operation time matrix (T') of the field blockN×1T of whichj' to (a) is,
where j is 1,2,3, …, N. Define the operation sequence matrix of the field as (SE)N×1,j=1,2,3,…,N,sejIs [1, N ]]Non-repeating whole within intervalA set of numbers. According to the operation Sequence (SE)N×1The time matrix after defining the sequence isIt is composed ofIn order to realize the purpose,
based on (SE)N×1Anddefining the operation time of the plant protection unmanned plane-field block as (TA)K×N’,taij′Indicates the time of j ' th operation of the ith aircraft, j ' is 1,2,3, …, N ' is max { sum [ x (i,:)]}。(TA)K×N’Each column of unassigned entries is completed with 0.
Considering the overlapping time of assembling batteries and medicine boxes simultaneously for a plurality of plant protection unmanned planes, the landing time matrix of the unmanned plane is defined as (T L)K×N’Takeoff time matrix is (TR)K×N’,tlij′And trij′Respectively as follows:
wherein, tsij′The time required for assembly to the next takeoff after landing of the aircraft, irrespective of the assembly time differences of the assemblers, may be set to a fixed length, ts, hereij′C, simplifying the computational complexity.
Seventh, GA-based field block operation sequencing optimization
To reduce multiple unmannedThe overall time of simultaneous machine operations requires determining the sequence of operations (SE) of a fieldN×1. The invention optimizes the operation sequence by GA algorithm.
7.1 evaluation of adaptive value for field work
The invention combines the operation Sequence (SE)N×1As chromosome, in determining (SE)N×1Under the condition, the final landing time tl of different plant protection unmanned aerial vehicles can be determinedKN’Defining an fitness value selection function FV for the chromosome2The definition is that,
wherein g (SEt) max [ tl [ ]KN′(SEt,c)](ii) a t represents the position of the chromosome in the population; wt is a time adaptive value parameter (min) with a value range of [300,600]。
7.2 genetic operator definition
a. Chromosome selection operations
Selecting initial population by roulette method, if fitness of a certain individual t is ftThe probability that the individual is selected can be expressed as Pt:
Where N represents the chromosome size.
b. Chromosome crossing operations
With reference to fig. 5 and fig. 6, a partial matching Crossover method (PMX) is used to perform Crossover operation of chromosomes: randomly generating two cross points and determining a matching section; and generating two sub-individuals according to the mapping relation between the matching segment and the non-matching segment in the parent. For example, there are two parents, and the repeated numbers are marked with a symbol, assuming that the matching segment generated randomly is [4,6 ]. Based on the intermediate mapping, the number of repeated portions may be determined, the first occurrence of font a1 may be assigned an initial value of 1, the occurrence of intersections of 1 originating from 6, and 6 participating in intersections of which the initial value is 5, 5 no longer participating in intersections, and the first occurrence of font a1 may be assigned a value of 5; the second occurrence of a1 is assigned a value of 4 in the same way. The assignment at child B1 is done in the same way.
c. Operation of chromosomal variation
And (3) randomly selecting N (>1) positions, and interchanging the corresponding genes. For chromosome A, the selected gene positions are 3, 5 and 6, and the compiled chromosome is shown in FIG. 7.
Eighthly, GA-based chromosome renewal process
The conditions for stopping genetic calculation of the chromosome in the invention are ① iteration to max (Gen) times, the algorithm is terminated, ② meets the requirement of an adaptive value, the algorithm is terminated, ③ the optimal chromosome is continuously maintained for 10 generations, and the algorithm is terminated, and the GA-based chromosome updating process is shown in figure 8.
In the invention, the population specification P is 20-100; the iteration number max (Gen) of genetic evolution is 100-500; cross probability pc0.4-0.99; probability of variation pm=0.0001~0.1。
Detailed description of the invention
With reference to fig. 9, based on the foregoing allocation method, the present invention further provides a plant protection unmanned aerial vehicle cluster operation task allocation device, where the allocation device includes an ARM embedded system, and a GPS positioning system, a temperature and humidity sensor, a data transmission module, a power supply module, an L CD touch screen, and a cache module, which are respectively connected to the ARM embedded system.
The temperature and humidity sensor is used for detecting the real-time temperature and the real-time humidity of the operation area and feeding back the detection result to the ARM embedded system.
The L CD touch screen is used for inputting field parameters and plant protection unmanned aerial vehicle parameters to the ARM embedded system, and the field parameters comprise boundary coordinates of each field.
The ARM embedded system combines the input field parameters, the plant protection unmanned aerial vehicle parameters and the environment information fed back by the temperature and humidity sensor, adopts the distribution method, jointly calculates all plant protection unmanned aerial vehicle cooperative task distribution models and the corresponding scheduling schemes of the plant protection unmanned aerial vehicle cluster operation, and sends the scheduling schemes to the flight control system of each plant protection unmanned aerial vehicle through the data transmission module, so that the flight control system controls machines and tools to perform cluster operation according to the received instructions, and the input parameters and the corresponding scheduling schemes are stored in the cache module.
The power supply module is used for providing electric energy required by normal work of the ARM embedded system, the GPS positioning system, the temperature and humidity sensor, the data transmission module, the power supply module, the L CD touch screen and the buffer module.
The GPS is used for acquiring the position information of each plant protection unmanned aerial vehicle and displaying the acquired position information to a user through an L CD touch screen.
Detailed description of the preferred embodiment
In order to determine the optimization performance of the multi-plant protection unmanned aerial vehicle operation planning based on the MOSF L A + GA collaborative task allocation algorithm, the target multi-field block with the total area of 330 mu is selected by the invention as shown in the following figures 10 and 11.
The maximum cruising time of the selected plant protection unmanned aerial vehicle is 25min, the drug loading is 15L, the plant protection unmanned aerial vehicle is divided again according to the multiple fields in the figure 10 based on the unmanned aerial vehicle, the division result is shown in figure 11, the area of the divided 22 fields is less than or equal to 15 mu, so that the plant protection unmanned aerial vehicle can complete spraying operation once for one frame, the rising/landing area of the unmanned aerial vehicle is positioned near the center of the fields, the area simultaneously meets the requirements of being convenient to take off, transporting liquid medicine batteries and the like, when multiple aircrafts land at the same time, the maximum distance between the aircrafts is 100m, so that the assembling work of the liquid medicine and the batteries is ensured to be completed within the specified time, in the invention, the position of 100m at the center of the fields is selected as the rising and landing assembling point of the aircrafts, the distance between the adjacent landing points is 25m, the assembling time c is 6min, and the.
In order to compare the weighted distance and the total flight time of the airplane entering/exiting the field under the operation of the unmanned aerial vehicle under the normal condition, the invention compares the following two traditional multi-machine operation modes:
in the conventional multi-aircraft operation mode ①, the aircraft is not sure of the field of operation, but the aircraft will perform the next field of operation after completing the previous field of operation.
Conventional multi-aircraft operation mode ②. unmanned aircraft have defined fields, each with a fixed sequence of operations.
The results of two conventional multi-machine operation modes are as follows:
compared with the traditional mode, the method can save 25min at most by adopting the MOSFA L + GA algorithm, and the working efficiency and the total time of the unmanned aerial vehicle cluster are shown in figure 12.
In this disclosure, aspects of the present invention are described with reference to the accompanying drawings, in which a number of illustrative embodiments are shown. Embodiments of the present disclosure are not necessarily defined to include all aspects of the invention. It should be appreciated that the various concepts and embodiments described above, as well as those described in greater detail below, may be implemented in any of numerous ways, as the disclosed concepts and embodiments are not limited to any one implementation. In addition, some aspects of the present disclosure may be used alone, or in any suitable combination with other aspects of the present disclosure.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to be limited thereto. Those skilled in the art can make various changes and modifications without departing from the spirit and scope of the invention. Therefore, the protection scope of the present invention should be determined by the appended claims.
Claims (10)
1. A plant protection unmanned aerial vehicle cluster operation task distribution method is characterized by comprising the following steps:
s1, dividing the operation field into N sub-fields according to the input field parameters and the plant protection unmanned aerial vehicle parameters, wherein the operation requirements of each sub-field are the same, and only corresponding to one unmanned aerial vehicle, the set of N sub-fields is defined as Hj,j=1,2,3,…,N;
S2, distributing all the sub-fields to all the plant protection unmanned aerial vehicles participating in the operation according to a preset distribution principle, wherein the preset distribution principle is as follows: a single plant protection unmanned aerial vehicle is allocated with at least one sub-field block, each sub-field block can be allocated only once, and all the sub-field blocks are allocated;
s3, setting the number of unmanned planes, the number of field blocks K and the number of sub-field blocks in the plant protection unmanned planeQuantity N, and the operation time matrix (T)K×NAnd the distance matrix of the in/out field is (L)1)K×N、(L2)K×NBased on the operation, an MOSE L A algorithm is adopted to find an optimal operation distribution matrix (X)K×N;
S4, distributing matrix (X) according to the optimal operationK×NAnd the aircraft operation time matrix (T)K×NDetermining an operation time matrix (T') of the sub-field blockN×1(ii) a Defining a sequence of operations matrix (SE) for the sub-fields in combination with predetermined assembly rulesN×1The preset assembly rule is that after each unmanned aircraft lands, assembly tasks including battery replacement and medicine chest replacement are carried out, but only one unmanned aircraft executes the assembly tasks at the same time;
s5, adopting genetic algorithm to work sequence matrix (SE) of the sub-farmland blocksN×1And (6) optimizing.
2. The plant protection unmanned aerial vehicle cluster operation task allocation method according to claim 1, wherein in step S1, the process of dividing the operation field into N sub-fields according to the input field parameters and plant protection unmanned aerial vehicle parameters comprises the following steps:
s11, putting the fields with the same operation requirements into the same set, wherein the number of the fields is M;
s12, further field processing is carried out on all the fields in the M sets according to the one-time lifting/landing operation area of the plant protection unmanned aerial vehicle, and the field processing comprises the operations of adjacent small field merging and large field dividing;
s13, extracting all the sub-field blocks in the M sets, setting the number of the sub-field blocks to be N, and defining the N sub-field block sets to be Hj,j=1,2,3,…,N。
3. The method for assigning task tasks to cluster plant protection unmanned aerial vehicle as claimed in claim 1, wherein in step S1 and step S2, the assigning all sub-fields to all plant protection unmanned aerial vehicles participating in operation according to the preset assignment rule is:
defining the operation distribution matrix of plant unmanned plane as (X)K×NWherein a parameter x is assignedijThe assignment is defined as:
where j is 1,2,3, …, N.
4. The plant protection unmanned aerial vehicle cluster task allocation method of claim 3, wherein in step S1 and step S3, the MOSE L A algorithm is used to find the optimal task allocation matrix (X)K×NComprises the following steps:
s31, inputting population scale S, sub-population number S, individual frog num of the population, total population iteration number D and population iteration number I;
s32, initializing the frog group P ═ X1,X2,…,XS};
S33, calculating adaptive values, arranging the adaptive values in a descending order, and dividing the group;
s34, intra-population evolution: each sub-population makes a first jump according to the following formula:
wherein: xbAnd XwRespectively, the optimal individual and the worst individual of the sub-group;
s35, merging of first population: and (4) performing descending arrangement according to the adaptive value, and combining the following formula, and performing second jump on the whole group:
wherein: xgAnd XwRespectively, the optimal individual and the worst individual of the whole population;
s36, second population merging: and (4) performing descending arrangement according to the adaptive value, and performing third jump on the weakest individuals of the population by combining the following formula:
(Xw)′=rand(X),X∈QM×N
wherein: rand (X) represents randomly generating a decision matrix in the feasible solution space;
s37, judging whether the current cycle times are less than the total group iteration times D, if so, turning to the step S33, otherwise, outputting the optimal individual X1;
Wherein, the related parameters of population update are defined as follows:
let matrix X1∈QM×N、X2∈QM×NWherein Q isM×NFor a feasible solution space, a matrix is definedRepresents X2And X1The matrix of the differences between the two,indicating an ectopic or operational point between matrices;
the meaning of z ═ cut (z, k) is: intercepting the front k columns of the one-dimensional array z and storing the front k columns into a new one-dimensional array z';
(X1) ' -F (B, Φ) means: updating an approximation matrix B for the matrix X according to the difference matrix B, counting element sums of all columns of B, if the element sums are larger than 0, storing the column number into a one-dimensional array z, and defining an operator phi as follows:
φ=~(B,X1,cut(z,ri(n)))
wherein: n is the length of the array z, ri (n) denotes the interval [1, n ]]Any integer within; "-" indicates that X is assigned to the column number ri (n) designated by z1Element correction in the column is consistent with the arrangement of the B element;
the meaning of (A) is: for matrix X1According to the operatorPerforming mutation operation, operatorIs defined as:
in the formula: switch (r)1,r2) Representing r in a switching matrix1Row and r2Line, swicth (c)1,c2) Representing c in the switching matrix1Column sum c2And (4) columns.
5. The plant protection unmanned aerial vehicle cluster task allocation method of claim 4, wherein in step S1, in step S33, the process of performing cluster division comprises the steps of:
s331, defining an adaptive value function FV of an individual t1As follows below, the following description will be given,
s332, taking the adaptive value of the individual as a reference for dividing the group according to FV1 tIn descending order, the individuals in the ith subgroup are:
{FV1 k|k∈[(i-1)×num+1,i×num]}
wherein i is 1,2,3, …, s.
6. The method of claim 4, wherein in step S1 and step S34, the process of evolving inside the population comprises the following steps:
s341, jumping to the optimal individual of the sub-population according to the following formula:
s342, introducing a mutation operator for detecting whether the plant protection unmanned aerial vehicle which is not distributed to the field exists, and performing mutation detection and processing by adopting the mutation operator:
(Xw)″=F((Xw)′,switch([r0,rm],c))
in the formula: r is0Representation matrix (X)w) ' line element and line number 0; r ismRepresentation matrix (X)w) ' line element sum is the largest line number; c represents rmAny column number with a median of 1; switch ([ r ]0,rm]And c) is represented in matrix (X)w) In column c of `, r is exchanged0And rmWhen r is an element of0When the value of (A) is not unique, a jump formula corresponding to the first jump is adopted for carrying out multiple variations;
s343, updating the optimal and worst individuals;
and S344, repeating the steps S341 to S343 until the cycle number reaches I times.
7. The method for assigning task tasks to cluster unmanned aerial vehicle for plant protection as claimed in claim 4, wherein in step S1, step S35, the process of merging the first population comprises the following steps:
s351, jumping to the population-optimized individual using the following formula:
s352, introducing a mutation operator for detecting whether the plant protection unmanned aerial vehicle which is not distributed to the field exists, and performing mutation detection and processing by adopting the mutation operator:
(Xw)″=F((Xw)′,switch([r0,rm],c))
in the formula: r is0Representation matrix (X)w) ' line element and line number 0; r ismRepresentation matrix (X)w) ' line element sum is the largest line number; c represents rmAny column number with a median of 1; switch ([ r ]0,rm]And c) is represented in matrix (X)w) In column c of `, r is exchanged0And rmWhen r is an element of0When the value of (A) is not unique, a jump formula corresponding to the first jump is adopted for carrying out multiple variations;
and S353, updating the optimal and worst individuals.
8. The plant protection unmanned aerial vehicle cluster task allocation method of claim 4, wherein in step S4, the operation sequence matrix (SE) defining the sub-fieldsN×1Comprises the following steps:
s41, according to the airplane operation time matrix, the time matrix is (T)K×NDetermining a time matrix (T') of the operation of the fieldN×1,tj' is:
wherein j is 1,2,3, …, N;
s42, defining the operation sequence matrix of the field block as (SE)N×1,sejIs [1, N ]]A non-repeating set of integers within the interval;
s42, according to the operation Sequence (SE)N×1The time matrix after defining the sequence isIt is composed ofComprises the following steps:
s43, based on (SE)N×1Anddefining the operation time of the plant protection unmanned plane-field block as (TA)K×N’,taij′Indicates the time of j ' th operation of the ith aircraft, j ' is 1,2,3, …, N ' is max { sum [ x (i,:)]};(TA)K×N’Each column of the non-assigned items is completed by 0;
s44, considering the overlapping time of assembling batteries and medicine boxes simultaneously on a plurality of plant protection unmanned planes, defining the landing time matrix of the unmanned planes as (T L)K×N’Takeoff time matrix is (TR)K×N’,tlij′And trij′Respectively as follows:
wherein, tsij′The time required for the aircraft to fit to the next takeoff after landing, here set to a fixed length, tsij′=c。
9. The method for assigning task tasks to cluster unmanned aerial vehicle for plant protection as claimed in claim 8, wherein in step S5, the application of genetic algorithm to operation sequence matrix (SE) of sub-field blockN×1The process of performing the optimization includes the steps of:
s51, evaluating the adaptive value of the sub-field operation:
order of operation (SE)N×1As chromosome, in determining (SE)N×1Under the condition of (3), determining the final landing time of different plant protection unmanned aerial vehicles to be tlKN’Defining an fitness value selection function FV of a chromosome2Comprises the following steps:
wherein, g (SE)t)=max[tlKN′(SEt,c)](ii) a t represents the position of the chromosome in the population; wt is a time adaptive value parameter, the unit is min, and the value range is [300,600,];
s52, genetic operator definition:
a. chromosome selection operations
Selecting initial population by roulette method, if fitness of a certain individual t is ftThen the probability that the individual is selected is denoted as Pt:
Wherein N represents the chromosome scale;
b. chromosome crossing operations
Carrying out chromosome crossing operation by adopting a partial matching crossing method: randomly generating two cross points and determining a matching section; generating two sub-individuals according to the mapping relation between the matching segments and the non-matching segments in the parent;
c. operation of chromosomal variation
Randomly selecting N positions, and interchanging genes corresponding to the N positions, wherein N is more than 1;
s53, updating chromosome based on GA algorithm to obtain optimal operation sequence matrix (SE) of sub-fieldN×1The method comprises the following substeps:
s531, inputting population size N, maximum iteration times max (Gen) of evolution, and cross rate pcThe rate of variation pm;
S532, generating an initial population { SE }PThe total population size is P;
s533, calculating the adaptive value of the population individuals;
s534, determining the following conditions for stopping genetic calculation: (1) the iteration number reaches max (Gen), (2) the requirement of an adaptive value is met, (3) the optimal chromosome is continuously kept for 10 generations; if any one of the conditions is satisfied, the iteration process is terminated, and the optimal individual { SE } is outputmax(f)Otherwise, go to step S535;
s535, selecting and calculating other chromosomes except the optimal chromosome;
s536, performing cross operation on other chromosomes except the optimal chromosome;
s537, performing mutation operation on other chromosomes except the optimal chromosome, and returning to step S533.
10. The plant protection unmanned aerial vehicle cluster operation task distribution device is characterized by comprising an ARM embedded system, and a GPS positioning system, a temperature and humidity sensor, a data transmission module, a power supply module, an L CD touch screen and a cache module which are respectively connected with the ARM embedded system;
the temperature and humidity sensor is used for detecting the real-time temperature and the real-time humidity of the operation area and feeding back the detection result to the ARM embedded system;
the L CD touch screen is used for inputting field parameters and plant protection unmanned aerial vehicle parameters to the ARM embedded system, and the field parameters comprise boundary coordinates of each field;
the ARM embedded system combines input field parameters, plant protection unmanned aerial vehicle parameters and environment information fed back by a temperature and humidity sensor, adopts the distribution method as any one of claims 1 to 9, jointly calculates all plant protection unmanned aerial vehicle cooperative task distribution models and corresponding scheduling schemes of plant protection unmanned aerial vehicle cluster operation, sends the scheduling schemes to the flight control system of each plant protection unmanned aerial vehicle through a data transmission module, and stores the input parameters and the corresponding scheduling schemes to a cache module;
the power supply module is used for providing electric energy required by normal work of the ARM embedded system, the GPS positioning system, the temperature and humidity sensor, the data transmission module, the power supply module, the L CD touch screen and the buffer module;
the GPS is used for acquiring the position information of each plant protection unmanned aerial vehicle and displaying the acquired position information to a user through an L CD touch screen.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010370842.1A CN111461464B (en) | 2020-05-06 | 2020-05-06 | Plant protection unmanned aerial vehicle cluster operation task allocation method and device |
AU2020416827A AU2020416827B2 (en) | 2020-05-06 | 2020-09-03 | Cluster job task assignment method and device for plant-protection unmanned aerial vehicles |
PCT/CN2020/113196 WO2021135336A1 (en) | 2020-05-06 | 2020-09-03 | Cluster job task assignment method and device for plant-protection unmanned aerial vehicles |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010370842.1A CN111461464B (en) | 2020-05-06 | 2020-05-06 | Plant protection unmanned aerial vehicle cluster operation task allocation method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111461464A true CN111461464A (en) | 2020-07-28 |
CN111461464B CN111461464B (en) | 2021-06-08 |
Family
ID=71679602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010370842.1A Active CN111461464B (en) | 2020-05-06 | 2020-05-06 | Plant protection unmanned aerial vehicle cluster operation task allocation method and device |
Country Status (3)
Country | Link |
---|---|
CN (1) | CN111461464B (en) |
AU (1) | AU2020416827B2 (en) |
WO (1) | WO2021135336A1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112799432A (en) * | 2021-04-08 | 2021-05-14 | 北京三快在线科技有限公司 | Obstacle avoidance control method and device for unmanned aerial vehicle, storage medium and electronic equipment |
CN112949292A (en) * | 2021-01-21 | 2021-06-11 | 中国人民解放军61540部队 | Cluster unmanned aerial vehicle return data processing method and device, equipment and storage medium |
WO2021135336A1 (en) * | 2020-05-06 | 2021-07-08 | 农业农村部南京农业机械化研究所 | Cluster job task assignment method and device for plant-protection unmanned aerial vehicles |
CN114021914A (en) * | 2021-10-22 | 2022-02-08 | 北京市农林科学院信息技术研究中心 | Unmanned aerial vehicle cluster flight control method and device |
CN116050779A (en) * | 2023-01-16 | 2023-05-02 | 农业农村部南京农业机械化研究所 | Plant protection unmanned aerial vehicle dynamic scheduling method based on train-dimensional simulated annealing algorithm |
CN116542468A (en) * | 2023-05-06 | 2023-08-04 | 中国人民解放军32370部队 | Unmanned aerial vehicle cluster task planning method |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114047777A (en) * | 2021-07-31 | 2022-02-15 | 海南微飞信息科技有限公司 | Mass operation scheduling method for plant protection unmanned aerial vehicle |
CN114020022B (en) * | 2021-11-04 | 2023-08-18 | 中国人民解放军陆军工程大学 | Heterogeneous unmanned aerial vehicle collaborative hit task planning method and device |
CN114020034A (en) * | 2021-11-26 | 2022-02-08 | 哈尔滨工程大学 | Task slice-based heterogeneous multi-unmanned aerial vehicle cooperative task allocation method |
CN114326827B (en) * | 2022-01-12 | 2023-06-09 | 北方工业大学 | Unmanned aerial vehicle cluster multitasking dynamic allocation method and system |
CN114594794B (en) * | 2022-03-08 | 2023-05-23 | 大连理工大学 | Multi-machine collaborative task planning method considering subsystem execution capacity |
CN115379386B (en) * | 2022-07-18 | 2024-05-14 | 深圳市人工智能与机器人研究院 | Unmanned plane cluster electromagnetic spectrum monitoring method and device based on wolf's swarm algorithm |
CN115309191B (en) * | 2022-09-21 | 2024-07-19 | 中国人民解放军国防科技大学 | EMARL unmanned aerial vehicle clustering method and device based on competition cooperation mechanism |
CN115454614B (en) * | 2022-11-14 | 2023-03-24 | 安徽大学 | Intelligent scheduling method for robot cluster energy supply |
CN115879603B (en) * | 2022-11-17 | 2024-05-14 | 武汉大学 | Multi-target-point-oriented multi-unmanned aerial vehicle collaborative data acquisition method and device |
CN116245346B (en) * | 2023-05-12 | 2023-07-28 | 深圳大学 | Multi-unmanned aerial vehicle task allocation method based on multiple swarm genetic algorithms and local search |
CN116483128B (en) * | 2023-06-19 | 2023-10-27 | 湖南林科达信息科技有限公司 | Unmanned aerial vehicle multitasking load device conversion method and system |
CN118195277B (en) * | 2024-05-16 | 2024-07-23 | 中国人民解放军海军航空大学 | Helicopter group offshore platform operation integrated scheduling method |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597850A (en) * | 2016-12-16 | 2017-04-26 | 新疆疆天航空科技有限公司 | Plant protection unmanned plane formation object distribution method based on chaotic leapfrog |
CN106873629A (en) * | 2017-04-14 | 2017-06-20 | 合肥工业大学 | Unmanned plane aviation job task distribution method and device |
CN106997209A (en) * | 2016-01-25 | 2017-08-01 | 深圳市鼎创旭飞科技有限公司 | Plant protection unmanned plane sprays operational method and system |
CN107037827A (en) * | 2017-04-14 | 2017-08-11 | 合肥工业大学 | Unmanned plane aviation job task is distributed and trajectory planning combined optimization method and device |
US10140875B1 (en) * | 2017-05-27 | 2018-11-27 | Hefei University Of Technology | Method and apparatus for joint optimization of multi-UAV task assignment and path planning |
CN109760846A (en) * | 2017-11-09 | 2019-05-17 | 湖南农业大学 | A kind of plant protection drone field automatically supplying apparatus and method |
CN109814597A (en) * | 2019-02-03 | 2019-05-28 | 唐山坤翼创新科技有限公司 | The control method of concentrating type plant protection drone control system |
CN110288075A (en) * | 2019-06-26 | 2019-09-27 | 电子科技大学 | A kind of feature selection approach based on improvement shuffled frog leaping algorithm |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SG11201806440WA (en) * | 2016-01-29 | 2018-08-30 | Garuda Robotics Pte Ltd | System and method for controlling an unmanned vehicle and releasing a payload from the same |
CN106054917A (en) * | 2016-05-27 | 2016-10-26 | 广州极飞电子科技有限公司 | Unmanned aerial vehicle flight control method and device, and remote controller |
CN106547276B (en) * | 2016-10-19 | 2019-07-26 | 上海圣尧智能科技有限公司 | Automatic spraying rectangular-ambulatory-plane paths planning method and fog machine spraying operation method |
CN107633105B (en) * | 2017-07-12 | 2020-09-08 | 西北工业大学 | Improved hybrid frog-leaping algorithm-based quad-rotor unmanned aerial vehicle parameter identification method |
CN111461464B (en) * | 2020-05-06 | 2021-06-08 | 农业农村部南京农业机械化研究所 | Plant protection unmanned aerial vehicle cluster operation task allocation method and device |
-
2020
- 2020-05-06 CN CN202010370842.1A patent/CN111461464B/en active Active
- 2020-09-03 WO PCT/CN2020/113196 patent/WO2021135336A1/en active Application Filing
- 2020-09-03 AU AU2020416827A patent/AU2020416827B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106997209A (en) * | 2016-01-25 | 2017-08-01 | 深圳市鼎创旭飞科技有限公司 | Plant protection unmanned plane sprays operational method and system |
CN106597850A (en) * | 2016-12-16 | 2017-04-26 | 新疆疆天航空科技有限公司 | Plant protection unmanned plane formation object distribution method based on chaotic leapfrog |
CN106873629A (en) * | 2017-04-14 | 2017-06-20 | 合肥工业大学 | Unmanned plane aviation job task distribution method and device |
CN107037827A (en) * | 2017-04-14 | 2017-08-11 | 合肥工业大学 | Unmanned plane aviation job task is distributed and trajectory planning combined optimization method and device |
US10140875B1 (en) * | 2017-05-27 | 2018-11-27 | Hefei University Of Technology | Method and apparatus for joint optimization of multi-UAV task assignment and path planning |
CN109760846A (en) * | 2017-11-09 | 2019-05-17 | 湖南农业大学 | A kind of plant protection drone field automatically supplying apparatus and method |
CN109814597A (en) * | 2019-02-03 | 2019-05-28 | 唐山坤翼创新科技有限公司 | The control method of concentrating type plant protection drone control system |
CN110288075A (en) * | 2019-06-26 | 2019-09-27 | 电子科技大学 | A kind of feature selection approach based on improvement shuffled frog leaping algorithm |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021135336A1 (en) * | 2020-05-06 | 2021-07-08 | 农业农村部南京农业机械化研究所 | Cluster job task assignment method and device for plant-protection unmanned aerial vehicles |
CN112949292A (en) * | 2021-01-21 | 2021-06-11 | 中国人民解放军61540部队 | Cluster unmanned aerial vehicle return data processing method and device, equipment and storage medium |
CN112949292B (en) * | 2021-01-21 | 2024-04-05 | 中国人民解放军61540部队 | Method, device, equipment and storage medium for processing return data of cluster unmanned aerial vehicle |
CN112799432A (en) * | 2021-04-08 | 2021-05-14 | 北京三快在线科技有限公司 | Obstacle avoidance control method and device for unmanned aerial vehicle, storage medium and electronic equipment |
CN114021914A (en) * | 2021-10-22 | 2022-02-08 | 北京市农林科学院信息技术研究中心 | Unmanned aerial vehicle cluster flight control method and device |
CN116050779A (en) * | 2023-01-16 | 2023-05-02 | 农业农村部南京农业机械化研究所 | Plant protection unmanned aerial vehicle dynamic scheduling method based on train-dimensional simulated annealing algorithm |
CN116050779B (en) * | 2023-01-16 | 2024-01-30 | 农业农村部南京农业机械化研究所 | Plant protection unmanned aerial vehicle dynamic scheduling method based on train-dimensional simulated annealing algorithm |
CN116542468A (en) * | 2023-05-06 | 2023-08-04 | 中国人民解放军32370部队 | Unmanned aerial vehicle cluster task planning method |
CN116542468B (en) * | 2023-05-06 | 2023-10-20 | 中国人民解放军32370部队 | Unmanned aerial vehicle cluster task planning method |
Also Published As
Publication number | Publication date |
---|---|
WO2021135336A1 (en) | 2021-07-08 |
AU2020416827B2 (en) | 2022-07-07 |
CN111461464B (en) | 2021-06-08 |
AU2020416827A1 (en) | 2021-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111461464B (en) | Plant protection unmanned aerial vehicle cluster operation task allocation method and device | |
CN111862579B (en) | Taxi scheduling method and system based on deep reinforcement learning | |
CN107977743B (en) | Multi-unmanned aerial vehicle cooperative task allocation method and device | |
CN112766813B (en) | Space-day collaborative observation complex task scheduling method | |
CN107886201B (en) | Multi-objective optimization method and device for multi-unmanned aerial vehicle task allocation | |
CN113009934A (en) | Multi-unmanned aerial vehicle task dynamic allocation method based on improved particle swarm optimization | |
CN110058573B (en) | Huff and puff flexible intelligent assembly logistics path planning platform | |
CN107330588B (en) | Task planning method for cooperative reconnaissance of multi-base heterogeneous unmanned aerial vehicle | |
CN109102124B (en) | Dynamic multi-target multi-path induction method and system based on decomposition and storage medium | |
CN111722643A (en) | Unmanned aerial vehicle cluster dynamic task allocation method imitating wolf colony cooperative hunting mechanism | |
CN108038306A (en) | A kind of method of power boiler burning modeling and optimization towards magnanimity high dimensional data | |
CN106843997A (en) | A kind of parallel virtual machine polymerization based on Spark with optimization MBBO algorithms | |
CN113485409B (en) | Geographic fairness-oriented unmanned aerial vehicle path planning and distribution method and system | |
CN110442143A (en) | A kind of unmanned plane situation data clustering method based on combination multiple target dove group's optimization | |
CN111121784A (en) | Unmanned reconnaissance aircraft route planning method | |
CN115202400A (en) | Unmanned aerial vehicle cluster task planning method based on self-adaptive penalty TAEA | |
CN113689122A (en) | Multi-target combination scheduling method considering transportation equipment in cloud-edge collaborative environment | |
CN113792933A (en) | Workshop layout optimization method and device, electronic equipment and storage medium | |
CN113821972A (en) | Multi-robot collaborative assembly line balancing method considering energy efficiency | |
CN114021914B (en) | Unmanned aerial vehicle cluster flight protection scheduling method and device | |
CN109889366A (en) | Network flow accrual accounting, analysis method and system | |
CN115826614A (en) | Multi-unmanned aerial vehicle energy support task allocation method | |
Zhun et al. | An improved pigeon-inspired optimization algorithm for solving dynamic facility layout problem with uncertain demand | |
CN113487264A (en) | Logistics distribution method and system based on heterogeneous multi-unmanned aerial vehicles | |
CN114417539A (en) | Mine ventilation network graph automatic drawing method based on improved differential evolution algorithm |
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 |