Next Article in Journal
Active Vibration Control and Parameter Optimization of Genetic Algorithm for Partially Damped Composites Beams
Previous Article in Journal
Development of a Vertical Submerging and Emerging Bat-Ray-Inspired Underwater Vehicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid Multi-Objective Chameleon Optimization Algorithm Based on Multi-Strategy Fusion and Its Applications

School of Intelligent Manufacturing and Electronic Engineering, Wenzhou University of Technology, Wenzhou 325035, China
*
Author to whom correspondence should be addressed.
Biomimetics 2024, 9(10), 583; https://doi.org/10.3390/biomimetics9100583
Submission received: 9 August 2024 / Revised: 23 September 2024 / Accepted: 23 September 2024 / Published: 25 September 2024

Abstract

:
Aiming at the problems of chameleon swarm algorithm (CSA), such as slow convergence speed, poor robustness, and ease of falling into the local optimum, a multi-strategy improved chameleon optimization algorithm (ICSA) is herein proposed. Firstly, logistic mapping was introduced to initialize the chameleon population to improve the diversity of the initial population. Secondly, in the prey-search stage, the sub-population spiral search strategy was introduced to improve the global search ability and optimization accuracy of the algorithm. Then, considering the blindness of chameleon’s eye turning to find prey, the Lévy flight strategy with cosine adaptive weight was combined with greed strategy to enhance the guidance of random exploration in the eyes’ rotation stage. Finally, a nonlinear varying weight was introduced to update the chameleon position in the prey-capture stage, and the refraction reverse-learning strategy was used to improve the population activity in the later stage so as to improve the ability of the algorithm to jump out of the local optimum. Eighteen functions in the CEC2005 benchmark test set were selected as an experimental test set, and the performance of ICSA was tested and compared with five other swarm intelligent optimization algorithms. The analysis of the experimental results of 30 independent runs showed that ICSA has stronger convergence performance and optimization ability. Finally, ICSA was applied to the UAV path-planning problem. The simulation results showed that compared with other algorithms, the paths generated by ICSA in different terrain scenarios are shorter and more stable.

Graphical Abstract

1. Introduction

In recent years, with the increasing scale and dimension of engineering optimization problems, it is difficult for traditional optimization algorithms to give the optimal solution [1,2]. High-dimensional function optimization problems have important application value in theoretical research and practical production processes. Many engineering problems require finding optimal solutions under various constraints, such as research on logistics distribution path optimization, high-speed rail dynamic scheduling method problems, and multi-class twin support vector machine parameter optimization problems [3,4,5]. After they are mathematically modeled, they are essentially function optimization problems. High-dimensional complex function optimization usually refers to function optimization problems with dimensions exceeding 100 dimensions [6,7]. As the dimensions of optimization problems increase, the search space and the number of local optima grow exponentially along with the computational complexity, making these problems increasingly difficult to solve [8,9,10]. The swarm intelligence optimization algorithm, inspired by nature, is widely used in solving complex optimization problems such as workshop scheduling and path planning because of its simple principle, easy implementation, and strong global search ability [11,12,13]. In industrial production, process parameters have a crucial impact on product quality, production quality, and production efficiency. Swarm intelligence optimization algorithms can establish a relationship model between parameters and production results by analyzing a large amount of production data in order to find the optimal combination of process parameters. Swarm intelligence optimization algorithms can analyze the historical operating data of devices and explore potential patterns between device operating status and faults and thus achieve the prediction of device faults. Swarm intelligence optimization algorithms can optimize production scheduling and resource allocation based on factors such as production tasks, equipment conditions, and personnel arrangements. With the deepening of research, swarm intelligence optimization algorithms are constantly emerging, such as the typical particle swarm optimization algorithm (PSO), ant colony optimization algorithm (ACO), dung beetle optimization algorithm (DBO), gorilla troops optimization algorithm (GTO), butterfly optimization algorithm (BOA), and so on [14,15].
The chameleon swarm algorithm (CSA) [16] is a swarm intelligence optimization algorithm proposed by Malik in 2021. It was inspired by the effective hunting behavior of chameleons in deserts and forests. Chameleons are a unique and highly specialized evolutionary branch with a wide range of species and are known for their ability to change color to blend in with their surroundings. The chameleon has excellent eyesight, and its eyes can independently rotate 360 degrees [17]. In the process of predation, its sticky tongue chases food and forms a small suction cup, which can catch food at a speed of 2590 m/s [18]. The main difference between the CSA algorithm and other biological behavior algorithms is that chameleons’ eyes can rotate 360°, which increases the probability of finding prey, and they can use their high-speed, sticky tongue to catch prey [19]. The prey-capture efficiency is very high, so the solution effect is better.
At present, many scholars have studied the CSA algorithm and applied it in different fields, such as economic load distribution optimization of power system, energy saving optimization of wireless network, energy optimization management of microgrids, path planning, etc. [20]. These studies have checked the flexibility and reliability of the algorithm to some extent.
The remainder of this paper is organized as follows: Section 2 reviews the related work on the CSA algorithm. Section 3 details the standard CSA algorithm. Section 4 introduces the proposed multi-strategy improved chameleon optimization algorithm (ICSA). Section 5 presents the experimental results and analysis of this algorithm on the CEC2005 test set. Section 6 validates the effectiveness of the ICSA algorithm in UAV 3D path planning. Finally, Section 7 concludes the study.

2. Related Work

Similar to other swarm intelligence optimization algorithms, the CSA is prone to falling into the local optimum and has weak global search ability. In order to enhance the optimization ability of CSA, many scholars have made many improvements to the CSA.
Anitha et al. proposed a modified gray wolf-based chameleon swarm (MGW-CS) algorithm, which enhances the global search capability of the CSA by introducing the leadership hierarchy and predation strategy of grey wolf optimization [21]. Braik et al. proposed an enhanced chameleon swarm algorithm (ECSA) by integrating roulette wheel selection, Lévy flight methods, and an add-on position updating strategy to develop a further balance between exploration and exploitation conducts [22]. Salawudeen, Ahmed et al. designed two novel variants of CSA called quasi-oppositional CSA (QCSA) and quasi-Lévy-oppositional CSA (QLCSA). In QCSA, a quasi-oppositional-based learning operator is used to improve the diversification and intensification of CSA. The QLCSA incorporates the Lévy flight operator into the QCSA to avoid stagnation and local minimal entrapment [23]. Hu, Gang et al. proposed an enhanced hybrid CSA called CCECSA. It introduces the crisscross optimization algorithm to avoid premature convergence. It incorporates an elite guidance mechanism to accelerate the convergence, adopts a competitive substitution mechanism to replace the worst individual, and implements an interference strategy to prevent the algorithm from getting trapped in local optima [24].
The aforementioned improved algorithms enhance the optimization capability of the standard CSA to varying degrees, yet issues such as premature convergence and susceptibility to local optima persist in some complex multimodal problems. To further enhance the performance of the CSA, this paper proposes a multi-strategy improved chameleon optimization algorithm (ICSA), with improvements in the following four aspects: (1) initializing the chameleon population using logistic mapping; (2) introducing a sub-population spiral search strategy during the prey-search phase; (3) updating chameleon positions during eyes’ rotation with a combination of the Lévy flight strategy with cosine adaptive factor and a greedy selection strategy; and (4) introducing nonlinear varying weight and a refraction reverse-learning strategy to stimulate the activity of the chameleon population during the prey-capture phase. Finally, the ICSA is applied to 3D UAV path planning, verifying the effectiveness of the algorithm in solving complex optimization problems.

3. Standard Chameleon Swarm Algorithm

In the standard CSA, chameleon food hunting is mainly divided into three main stages: searching for prey, discovering prey through eyes’ rotation, and hunting prey [25,26].

3.1. Search for Prey

Like other creatures in nature, chameleons continuously adjust their positions during roaming prey searches, using previous positions and social experiences to track and discover prey [27]. The chameleons’ movement behavior during searching can be characterized based on the updating strategy of position, as in Equation (1):
x t + 1 i , j = x t i , j + p 1 r 2 P t i , j G t j + p 2 r 1 G t j x t i , j r P p x t i , j + μ r 3 u b l b + l b sgn r a n d 0.5     r < P p
where x t i , j presents the current position of the i-th chameleon in the j-th dimension at iteration t. x t + 1 i , j is the new position of chameleon. P t i , j is the individual best position of the i-th chameleon in the j-th dimension at iteration loop t. G t j is the global best position of a chameleon by any chameleon in the t-th iteration.
Pp represents the probability of perceiving prey. p1, p2 are positive coefficients controlling the algorithm’s development capability, and r1, r2, r3, r, and rand are random numbers generated uniformly in the range of [0, 1] [28]. sgn(rand − 0.5) indicates the chameleon’s rotation direction, taking a value of +1 or −1. If rPp, the chameleon can alter its position based on the prey observed in the search space [29]. Similarly, when r < Pp, the chameleon will randomly explore the search space in various directions and regions to locate prey, increasing the likelihood of detecting nearby target prey. μ is a function of iterations parameter that decreases with the number of iteration, as in Equation (2) [30]:
μ = e k t / T m a x 3
The original study [16] also performed a sensitivity analysis of the algorithm’s parameters, confirming that the algorithm exhibits optimal search performance when k = 3.5. Here, Tmax represents the maximum number of iterations, and t denotes the current iteration number.

3.2. Chameleon Eyes’ Rotation

Chameleons have a unique advantage over other organisms: the ability to independently rotate their eyes 360° to precisely locate prey. This process can be outlined in four steps [31]:
  • The original position of the chameleon is the center of gravity (i.e., the origin);
  • Discover the rotation matrix that identifies the position of prey;
  • Update the chameleon’s position using the matrix of rotations at the center of gravity;
  • Translate the chameleons back to the original position.
The aforementioned steps can be represented by the following equations [32]:
x t + 1 i = m x t i x ¯ t i + x ¯ t i
where x t i presents the current position of the i-th chameleon at iteration t. x ¯ t i represents the center of the current position of the chameleon before rotation, which is also the average location of the individual in each dimension before position rotation. x t + 1 i represents the new position of the i-th chameleon after rotation. m is a rotation matrix, so m x t i x ¯ t i represents the rotating centered coordinates of the chameleon in the search space. m is defined as shown below [33]:
m = R θ , V z 1 , z 2
where z1 and z2 are two orthogonal vectors in an n-dimensional search space where the size of each vector is d. R denotes the rotation matrix on each axis [34].
θ = r · sgn r a n d 0.5 × π
where θ denotes the random rotation angle of the chameleon’s eyes. r is a random number in the range of [0, 1], and sgn(rand − 0.5) is either 1 or −1. Thus, the θ range is [−π, π].

3.3. Hunting Prey

The chameleon uses its tongue to attack prey. Once the tongue makes contact with the prey, it rapidly adheres to the surface of the prey, creating a small suction cup-like effect for effective capture [35]. Its position is updated slightly, as it can drop its tongue to double its length. The speed of the chameleon’s tongue when it is extended toward prey can be mathematically modeled according to Equation (6) [36]:
v t + 1 i , j = λ t v t i , j + c 1 r 1 G t j x t i , j + c 2 r 2 P t i , j x t i , j λ t = 1 t / T max ρ t / T max
where v t i , j presents the current velocity of the i-th chameleon in the j-th dimension at iteration t, and v t + 1 i , j is new velocity of chameleon. r1 and r2 are random numbers in the range [0, 1]. c1 and c2 are two positive constants that control the influence of G t j and P t i , j on the chameleon’s speed. λ t is the inertia weight that is linearly decreased with the iterative generations, and ρ is equal to 1 [37].
Equation (7) was proposed to update the position of the chameleon:
x t + 1 i , j = x t i , j + v t i , j 2 v t 1 i , j 2 2 a
where v t i , j represents the current speed, v t 1 i , j denotes the speed from the previous iteration, and x t i , j indicates the current position of the chameleon. a represents acceleration rate of the chameleon’s tongue projection, and it is defined as given in Equation (8):
a = 2590 × 1 e log t
where t denotes the current iteration number.

4. Multi-Strategy of the Chameleon Optimization Algorithm

4.1. Logistic Chaotic Map Initialization

The population initialization method has certain influence on the optimization effect and robustness of swarm intelligence algorithm. In the CSA algorithm, the initial position of chameleon swarm is randomly distributed, which may lead to uneven distribution or excessive concentration of the population and then lead to premature convergence of the algorithm falling into the local optimum or the convergence speed slowing down, which affects the final optimization effect and reduces the robustness of the algorithm [38].
Therefore, introducing appropriate chaotic sequence initialization can enrich the diversity of the initial population and reduce its above-mentioned correlation possibilities. At present, commonly used chaos optimization methods include logistic mapping, Singer mapping, Bernoulli mapping, Gauss mapping, etc. Among them, logistic mapping features a heterogeneous distribution, which helps prevent the initial population from being either too tightly clustered or too dispersed. Therefore, this paper selects logistic mapping to initialize the population, and the specific operation is shown in Equation (9):
x i = l b j + L j × u b j l b j i = 1 , 2 , , N , j = 1 , 2 , , d i m
where xi is the initial position of the i-th chameleon among N chameleon individuals. ubj and lbj are the lower and upper bounds of the search space, respectively. Lj is the generated logistic chaotic sequence. The definition of the logistic mapping is shown in Equation (10):
L j + 1 = r L j 1 L j
where L j 0 , 1 . r is the control parameter, and r is set to 0.3 in this paper.

4.2. Sub-Population Spiral Search Strategy

In the standard CSA, during the prey-search phase, chameleons have a 90% probability of moving towards the individual best position and the global best position. This local search method, however, is prone to getting trapped in the local optimum. Chameleons have a 10% probability of randomly exploring different directions and areas within the search space. This is a global search method, but it lacks guidance.
In order to balance the local and global search ability of the population at this stage, this paper firstly divides chameleon individuals into class A and class B according to their fitness values. Class A is chameleon individuals whose fitness value is lower than the average fitness value of the population, and class B is chameleon individuals whose fitness value is higher than the average fitness value of the population.
For a class A chameleon, the position is updated by using the above formula in Equation (11) in order to guide it to quickly search for prey by the individual best position and the global best position and speed up the convergence of the algorithm. For a class B chameleon, a new position update method is proposed. That is to say, on the basis of retaining a certain search randomness, guided by the global best position and combined with the variable spiral search strategy, more local search paths are developed for class B chameleons, thus improving the ability of the algorithm to jump out of the local optimum and find the global optimum. In summary, the position update formula of a chameleon in the prey-search stage is shown in Equation (11):
x t + 1 i , j = x t i , j + p 1 r 2 P t i , j G t j + p 2 r 1 G t j x t i , j , f i t t i < m e a n f i t t x t + 1 i , j = G t j + e z l · cos 2 π l · G t j x t i , j · sgn rand 0.5 , f i t t i m e a n f i t t z = e k · cos π · 1 t / T
where f i t t i is the fitness value of i-th chameleon at iteration t, meanfitt is the average fitness value of the population at iteration t, l is the uniform random number of [−1, 1], and k is the variation coefficient, which is set to 5 in this paper.

4.3. Lévy Flight Strategy with Cosine Adaptive Factor and Greedy Selection Strategy

4.3.1. Lévy Flight Strategy with Cosine Adaptive Factor

In the stage of chameleons’ eyes turning to find prey, the chameleon changes the search direction and search area to approach the center of gravity with the help of the current-position center of gravity and rotation matrix [39]. While this process can enhance the global search capability of the chameleon algorithm, it introduces a high degree of randomness. This randomness may cause the algorithm to inadequately cover the search space during early iterations, potentially missing the optimal solution. Additionally, it may lack adaptability in later stages, resulting in suboptimal search performance [40].
To solve this problem, this paper introduces the Lévy flight strategy with cosine adaptive factor. The formula for calculating the cosine adaptive factor is as follows:
c t = 0.075 1 + cos π t T max
where t denotes the current iteration number. Tmax represents the maximum number of iterations. The variation curve of cosine adaptive factor ct with the number of iterations is shown in Figure 1.
Lévy flight is a random flight with alternating near and far search ranges, and the probability distribution of its step size is heavy-tailed distribution. When defining walking in a space larger than one dimension, the walking step size is a random direction of each peer. The expression for the Lévy step length is as follows:
L evy β = u v 1 β
where u ~ N 0 , σ u 2 , v ~ N 0 , σ v 2 . σ u , σ v are defined as follows:
σ u = Γ 1 + β sin π β / 2 2 β 1 / 2 Γ 1 + β / 2 β σ v = 1
where Γ represents the standard gamma function. β is a constant in the range of [0, 2], and it is set to 1.5 in this paper.
In summary, after updating its position through eye rotation, the chameleon uses the Lévy flight strategy with cosine adaptive factor to expand its search area. As shown in Figure 1, the cosine adaptive factor ct stabilizes at a higher value during early iterations, increasing the probability of long-step movements, which helps expand the global search range. In later iterations, ct stabilizes at a lower value, increasing the probability of short-step movements, which allows for more precise local searches and enhances optimization accuracy. Furthermore, considering the guiding influence of the group’s optimal position, the position update formula after incorporating the Lévy flight strategy with cosine adaptive factor during the eye-rotation stage is presented in Equation (15):
x l t i , j = x t i , j + c t · G t j x t i , j · L evy β t i , j

4.3.2. Greedy Selection Strategy

Due to the random characteristics of Lévy flight, the updated chameleon position does not necessarily have better fitness. Thus, the greedy strategy is adopted to choose whether to update the new position, and the formula is shown in Equation (16):
x t i , j = x t i , j , f i t n e s s x t i , j f i t n e s s x l t i , j x l t i , j , f i t n e s s x t i , j < f i t n e s s x l t i , j
where f i t n e s s x t i , j and f i t n e s s x l t i , j represent the fitness value of the chameleon’s eye-rotation position update and the fitness value calculated by the Lévy flight strategy update position with cosine adaptive factor, respectively.

4.4. Nonlinear Varying Weight and Refraction Reverse-Learning Strategy

4.4.1. Nonlinear Varying Weight

As seen in Equation (6), the speed update during the chameleon’s prey-hunting stage is similar to that of particle swarm optimization (PSO), and c1 and c2 are the control factors of individual optimal position and global optimal position for velocity update. In the early stage of algorithm iteration, the chameleon position is far away from the optimal position. In order to achieve faster convergence speed, it is expected to increase c1 and c2. In the later iteration stage of the algorithm, the chameleon position is relatively close to the optimal position. In order to ensure the convergence accuracy and avoid falling into the local optimum, it is expected to reduce c1 and c2. In this paper, a nonlinear varying weight is adopted, as shown in Equation (17):
ω t = 1 t T max 2 t T max
where t denotes the current iteration number. Tmax represents the maximum number of iterations. The variation curve of nonlinear varying weight ωt with the number of iterations is shown in Figure 2.
The speed update formula is changed as follows:
v t + 1 i , j = λ t v t i , j + ω t r 1 G t j x t i , j + ω t r 2 P t i , j x t i , j

4.4.2. Refraction Reverse-Learning Strategy

Similar to many swarm intelligence optimization algorithms, the CSA can easily fall into the local optimum in the later iteration stage because of its high population similarity and low activity. To address this issue, this paper introduces a refraction reverse-learning strategy following the chameleon’s prey-hunting position update. The refraction reverse-learning strategy is based on the refraction principle of light, and its principle is shown in Figure 3 below.
In Figure 3, O is the midpoint of ub and lb, which are the upper and lower limits of the interval, respectively. The refracted ray OX*′ is formed by the refraction of the incident ray XO by boundary mass, and the heights are h′ and h, respectively.
Assuming that the refractive index of free space and dielectric space is the same, then θ1 = θ2. X* and X, respectively, represent the projections of OX*′ and XO on the x-axis, respectively. At this time, the expression of the refraction imaging principle can be described as follows:
l b + u b 2 X X * l b + u b 2 = h h
Let k = h h . k denotes the dynamic scaling factor, and its value is equal to t, which changes linearly with the current number of iterations; then, the expression of X* is as follows:
X * = l b + u b 2 + l b + u b 2 t X t
Through the refraction reverse-learning strategy, the search radius is dynamically adjusted. In the early stage of algorithm iteration, potential feasible solutions are developed on a large scale to increase the possibility of global search. In the later stage of algorithm iteration, high-quality solutions are developed in a small-scale and fine way to improve population activity, which will improve the optimization accuracy and convergence speed to a certain extent. However, the updated chameleon position does not necessarily have better fitness, so it is necessary to adopt greedy strategy according to Equation (16) to choose whether to update the new position.

4.5. Pseudo-Code of the Multi-Strategy Improved Chameleon Algorithm

To sum up, the multi-strategy improved chameleon algorithm (ICSA) (Algorithm 1) proposed in this paper can increase the group diversity to a certain extent, balance the local search and global search capabilities of the algorithm, and improve the optimization ability of the algorithm. The pseudo-code of the ICSA is as follows:
Algorithm 1: ICSA pseudo-code
Begin
1.  
Setting parameters: N ← population size, dim ← spatial dimension,
         Tmax ← maximum number of iterations,
         lb ← lower limit of the search space,
         ub ← upper limit of the search space.
2.  
Initialize the population X by logistic chaotic mapping using Equation (9).
3.  
Calculate the individual fitness to determine the individual best position P and the global best position G.
4.  
While (t ≤ Tmax).
5.  
Calculate the mean individual fitness meanfitt.
6.  
for i = 1 to N
7.  
Update the chameleon’s position based on Equation (11).
8.  
end for
9.  
Compute the cosine adaptive factor ct using Equation (12).
10.
for i = 1 to N
11.
Update the chameleon’s position according to Equations (3)–(5).
12.
Calculate the chameleon’s random position using the Lévy flight strategy with the
  cosine adaptive factor using Equation (15).
13.
Determine whether to update the chameleon’s position to the random position based
  on the greedy selection using Equation (16).
14.
end for
15.
Calculate the nonlinear varying weight ωt according to Equation (17).
16.
for i = 1 to N
17.
Update the chameleon’s velocity using Equation (18).
18.
Update the chameleon’s position using Equation (7).
19.
Calculate the chameleon’s expanded position by incorporating the refraction reverse
  learning strategy as defined by Equation (20).
20.
Determine whether to update the chameleon’s position to the expanded position based
  on the greedy selection using Equation (16).
21.
end for
22.
Calculate the individual fitness and update the individual best position Pt and the global
 best position Gt.
23.
t = t + 1.
24.
end while
25.
Output the fitness array, the global best position, and the optimal fitness value.
End

4.6. Flowchart of the Multi-Strategy Improved Chameleon Algorithm

The flowchart of the multi-strategy improved chameleon algorithm (ICSA) is shown in Figure 4.

5. Experimental Results and Simulation Analysis

5.1. Experimental Environment Settings

To verify the performance of the ICSA algorithm, the experimental environment was configured as follows: The operating system was Windows 10 (64-bit), the processor was an Intel(R) Core(TM) i7-10875H with a main frequency of 2.3 GHz, the computer memory was 16 GB, and the simulation software used was MATLAB R2020b.

5.2. Comparison Algorithms and Parameter Settings

To validate the performance of the ICSA algorithm, we selected the classic CSA algorithm and four other swarm intelligence optimization algorithms that are either classic or have shown good performance in recent years for comparative analysis. These algorithms include the dung beetle optimizer (DBO) [41], the subtraction-average-based optimizer (SABO) [42], the gorilla troop optimization (GTO) [43], and the ant colony optimization (ACO) [44].
CEC2005 was chosen as the benchmark function. To ensure the reliability of the test results, the population size for all algorithms was set to 30 and the number of iterations to 1000. The specific parameter settings for each algorithm are detailed in Table 1.

5.3. Test Functions

To fairly evaluate the effectiveness of the ICSA algorithm, we selected 18 test functions from the CEC2005 test set for validation. These included six unimodal functions (F1(x)~F6(x)), five multimodal functions (F7(x)~F11(x)), and seven fixed-dimension functions (F12(x)~F18(x)). Unimodal functions were used to evaluate the algorithm’s local search capability, multimodal functions to evaluate its global search capability, and fixed-dimension functions to assess the balance between exploration and exploitation. Detailed information about the benchmark test functions is presented in Table 2.

5.4. Optimization Results and Analysis of Test Functions

Using the experimental environment and parameter settings described above, each algorithm was run independently 30 times on the 18 test functions. The results are presented in Table 3, including metrics such as average value, standard deviation, best value, and worst value.
To eliminate randomness, the average value of the experimental results was used to measure the optimization accuracy of the algorithms. By using the mean value as a comparison metric, it is evident that the ICSA algorithm achieves better solutions in the test functions compared to the standard CSA algorithm, indicating that the improvements enhance the algorithm’s optimization accuracy. Compared to other optimization algorithms, ICSA consistently finds the optimal or near-optimal solutions, demonstrating its strong optimization capability, with GTO being the next best performer.
The standard deviation of the experimental results measures the degree of deviation in the algorithm’s results and is a key indicator of algorithm stability. Using the standard deviation as a comparison metric, it is clear that ICSA exhibits stable and superior optimization performance compared to both the standard CSA algorithm and other optimization algorithms.
In summary, within the specified number of iterations, the optimization accuracy and robustness of ICSA show significant advantages.

5.5. Convergence Curve Analysis

The convergence curve visually demonstrates the convergence speed and optimization accuracy of the algorithms. Figure 5 shows the convergence curves for each algorithm obtained simultaneously in the aforementioned experiments. The horizontal axis of each curve represents the number of iterations, and the vertical axis represents the fitness value.
It is evident that, compared to the standard CSA algorithm, ICSA significantly accelerates convergence and enhances optimization accuracy. Compared to other optimization algorithms, ICSA also demonstrates substantial advantages in both convergence speed and optimization accuracy. In summary, within the specified number of iterations, ICSA exhibits significant advantages in convergence speed and optimization accuracy.

6. Application and Analysis of ICSA in UAV Path Planning

To validate the effectiveness of the ICSA algorithm in addressing real-world optimization problems, this paper applies ICSA to the UAV 3D path planning. The objective of UAV 3D path planning is to find a safe and efficient route from the starting position to the target position under specific constraints [45].

6.1. Terrain Obstacle Model

For the UAV path-planning problem, terrain obstacle modeling is the first step. This study selected a three-dimensional space of 1000 m × 1000 m × 100 m, used an exponential function to simulate mountain peaks as the terrain, and modeled obstacle threat areas using cylinders [46]. The mathematical model for the mountain peak terrain is as follows:
z x , y = i = 1 n h i exp x x i x s i 2 y y i y s i 2
where z represents the mountain peak terrain data. n denotes the number of mountain peaks. hi is the height of the i-th mountain peak. [xi, yi] is the center coordinate of the i-th mountain peak, and xsi and ysi are the descending slopes of the i-th mountain peak along the x-axis and y-axis, respectively.

6.2. Cost Function

6.2.1. Path Length Cost

Path length is a crucial factor in UAV path planning. A shorter path length indicates a shorter time required for the UAV to complete the mission, thereby reflecting higher path-planning quality [47]. After determining the starting and ending points of the path, M track points are selected between these points. Using the cubic spline interpolation method, the planned path is then discretized into a path with N nodes, where each node Xi corresponds to a spatial coordinate (xi, yi, zi). The path length cost is the sum of the Euclidean distances between adjacent nodes in the path, calculated as follows:
F 1 X = i = 1 N 1 x i + 1 x i 2 + y i + 1 y i 2 + z i + 1 z i 2

6.2.2. Flight Altitude Cost

Choosing the correct flight altitude is also critical for UAV path planning. A reasonable altitude can reduce collision risk and enhance safety while also affecting flight speed and energy consumption, thus influencing flight efficiency [48].
Ground obstacle threat areas are modeled as cylinders, with mountain peaks representing the terrain. Therefore, effective path planning should ensure the flight path avoids intersecting threat areas and that the flight altitude at any node is above the corresponding terrain height. If the flight path is valid, the flight altitude cost is calculated using the formula in the upper part of Equation (23). Otherwise, the flight altitude cost is set to 1032. The total flight altitude cost calculation is as follows:
F 2 X = n = 1 N z n z ¯ , IsPathOk = true 10 32 , IsPathOk = false
where z n represents the altitude of each path node above ground level, and z ¯ denotes the average altitude of the respective path nodes.

6.2.3. Turning Angle Cost

The UAV’s turning angle directly affects its stability and maneuverability. During flight, the turning angle should not exceed the predefined maximum turning angle. The formula for angle change during UAV flight is as follows:
θ i = arccos α i β i α i β i
where αi represents the difference in spatial coordinates between the path node Xi+1 (xi+1, yi+1, zi+1) and the path node Xi (xi, yi, zi). βi represents the difference in spatial coordinates between the path node Xi+2 (xi+2, yi+2, zi+2) and the path node Xi+1 (xi+1, yi+1, zi+1). The formula for calculating the turning angle cost is as follows:
F 3 X = i = 1 N 1 θ i + 1 θ i

6.2.4. Total Cost Function

By integrating the three path optimization-related cost functions discussed above, the total cost function is defined as the weighted sum of the following components:
F X = b 1 F 1 X + b 2 F 2 X + b 3 F 3 X
The weights b1, b2, and b3 are coefficients that satisfy b1 + b2 + b3 = 1, and they can be adjusted according to different scenarios. Under the condition of safe obstacle avoidance, a smaller total cost function value indicates smoother and more efficient path planning. In this paper, the weights b1, b2, and b3 are set to 0.4, 0.4, and 0.2, respectively.

6.3. UAV Path-Planning Simulation

In this study, two different threat obstacle terrains were selected. Terrain 1 contains seven mountain peaks and one cylindrical obstacle threat area, while terrain 2 contains four mountain peaks and one cylindrical obstacle threat area. The specific parameter settings for the peaks and cylindrical obstacle threat areas are provided in Table 4 and Table 5.
The previously proposed ICSA was applied alongside five other swarm intelligence optimization algorithms for path-planning simulations and comparative analysis in these terrains. Each algorithm was tested independently 20 times, with the population size set at 100 and the maximum number of iterations set at 1000. The starting coordinates were (0, 0, 20) and the ending coordinates were (977, 970, 60). The path-planning simulation results for the two different terrains are shown in Figure 6.
As shown in Figure 6, terrain 1 is more complex, while terrain 2 is relatively simpler. In both terrains, the ICSA algorithm demonstrates superior path planning compared to the standard CSA algorithm, resulting in shorter flight distances. Figure 7 indicates that the ICSA algorithm achieves a lower total cost, faster convergence speed, and better path-planning performance.
Table 6 provides a summary of the UAV path-planning results in the two different terrains, including metrics such as the average value, standard deviation, best value, and worst value.
Using the average value as a comparison metric, it is evident that the ICSA achieves a smaller total cost function value compared to the standard CSA algorithm, demonstrating that the improved strategy can effectively enhance the algorithm’s optimization accuracy. When compared to other optimization algorithms, the ICSA consistently finds the optimal solution or gets closer to it, indicating significant advantages, with the next best performance observed in the GTO algorithm.
The standard deviation of the experimental results primarily measures the deviation of the algorithm’s results, serving as a key indicator of algorithm stability. Using standard deviation as a comparison metric, the ICSA exhibits stable and excellent optimization performance relative to both the standard CSA algorithm and other optimization algorithms. Overall, the ICSA algorithm demonstrates superior performance in solving the 3D path-planning problem.

7. Conclusions

The standard CSA algorithm faces issues such as slow convergence speed, poor robustness, and a tendency to get trapped in local optima, indicating significant room for performance improvement in practical applications. Thus, this paper introduces a multi-strategy improved chameleon optimization algorithm (ICSA), incorporating strategies such as logistic mapping, sub-population spiral search strategy, Lévy flight with cosine adaptive factor, greedy strategy, and refraction reverse-learning strategy.
Simulation experiments using the CEC2005 standard test functions demonstrated that the ICSA algorithm offers faster convergence speed, higher optimization accuracy, and greater robustness compared to five other swarm intelligence optimization algorithms. Finally, the ICSA was applied to the UAV path-planning problem, accounting for path length cost, altitude cost, and turning angle cost. The ICSA shows superior performance in terms of path-planning length, convergence speed, and stability, effectively achieving 3D path planning across different terrains.
Like other intelligent optimization algorithms, ICSA has strong adaptability. It can also be verified through experiments that ICSA shows good performance in classic engineering problems such as pressure vessel design, three-bar truss design, and tension-compression spring design. Going forward, the performance of the ICSA algorithm could be further enhanced based on the current research. It is believed that ICSA will have certain application value in various industrial fields.

Author Contributions

Conceptualization, Y.C.; writing—review and editing, Y.C.; validation, L.C.; data curation, L.C.; writing—original draft preparation, Y.Y. methodology, Y.C.; software, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Natural Science Foundation of Zhejiang Province under Grant LY23F010002 and the Industrial Science and Technology Project of Yueqing City under Grant 2022G007.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon request. There are no restrictions on data availability.

Conflicts of Interest

The authors declare no competing interests.

References

  1. Tang, J.; Duan, H.; Lao, S. Swarm intelligence algorithms for multiple unmanned aerial vehicles collaboration: A comprehensive review. Artif. Intell. Rev. 2023, 56, 4295–4327. [Google Scholar] [CrossRef]
  2. Tang, J.; Liu, G.; Pan, Q. A review on representative swarm intelligence algorithms for solving optimization problems: Applications and trends. IEEE/CAA J. Autom. Sin. 2021, 8, 1627–1643. [Google Scholar] [CrossRef]
  3. Madani, A.; Engelbrecht, A.; Ombuki-Berman, B. Cooperative coevolutionary multi-guide particle swarm optimization algorithm for large-scale multi-objective optimization problems. Swarm Evol. Comput. 2023, 78, 101262. [Google Scholar] [CrossRef]
  4. Nourani, B.; Arvanaghi, H.; Pourhosseini, F.A.; Javidnia, M.; Abraham, J. Enhanced Support Vector Machine with Particle Swarm Optimization and Genetic Algorithm for Estimating Discharge Coefficients of Circular-Crested Oblique Weirs. Iran. J. Sci. Technol. Trans. Civ. Eng. 2023, 47, 3185–3198. [Google Scholar] [CrossRef]
  5. Chen, C.; Cao, L.; Chen, Y.; Chen, B.; Yue, Y. A comprehensive survey of convergence analysis of beetle antennae search algorithm and its applications. Artif. Intell. Rev. 2024, 57, 141. [Google Scholar] [CrossRef]
  6. Cao, L.; Chen, H.; Chen, Y.; Yue, Y.; Zhang, X. Bio-Inspired Swarm Intelligence Optimization Algorithm-Aided Hybrid TDOA/AOA-Based Localization. Biomimetics 2023, 8, 186. [Google Scholar] [CrossRef]
  7. Yue, Y.; Cao, L.; Chen, H.; Chen, Y.; Su, Z. Towards an Optimal KELM Using the PSO-BOA Optimization Strategy with Applications in Data Classification. Biomimetics 2023, 8, 306. [Google Scholar] [CrossRef]
  8. Chen, B.; Cao, L.; Chen, C.; Chen, Y.; Yue, Y. A comprehensive survey on the chicken swarm optimization algorithm and its applications: State-of-the-art and research challenges. Artif. Intell. Rev. 2024, 57, 170. [Google Scholar] [CrossRef]
  9. Yue, Y.; Cao, L.; Zhang, Y. Novel WSN Coverage Optimization Strategy Via Monarch Butterfly Algorithm and Particle Swarm Optimization. Wirel. Pers. Commun. 2024, 135, 2255–2280. [Google Scholar] [CrossRef]
  10. Jiang, S.; Yue, Y.; Chen, C.; Chen, Y.; Cao, L. A Multi-Objective Optimization Problem Solving Method Based on Improved Golden Jackal Optimization Algorithm and Its Application. Biomimetics 2024, 9, 270. [Google Scholar] [CrossRef]
  11. Wang, S.; Cao, L.; Chen, Y.; Chen, C.; Yue, Y.; Zhu, W. Gorilla optimization algorithm combining sine cosine and cauchy variations and its engineering applications. Sci. Rep. 2024, 14, 1–20. [Google Scholar] [CrossRef]
  12. Jaafari, A.; Panahi, M.; Mafi-Gholami, D.; Rahmati, O.; Shahabi, H.; Shirzadi, A.; Lee, S.; Bui, D.T.; Pradhan, B. Swarm intelligence optimization of the group method of data handling using the cuckoo search and whale optimization algorithms to model and predict landslides. Appl. Soft Comput. 2022, 116, 108254. [Google Scholar] [CrossRef]
  13. Yue, Y.; Cao, L.; Lu, D.; Hu, Z.; Xu, M.; Wang, S.; Li, B.; Ding, H. Review and empirical analysis of sparrow search algorithm. Artif. Intell. Rev. 2023, 56, 10867–10919. [Google Scholar] [CrossRef]
  14. Tang, W.; Cao, L.; Chen, Y.; Chen, B.; Yue, Y. Solving Engineering Optimization Problems Based on Multi-Strategy Particle Swarm Optimization Hybrid Dandelion Optimization Algorithm. Biomimetics 2024, 9, 298. [Google Scholar] [CrossRef] [PubMed]
  15. Yue, Y.; You, H.; Wang, S.; Cao, L. Improved whale optimization algorithm and its application in heterogeneous wireless sensor networks. Int. J. Distrib. Sens. Netw. 2021, 17, 15501477211018140. [Google Scholar] [CrossRef]
  16. Braik, M.S. Chameleon Swarm Algorithm: A bio-inspired optimizer for solving engineering design problems. Expert Syst. Appl. 2021, 174, 114685. [Google Scholar] [CrossRef]
  17. Said, M.; El-Rifaie, A.M.; Tolba, M.A.; Houssein, E.H.; Deb, S. An efficient chameleon swarm algorithm for economic load dispatch problem. Mathematics 2021, 9, 2770. [Google Scholar] [CrossRef]
  18. Sankara Srinivasa Rao, I.; Rajalakshmi, N.R. Energy efficient optimal parent selection with Chameleon Swarm algorithm and non-overlap backoff algorithm based dynamic slot allocation in wireless body area network. Wirel. Netw. 2024, 30, 2219–2234. [Google Scholar] [CrossRef]
  19. Wang, J.; Lv, M.; Li, Z. Multivariate selection-combination short-term wind speed forecasting system based on convolution-recurrent network and multi-objective chameleon swarm algorithm. Expert Syst. Appl. 2023, 214, 119129. [Google Scholar] [CrossRef]
  20. Ballous, K.A.; AlShabi, M.; Nassif, A.B.; Bettayeb, M. Path planning and obstacle avoidance utilizing chameleon swarm algorithm. Unmanned Syst. Technol. 2023, 12549, 68–74. [Google Scholar]
  21. Zhou, J.; Xu, Z. Optimal sizing design and integrated cost-benefit assessment of stand-alone microgrid system with different energy storage employing chameleon swarm algorithm: A rural case in Northeast China. Renew. Energy 2023, 202, 1110–1137. [Google Scholar] [CrossRef]
  22. Zahraoui, Y.; Alhamrouni, I.; Mekhilef, S. A Chameleon algorithm for solving economic dispatch problem in microgrid system. Bull. Electr. Eng. Inform. 2023, 12, 1982–1992. [Google Scholar] [CrossRef]
  23. Salawudeen, A.T.; Moritz, M.; Jahn, I.; Johnson, O.; Monti, A. Enhanced Chameleon Swarm Algorithms for Nested Identical Control of Load Frequency in Autonomous Microgrid. IEEE Access 2024, 12, 42544–42571. [Google Scholar] [CrossRef]
  24. Hu, G.; Yang, R.; Wei, G. Hybrid chameleon swarm algorithm with multi-strategy: A case study of degree reduction for disk Wang–Ball curves. Math. Comput. Simul. 2023, 206, 709–769. [Google Scholar] [CrossRef]
  25. Braik, M.S.; Hammouri, A.I.; Awadallah, M.A.; Al-Betar, M.A.; Khtatneh, K. An improved hybrid chameleon swarm algorithm for feature selection in medical diagnosis. Biomed. Signal Process. Control. 2023, 85, 105073. [Google Scholar] [CrossRef]
  26. Rizk-Allah, R.M.; Hassanien, A.E.; Snášel, V. A hybrid chameleon swarm algorithm with superiority of feasible solutions for optimal combined heat and power economic dispatch problem. Energy 2022, 254, 124340. [Google Scholar] [CrossRef]
  27. Elaziz, M.A.; Ahmadein, M.; Ataya, S.; Alsaleh, N.; Forestiero, A.; Elsheikh, A.H. A quantum-based chameleon swarm for feature selection. Mathematics 2022, 10, 3606. [Google Scholar] [CrossRef]
  28. Liu, G.; Yuan, J.; Lin, K.P.; Miao, Y.; Li, R. Developing an improved chameleon swarm algorithm for combined cooling, heating and power micro-grid system. Expert Syst. Appl. 2024, 237, 121540. [Google Scholar] [CrossRef]
  29. Mostafa, R.R.; Houssein, E.H.; Hussien, A.G.; Singh, B.; Emam, M.M. An enhanced chameleon swarm algorithm for global optimization and multi-level thresholding medical image segmentation. Neural Comput. Appl. 2024, 36, 8775–8823. [Google Scholar] [CrossRef]
  30. Dinh, P.H. Medical image fusion based on enhanced three-layer image decomposition and chameleon swarm algorithm. Biomed. Signal Process. Control. 2023, 84, 104740. [Google Scholar] [CrossRef]
  31. Braik, M.S. Modified chameleon swarm algorithm for brightness and contrast enhancement of satellite images. Multimed. Tools Appl. 2024, 83, 26819–26870. [Google Scholar] [CrossRef]
  32. Zhou, W.; Wang, L.; Han, X.; Wang, Y.; Zhang, Y.; Jia, Z. Adaptive density spatial clustering method fusing chameleon swarm algorithm. Entropy 2023, 25, 782. [Google Scholar] [CrossRef] [PubMed]
  33. Uluslu, A. Chameleon Swarm Algorithm assisted optimization of U-slot patch antenna for quad-band applications. IEEE Access 2022, 10, 74152–74163. [Google Scholar] [CrossRef]
  34. Braik, M.S.; Awadallah, M.A.; Al-Betar, M.A.; Hammouri, A.I.; Zitar, R.A. A non-convex economic load dispatch problem using chameleon swarm algorithm with roulette wheel and levy flight methods. Appl. Intell. 2023, 53, 17508–17547. [Google Scholar] [CrossRef]
  35. Anitha, S.; Saravanan, S.; Chandrasekar, A. A modified gray wolf-based chameleon swarm algorithm for minimizing energy consumption and enabling secure communication in wireless sensor network. Concurr. Comput. Pract. Exp. 2022, 34, e7295. [Google Scholar] [CrossRef]
  36. Ren, H.; Hou, X.; Jia, Z.; Mashhadi, A. A new optimal energy management strategy of microgrids using chaotic map-based chameleon swarm algorithm. IET Renew. Power Gener. 2023, 17, 1978–1992. [Google Scholar] [CrossRef]
  37. Rizk-Allah, R.M.; El-Hameed, M.A.; El-Fergany, A.A. Model parameters extraction of solid oxide fuel cells based on semi-empirical and memory-based chameleon swarm algorithm. Int. J. Energy Res. 2021, 45, 21435–21450. [Google Scholar] [CrossRef]
  38. Tian, D. Particle swarm optimization with chaos-based initialization for numerical optimization. Intell. Autom. Soft Comput. 2017, 3, 1–12. [Google Scholar] [CrossRef]
  39. Houssein, E.H.; Saad, M.R.; Hashim, F.A.; Shaban, H.; Hassaballah, M. Lévy flight distribution: A new metaheuristic algorithm for solving engineering optimization problems. Eng. Appl. Artif. Intell. 2020, 94, 103731. [Google Scholar] [CrossRef]
  40. Li, J.; An, Q.; Lei, H.; Deng, Q.; Wang, G.G. Survey of lévy flight-based metaheuristics for optimization. Mathematics 2022, 10, 2785. [Google Scholar] [CrossRef]
  41. Xue, J.; Shen, B. Dung beetle optimizer: A new meta-heuristic algorithm for global optimization. J. Supercomput. 2023, 79, 7305–7336. [Google Scholar] [CrossRef]
  42. Yang, T.; Sun, X.; Yang, H.; Liu, Y.; Zhao, H.; Dong, Z.; Mu, S. Integrated thermal error modeling and compensation of machine tool feed system using subtraction-average-based optimizer-based CNN-GRU neural network. Int. J. Adv. Manuf. Technol. 2024, 131, 6075–6089. [Google Scholar] [CrossRef]
  43. Abdollahzadeh, B.; Soleimanian Gharehchopogh, F.; Mirjalili, S. Artificial gorilla troops optimizer: A new nature-inspired metaheuristic algorithm for global optimization problems. Int. J. Intell. Syst. 2021, 36, 5887–5958. [Google Scholar] [CrossRef]
  44. Zhou, X.; Ma, H.; Gu, J.; Chen, H.; Deng, W. Parameter adaptation-based ant colony optimization with dynamic hybrid mechanism. Eng. Appl. Artif. Intell. 2022, 114, 105139. [Google Scholar] [CrossRef]
  45. Aggarwal, S.; Kumar, N. Path planning techniques for unmanned aerial vehicles: A review, solutions, and challenges. Comput. Commun. 2020, 149, 270–299. [Google Scholar] [CrossRef]
  46. Puente-Castro, A.; Rivero, D.; Pazos, A.; Fernandez-Blanco, E. A review of artificial intelligence applied to path planning in UAV swarms. Neural Comput. Appl. 2022, 34, 153–170. [Google Scholar] [CrossRef]
  47. Jones, M.; Djahel, S.; Welsh, K. Path-planning for unmanned aerial vehicles with environment complexity considerations: A survey. ACM Comput. Surv. 2023, 55, 234. [Google Scholar] [CrossRef]
  48. Wang, X.; Gursoy, M.C.; Erpek, T.; Sagduyu, Y.E. Learning-based UAV path planning for data collection with integrated collision avoidance. IEEE Internet Things J. 2022, 9, 16663–16676. [Google Scholar] [CrossRef]
Figure 1. The variation curve of cosine adaptive factor ct.
Figure 1. The variation curve of cosine adaptive factor ct.
Biomimetics 09 00583 g001
Figure 2. The variation curve of nonlinear varying weight ωt.
Figure 2. The variation curve of nonlinear varying weight ωt.
Biomimetics 09 00583 g002
Figure 3. Schematic diagram of refraction reverse learning.
Figure 3. Schematic diagram of refraction reverse learning.
Biomimetics 09 00583 g003
Figure 4. The flowchart of ICSA algorithm.
Figure 4. The flowchart of ICSA algorithm.
Biomimetics 09 00583 g004
Figure 5. Convergence curves of test functions (N = 30).
Figure 5. Convergence curves of test functions (N = 30).
Biomimetics 09 00583 g005aBiomimetics 09 00583 g005bBiomimetics 09 00583 g005c
Figure 6. UAV 3D path planning (N = 100).
Figure 6. UAV 3D path planning (N = 100).
Biomimetics 09 00583 g006
Figure 7. Total cost function convergence curve.
Figure 7. Total cost function convergence curve.
Biomimetics 09 00583 g007
Table 1. Algorithm parameter settings.
Table 1. Algorithm parameter settings.
AlgorithmMain Parameter Configurations
ICSAρ = 1, r = 0.3, β = 1.5
CSAPp = 0.1, k = 3.5, ρ = 1
DBOB = 0.3, k = 0.1, S = 0.5
SABO-
GTOp = 0.03, β = 3, ω = 0.8
Table 2. Test functions.
Table 2. Test functions.
Function TypeFunction ExpressionDimension Search Range Optimal Value
Unimodal Functions F 1 x = i = 1 n x i 2 30[−100, 100]0
F 2 x = i = 1 n x i + i = 1 n x i 30[−10, 10]0
F 3 x = i = 1 n j 1 i x j 2 30[−100, 100]0
F 4 x = max i x i , 1 i n 30[−10, 10]0
F 5 x = i = 1 n 1 100 x i + 1 x i 2 2 + x i 1 2 30[−30, 30]0
F 6 x = i = 1 n x i + 0.5 2 30[−100, 100]0
Multimodal Functions F 7 x = i = 1 n x 1 2 10 cos 2 π x i + 10 30[−5.12, 5.12]0
F 8 x = 20 exp 0.2 1 n i = 1 n x i 2 exp 1 n i = 1 n cos 2 π x i + 20 + e 30[−32, 32]0
F 9 x = 1 4000 i = 1 n x i 2 i = 1 n cos x i i + 1 30[−600, 600]0
F 10 x = π n 10 sin π y 1 + i = 1 n 1 y i 1 2 1 + 10 sin 2 π y i + 1   + y n 1 2 + i = 1 n u x i , 10 , 100 , 4
y i = 1 + x i + 1 4 u x i , a , k , m = k x i a m , x i > a 0 , a < x i < a k x i a m , x i < a
30[−50, 50]0
F 11 x = 0.1 sin 2 3 π x 1 + i = 1 n x i 1 2 1 + sin 2 3 π x i + 1 · x n 1 2 1 + sin 2 2 π x n + i = 1 n u x i , 5 , 100 , 4   30[−50, 50]0
Fixed-Dimension Functions F 12 x = 1 500 + j = 1 25 1 j + i = 1 2 x i a i j 6 1 2[−65, 65]1
F 13 x = i = 1 11 a i x 1 b i 2 + b i x 2 b i 2 + b i x 3 + x 4 2 4[−5, 5]0.0003
F 14 x = x 2 5.1 4 π 2 x 1 2 + 5 π x 1 6 2 + 10 1 1 8 π cos x 1 + 10 2[−5, 5]0.398
F 15 x = 1 + x 1 + x 2 + 1 2 19 14 x 1 + 3 x 1 2 14 x 2 + 6 x 1 x 2 + 3 x 2 2     × 30 + 2 x 1 3 x 2 2 18 32 x 1 + 12 x 1 2 + 48 x 2 36 x 1 x 2 + 27 x 2 2 2[−2,2]3
F 16 x = i = 1 4 c i exp j = 1 6 a i j x j p i j 2 6[0, 1]−3.32
F 17 x = i = 1 5 X a i X a i T + c i 1 4[0, 10]−10.1532
F 18 x = i = 1 7 X a i X a i T + c i 1 4[0, 10]−10.4028
Table 3. Optimization results of test functions (N = 30).
Table 3. Optimization results of test functions (N = 30).
Test FunctionIndicatorICSACSADBOSABOGTOACO
F 1 x Avg03.92351 × 10−373.5695 × 10−211008.28236 × 10−53
Std08.74709 × 10−370002.17691 × 10−52
Best01.53927 × 10−391.9754 × 10−304003.08683 × 10−56
Worst04.05902 × 10−361.0708 × 10−209001.04369 × 10−51
Ranking142113
F 2 x Avg1.3023 × 10−2560.000751462.436 × 10−1101.4396 × 10−21909.08987 × 10−32
Std00.0039403461.3342 × 10−109009.92541 × 10−32
Best1.7608 × 10−2773.57338 × 10−162.0974 × 10−1612.3782 × 10−22306.93638 × 10−33
Worst3.9048 × 10−2550.0216033667.308 × 10−1092.1157 × 10−21805.00026 × 10−31
Ranking264315
F 3 x Avg07.92906 × 10−269.4217 × 10−1816.8492 × 10−18402.56123 × 10−08
Std03.36937 × 10−250005.22913 × 10−08
Best06.50681 × 10−311.6863 × 10−2963.6308 × 10−23701.9426 × 10−10
Worst01.80752 × 10−242.8265 × 10−1792.0547 × 10−18202.60581 × 10−07
Ranking143215
F 4 x Avg2.8539 × 10−2590.0999372362.7047 × 10−1031.9581 × 10−17004.76818 × 10−14
Std00.2947008291.4814 × 10−102001.39493 × 10−13
Best4.715 × 10−2803.21841 × 10−105.5798 × 10−1603.0862 × 10−17401.84925 × 10−15
Worst8.5334 × 10−2581.4154561268.1139 × 10−1023.1922 × 10−16907.76042 × 10−13
Ranking264315
F 5 x Avg0.01112613412.676250064.1196949937.2734225330.0059783633.593146727
Std0.02859598425.742756840.2800244480.5854059570.026227530.963212705
Best5.46113 × 10−116.3417737943.5172003596.265471372.02995 × 10−110.005562942
Worst0.120001233148.93901024.7971451538.7069966610.1424920674.118201494
Ranking264513
F 6 x Avg007.60789 × 10−300.0446418871.11302 × 10−260
Std001.96701 × 10−290.086786593.01425 × 10−260
Best003.08149 × 10−330.0002472414.93038 × 10−320
Worst008.81829 × 10−290.2673694841.22077 × 10−250
Ranking112431
F 7 x Avg06.8983798891.2271148270030.17059974
Std03.3146958644.758636727004.688056919
Best01.98991812800017.87799085
Worst016.9142903521.889074070038.05089992
Ranking132114
F 8 x Avg8.88178 × 10−168.8537898398.88178 × 10−164.20404 × 10−158.88178 × 10−164.44089 × 10−15
Std09.86513438609.01352 × 10−1600
Best8.88178 × 10−164.44089 × 10−158.88178 × 10−168.88178 × 10−168.88178 × 10−164.44089 × 10−15
Worst8.88178 × 10−1619.966106898.88178 × 10−164.44089 × 10−158.88178 × 10−164.44089 × 10−15
Ranking141213
F 9 x Avg00.155762810.015919282000.366971328
Std00.0917229320.048064635000.133866692
Best00.027093650000.160998083
Worst00.4086782610.219083352000.68081313
Ranking132114
F 10 x Avg4.71163 × 10−329.45344 × 10−074.51664 × 10−190.0283532242.75368 × 10−264.71486 × 10−32
Std1.67022 × 10−475.17786 × 10−062.39465 × 10−180.0327991161.39317 × 10−251.76746 × 10−34
Best4.71163 × 10−321.62026 × 10−264.83264 × 10−320.0007995255.00206 × 10−324.71163 × 10−32
Worst4.71163 × 10−322.83603 × 10−051.3124 × 10−170.1745226187.6356 × 10−254.80844 × 10−32
Ranking154632
F 11 x Avg1.34978 × 10−325.42618 × 10−050.0280776210.1361360580.0017174041.98295 × 10−32
Std5.5674 × 10−480.0001552330.0541851420.1338857710.00760063.468 × 10−32
Best1.34978 × 10−322.86319 × 10−291.5963 × 10−320.0019682162.95216 × 10−321.34978 × 10−32
Worst1.34978 × 10−320.000615520.1962538080.4212688510.0405347562.03448 × 10−31
Ranking135642
F 12 x Avg0.9980038380.9980038381.2946386354.3383419630.9980038381.031138073
Std04.12326 × 10−171.0093826353.60295764800.181483682
Best0.9980038380.9980038380.9980038380.9996218130.9980038380.998003838
Worst0.9980038380.9980038385.92884512513.218157880.9980038381.9920309
Ranking113412
F 13 x Avg0.0003074860.0007436350.0005551780.0008419770.0003380090.001110708
Std1.9311 × 10−190.0002543040.0002626390.0010247630.0001671810.000280454
Best0.0003074860.0003075170.0003074860.0003142970.0003074860.000827985
Worst0.0003074860.0012232780.0012481670.004828560.0012231730.001873158
Ranking143526
F 14 x Avg0.3978873580.3978873580.3978873580.4863772680.3978873580.397887358
Std02.17226 × 10−1500.15758542700
Best0.3978873580.3978873580.3978873580.3979173330.3978873580.397887358
Worst0.3978873580.3978873580.3978873581.035954750.3978873580.397887358
Ranking111211
F 15 x Avg3335.10570114133
Std1.91632 × 10−151.1057 × 10−131.91455 × 10−155.1705607662.13615 × 10−151.35253 × 10−15
Best3333.00069794333
Worst33328.9715801233
Ranking111211
F 16 x Avg−3.310105859−3.278107229−3.210739931−3.288603505−3.274437923−3.278401027
Std0.0362776890.0586721770.0910773170.0498106690.0592412180.058273385
Best−3.321995172−3.321995172−3.321995172−3.321736945−3.321995172−3.321995172
Worst−3.20310205−3.19804659−3.020699112−3.10979718−3.20310205−3.20310205
Ranking146253
F 17 x Avg−10.15319968−9.311128423−6.766641721−4.943759086−10.15319968−5.772228226
Std6.73596 × 10−151.9151177162.4353756890.3548304296.50588 × 10−153.651680149
Best−10.15319968−10.15319968−10.15319968−5.054837209−10.15319968−10.15319968
Worst−10.15319968−5.10077214−4.851276173−3.2791667−10.15319968−2.630471668
Ranking123514
F 18 x Avg−10.40294057−10.04996102−7.81624257−5.00084003−10.40294057−8.639825919
Std8.07992 × 10−161.3433176012.8471958960.5188370429.32988 × 10−162.991783299
Best−10.40294057−10.40294057−10.40294057−6.509305951−10.40294057−10.40294057
Worst−10.40294057−5.087671825−2.765897328−3.224668878−10.40294057−2.765897328
Ranking124513
Average Ranking1.173.3333.281.673.17
Final Ranking163524
Table 4. The mountain peaks information in terrain 1 and terrain 2.
Table 4. The mountain peaks information in terrain 1 and terrain 2.
TerrainMountain PeakCenter CoordinateHeight The Descending Slope of the Peaks Along the x-AxisThe Descending Slope of the Peaks Along the y-Axis
Terrain 1Peak 1[300, 300]559090
Peak 2[500, 500]60105105
Peak 3[800, 800]806060
Peak 4[250, 700]706060
Peak 5[750, 225]657575
Peak 6[550, 750]559090
Peak 7[850, 600]756060
Terrain 2Peak 1[750, 250]559090
Peak 2[750, 600]75120120
Peak 3[250, 450]807575
Peak 4[450, 800]709090
Table 5. The cylindrical obstacle threat area information in Terrain 1 and Terrain 2.
Table 5. The cylindrical obstacle threat area information in Terrain 1 and Terrain 2.
TerrainObstacle Threat AreaThe Center Coordinates of the Threat AreaThe Radius of the Threat Area
Terrain 1Threat 1[250, 700]150
Terrain 2Threat 2[500, 500]150
Table 6. Path-planning results in different terrains.
Table 6. Path-planning results in different terrains.
TerrainIndicatorICSACSADBOSABOGTOACO
Terrain 1Avg809.1455837.2897835.7821932.447809.3106827.3973
Std8.227719.439729.1549170.86129.289510.8471
Best793.578816.3347794.6715835.8318793.4277798.1731
Worst816.6946876.5999896.35911536.7116819.239855.612
Ranking154623
Terrain 2Avg819.0577824.8078845.81591053.9564824.4381831.7481
Std1.25463.873318.6131151.935214.131520.9946
Best818.1939818.2718818.4718944.4681818.1781818.5506
Worst821.6343833.3729886.51011655.4466865.5335867.3063
Ranking135624
Average Ranking144.5623.5
Final Ranking145623
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

Chen, Y.; Cao, L.; Yue, Y. Hybrid Multi-Objective Chameleon Optimization Algorithm Based on Multi-Strategy Fusion and Its Applications. Biomimetics 2024, 9, 583. https://doi.org/10.3390/biomimetics9100583

AMA Style

Chen Y, Cao L, Yue Y. Hybrid Multi-Objective Chameleon Optimization Algorithm Based on Multi-Strategy Fusion and Its Applications. Biomimetics. 2024; 9(10):583. https://doi.org/10.3390/biomimetics9100583

Chicago/Turabian Style

Chen, Yaodan, Li Cao, and Yinggao Yue. 2024. "Hybrid Multi-Objective Chameleon Optimization Algorithm Based on Multi-Strategy Fusion and Its Applications" Biomimetics 9, no. 10: 583. https://doi.org/10.3390/biomimetics9100583

APA Style

Chen, Y., Cao, L., & Yue, Y. (2024). Hybrid Multi-Objective Chameleon Optimization Algorithm Based on Multi-Strategy Fusion and Its Applications. Biomimetics, 9(10), 583. https://doi.org/10.3390/biomimetics9100583

Article Metrics

Back to TopTop