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.
Category | Author | Algorithm | Year |
---|
DE | Fu [25] | DE + QPSO | 2013 |
| Yu [26] | DE + GWO | 2013 |
| Adhikari [27] | Adaptive Fuzzy DE algorithm | 2017 |
| Ali [28] | DE + MMACO | 2023 |
RRT | Fathy S. Elkazzaz [29] | DE + RRT | 2017 |
| Farzad Kiani [30] | GWO + RRT | 2021 |
| Fan Yang [31] | ACO + RRT | 2022 |
| MF Aslan [32] | PSO + GDRRT* | 2023 |
CCEA | Wu [33] | CCEA | 2023 |
| Shao [34] | CCEA + DCPSO | 2019 |
| Kesireddy [35] | CCEA + NSGA-II∖NSGA-III | 2019 |
Others | Fang [23] | APF | 2023 |
| Goel [19] | FA | 2018 |
| Chen [11] | ACO + GA | 2021 |
| Majumder [24] | MOTLBO | 2021 |
| Shao [13] | PSO | 2020 |
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 . 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 in the coordinate system. Then, the path of UAV m can be represented as , where represents the i-th waypoint of the m-th UAV. Our main objective is to find a set of paths 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:
where
represents the vector between two waypoints
and
, and
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
when the
m-th UAV passes through obstacles can be expressed as follows:
where
Z is the total number of no-fly zones, and
is the threat cost of the
z-th no-fly zone to the path segment
. The calculation method is as follows:
where
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
, and the maximum flight altitude be
. The flight altitude of the
m-th UAV should be between the two, and the cost associated with the path and altitude
is calculated as follows:
where
and
are the maximum and minimum values of the altitude limit, respectively. The function
represents the altitude at coordinates
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:
where
and
represent the projections of the path segments
and
on the xy plane, respectively.
Thus, the smoothness cost
can be expressed as
where
and
are the cost functions for the turning angle and climb angle, respectively:
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
where
is the given position of the
l-th combat platform,
, and
C is the number of combat platforms. Let the
l-th radar have a radius
of a circle. Then, the cost calculation formula between the UAV and the radar is
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
The missile kill cost of the
m-th UAV for the entire path is calculated as
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
, the difference in their arrival times (
and
) is checked to see if it is less than the safe time interval
to avoid collision risks. The arrival times
and
can be obtained by dividing the path by the speed. Thus, the total collision cost can be expressed as
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:
where
represents the weight coefficient of the
i-th constraint or cost.
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.