Next Article in Journal
An Improved YOLOv7 Model for Surface Damage Detection on Wind Turbine Blades Based on Low-Quality UAV Images
Previous Article in Journal
Cooperative Drone Transportation of a Cable-Suspended Load: Dynamics and Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-UAV Path Planning Based on Cooperative Co-Evolutionary Algorithms with Adaptive Decision Variable Selection

College of Systems Engineering, National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Drones 2024, 8(9), 435; https://doi.org/10.3390/drones8090435
Submission received: 12 July 2024 / Revised: 16 August 2024 / Accepted: 23 August 2024 / Published: 26 August 2024
(This article belongs to the Section Drone Communications)

Abstract

:
When dealing with UAV path planning problems, evolutionary algorithms demonstrate strong flexibility and global search capabilities. However, as the number of UAVs increases, the scale of the path planning problem grows exponentially, leading to a significant rise in computational complexity. The Cooperative Co-Evolutionary Algorithm (CCEA) effectively addresses this issue through its divide-and-conquer strategy. Nonetheless, the CCEA needs to find a balance between computational efficiency and algorithmic performance while also resolving convergence difficulties arising from the increased number of decision variables. Moreover, the complex interrelationships between the decision variables of each UAV add to the challenge of selecting appropriate decision variables. To tackle this problem, we propose a novel collaborative algorithm called CCEA-ADVS. This algorithm reduces the difficulty of the problem by decomposing high-dimensional variables into sub-variables for collaborative optimization. To improve the efficiency of decision variable selection in the collaborative algorithm and to accelerate the convergence speed, an adaptive decision variable selection strategy is introduced. This strategy selects decision variables according to the order of solving single-UAV constraints and multi-UAV constraints, reducing the cost of the optimization objective. Furthermore, to improve computational efficiency, a two-stage evolutionary optimization process from coarse to fine is adopted. Specifically, the Adaptive Differential Evolution with Optional External Archive algorithm (JADE) is first used to optimize the waypoints of the UAVs to generate a basic path, and then, the Dubins algorithm is combined to optimize the trajectory, yielding the final flight path. The experimental results show that in four different scenarios involving 40 UAVs, the CCEA-ADVS algorithm significantly outperforms the Grey Wolf Optimizer (GWO), Particle Swarm Optimization (PSO), Artificial Bee Colony (ABC), and JADE algorithms in terms of path performance, running time, computational efficiency, and convergence speed. In addition, in large-scale experiments involving 500 UAVs, the algorithm also demonstrates good adaptability, stability, and scalability.

1. Introduction

As the technology of unmanned aerial vehicles (UAVs) matures, their application in fields such as military reconnaissance, disaster relief, and environmental monitoring is becoming increasingly widespread [1]. To achieve a higher level of autonomous flight, efficient multi-UAV path planning algorithms are essential. Multi-UAV path planning involves complex issues such as coordinating movements, avoiding obstacles, and optimizing paths for multiple UAVs. Through collaborative path planning, this approach can better meet the specific requirements of various missions, ensuring UAVs can safely and effectively adapt to changing environments and quickly generate new paths. More importantly, as the number of UAVs increases, the complexity of path planning also increases accordingly. Effective multi-UAV path planning algorithms can better adapt to the growing scale of UAV groups. In Table 1, we present a comprehensive overview of representative evolutionary algorithms. Additionally, these algorithms are summarized in terms of the algorithm’s authors, names, establishment times, and other aspects.
The UAV path planning problem is highly complex. Although researchers have proposed traditional heuristic methods such as mixed-integer linear programming [2,3,4], A* [5,6], RRT [7], and nonlinear programming [8], these methods have certain limitations when addressing this problem. They typically require discretization of the solution space and, in some cases, linearization of the model to simplify the problem. However, such processing leads to a lack of flexibility in these methods, making it difficult to deal with additional constraints or objectives. In contrast, metaheuristic algorithms demonstrate higher flexibility and effectiveness. These methods avoid the need to construct a time-consuming configuration space and can easily handle various constraints, such as threat areas, flight altitude restrictions, and turning radius requirements, as well as objectives. They also possess strong global search capabilities, which help to avoid falling into local optimal solutions. Current metaheuristic algorithms include Differential Evolution (DE) [9,10], Ant Colony Optimization (ACO) [11,12], Particle Swarm Optimization (PSO) [13,14,15], Genetic Algorithm (GA) [16,17], Artificial Bee Colony (ABC) [18], Firefly Algorithm (FA) [19,20], Grey Wolf Optimizer (GWO) [21,22], Artificial Potential Field (APF) [23], and Teaching–Learning-Based Optimization (TLBO) [24].
Table 1. Related work on evolutionary algorithms.
Table 1. Related work on evolutionary algorithms.
CategoryAuthorAlgorithmYear
DEFu [25]DE + QPSO2013
Yu [26]DE + GWO2013
Adhikari [27]Adaptive Fuzzy DE algorithm2017
Ali [28]DE + MMACO2023
RRTFathy S. Elkazzaz [29]DE + RRT2017
Farzad Kiani [30]GWO + RRT2021
Fan Yang [31]ACO + RRT2022
MF Aslan [32]PSO + GDRRT*2023
CCEAWu [33]CCEA2023
Shao [34]CCEA + DCPSO2019
Kesireddy [35]CCEA + NSGA-II∖NSGA-III2019
OthersFang [23]APF2023
Goel [19]FA2018
Chen [11]ACO + GA2021
Majumder [24]MOTLBO2021
Shao [13]PSO2020
The Differential Evolution (DE) algorithm, among other metaheuristic algorithms, has shown excellent performance in multi-UAV path planning. The DE algorithm guides the search through the differential vector between individuals, enabling it to escape local optima more effectively and enhance global search capabilities. Additionally, its mutation strategy helps maintain population diversity throughout the evolutionary process. In recent years, numerous researchers have attempted to address path planning problems using the DE algorithm and its variants. Fu et al. [25] proposed an algorithm that hybridizes DE with Quantum-Based Particle Swarm Optimization (QPSO), termed DEQPSO. Yu and Wang [26] proposed a Hybrid Group Search Optimizer (GSO) integrated with DE, utilizing Grey Wolf Optimization (GWO) to update flight paths through searching angles and distances, with DE employed to refine feasible UAV paths within the search area. Adhikari and Kim [27] proposed an Adaptive Fuzzy DE algorithm for UAV 3D path planning. In 2023, Ali et al. [28] proposed an innovative metaheuristic algorithm, a hybrid of Differential Evolution (DE) and Maximum–Minimum Ant Colony Optimization (MMACO). This new algorithm aims to improve upon the traditional DE algorithm by mitigating the blank areas that may arise during the optimization process of classical Ant Colony Optimization (ACO) and its enhanced version MMACO. Furthermore, the RRT algorithm has a distinct advantage in terms of solution speed due to its ability to quickly generate local optima. Many researchers have conducted studies on the RRT algorithm, resulting in numerous variants that integrate with metaheuristic algorithms. MF Aslan et al. [32] combined PSO with GDRRT* to propose a PSO-GDRRT* algorithm for rapidly finding the optimal path. Fan Yang et al. [31] addressed the local optimum issue of RRT by integrating Ant Colony Optimization with RRT, enabling the algorithm to achieve an ideal path plan. Farzad Kiani et al. [30] combined the RRT algorithm with three Grey Wolf Optimizers and their variants, leveraging the characteristics of RRT sampling and metaheuristic algorithms to identify the optimal path. F. S. Elkazzaz et al. [29] introduced a hybrid algorithm based on RRT* and Differential Evolution to solve the optimized path planning problem.
One notable characteristic of multi-UAV path planning is its high variable dimensionality. For instance, when dealing with 100 UAVs, if each UAV has 10 waypoints and each waypoint is represented by three variables, the total number of variables would reach 3000. To address this issue, an effective approach is to employ the Cooperative Co-Evolutionary Algorithm (CCEA). This algorithm effectively reduces the dimensionality of the problem by decomposing the complex problem into multiple sub-problems, making the search space more manageable and optimized. Moreover, the CCEA leverages parallel processing mechanisms to handle each sub-problem simultaneously, thereby significantly improving computational efficiency. Additionally, its unique co-evolutionary strategy balances local and global searches, not only enhancing the algorithm’s global optimization capability but also improving its ability to escape local optima. Consequently, many researchers have utilized the CCEA to solve path planning problems. For example, Shao et al. [34], building upon the foundation of co-evolution, proposed a Distributed Cooperative Particle Swarm Optimization (DCPSO) algorithm with an elitist preservation strategy. Kesireddy et al. [35] combined CCEAs with Non-Dominated Sorting Genetic Algorithm II (NSGA-II) and NSGA-III to address multi-objective problems for multiple agents. Wu et al. [33] proposed a multi-UAV collaborative path planning algorithm based on co-evolutionary optimization, designing two information sharing strategies to cope with the combined path search among multiple UAVs.
However, the application of the Cooperative Co-Evolutionary Algorithm (CCEA) for solving multi-UAV path planning still faces the following challenges: 1. The multi-UAV path planning problem is constrained by various conditions, and there is a complex interrelationship among the decision variables of each UAV, increasing the difficulty of selecting decision variables. 2. With the increase in the number of UAVs, the scale of the path planning problem grows exponentially, leading to a significant increase in computational complexity. In the practical application of multi-UAV path planning, it is necessary to enhance the computational speed of the algorithms on one hand, and on the other hand, to ensure that the algorithms can find effective flight routes. The CCEA needs to address the challenge of improving computational efficiency while maintaining the performance of the algorithm.
To address the first challenge, we propose an adaptive decision variable selection method guided by optimization objectives. This method fully leverages the collaborative information between UAVs to optimize the selection and execution of paths, enabling the algorithm to converge more rapidly to the optimal or near-optimal solution, thereby reducing computational time and the number of iterations. To tackle the second challenge, we propose a two-phase strategy. Initially, we plan a feasible solution for the path, and subsequently, based on this feasible solution, we refine the planning. This effectively reduces the number of evaluations the algorithm must perform and reduces the computational complexity of the algorithm. By integrating the advantages of the Cooperative Co-Evolutionary Algorithm (CCEA) and the Dubins path planning algorithm, we propose a Cooperative Co-Evolutionary Algorithm for UAV path planning problems, named CCEA-ADVS. Our contributions in this work are as follows:
  • To enhance the efficiency of decision variable selection and to accelerate the convergence of the algorithm, an adaptive decision variable selection strategy is proposed. This strategy first addresses the single-UAV constraint problem, then the multi-UAV constraint problem, and finally, the path planning problem, adaptively selecting decision variables based on the current population state.
  • To reduce the computational complexity and cost of the algorithm, a two-stage strategy is proposed. In the first stage, the algorithm generates a rough path for the UAV from the starting point to the destination. This path may not be optimal, but it provides a basic framework that can be refined in subsequent stages. The goal of the first stage is to quickly generate a feasible path rather than an optimal solution, thereby reducing computational load. The second stage involves Dubins path planning [36,37]. In this stage, the algorithm refines the path generated in the first stage by adjusting points and curves along the path, making the trajectory smoother and more precise, and meeting the actual flight requirements of the UAV. This kind of optimization usually requires more computational resources, but since a basic path is already in place, the optimization can be more targeted, thus reducing unnecessary computations. By adopting this coarse-to-fine approach, an effective balance between computational cost and path quality is achieved. A reasonable path is quickly generated in the first stage, and then, refined in the second stage, reducing the number of evaluations needed in the process of finding the optimal solution, thereby reducing the overall computational complexity.
  • To demonstrate the advantages of the proposed algorithm, we conducted comparative experiments and algorithm analysis. We analyzed the encoding method, adaptive decision variable strategy, evolutionary algorithm, and two-stage strategy in the proposed algorithm, all of which showed good performance. Comparative experiments showed that the CCEA-ADVS algorithm significantly outperformed the GWO, PSO, ABC, and JADE algorithms in terms of path performance, running time, computational efficiency, and convergence speed. Moreover, in large-scale experiments involving 500 UAVs, the algorithm also demonstrated good adaptability, stability, and scalability.
The remainder of this paper is organized as follows: Section 2 designs the cost function for the multi-UAV path planning problem. Section 3 explains the basic principles of the CCEA and Dubins algorithms. Section 4 details the proposed CCEA-ADVS algorithm. Section 5 presents the comparative experiments with the algorithms, as well as an analysis of the algorithmic schemes. Finally, Section 6 and Section 7 are the discussion and conclusions, respectively.

2. Problem Description

Path planning for multiple UAVs is an extremely challenging problem, which not only requires consideration of the UAVs’ mission types, flight areas, and the UAVs’ own performance, but also the impact of obstacles on the path during flight and the issue of mutual collisions between UAVs. In this paper, we aim to develop an algorithm to generate a series of optimal or near-optimal flight paths from the starting point to the destination.
Assume the flight environment of the UAV is defined as E, the number of UAVs is M, and the velocity of each UAV is defined as v uav . In this paper, we conduct planning at both the path point and trajectory levels. To evaluate the scores of the results, we need to discretize the path or trajectory into N points, where each point has the value of ( x i , y i , z i ) in the coordinate system. Then, the path of UAV m can be represented as P m = { P 1 m , P 2 m , , P N m m } , where P i m represents the i-th waypoint of the m-th UAV. Our main objective is to find a set of paths P = { P 1 , P 2 , , P M } that minimizes the total cost of all paths in the set.
To transform this problem into a global optimization problem, thereby increasing the likelihood of finding a global optimal solution or an approximate global optimal solution, we adopt the method of converting all constraints into a cost function. The cost function is configured in accordance with the references [38,39], Among them, path cost, threat cost, altitude cost, smoothing cost, and radar missile cost belong to the single-UAV constraint, and multi-UAV collision cost belongs to the multi-UAV constraint. Six cost constraints are set up as follows:
  • Path cost: The path cost is calculated by summing the Euclidean distances between each waypoint and dividing by the Euclidean distance between the start and end points (minimum distance). The path cost for UAV m is calculated as follows:
    F 1 m = i = 1 N m 1 P i m P i + 1 m P 1 m P N m m
    P i m P i + 1 m = ( x i + 1 m x i m ) 2 + ( y i + 1 m y i m ) 2 + ( z i + 1 m z i m ) 2
    where P i m P i + 1 m represents the vector between two waypoints P i m and P i + 1 m , and ( x i m , y i m , z i m ) are the coordinates of the i-th waypoint for the m-th UAV.
  • Flight-prohibited zone cost: During actual flight, the path of the UAV must not overlap with the location of prohibited flight areas. The cost F 2 m when the m-th UAV passes through obstacles can be expressed as follows:
    F 2 m = i = 1 N m z = 1 Z T z ( P i m P i + 1 m )
    where Z is the total number of no-fly zones, and T z is the threat cost of the z-th no-fly zone to the path segment P i m P i + 1 m . The calculation method is as follows:
    T z ( P i ) = 0 , ( x i m , y i m , z i m ) I F z , otherwise
    where I F z is the z-th no-fly zone. In this paper, two shapes of no-fly zones are considered: rectangular and circular.
  • Altitude cost: In real life, UAVs often need to fly at low altitudes, so the impact of terrain on the safety of UAV flight must also be considered by adding altitude constraints. Let the minimum flight altitude of the UAV be h min , and the maximum flight altitude be h max . The flight altitude of the m-th UAV should be between the two, and the cost associated with the path and altitude F 3 m is calculated as follows:
    F 3 m = 1 N m i = 1 N m H i m
    H i m = h i m , h min h i m h max , otherwise
    h i m = z i m map ( x i m , y i m )
    where h max and h min are the maximum and minimum values of the altitude limit, respectively. The function map ( x i m , y i m ) represents the altitude at coordinates ( x i m , y i m ) on the map.
  • Smoothness cost: Smoothness constraints are an important indicator to evaluate whether the path generated by the UAV is feasible. UAVs have maximum angles when changing direction and climbing, so the limits of the UAV’s turning rate and climb rate must be considered when generating the path. The formulas for calculating the turning angle and climb angle of the m-th UAV are as follows:  
    ϕ i m = arctan P i m P i + 1 m × P i + 1 m P i + 2 m P i m P i + 1 m · P i + 1 m P i + 2 m
    φ i m = arctan z i + 1 m z i m P i m P i + 1 m
    where P i m P i + 1 m and P i + 1 m P i + 2 m represent the projections of the path segments P i m P i + 1 m and P i + 1 m P i + 2 m on the xy plane, respectively.
    Thus, the smoothness cost F 4 m can be expressed as
    F 4 m = i = 1 N m 1   c t u r n i m + i = 1 N m 1   c s l o p e i m
    where c t u r n i m and c s l o p e i m are the cost functions for the turning angle and climb angle, respectively:
    c t u r n i m = 0 , | ϕ i m | < ϕ max , o t h e r w i s e ,   c s l o p e i m = 0 , | φ i m | < φ max , o t h e r w i s e
  • Radar detection and missile kill cost: When planning UAV paths, the potential impact of combat platforms such as radar missiles on the UAV flight path must be considered. In this paper, the radars and missiles in different combat platforms are set to the same coordinates as the combat platform. The formula for calculating the distance between the path point and the l-th combat platform for the m-th UAV is
    d i s i , l m = ( x i m m x l ) 2 + ( y i m m y l ) 2 + ( z i m m z l ) 2
    where ( m x l , m y l , m z l ) is the given position of the l-th combat platform, l = ( 1 , 2 , , C ) , and C is the number of combat platforms. Let the l-th radar have a radius R R D max l of a circle. Then, the cost calculation formula between the UAV and the radar is
    R D i , l m = δ 4 ( d i s i , l m ) 4 d i s i , l m R R D max l 0 o t h e r w i s e
    If the UAV is within the range of the enemy missile, it faces risks. A path with a lower kill risk is safer than one with a higher probability. For each path point, when the UAV is within the maximum risk distance area defined by the l-th enemy missile, it poses a certain kill risk to the UAV, represented as
    R K i , l m = d i s i , l m R R K m a x l 0 otherwise
    The missile kill cost of the m-th UAV for the entire path is calculated as   
    F 5 m = i = 1 N m l = 1 C ( R K i , l m + R D i , l m )
  • Multi-UAV collision cost: Assume that UAVs have a constant speed flight. To ensure safety, each point of the path of the u-th UAV is compared with every point of the v paths of the remaining UAVs. If the distance between two points is smaller than the safety threshold d m i n , the difference in their arrival times ( t i u and t j v ) is checked to see if it is less than the safe time interval t m i n to avoid collision risks. The arrival times t i u and t j v can be obtained by dividing the path by the speed. Thus, the total collision cost can be expressed as
    F 6 = u = 1 M v = u + 1 M i = 1 N u j = 1 N v c o l i d e i j u v
    c o l i d e i j u v = , if d i j u v < d min and | t i u t j v | < t min 0 , otherwise
    d i j u v = ( x i u x j v ) 2 + ( y i u y j v ) 2 + ( z i u z j v ) 2
  • Total cost: The constrained costs for UAVs are transformed into cost functions, quantifying the quality of UAV paths by the total cost of the cost functions. The total cost function is shown as follows:
    F = m = 1 M i = 1 5 α i F i m + α 6 F 6
    where α i represents the weight coefficient of the i-th constraint or cost.

3. Preliminary Knowledge of CCEA-ADVS Algorithm

The heuristic co-evolutionary path planning method based on CCEA-ADVS employs the CCEA algorithm, JADE, and Dubins path planning to offer an efficient solution. By applying the co-evolutionary principle inherent to CCEA, the multi-UAV path planning problem is decomposed into individual sub-problems. JADE’s adaptive parameter adjustment is integrated to optimize and solve each. Using the preliminary paths from CCEA and JADE optimizations, the Dubins path planning method is then applied to generate UAV paths that are compatible with actual flight requirements.

3.1. Cooperative Co-Evolutionary Algorithm (CCEA)

CCEA, an optimization algorithm that employs co-evolutionary strategies [40,41], decomposes complex optimization problems into several sub-problems. These sub-problems are addressed through the co-evolution of multiple subpopulations within an evolutionary algorithm (EA) framework. The individual solutions are then integrated to form a comprehensive solution to the overarching problem.
The CCEA algorithm framework is presented in Algorithm 1. This framework applies local search and fitness sharing among individuals while ensuring adherence to constraint conditions to explore the solution space and escape local optima. In this context, EAStep ( K i , b, N p o p ) represents the evolutionary algorithm function for generating and evaluating the offspring within subpopulations. SelectBest ( K i , O i ) denotes the selection process for the next generation of parent subpopulations.
Algorithm 1 Cooperative Co-Evolutionary Algorithm
Require: The number of variables n, Subpopulation size N p o p .
Ensure: Solution P
 1: Problem decomposition, n-dimensional decision vector is divided into n 1-dimensional subcomponents { x 1 },{ x 2 },…,{ x n }.
 2: Initialize the context vector b { b 1 , b 2 , . . . , b n } . b i is a solution for the i-th subpopulation.
 3: set F b e s t f ( b ) ▹ Evaluate the fitness
 4: for i← 1 to n do 
 5:       Initialize the i-th subpopulation K i K i { K i , 1 , . . . , K i , N p o p }
 6:       for j ← 1 to N p o p  do 
 7:             F K i j f ( b 1 , . . . b i 1 , K i , j , b i + 1 , . . . , b n ) ▹ Evaluate each individual
 8:       end for 
 9: end for 
10: while stopping criterion is not met do 
11:       for i ← 1 to n do 
12:              O i , F O i EAStep ( K i , b, N p o p )▹ generate and evalute offspring sub-populations
13:              j the id of the best individual in O i
14:             if  F O i j < F b e s t  then 
15:                  F b e s t F O i j
16:                  b i O i j
17:             end if 
18:              K i SelectBest ( K i , O i )
19:       end for 
20: end while 
21: function EAStep( K i , b , N p o p )
22:          O i GeneticOperators ( K i ) ▹ Generate offspring sub-population
23:    for j ← 1 to N p o p  do 
24:          F O i j f ( b 1 , . . . b i 1 , O i , j , b i + 1 , . . . , b n )
25:       end for 
26:       return  O i , F O i
27: end function 

3.2. Dubins Algorithm

The Dubins path planning algorithm combines simple linear motions, including straight segments, left turns, and right turns, aiming to generate the shortest and most efficient path [36]. The core idea of this algorithm is to accurately determine the shortest path connecting two points in three-dimensional space while meeting the UAV’s curvature constraints and specific entry and exit directions. In essence, a Dubins path represents the shortest feasible path connecting two points in three-dimensional space with specific directions under a given maximum curvature constraint. This path typically consists of a sequence of CSC (circular arc–straight segment–circular arc) or CCC (three consecutive circular arcs) [37].
The Dubins algorithm’s advantage lies in its rigorous geometric principles, enabling the precise calculation of the shortest path within the physical constraints of flight and avoiding the computational errors and approximate solutions common to other algorithms. It offers a solution for UAVs that meets the physical constraints of flight while achieving the shortest flight distance.

4. Proposed Method

4.1. Algorithm Framework

In response to the issues with decision variable selection in the CCEA algorithm and the computational complexity of multi-UAV path planning, we propose the CCEA-ADVS algorithm. To reduce computational complexity and generate flight trajectories that meet practical requirements, the CCEA-ADVS algorithm employs a two-phase strategy. This approach transitions from coarse to fine, dividing the optimization into two phases. In the first phase, we conduct a rough search, using the JADE algorithm to quickly identify a feasible initial path solution. While this solution is not optimal, it serves as a solid foundation for the subsequent phase. In the second phase, we build on the initial path and conduct a more refined search. Applying the Dubins algorithm to refine the flight paths, we identify the optimal route that satisfies all constraints and closely aligns with the flight trajectories of UAVs in the real world. In order to improve the efficiency of decision variable selection in cooperative algorithms and accelerate the convergence of the algorithm, we introduce an adaptive decision variable selection strategy. The main process is shown in Figure 1.

4.1.1. Individual Encoding Scheme

The choice of encoding is crucial for evolutionary algorithms, as it directly affects the efficiency of the algorithm, the representation of individuals, and the design of operations for the associated evolutionary algorithms.
In this paper, the encoding of an individual is represented by a list of 3D waypoints with a length of 3 N W P . During path planning, two waypoints are directly connected. In trajectory planning, two waypoints are connected using Dubins curves.
In existing evolutionary algorithms, these waypoints are typically represented in 3D coordinates using encoding methods such as Cartesian coordinate systems [13], polar coordinate systems [42], spherical vectors [43], and rotated Cartesian coordinate systems [39]. Among these, the rotated Cartesian coordinate system has shown better results and is simple to implement; therefore, we choose the rotated coordinate encoding method to encode the path solutions. For any waypoint ( x i , y i , z i ) , i = 1 , 2 . . . , N w , in a rotated coordinate system, its codification in the Cartesian coordinate system ( x i , y i , z i ) is mathematically defined as
x i = x 1 + cos ( φ r o t + ϕ r o t ) · x i 2 + y i 2 y i = y 1 + sin ( φ r o t + ϕ r o t ) · x i 2 + y i 2 z i = z i
where φ r o t is the angle included by the direction from the start to the waypoint and x -axis, and ϕ r o t is the angle between the x -axis and x-axis.

4.1.2. Multi-UAV Path Planning Problem Decomposition

When dealing with multi-UAV problems, the originally complex multi-UAV path planning problem is transformed into a series of more easily solvable single-UAV path planning problems, as shown in Figure 2. Assuming there are M UAVs, this results in M × 3 N W P decision variables. If the decision variable b is divided into M groups according to 3 N W P , then b = { b 1 , b 2 , . . . , b M } , where each b i contains 3 × N W P decision variables { b ( i 1 ) × 3 N W P + 1 , . . . , b i × 3 N W P } . In this way, the path planning problem for each UAV can be carried out independently while ensuring that collision avoidance requirements between all UAVs are met within the overall mission. Although each UAV’s path planning problem only considers collision avoidance constraints directly related to that UAV, other global constraints, such as flight regulations and environmental limitations, still apply to the planning process of each UAV and remain unchanged. Let a single UAV be denoted as m. The collision cost for UAV m is as follows:
F 6 d m = v m i = 1 N m j = 1 N v c o l i d e i j
Then, the calculation of the sub-problem cost of UAV m becomes
F m = c = 1 5 α c F c m + α 6 F 6 d m

4.1.3. CCEA-ADVS Algorithm

Algorithm 2 provides a detailed description of the specific steps of the CCEA-ADVS algorithm. It begins with problem decomposition and the initialization of the environmental vector b and the initial subpopulation (lines 1–6). It then proceeds with the decision variable selection strategy (line 8) and a two-phase optimization strategy. The decision variable selection is represented by the function DecisionVariableSelect ( M , K , ω ) , detailed in Algorithm 3. Here, Ord denotes the order set of decision variable selection, and I represents the set of decision variables that require updating. Subsequently, the decision variables (UAV paths) are updated according to the selection order. The first phase optimizes the flight path through the JADE algorithm (line 11), and the second phase refines the flight path through the Dubins algorithm (line 13), followed by an update of the environmental variables b and the selection of the next parent subpopulation (lines 15–22). Finally, the algorithm outputs a path planning solution that effectively selects decision variables during the planning process, improves the computational efficiency of the algorithm, and ensures the efficient and safe operation of the UAVs.
Algorithm 2 CCEA-ADVS algorithm
Require: Number of UAVs M, the number of waypoints N W P , Subpopulation size N p o p .
Ensure: Solution P
 1: Problem decomposition, the decision vector is divided into M  3 N W P -dimensional subcomponents { x 1 ,…, x 3 N W P },…,{ x ( M 1 ) 3 N W P + 1 ,…, x M 3 N W P }.
 2: Initialize the context vector b { b 1 , b 2 , . . . , b M } . b i is a solution for the i-th subpopulation.
 3: Assign strategy weights ω { 1 , 1 , 1 } .
 4: set F b e s t f ( b ) ▹ Evaluate the fitness
 5: for i← 1 to M do
 6:        Initialize the i-th subpopulation K i
 7:        for j ← 1 to N p o p  do
 8:             F K i j f ( b 1 , . . . b i 1 , K i , j , b i + 1 , . . . , b n )
 9:        end for
10: end for
11: for iter = 1 to N i t e r 1 + N i t e r 2  do
12:        F o l d F b e s t
13:        O r d , ω , S t r a t e g y DecisionVariableSelect ( M , K , ω ) ▹ Algorithm 3
14:       for i O r d 1 to O r d e n d  do
15:            if iter  N i t e r 1  then▹ Maximum Number of Phase 1 Iterations
16:                O i , F O i J A D E ( K i , b , N p o p ) ▹ Phase 1, Algorithm 4
17:            else
18:                O i , F O i JADEWithDubins ( K i , b , N p o p ) ▹ Phase 2
19:            end if
20:             j the id of the best individual in O i
21:            if  F O i j < F b e s t  then
22:                F b e s t F O i j , b i O i j ▹ Update context vector
23:            end if
24:             K i SelectBest ( K i , O i )
25:       end for
26:       if  F b e s t < F o l d  then
27:             ω s t r a t e g y min ( ( 1 ρ w ) ω s t r a t e g y + ρ w 0.5 , 10 ) ▹ Update weights
28:       else
29:             ω s t r a t e g y max ( ( 1 ρ w ) ω s t r a t e g y ρ w 0.1 , 0.1 ) ▹ Update weights
30:       end if
31: end for
32: P g
33: return P

4.2. Adaptive Decision Variable Selection Strategy

The selection of decision variables is a crucial step in multi-UAV path planning. Multi-UAV path planning in complex environments is inherently complex, further complicated by the need for UAV collaboration, requiring decision variables to balance individual UAV path optimization with overall collaborative objectives. In cooperative co-evolutionary path planning, most algorithms use a round robin or random method to select decision variables, which do not fully utilize the structural information of the problem and may lack specificity and optimization when selecting variables. To optimize the target-guided decision variable selection method, we proposes an adaptive decision variable selection strategy. The strategy first addresses single-UAV constraints, ensuring that the flight path of an individual UAV meets all basic flight constraints, such as avoiding obstacles, and complying with speed and altitude limits. Then, it addresses multi-UAV constraints, considering mutual constraints between multiple UAVs on the basis of single-UAV constraints, such as avoiding mutual collisions and coordinating flight paths. Finally, on the premise of meeting both single- and multi-UAV constraints, the path is further optimized to complete the path planning. The decision variable selection strategy is shown in Algorithm 3, with the following steps.
Algorithm 3 Adaptive Decision Variables Selection Strategy
Require: Number of UAVs M, context vector b, strategy weights ω
Ensure: Decision variable set O r d , Strategy weights ω , Selected Strategy S t r a t e g y
 1: if constraints for single UAV are not met for b then 
 2:      O r d UAVs ids that do not meet the single-UAV constraints in b.
 3: else 
 4:     if exist Multi-UAV collision in b then 
 5:         Randomly assign strategies S t r a t e g y according to weights ω .
 6:     else 
 7:          S t r a t e g y 3
 8:     end if
 9:      n s e l 8
10:      c o l ( 1 , . . . , M ) Count the number of collisions for each UAV
11:     if  S t r a t e g y =1 then 
12:          O r d GreedyTopK ( c o l , n s e l )
13:     else if  S t r a t e g y =2 then 
14:         Initializes decision variable set O r d
15:          P r o b ( 1 , . . . , M ) VariableProbabilities ( c o l )
16:         Randomly select n s e l different UAVs based on the P r o b .
17:     else if  S t r a t e g y = 3 then 
18:          O r d Randperm ( M , n s e l )
19:     end if
20: end if
(1) Solving the constraints of a single UAV (lines 1–3): Each UAV or sub-problem can independently initiate with a solution, thereby taking turns. This solution ensures that each UAV has a feasible solution that is not infinite (a solution is considered feasible when the fitness costs of F 2 m to F 5 m are not infinite), and it also guarantees that each UAV or sub-problem has the opportunity to engage in the optimization process, establishing an initial reference point for subsequent cooperative evolution.
(2) Solving the constraints of multiple UAVs (lines 4–18): During the decision variable selection phase, autonomously select the appropriate decision variable selection strategy based on the operator weight ω . The following are the strategies for decision variable selection:
  • Strategy 1 (greedy strategy): Optimizes the UAV with the most conflicts. The function GreedyTopK(K, n s e l ) represents the greedy function, which returns the top n s e l UAVs with the most conflicts.
  • Strategy 2 (roulette wheel selection strategy): The VariableProbabilities(K) function assigns selection probabilities based on the number of UAV conflicts, and the roulette wheel is used to select n s e l sub-problems according to the selection probabilities. This strategy combines randomness and conflict information, which helps to discover new optimization opportunities.
  • Strategy 3 (random selection strategy): Represented by the Randperm(M, n s e l ) function, it randomly selects n s e l sub-problems for decision variable selection. This strategy is simple and easy to implement but may lack specificity, leading to slow convergence and wasting a lot of computational resources.
(3) Reducing the cost of the optimization objective (lines 7 and 17): At this point, a feasible solution has been found, and variables can be selected in a completely random manner.

4.3. Adaptive Differential Evolution with Optional External Archive (JADE)

The JADE algorithm is an improved version of the Differential Evolution (DE) algorithm, automatically adjusting control parameters to enhance robustness. It uses historical data to inform evolutionary direction, effectively maintaining population diversity and improving convergence performance.
The JADE algorithm follows the basic steps of differential evolution, including initialization, mutation, crossover, selection, and adaptive parameter control, as shown in Algorithm 4. During mutation, new candidate solutions are generated using the ‘DE/current-to-pbest’ strategy. During crossover, individuals are produced using the crossover probability CR. During selection, the fittest individuals are selected based on fitness. Here, C R i Normrnd( μ C R , 0.1 ) indicates that C R i follows a normal distribution with mean μ C R and standard deviation 0.1, and F i Cauchyrnd( μ F , 0.1 ) indicates that F i follows a Cauchy distribution with location parameter μ F and scale parameter 0.1. Randint(1,3) randomly selects an integer between 1 and 3.
Additionally, within the CCEA framework, the JADE algorithm does not require convergence in a single call, thus a smaller number of iterations, N E A , is employed. In this study, we set N E A to 20.
Algorithm 4 JADE algorithm framework
Require: Subpopulation K i , context vector b, Subpopulation size N p o p .
Ensure: Solution O i
 1: set parameters μ C R , μ F , c , A , Iterations N E A
 2: x a l l i , 1 K i x a l l i , 1 represents current population of the first iteration
 3: for g← 1 to N E A  do 
 4:       S F ; S C R
 5:      for j← 1 to N p o p  do
 6:           C R i Normrnd ( μ C R , 0.1 ) , F i Cauchyrnd ( μ F , 0.1 ) , J r a n d randint ( 1 , 3 )
 7:          Randomly choose x b e s t i , g as one of the 100 ρ % best vectors.
 8:          Randomly choose x r 1 i , g x j i , g from x a l l i , g .
 9:          Randomly choose x r 2 i , g x j i , g , x r 2 i , g x r 1 i , g from x a l l i , g A .
10:           v j i , g x j i , g + F i · ( x b e s t i , g x j i , g ) + F i · ( x r 1 i , g x r 2 i , g )
11:          for k ← 1 to 3 do 
12:               if  k = J r a n d or rand(0,1)< C R i  then
13:                    u j , k i , g v j , k i , g
14:               else
15:                    u j , k i , g x j , k i , g
16:               end if
17:          end for
18:          if  F x j i , g F u j i , g  then
19:               x j i , g + 1 x j i , g
20:          else
21:               x j i , g + 1 u j i , g , S C R C R i , S F F i , A x j i , g
22:          end if
23:     end for
24:     Randomly remove solutions from A so that | A | N p o p
25:      μ C R ( 1 c ) · μ C R + c · m e a n ( S C R ) , μ F ( 1 c ) · μ F + c · K S F K 2 K S F K
26: end for
27: O i SelectBest ( x a l l i , g )

5. Experiment

5.1. Experiment Setup

This experiment is based on the simulation of real-world scenarios. The simulation experiments were conducted on a Lenovo laptop equipped with an Intel 2.5 GHz processor and 8 GB of RAM, according to the algorithm design. The experimental scenarios include four three-dimensional terrain maps generated from digital elevation model (DEM) data. We assume known position information for obstacles and combat platforms, which, with known center coordinates, radius, and height, are represented on the 3D maps. To thoroughly test the algorithm’s performance, the fourth scenario included 20 combat platforms and 7 no-fly zones, while the others included 5 combat platforms and 7 no-fly zones each. The four experimental scenarios are shown in Figure 3. We set the UAV speed, v u a v , to 100 km/h. The Dubins path’s turning radius is 3.8 km. The safety distance S 1 between the UAV and the edge of the threat is set to 15 km. The safety distance S 2 between UAVs is set to 10 km. The no-fly zones’ length and width are both set to 100 km. The radar radius is set to 50 km, and the missile radius is set to 15 km.
In the absence of specific designations, the parameters related to the algorithms in this paper are shown in Table 2, with other parameters set according to Ref. [44].

5.2. Comparative Experiment with Existing Algorithms

5.2.1. Runtime and Quality

The comparative experiments of the algorithm were conducted with four algorithms—PSO [13,14,15], JADE [44,45], ABC [18], and GWO [21,22]—each compared across four scenarios, and all four methods used the priority planning method [46] to address UAV collision issues. The experimental setup in this section involved 40 UAVs, with each algorithm running 20 times and the results taken as the average. Table 3 lists the parameter settings for the comparative algorithms, which were configured according to the aforementioned references.
Figure 4 illustrates the convergence curves of the five algorithms in four scenarios. It can be observed that the CCEA-ADVS algorithm demonstrates a faster convergence rate and the lowest cost consumption compared to the other four algorithms. Among them, the ABC algorithm performs the worst, with the slowest convergence speed and the poorest convergence precision. The ABC algorithm is greatly influenced by parameters and has high requirements for the continuity of the search space. When facing complex high-dimensional optimization problems, it takes a longer time to find better solutions and is prone to falling into local optima, with relatively weak global optimization capabilities. Although the GWO and PSO algorithms also have strong global search capabilities, their search processes rely more on guidance by the global optimum. In complex or high-dimensional search spaces, as the scale of UAVs and environmental complexity increase, this reliance leads to the algorithms being more likely to fall into local optima in later stages. The JADE algorithm, compared to GWO and PSO, has a slower convergence speed but results in lower cost consumption. The JADE algorithm, through its adaptive parameter control and “DE/current-to-pbest” mutation strategy, ensures the reliability of the algorithm while obtaining lower-cost paths, but this also leads to a relatively slower convergence rate of the algorithm.
Figure 5 shows a boxplot of the optimal costs of the five algorithms over 20 independent iterations in four scenarios, with statistical results presented in Table 4. It can be seen that the CCEA-ADVS algorithm performs very well in the four scenarios, obtaining the lowest path cost and consuming the least time compared to the other four algorithms. The smaller variance coefficient of the CCEA-ADVS algorithm proves its superiority in terms of stability. Less running time and lower path cost indicate that the algorithm’s planned path has higher performance in complex flight environments. In conjunction with the experiments in Section 5.3, the adaptive decision variable selection strategy adopted by the CCEA-ADVS algorithm improves the efficiency and cost effectiveness of path optimization, enabling the algorithm to obtain the lowest-cost flight path. The two-stage strategy greatly reduces the algorithm’s running time and computational complexity. In addition, the rotated coordinate encoding method adopted by the algorithm also has certain optimization effects on cost and running time. In summary, the CCEA-ADVS algorithm effectively enhances the computational efficiency and decision variable selection efficiency of the CCEA algorithm. The algorithm can more rapidly determine a flight path that is less costly and more secure and efficient.

5.2.2. Scalability

When the number of UAVs reaches the level of hundreds, the comparative algorithms in the experiment (PSO, GWO, ABC, JADE) can no longer meet the path planning requirements. In order to test whether the CCEA-ADVS algorithm meets the requirements for large-scale UAV path planning, we compared the two situations: one without using the Dubins algorithm and the other with the Dubins algorithm. The first group of parameters was set to N i t e r 1 = 40 , N i t e r 2 = 0 ; and the second group of parameters was set to N i t e r 1 = 20 , N i t e r 2 = 20 , with all other parameters set according to Table 2.
The scalability and cost experimental results are depicted in Figure 6, with corresponding data detailed in Table 5. It can be seen from the figure that using Dubins planning (i.e., employing a two-phase strategy) yields a lower cost, and as the number of UAVs increases, the gap in cost widens. The time results for the algorithm are depicted in Figure 7, indicating that not using Dubins planning (i.e., conducting only the first-phase strategy) has a shorter running time, and thus, can obtain the UAV flight path more quickly compared to using only Dubins or employing a two-phase strategy. In summary, the CCEA-ADVS algorithm can meet the needs of large-scale UAV path planning and has excellent scalability. Furthermore, the scalability experiments demonstrate that the two-phase strategy combines the advantages of fast operation from the first phase and low path cost from the second phase, enabling the acquisition of an optimal path in a shorter amount of time.

5.2.3. Time Complexity Analysis

This section analyzes the time complexity of the proposed CCEA-ADVS algorithm. Assuming the complexity of calculating fitness is O ( 1 ) .
Firstly, the time complexity of the CCEA based on the JADE evolutionary algorithm is discussed. The basic steps of the algorithm are shown in Algorithm 1. The time complexity of each iteration of the algorithm = time complexity of initializing the population + time complexity of the co-evolutionary process. Time complexity for initializing the population: The time complexity for initializing the population of each sub-problem is O ( M N W P N p o p ) . Time complexity of co-evolutionary process: O ( M N W P N p o p ) . Thus, the total time complexity is O ( M N W P N p o p ) + O ( N i t e r M N W P N p o p ) . Since N i t e r is much greater than 1, we can simplify this to O ( N i t e r M N W P N p o p ) .
In the CCEA-ADVS algorithm, the time complexity of the adaptive decision variable strategy is as shown in Algorithm 3, time complexity = solving the single-UAV constraint + solving the multi-UAV constraints + reducing the cost of the optimization objective time complexity. The time complexity is then O ( M N W P N p o p ) + O ( n s e l N W P N p o p ) + O ( n s e l N W P N p o p ) . Since in the CCEA-ADVS algorithm, the execution times for solving the single-UAV constraints are much less than the latter two, we can simplify this to O ( n s e l N W P N p o p ) . The two-stage strategy is as shown in Algorithm 2. The time complexity of the CCEA-ADVS algorithm is O ( M N W P N p o p ) + O ( N i t e r 1 n s e l N W P N p o p ) + O ( N i t e r 2 n s e l N W P N p o p ) . Since N i t e r 1 n s e l is much greater than M, it can be simplified to O ( N i t e r n s e l N W P N p o p ) . Because N i t e r n s e l < N i t e r M , the proposed CCEA-ADVS algorithm enhances the convergence speed, making it more efficient.

5.3. Algorithm Analysis

5.3.1. Encoding Method Analysis

To validate that the rotated Cartesian vector (ROC) [39] encoding method is more suitable for the CCEA-ADVS algorithm to perform multi-UAV path planning, we compared it with the Cartesian (CAR) [13], relative polar (REP) [42], and relative spherical vector (RES) [43] encoding methods in four scenarios with 10 UAVs each.
The experimental results for the four coding methods are presented in Table 6. In terms of algorithmic cost, the rotated coordinate coding method consistently outperforms the other three, achieving the lowest cost and the shortest time consumption across all four scenarios. Among the methods evaluated, the Cartesian encoding approach exhibits the poorest performance. With an increasing number of UAVs, it incurs the highest costs and time consumption.
In summary, the rotated coordinate coding method is superior to the other three, effectively reducing the errors associated with coordinate transformation and more naturally addressing the relative motion among multiple UAVs. By employing the rotated coordinate coding, a superior path can be obtained at a reduced cost and with minimized time expenditure.

5.3.2. Adaptive Decision Variable Strategy Analysis

To verify the superior performance of the proposed adaptive decision variable strategy, an analysis of the decision variable strategy is presented below. The strategy is compared with the greedy strategy (selecting only the UAV with the most conflicts with other UAVs) and the random selection strategy (randomly choosing UAVs in conflict). The comparative experiments were conducted in four scenarios, with each scenario run 20 times, and the results are taken as the average values.
The experimental results, as detailed in Table 7, indicate that the adaptive decision variable strategy results in the lowest cost compared to the greedy and random strategies. The random strategy exhibits the least favorable performance, characterized by the longest execution times and the highest cost expenditure. In comparison to the greedy strategy, the adaptive decision variable strategy more effectively evades local optima, leading to reduced cost consumption and expedited running times. The experiments confirm that the adaptive decision variable strategy enhances the efficiency and cost-effectiveness of path optimization, enabling the generation of more rapid and superior flight paths for UAVs.

5.3.3. Evolutionary Algorithm Analysis

To verify the effectiveness of the JADE evolutionary algorithm used in the path planning for multiple UAVs, this section conducts a comparative analysis of the performance of the JADE evolutionary algorithm. We compared JADE with PSO, GA, and DE using 10 UAVs across four scenarios, with each set of experiments conducted 20 times and the average values taken. The algorithm parameter settings can be found in Table 2.
The relevant experimental data of the evolutionary algorithms in the four scenarios are shown in Table 8. It can be observed that the JADE algorithm outperforms the other three algorithms, capable of obtaining flight paths with lower costs. Additionally, the time consumed by the algorithm is less than that of the PSO and GA algorithms, and slightly more than that of the DE algorithm. Therefore, selecting the JADE algorithm can yield a more optimal path that takes into account both cost and time.

5.3.4. Two-Phase Strategy Analysis

To validate that the two-stage strategy can effectively reduce the computational complexity of the algorithm, this section conducts an experimental analysis of the strategy, exploring its performance. We deployed 10 UAVs across four scenarios, dividing the experiments into three groups: one conducting only the first stage of the algorithm (JADE), another conducting only the second stage (JADEWithDubins), and the third using the CCEA-ADVS algorithm. Each group was tested 20 times, and the results were averaged. The parameter settings for each group were as follows: the first group was set to N i t e r 1 = 100 , N i t e r 2 = 0 . The second group was set to N i t e r 1 = 0 , N i t e r 2 = 100 . The third group was set to N i t e r 1 = 50 , N i t e r 2 = 50 . All other parameters were set according to Table 2.
The experimental results are presented in Table 9. The data indicate that the two-stage strategy employed by our algorithm exhibits superior performance across four different scenarios, incurring the lowest cost consumption. In terms of time efficiency, the two-stage strategy requires significantly less time compared to using only the Dubins path planning strategy, yet is slightly higher than the time taken by the strategy that executes only phase 1. Table 10 presents a comparative analysis of the rough path for phase 1 and the optimized Dubins path for phase 2 within the third group of UAVs, incorporating the scalability experiments detailed in Section 5.2.2. The two-stage strategy combines the fast running speed of phase 1 with the low path cost of phase 2, significantly reducing both running time and computational complexity. In conclusion, the two-stage strategy can achieve superior performance and expedited flight paths.

5.3.5. Visual Display of Planned Paths

In this paper, 10 UAVs were randomly selected for visual presentation. The algorithm’s planning results in a two-dimensional view are shown in Figure 8.
It can be observed that by considering the airspace environment, flight conditions, and the UAVs’ own performance parameters, the algorithm can effectively reduce the risk of collisions with obstacles, radar and missile threats, and other UAVs, planning a flight path that is both safe and efficient. Furthermore, some UAVs cross the radar area due to the influence of their turning radii in Figure 8a,b,d. The algorithm significantly enhances both the safety and success rate of UAV missions, improves flight efficiency, and reduces unnecessary energy consumption and time costs. In summary, the algorithm demonstrates outstanding performance and practicality in multi-UAV path planning and is capable of conducting path planning in complex flight environments.

6. Discussion

This study provides innovative ideas for further improving the performance of the algorithm through an in-depth discussion of the decision variable subsets among CCEA collaborators, as well as a two-stage hierarchical optimization method. However, the obstacles and threats involved in the currently proposed methods are regular and static, which has certain limitations in practical applications. To more closely simulate real flight conditions and real-world flight environments, future work needs to be expanded in multiple areas. One area is the issue of obstacles and threats such as radar missiles with regular shapes. In reality, the shapes of obstacles are diverse, and threats like radar missiles are not just simple circles; everything still needs to be considered from various aspects based on actual conditions. Another area that needs to be explored is the algorithm’s adaptability to dynamic environments. In practical applications, the environment in which UAVs perform tasks is often dynamic, including sudden obstacles, weather changes, and dynamic threats such as weaponry. Static obstacles and threats often cannot meet actual needs. Enhancing the algorithm’s perception and adaptability to dynamic environments, enabling UAVs to adjust their path planning in a timely manner when the environment changes, is one of the key points of future research work.
In summary, although this study has made certain progress in the application of the CCEA algorithm, further exploration is still needed in the modeling of obstacles and threats in complex environments, as well as in the adaptability and response capabilities to dynamic environments. These research directions are crucial for promoting the comprehensive development of UAV technology, expanding its application prospects, and enhancing overall efficiency, and will lay a solid theoretical and technical foundation for the future application of the UAV field.

7. Conclusions

Addressing the large-scale three-dimensional UAV path planning problem, this paper proposes a CCEA-ADVS-based Cooperative Co-Evolutionary Algorithm. In the proposed method, the multi-UAV path planning problem is transformed into single-UAV path planning problems, and a cost function for individual UAVs is designed. Secondly, to enhance computational efficiency and reduce the scale of problem solving, a two-stage planning strategy is devised. This strategy progresses from coarse to fine, reducing the number of algorithmic evaluations in each phase by focusing on different optimization priorities, thereby increasing the operational speed and computational efficiency. For the conflict optimization of multiple UAVs, an adaptive decision variable selection strategy is proposed, which effectively combines environmental characteristics and UAV path planning features to improve path optimization efficiency and cost effectiveness. Finally, comparative experimental results show that the CCEA-ADVS algorithm outperforms other algorithms in both convergence accuracy and speed. In larger-scale experiments involving 500 UAVs, the proposed algorithm also demonstrates excellent processing capabilities and good scalability. The analysis of the algorithm indicates that the encoding method, evolutionary algorithm, and the adaptive decision variable strategy and two-stage strategy employed all offer better performance and higher efficiency. The CCEA-ADVS algorithm can effectively reduce the computational complexity of the algorithm, improve running speed, and plan safe, efficient, and low-cost flight paths.

Author Contributions

Conceptualization, Q.M., T.Y. and Q.Q.; methodology, Q.M., T.Y. and K.C.; software, Q.M., T.Y. and K.C.; formal analysis, Q.M. and Q.Q.; investigation, Q.M.; resources, Q.M. and Q.Q.; data curation, Q.M.; writing—original draft preparation, Q.M.; writing—review and editing, Q.M.; visualization, Q.M.; supervision, T.Y.; project administration, T.Y.; funding acquisition, Q.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Hunan Provincial Department of Education Scientific Research Outstanding (22B0222).

Data Availability Statement

The data used to support this study have not been made available.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Qu, Q.; Meng, Q.; Chen, K.; Zhang, L. Research and Implementation of Ship/Unmanned Vehicle Collaborative Warfare Simulation System. In Proceedings of the 2023 IEEE 3rd International Conference on Data Science and Computer Application (ICDSCA), Dalian, China, 27–29 October 2023; pp. 275–279. [Google Scholar]
  2. Bellingham, J.S. Coordination and Control of UAV Fleets Using Mixed-Integer Linear Programming. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2002. [Google Scholar]
  3. Kuwata, Y.; How, J. Three dimensional receding horizon control for UAVs. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Providence, RI, USA, 16–19 August 2004; p. 5144. [Google Scholar]
  4. Richards, A.; How, J.P. Aircraft trajectory planning with collision avoidance using mixed integer linear programming. In Proceedings of the 2002 American Control Conference (IEEE Cat. No. CH37301), Anchorage, AK, USA, 8–10 May 2002; Volume 3, pp. 1936–1941. [Google Scholar]
  5. Melchior, P.; Orsoni, B.; Lavialle, O.; Poty, A.; Oustaloup, A. Consideration of obstacle danger level in path planning using A* and fast-marching optimisation: Comparative study. Signal Process. 2003, 83, 2387–2396. [Google Scholar] [CrossRef]
  6. Szczerba, R.J.; Galkowski, P.; Glicktein, I.S.; Ternullo, N. Robust algorithm for real-time route planning. IEEE Trans. Aerosp. Electron. Syst. 2000, 36, 869–878. [Google Scholar] [CrossRef]
  7. Wang, W.; Deng, H.; Wu, X. Path planning of loaded pin-jointed bar mechanisms using Rapidly-exploring Random Tree method. Comput. Struct. 2018, 209, 65–73. [Google Scholar] [CrossRef]
  8. Raghunatan, A.; GV, S.D.; Biegler, L.; Samad, T. Dynamic Optimization Strategies for 3D Conflict Resolution of Multiple Aircraft. AIAA J. Guid. Control. Dyn. 2003, 27, 586–594. [Google Scholar] [CrossRef]
  9. Zhang, X.; Duan, H. An improved constrained differential evolution algorithm for unmanned aerial vehicle global route planning. Appl. Soft Comput. 2015, 26, 270–284. [Google Scholar] [CrossRef]
  10. Chakraborty, J.; Konar, A.; Jain, L.C.; Chakraborty, U.K. Cooperative multi-robot path planning using differential evolution. J. Intell. Fuzzy Syst. 2009, 20, 13–27. [Google Scholar] [CrossRef]
  11. Chen, Y.; Chen, M.; Chen, Z.; Cheng, L.; Yang, Y.; Li, H. Delivery path planning of heterogeneous robot system under road network constraints. Comput. Electr. Eng. 2021, 92, 107197. [Google Scholar] [CrossRef]
  12. Wu, H.; Li, H.; Xiao, R.; Liu, J. Modeling and simulation of dynamic ant colony’s labor division for task allocation of UAV swarm. Phys. A Stat. Mech. Its Appl. 2018, 491, 127–141. [Google Scholar] [CrossRef]
  13. Shao, S.; Peng, Y.; He, C.; Du, Y. Efficient path planning for UAV formation via comprehensively improved particle swarm optimization. ISA Trans. 2020, 97, 415–430. [Google Scholar] [CrossRef]
  14. Zhang, Y.; Gong, D.W.; Zhang, J.H. Robot path planning in uncertain environment using multi-objective particle swarm optimization. Neurocomputing 2013, 103, 172–185. [Google Scholar] [CrossRef]
  15. Qiaorong, Z.; Guochang, G. Path planning based on improved binary particle swarm optimization algorithm. In Proceedings of the 2008 IEEE Conference on Robotics, Automation and Mechatronics, Chengdu, China, 21–24 September 2008; pp. 462–466. [Google Scholar]
  16. Ye, F.; Chen, J.; Tian, Y.; Jiang, T. Cooperative multiple task assignment of heterogeneous UAVs using a modified genetic algorithm with multi-type-gene chromosome encoding strategy. J. Intell. Robot. Syst. 2020, 100, 615–627. [Google Scholar] [CrossRef]
  17. Shorakaei, H.; Vahdani, M.; Imani, B.; Gholami, A. Optimal cooperative path planning of unmanned aerial vehicles by a parallel genetic algorithm. Robotica 2016, 34, 823–836. [Google Scholar] [CrossRef]
  18. Xue, Y.; Jiang, J.; Zhao, B.; Ma, T. A self-adaptive artificial bee colony algorithm based on global best for global optimization. Soft Comput. 2018, 22, 2935–2952. [Google Scholar] [CrossRef]
  19. Goel, U.; Varshney, S.; Jain, A.; Maheshwari, S.; Shukla, A. Three dimensional path planning for UAVs in dynamic environment using glow-worm swarm optimization. Procedia Comput. Sci. 2018, 133, 230–239. [Google Scholar] [CrossRef]
  20. Aljarah, I.; Ludwig, S.A. A new clustering approach based on glowworm swarm optimization. In Proceedings of the 2013 IEEE Congress on Evolutionary Computation, Cancun, Mexico, 20–23 June 2013; pp. 2642–2649. [Google Scholar]
  21. Dewangan, R.K.; Saxena, P. Three-dimensional route planning for multiple unmanned aerial vehicles using Salp Swarm Algorithm. J. Exp. Theor. Artif. Intell. 2023, 35, 1059–1078. [Google Scholar] [CrossRef]
  22. Meidani, K.; Hemmasian, A.; Mirjalili, S.; Farimani, A.B. Adaptive grey wolf optimizer. Neural Comput. Appl. 2022, 34, 7711–7731. [Google Scholar] [CrossRef]
  23. Fang, Y.; Yao, Y.; Zhu, F.; Chen, K. Piecewise-potential-field-based path planning method for fixed-wing UAV formation. Sci. Rep. 2023, 13, 2234. [Google Scholar] [CrossRef]
  24. Majumder, A.; Majumder, A.; Bhaumik, R. Teaching–learning-based optimization algorithm for path planning and task allocation in multi-robot plant inspection system. Arab. J. Sci. Eng. 2021, 46, 8999–9021. [Google Scholar] [CrossRef]
  25. Fu, Y.; Ding, M.; Zhou, C.; Hu, H. Route Planning for Unmanned Aerial Vehicle (UAV) on the Sea Using Hybrid Differential Evolution and Quantum-Behaved Particle Swarm Optimization. IEEE Trans. Syst. Man Cybern. Syst. 2013, 43, 1451–1465. [Google Scholar] [CrossRef]
  26. Yu, C.; Wang, Z. UAV path planning using GSO-DE algorithm. In Proceedings of the 2013 IEEE International Conference of IEEE Region 10 (TENCON 2013), Xi’an, China, 22–25 October 2013; pp. 1–4. [Google Scholar] [CrossRef]
  27. Adhikari, D.; Kim, E.; Reza, H. A fuzzy adaptive differential evolution for multi-objective 3D UAV path optimization. In Proceedings of the 2017 IEEE Congress on Evolutionary Computation (CEC), Donostia, Spain, 5–8 June 2017; pp. 2258–2265. [Google Scholar] [CrossRef]
  28. Ali, Z.A.; Zhangang, H.; Zhengru, D. Path planning of multiple UAVs using MMACO and DE algorithm in dynamic environment. Meas. Control 2023, 56, 459–469. [Google Scholar] [CrossRef]
  29. Elkazzaz, F.S.; Abozied, M.A.; Hu, C. Hybrid RRT/DE Algorithm for High Performance UCAV Path Planning. In Proceedings of the 2017 VI International Conference on Network, Communication and Computing, Kunming, China, 8–10 December 2017; pp. 235–242. [Google Scholar]
  30. Kiani, F.; Seyyedabbasi, A.; Aliyev, R.; Gulle, M.U.; Basyildiz, H.; Shah, M.A. Adapted-RRT: Novel hybrid method to solve three-dimensional path planning problem using sampling and metaheuristic-based algorithms. Neural Comput. Appl. 2021, 33, 15569–15599. [Google Scholar] [CrossRef]
  31. Yang, F.; Fang, X.; Gao, F.; Zhou, X.; Li, H.; Jin, H.; Song, Y. Obstacle avoidance path planning for UAV based on improved RRT algorithm. Discret. Dyn. Nat. Soc. 2022, 2022, 4544499. [Google Scholar] [CrossRef]
  32. Aslan, M.F.; Durdu, A.; Sabanci, K. Goal distance-based UAV path planning approach, path optimization and learning-based path estimation: GDRRT*, PSO-GDRRT* and BiLSTM-PSO-GDRRT. Appl. Soft Comput. 2023, 137, 110156. [Google Scholar] [CrossRef]
  33. Wu, Y.; Nie, M.; Ma, X.; Guo, Y.; Liu, X. Co-Evolutionary Algorithm-Based Multi-Unmanned Aerial Vehicle Cooperative Path Planning. Drones 2023, 7, 606. [Google Scholar] [CrossRef]
  34. Shao, Z.; Yan, F.; Zhou, Z.; Zhu, X. Path Planning for Multi-UAV Formation Rendezvous Based on Distributed Cooperative Particle Swarm Optimization. Appl. Sci. 2019, 9, 2621. [Google Scholar] [CrossRef]
  35. Kesireddy, A.; Shan, W.; Xu, H. Global optimal path planning for multi-agent flocking: A multi-objective optimization approach with NSGA-III. In Proceedings of the 2019 IEEE Symposium Series on Computational Intelligence (SSCI), Xiamen, China, 6–9 December 2019; pp. 64–71. [Google Scholar]
  36. Dubins, L. On Curves of Minimal Length with a Constraint on Average Curvature, and with Prescribed Initial and Terminal Positions and Tangents. Am. J. Math. 1957, 79, 497–516. [Google Scholar] [CrossRef]
  37. Song, X.; Hu, S. 2D path planning with dubins-path-based A* algorithm for a fixed-wing UAV. In Proceedings of the 2017 3rd IEEE International Conference on Control Science And Systems Engineering (ICCSSE), Beijing, China, 17–19 August 2017; pp. 69–73. [Google Scholar]
  38. Besada-Portas, E.; de la Torre, L.; de la Cruz, J.M.; de Andrés-Toro, B. Evolutionary trajectory planner for multiple UAVs in realistic scenarios. IEEE Trans. Robot. 2010, 26, 619–634. [Google Scholar] [CrossRef]
  39. Yang, P.; Tang, K.; Lozano, J.A.; Cao, X. Path planning for single unmanned aerial vehicle by separately evolving waypoints. IEEE Trans. Robot. 2015, 31, 1130–1146. [Google Scholar] [CrossRef]
  40. Potter, M.A.; De Jong, K.A. A cooperative coevolutionary approach to function optimization. In International Conference on Parallel Problem Solving from Nature; Springer: Berlin/Heidelberg, Germany, 1994; pp. 249–257. [Google Scholar]
  41. Ma, X.; Li, X.; Zhang, Q.; Tang, K.; Liang, Z.; Xie, W.; Zhu, Z. A survey on cooperative co-evolutionary algorithms. IEEE Trans. Evol. Comput. 2018, 23, 421–441. [Google Scholar] [CrossRef]
  42. Yang, P.; Tang, K.; Lozano, J.A. Estimation of distribution algorithms based unmanned aerial vehicle path planner using a new coordinate system. In Proceedings of the 2014 IEEE Congress on Evolutionary Computation (CEC), Beijing, China, 6–11 July 2014; pp. 1469–1476. [Google Scholar]
  43. Phung, M.D.; Ha, Q.P. Safety-enhanced UAV path planning with spherical vector-based particle swarm optimization. Appl. Soft Comput. 2021, 107, 107376. [Google Scholar] [CrossRef]
  44. Pan, J.S.; Liu, N.; Chu, S.C. A hybrid differential evolution algorithm and its application in unmanned combat aerial vehicle path planning. IEEE Access 2020, 8, 17691–17712. [Google Scholar] [CrossRef]
  45. Zhang, H.H.; Wang, T.; Peng, X. Underwater glider 3D path planning with adaptive segments and optimal motion parameters based on improved JADE algorithm. Ocean Eng. 2024, 299, 117377. [Google Scholar]
  46. Li, H.; Long, T.; Xu, G.; Wang, Y. Coupling-Degree-Based Heuristic Prioritized Planning Method for UAV Swarm Path Generation. In Proceedings of the 2019 Chinese Automation Congress (CAC), Hangzhou, China, 22–24 November 2019. [Google Scholar]
Figure 1. Main flow of CCEA-ADVS.
Figure 1. Main flow of CCEA-ADVS.
Drones 08 00435 g001
Figure 2. Problem decomposition.
Figure 2. Problem decomposition.
Drones 08 00435 g002
Figure 3. Top view in four scenarios. Blue represents missile radar; gray is the no-fly zone. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Figure 3. Top view in four scenarios. Blue represents missile radar; gray is the no-fly zone. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Drones 08 00435 g003
Figure 4. Convergence curves of the five algorithms. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Figure 4. Convergence curves of the five algorithms. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Drones 08 00435 g004
Figure 5. Box diagram of five algorithms run for 20 iterations. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Figure 5. Box diagram of five algorithms run for 20 iterations. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Drones 08 00435 g005
Figure 6. Fitness of the experiment.
Figure 6. Fitness of the experiment.
Drones 08 00435 g006
Figure 7. Runtimes of the experiment.
Figure 7. Runtimes of the experiment.
Drones 08 00435 g007
Figure 8. Flight path planned by the proposed algorithm. Blue represents missile radar; gray is the no-fly zone. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Figure 8. Flight path planned by the proposed algorithm. Blue represents missile radar; gray is the no-fly zone. (a) Scenario 1. (b) Scenario 2. (c) Scenario 3. (d) Scenario 4.
Drones 08 00435 g008
Table 2. CCEA-ADVS algorithm’s parameters.
Table 2. CCEA-ADVS algorithm’s parameters.
ParameterSetting
Maximum number of phase 1 iterations N i t e r 1 = 50
Maximum number of phase 2 iterations N i t e r 2 = 50
Population size N p o p = 100
Weight change rate of decision variable selection strategy ρ ω = 0.05
Internal iteration count of JADE N E A = 20
Zoom factor of JADE μ F = 0.5
Crossover probability of JADE μ C R = 0.5
Optimal individual density of JADE ρ = 0.06
Table 3. Parameter settings for the comparison algorithms.
Table 3. Parameter settings for the comparison algorithms.
AlgorithmParameter
JADE N i t e r = 100, N p o p = 100, F = 0.8, C R = 0.9, r e f r e s h = 10
PSO N i t e r = 100, N p o p = 100, w = 0.8, c 1 = 1.45, c 2 = 1.5
ABC N i t e r = 100, N p o p = 100, F o o d N u m b e r = 150, l i m i t = 100
GWO N i t e r = 100, N p o p = 100
Table 4. The comparison results of the five algorithms in terms of cost and runtime. The best results achieved among all algorithms are shown in bold.
Table 4. The comparison results of the five algorithms in terms of cost and runtime. The best results achieved among all algorithms are shown in bold.
Scenario PSOJADEABCGWOCCEA-ADVS
Scenario 1Best Cost871.1663.21312.7660615
Worst Cost1305.4849.71400.2733.3689
Mean Cost1204.1742.71358.0700.3647.6
Std Cost97.775.827.921.724.2
Runtime (s)985.5849.71028.2715.8431.3
Scenario 2Best Cost844.3642.1762.2810.3632
Worst Cost1050.67251479.8886.6664
Mean Cost945.5675.91111.3841.4643.1
Std Cost53.120.9269.524.58.6
Runtime (s)281.4379330.3300.1208.2
Scenario 3Best Cost690.5672.5694.2688.4663
Worst Cost740.2698.9987.2839.8686
Mean Cost709.8685.3813.7740.4672.1
Std Cost15.97.387.251.15.6
Runtime (s)570.9432.8452.2562.3375.4
Scenario 4Best Cost625.7605.1766.1684.6575.8
Worst Cost688.6625.9825.9742.8610.5
Mean Cost647.3616.2793.4709.6591.7
Std Cost19.85.418.921.513.0
Runtime (s)809.0733.1871.91091.0679.5
Table 5. Results of scalability analysis.
Table 5. Results of scalability analysis.
Number 40100200300400500
N iter
N i t e r 1 = 40, N i t e r 2 = 0Best Cost95224974949724610,51512,617
Time (s)155.2536.11186.53152.43752.16017
N i t e r 1 = 20, N i t e r 1 = 20Best Cost100526815172755411,09914,550
Time (s)96.2378.1891.32594.43510.75431.8
Table 6. The results of different encoding methods are compared. The best results achieved among all algorithms are shown in bold.
Table 6. The results of different encoding methods are compared. The best results achieved among all algorithms are shown in bold.
ScenarioNumberCost ( × 10 0 )Average Time (s)
ROCREPRESCARROCREPRESCAR
Scen 1510513317714256.865.353.871.1
1019221520922253.7103.392.1132.1
20420424489487248.1254.2266.2284.7
30623663670692320.6371.8438.6420.5
401045190616601881621.2642.7644.6676.3
Scen 2510510811114053.050.054.461.2
10200234208250113.6106.9102.4100.4
20778750808901269.7261.5245.5267.4
30611664665650620.3460.8512.4580.5
40810887892920610.6668.3815.6831.7
Scen 3510311511112059.453.755.652.5
10204220225224118.3112.7132.7131.2
20288311305291188.1214.2211.3204.6
30413443471478419.9279.5238.1345.2
40826882923865596.3597.7563.4632.4
Scen 45103110188167100.1102.196.7121.3
10225232246256198.6211.5200.4172.4
20426443454466453.5410.0453.8444.6
30640641686694871.5811.8875.7903.4
4086490211611421932.1978.41171.21241.2
Table 7. Comparison experiment of three decision variable selection methods. The best results achieved among all algorithms are shown in bold.
Table 7. Comparison experiment of three decision variable selection methods. The best results achieved among all algorithms are shown in bold.
Scenario ADVSGreedyRandom
Scenario 1Best Cost497.0501.7516.9
Worst Cost500.8506.7537.9
Mean Cost499.2504.8525.9
Std Cost1.21.16.9
Runtime (s)806.6849.71228.2
Scenario 2Best Cost466.2467.5471.0
Worst Cost467.7471.7480.1
Mean Cost467.0469.6475.5
Std Cost0.61.32.8
Runtime (s)571.6558.4774.8
Scenario 3Best Cost493.3494.0495.8
Worst Cost494.9495.7502.8
Mean Cost494.2494.9499.2
Std Cost0.50.62.1
Runtime (s)661.2686.4741.4
Scenario 4Best Cost581.1591.4598.4
Worst Cost587.0601.8608
Mean Cost584.1597.4603.7
Std Cost1.73.32.6
Runtime (s)702.1789.6812.6
Table 8. The comparison results of the four algorithms in terms of cost and runtime. The best results achieved among all algorithms are shown in bold.
Table 8. The comparison results of the four algorithms in terms of cost and runtime. The best results achieved among all algorithms are shown in bold.
Scenario PSOGADEJADE
Scenario 1Best Cost208213209201
Worst Cost219222213210
Mean Cost212.6217210.8205.4
Std Cost3.62.71.22.7
Runtime (s)98.5118.064.277.6
Scenario 2Best Cost203213205202
Worst Cost206222208204
Mean Cost204.7217.9206.4203.0
Std Cost1.02.71.00.7
Runtime (s)247.6215.6181.3206.6
Scenario 3Best Cost206216206204
Worst Cost217226208206
Mean Cost210.5221.5207.0205.0
Std Cost3.43.10.70.6
Runtime (s)230.9266.5190.8191.1
Scenario 4Best Cost206216212202
Worst Cost209218213206
Mean Cost207.3216.9212.5204.2
Std Cost0.90.60.31.2
Runtime (s)260.7293.9278.8253.7
Table 9. The comparison results of the two-stage strategies in terms of cost and runtime.
Table 9. The comparison results of the two-stage strategies in terms of cost and runtime.
Scenario N iter 1 = 100, N iter 2 = 0 N iter 1 = 0, N iter 2 = 100 N iter 1 = 50, N iter 2 = 50
Scenario 1Best Cost207207202
Worst Cost220210206
Mean Cost214208204
Std Cost4.81.11.3
Runtime (s)36.197.938.3
Scenario 2Best Cost204205203
Worst Cost206217204
Mean Cost205210203
Std Cost0.64.00.3
Runtime (s)36.8117.243.8
Scenario 3Best Cost205205202
Worst Cost208206203
Mean Cost206205202
Std Cost0.90.30.4
Runtime (s)38.682.149.2
Scenario 4Best Cost210211208
Worst Cost212213210
Mean Cost211212208
Std Cost0.70.50.6
Runtime (s)82.1160.5103.1
Table 10. A comparative analysis of the paths of phase 1 and phase 2 for the third group ( N i t e r 1 = 50, N i t e r 2 = 50).
Table 10. A comparative analysis of the paths of phase 1 and phase 2 for the third group ( N i t e r 1 = 50, N i t e r 2 = 50).
Scenario Scenario 1Scenario 2Scenario 3Scenario 4
Phase
Phase 1Path length235229235236
Time (s)13.214.819.744.8
Phase 2Path length204203202208
Time (s)38.343.849.2103.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Meng, Q.; Qu, Q.; Chen, K.; Yi, T. Multi-UAV Path Planning Based on Cooperative Co-Evolutionary Algorithms with Adaptive Decision Variable Selection. Drones 2024, 8, 435. https://doi.org/10.3390/drones8090435

AMA Style

Meng Q, Qu Q, Chen K, Yi T. Multi-UAV Path Planning Based on Cooperative Co-Evolutionary Algorithms with Adaptive Decision Variable Selection. Drones. 2024; 8(9):435. https://doi.org/10.3390/drones8090435

Chicago/Turabian Style

Meng, Qicheng, Qingjun Qu, Kai Chen, and Taihe Yi. 2024. "Multi-UAV Path Planning Based on Cooperative Co-Evolutionary Algorithms with Adaptive Decision Variable Selection" Drones 8, no. 9: 435. https://doi.org/10.3390/drones8090435

APA Style

Meng, Q., Qu, Q., Chen, K., & Yi, T. (2024). Multi-UAV Path Planning Based on Cooperative Co-Evolutionary Algorithms with Adaptive Decision Variable Selection. Drones, 8(9), 435. https://doi.org/10.3390/drones8090435

Article Metrics

Back to TopTop