Adaptive Operator Quantum-Behaved Pigeon-Inspired Optimization Algorithm with Application to UAV Path Planning

: Path planning of unmanned aerial vehicles (UAVs) in threatening and adversarial areas is a constrained nonlinear optimal problem which takes a great amount of static and dynamic constraints into account. Quantum-behaved pigeon-inspired optimization (QPIO) has been widely applied to such nonlinear problems. However, conventional QPIO is suffering low global convergence speed and local optimum. In order to solve the above problems, an improved QPIO algorithm, adaptive operator QPIO, is proposed in this paper. Firstly, a new initialization process based on logistic mapping method is introduced to generate the initial population of the pigeon-swarm. After that, to improve the performance of the map and compass operation, the factor parameter will be adaptively updated in each iteration, which can balance the ability between global and local search. In the ﬁnal landmark operation, the gradual decreasing pigeon population-updating strategy is introduced to prevent premature convergence and local optimum. Finally, the demonstration of the proposed algorithm on UAV path planning problem is presented, and the comparison result indicates that the performance of our algorithm is better than that of particle swarm optimization (PSO), pigeon-inspired optimization (PIO), and its variants, in terms of convergence and accuracy.


Introduction
Unmanned aerial vehicles (UAVs) have been widely and more extensively applied for various practical military and civilian missions, e.g., disaster circumstances, modern warfare, and weather monitoring in recent years, which benefit from their strong viability, low cost, and excellent maneuverability [1,2].Autonomous navigation and guidance are the most important requirements of the UAV system, of which the basic component is path planning.Generally, path planning aims not only to generate a trajectory to a target in avoiding obstacles or collisions (assuming reference flight conditions and providing maps of the environment), but also to optimize a given function under kinematic and/or dynamic constraints [3].The classical UAV path planning problem can be modeled as a nonlinear optimization problem aiming to find a flyable path from a given start position to a target position [4].There does not exist an algorithm that provides an exact analytic solution to such a problem.A series of methods have been proposed and widely applied to solve the optimization problem, such as simulated annealing algorithm [5], artificial potential field algorithm [6], A* algorithm [7], RRT algorithm [8], polynomial optimization method [9], and heuristic approaches [10].
In recent decades, inspired by the organized behavior of natural biological groups, numerous swarm-intelligence optimization algorithms have been proposed to be applied to UAV path planning problem [11,12].Notable examples include ant colony optimization algorithm (ACO) [13], particle swarm optimization algorithm (PSO) [14], fruit fly optimization algorithm (FOA) [15], and pigeon-inspired optimization algorithm (PIO) [16].Various merits, including simple structure, general problem adaptability, and rapid search rate, make swam-intelligence optimization algorithms a promising tool for solving UAV path planning problems, especially under varied and complex environments.
The ACO algorithm can imitate the behavior of pheromone-mediated information transferred among ant colonies to optimize path planning problems [17].However, ACO usually incurs a large number of calculations, and it is easy to fall into the local optimum.To address these problems, direction guidance and chaotic theory were applied to ACO for efficiency improvement [18,19].By adding the chaos disturbance factor and direction guiding factor into the ACO, the direction of the search is strengthened, and the ACO can jump out of the local optimum.An improved ACO with redefined heuristic information function and pheromone updating method was applied to multi-UAV path planning to generate initial path [20].With the redefined updating method, the probability of falling into local optimum and the global search ability are improved.
Besides, fruit fly behavior-inspired algorithm, named FOA, has been successfully applied to solve optimal path planning problems.Similar to the other swarm-intelligence algorithms, the conventional FOA still has the drawback that it easily falls into a local optimum.In order to improve the search efficiency and global search ability of the conventional FOA, Zhang et al. proposed a novel angle-encoded FOA with mutation adaptation mechanism (θ-MAFOA) and adopted the novel algorithm to the UAV path planning problem [21].With the mutation adaption mechanism, the ability of exploitation and exploration of FOA are balanced.
PSO is a random search algorithm and has been widely applied to UAV path planning [22].To avoid falling into the local optimum, integrating the PSO with other intelligence algorithms becomes even more efficient.Vincent made a comparison of parallel genetic algorithm and PSO for real-time UAV path planning to avoid poor solutions [23].The generated UAV path is composed of vertical helices, line segments, and circular arcs, and a multi-objective cost function was established to evaluate the characteristic of the path by using this planning method.A hybrid differential evolution (DE) and quantum-behaved particle swarm optimization (DEQPSO) were used to plan the route for UAV on the sea [24].Optimizing the design parameters is also another effective method for performance improvement.With the hybrid PSO strategy [25], the accuracy and the convergence rate of UAV path planning can be further increased by adding a contraction factor.Asma Ayari also proposed another improved PSO algorithm that satisfies Gaussian dynamic distribution and reduces the randomness of particles [26].
PIO is a novel swarm-intelligence optimization algorithm which is inspired by the behavior of pigeon homing with the help of landmark, map, and compass [27].The PIO algorithm is composed of two operators, i.e., map and compass operator and landmark operator.Since PIO has demonstrated its fast convergence rate and adaptability in nonlinear optimizations, it has been widely applied to UAV path planning problems [28], multiple spacecraft position problem [29], and multi-UAV cooperation control [30,31].In order to further improve the convergence rate and avoid the premature convergence problem, new biological behavior characteristics were introduced into PIO.Some predatory-prey pigeon-inspired optimization (PPPIO) algorithms were proposed and applied to UAV path planning in 3D static and dynamic environment [32,33].The most recent hybrid pigeon-inspired optimization with quantum theory (QPIO) was proposed to solve the premature convergence problem and continuous optimization problems [34].However, the existing PIO algorithm and its variants did not perform initial route optimization process.PPPIO, along with QPIO, failed to take account of global search and local search synchronously.Moreover, the path generated by QPIO is not smooth, and the existing problem of population decline too fast in landmark operator.
To address the above shortcomings, we introduce some adaptive operators into conventional QPIO to improve the capability of searching appropriate solutions.The chaotic strategy is introduced in our algorithm to generate some initial solutions with wider solution space coverage.Differing from other conventional swarm-intelligence algorithms, the map and compass factor R, varying with the number of iterations, is adopted to balance global search ability and local search ability in our algorithm.At the landmark operator, a new population updating equation is proposed to increase the search diversity.Simulation results and comparisons verify the feasibility and effectiveness of our proposed algorithm.
The rest of the paper is organized as follows: Section 2 describes the preliminaries, including the problem formulation of UAV path planning and the theory of QPIO.Section 3 shows the implementation procedure of UAV path planning using our proposed adaptive operator quantum-behaved pigeon-inspired optimization (AOQPIO).Then, we compare the quality of the path produced by PIO, PSO, QPIO, and AOQPIO in Section 4. Section 5 is the conclusion.

Problem Formulation
UAV needs to plan a safe and short flying path from starting position to target position.However, numerous obstacles and no-fly zones fulfill complex environments, such as urban, indoor, or forest area.It is crucial to take the constraints of nonlinear dynamics and unstructured environment into account.Further, the fly speed and altitude changing need to be also considered.In the paper, we model the above planning problems as a mixed objective optimization problem with constraints on UAV dynamic and threats, and assume that UAVs satisfy following dynamics.
u cos θ cos ψ + v(sin φ sin θ cos ψ − cos φ sin ψ) + w(sin φ sin ψ + cos φ sin θ cos ψ) u cos θ cos ψ + v(sin φ sin θ sin ψ + cos φ sin ψ) + w(− sin φ sin ψ + cos φ sin θ sin ψ) where η = θ φ ψ are Euler angles; and Ω = [ p q r ] are angular velocities; X = [ x y z ] and V = [ u v w ] are position and linear velocities respectively in the North-East-Down inertial reference.Further, the threats in the environment are simplified as cylinders that denoted by distinct vectors, including position coordinates (x k , y k , z k ), radius r k , and corresponding threat levels t k .Threat level t k is used to indicate the impact degree of threats on the UAV path.
By taking fly time, path length, and security requirements into account, the UAV path planning problem can be established as the following optimization problem: The cost function J is consisted of J l , J a , and J t , which represent path length cost, the altitude cost and the threat cost, respectively.The path length cost and altitude cost can be directly obtained by where n is the amount of path segments; l i and h i are the length and the average altitude above the sea level of the i-th path segment, respectively.To simplify the total threat cost model, an efficient approximation is adopted to the exact solution.In our work, threat cost, which can be obtained by five points, of each edge connects two discrete points along the edge [16].
where d 0.1,k is the length between the 1/10 point and the k-th threat center.We divide the line segment between two discrete points (x i−1 , y i−1 , z i−1 ), (x i , y i , z i ) into five segments of the same length, and select five of them.The threat cost of the line segment between two discrete points is approximately equal to the total cost of the selected five points.Threat computation model is shown in Figure 1.
where n is the amount of path segments; i l and i h are the length and the average altitude above the sea level of the i-th path segment, respectively.To simplify the total threat cost model, an efficient approximation is adopted to the exact solution.In our work, threat cost, which can be obtained by five points, of each edge connects two discrete points along the edge [16].
where 0.1,k d is the length between the 1/10 point and the k-th threat center.We divide the line segment between two discrete points x y z into five segments of the same length, and select five of them.The threat cost of the line segment between two discrete points is approximately equal to the total cost of the selected five points.Threat computation model is shown in Figure 1.
The image of threat computation model.

QPIO
QPIO, which consists of two operators-the quantum behavior map and compass operator and the landmark operator, is an emerging variant of PIO algorithm.The quantum behavior is introduced to PIO to increase the local search capacity as well as the randomness of the position, so that the QPIO algorithm has the ability to improve the optimization efficiency and avoid premature convergence problem.

Map and Compass Operator with Quantum Behavior
Each pigeon has its position  in every feasible position of pigeon, where 1, 2, 3 p N , and D refers to the number of feasible solutions and their dimensions, respectively.In the D dimensional solution space, the position and speed of individuals are updated as follows: ) where ( ) i V t and ( ) i X t are the speed and position of the pigeon i at iteration t ; g X represents the best position that the whole pigeon obtained by comparing the positions of all the pigeons until iteration is a randomly generated number; ( ) refers to the map and compass factor.Importantly, the map and compass factor R is a constant value that plays an important role in convergence of the algorithm.The smaller R becomes, the larger value of

QPIO
QPIO, which consists of two operators-the quantum behavior map and compass operator and the landmark operator, is an emerging variant of PIO algorithm.The quantum behavior is introduced to PIO to increase the local search capacity as well as the randomness of the position, so that the QPIO algorithm has the ability to improve the optimization efficiency and avoid premature convergence problem.

Map and Compass Operator with Quantum Behavior
Each pigeon has its position ] in every feasible position of pigeon, where i = 1, 2, 3 • • • N p .N p , and D refers to the number of feasible solutions and their dimensions, respectively.In the D dimensional solution space, the position and speed of individuals are updated as follows: where V i (t) and X i (t) are the speed and position of the pigeon i at iteration t; X g represents the best position that the whole pigeon obtained by comparing the positions of all the pigeons until iteration t − 1; rand ∈ (0, 1) is a randomly generated number; R ∈ (0, 1) refers to the map and compass factor.Importantly, the map and compass factor R is a constant value that plays an important role in convergence of the algorithm.The smaller R becomes, the larger value of e −Rt becomes.Therefore, pigeons inherit a greater speed, which is conducive to fast convergence and better global search ability; conversely, large R would lead to only local search.
In order to improve the optimization efficiency and avoid premature convergence problem, quantum-behaved theory is added to the conventional PIO ( 7)-( 8) [35,36].
where m_best i is the mean value of best pigeon positions, and X pi (t) is the local best solution that pigeon i has obtained before the iteration t; ε is the expansion-contraction coefficient: where q is an experimentally defined parameter; N c1max refers to the maximum number of generation; ϕ and m are random numbers in the given interval [0, 1].

Landmark Operator
In each iteration loop, the landmark operator will sort fitness values of the current individual, and round off each half after discarding the poor quality individuals.The position X c (t) of remaining pigeon center is used as a landmark and the reference direction of the flight.The update formula is as follows: In the Equations ( 13) and ( 14), N p is the population of pigeons and f itness(X i (t)) is the quality of the pigeon individual which is constructed by the cost function according to Equations ( 3)- (6).
Remark.For the minimum optimization problems such as the path planning problem, we can choose f itness(X i (t)) = 1  J+eps , where J is the optimization function described in Equation ( 3) and eps refers to the relative accuracy of the float point.

Method
In this section, the population initialization process of QPIO is optimized by logistic mapping and the map and compass factor R is modified to balance the global and local search ability.In the landmark operator, new population updating and pigeon position X i updating strategies are given.

Initialization Process
QPIO presents the best convergence rate among all PIO algorithms, but it is still suffering from a major disadvantage that is vulnerable to fall into local optimum [37].Specifically, when the initialization is away from the global optimum, this error will guide other pigeons to search in wrong spaces and even leading, eventually, to failure.Since most initialization of pigeon population is achieved by random distribution, their performance quality could not be ensured in conventional methods.Therefore, the chaotic strategy is introduced in our algorithm to generate some initial solutions that can obtain wider solution space coverage.
Chaos is a kind of seemingly random motion appearing in the decisive dynamic system.Chaotic motion has three distinct characteristics: (1) high sensitivity to initial values; (2) ergodicity of motion trajectories; (3) randomness.These properties of chaos can be used to initialize the pigeon with chaotic variables.Logistic mapping is a common method used in this paper to produce chaotic variables, and avoid premature convergence [38].In the initialization process, the logistic mapping strategy is used to make the pigeon position random and traversable, so that it has the opportunity to jump out of the local optimum and find the global optimum, thus effectively avoiding premature convergence.The D dimensions vector u 1 = (u 1,1 , u 1,2 . . .u 1,D ) ∈ [0, 1] D is generated randomly.Then, a chaotic sequence U with N vectors u 1 , u 2 , . . .u N is obtained according to the recurrence Equation ( 16): After the chaos process, the chaos sequence U is converted into the problem solution space [X min , X max ] D .The corresponding individual initialization formula can be converted as Initialization of pigeon by chaotic sequence is based on the generation of initial values, from which the initial population can be determined by selecting the best one.Although it does not change the stochastic nature of AOQPIO initialization, the chaotic process improves the diversity of population and the ergodicity of pigeon search by applying chaotic theory.

Adaptive Map and Compass Factor Strategy
The selection of parameters significantly has an impact on the result.Conventional methods of choosing appropriate parameters can be roughly divided into two categories [39]: parameter adjustment and parameter control.Parameter adjustment is used for parameters preset before the beginning, and the parameters that are kept constant throughout the execution of the algorithm.Parameter control means that parameters run with the initial value at the beginning of the algorithm.Then, they can be adaptively updated according to the operation of the algorithm.
QPIO is simple to be implemented as its parameter number is limited.However, it still poses a great challenge to adjusting the parameters in QPIO.As mentioned in a previous section, the turning of map and compass factor R can adjust the algorithm on better global or local searching capability.If R is set to some constant values, the algorithm is not capable of coordinating the global and local search synchronously during the optimization procedure.Thus, we aim to find a function with the following ideal characteristics.The initial value of R is small with gradual changes; as the independent variable increases, it rapidly reaches the preset values [40].
Apparently, R(τ) belongs to [0, 1/a]; τ ∈ [−10, 10] is determined by the number of iteration and its maximum; N c is the number of iteration; parameter b can adjust the rising speed.Specifically, the value of R is coincident with conventional QPIO algorithm when the value of a is set to 1.By using the dynamic parameters tuning strategy ( 18) and ( 19), R is small at the beginning and it gradually increases after fully global search.The local search is better completed accordingly.With the consideration on the value range of R and the balance of global-local search, we set a = 1 and b = 100 in this paper.

Adaptive Compression Factor Strategy
The number of pigeons decreases by half after each iteration in a conventional QPIO algorithm, according to Equation (14).However, if the pigeon number decreases too fast, only one pigeon can survive after a small amount of iterations eventually, which prevents globally searching for the algorithm.The optimization performance is reduced in the later searching stage.Besides, the landmark operator in the conventional QPIO do not consider the waypoints distribution, which will influence the smoothness of the planning path.
To address the above problem, we propose a new pigeon number updating strategy in the landmark operator.
where w ∈ [0, 1] is constant and N dec is a constant parameter in the initial.The new formula of pigeon position updating in the landmark operator is shown, as below.
where s ∈ (0, 1) is a constant and N c2max is the maximum number of generation after the landmark operation finished.

Procedures of AOQPIO
Differing from other conventional PIO algorithms, the process of the proposed method in this paper has mainly three improvements.Firstly, the chaotic strategy is introduced in our algorithm to generate some initial solutions with wider solution space coverage.Further, the map and compass factor R, varying with the number of iterations, is adopted to balance global search ability and local search ability.Last but not least, the new population updating equation is proposed to increase the search diversity at the landmark operator period.Overall, the complete procedure of AOQPIO algorithm can be presented as the following Algorithm 1: Set initial values for D, N p , N c1max , N c2max , q, w, N dec ; 2 Generate the chaos sequence U according to Equation ( 20); X g is output as the global optimal of the fitness function

AOQPIO-Based UAV Path Planning
In order to further illustrate AOQPIO, in this paper, to solve UAV path planning problem, the detailed steps of AOQPIO-based UAV path planning are described in this section.If the user-defined number of control points of the UAV path is D, then the path planner should determine D1 = 3 * D coordinate parameters, namely X 1 , X 2 , • • • , X D1 , where (X i , X D+i , X 2D+i ), i = 1, 2, • • • , D denote the three dimensional coordinates in solution space.The point (x i , x D+i , x 2D+i ) in solution space can be transformed into the point (x(i), y(i), z(i)) in the real flying space by the following transform equation [41]: where the (x start , y start , z start ) represents the coordinate of the given start point in the flying space, and the ξ is the angle that the x-axis in solution space counterclockwise rotate to parallel segment start and goal point.The implementation procedure of AOQPIO algorithm applied to UAV path planning can be described as follows: Step 1: Environmental modeling, initialize terrain and threat information, including the center coordinates of threats, radius, and threat level of the threat; Step 2: According the environment model, build the path planning optimization function on the basis of Equations ( 3)-( 6), and initialize the detailed information about the path planning task; Step 3: The population information and algorithm parameters are initialized, including the population size N p , the number of waypoints D, the solution dimension D1, the operation operator parameter, and the number of iteration N c1max and N c2max for two operators, and N c2max > N c1max ; Step 4: Initialize individual velocity and position information by chaos sequence method according to Equations ( 16) and (17).Compare the fitness value f itness(), which is defined based on the optimization function, of each pigeon and find the current best position of pigeons; Step 5: Operate the map and compass operator.Firstly, we update the parameter R according to Equations ( 18) and (19).Next, we update the velocity and position of each pigeon using quantum-behaved method.Then, we compare all the pigeons' fitness value and update the new global best position X p ; Step 6: If the number of iterations N c > N cmax1 , the iteration switches from the map compass operator to the landmark operator, otherwise, it returns to the step 5; Step 7: Operate the landmark operator.According to the adaptation value, update the number of pigeons by Equation (20), then calculate the center of the pigeon according to Equation ( 14) and adjust the position of each pigeon to fly to the center of the pigeon with adaptive compression factor strategy according to the Equations ( 21) and (22); Step 8: If the number of iterations N c > N cmax2 , the iteration terminates, and the result is output, otherwise, it returns to the step 7.
Step 9: Transform the best solution result into the waypoint in real flying space according to Equation (24).
To illustrate this procedure further, the pseudocode of this procedure is also given as the following Algorithm 2.

Results and Discussion
To validate the efficiency of AOQPIO algorithm on UAV path planning problem, the optimization results in 3D environment are presented.The 3D environment containing a complex undulating terrain as shown in Figure 2, can also be found in [42]: where α, β are constant parameters, and z represents the altitude of a certain point.

Results and Discussion
To validate the efficiency of AOQPIO algorithm on UAV path planning problem, the optimization results in 3D environment are presented.The 3D environment containing a complex undulating terrain as shown in Figure 2, can also be found in [42]: ( , ) sin( / 5+1) sin( / 5) cos * sin( * ) z x y x y x y x y where , α β are constant parameters, and z represents the altitude of a certain point.Compared with existing swarm-intelligence algorithms, i.e., PIO, QPIO, PSO, we utilized path distance, cost, and altitude variations to evaluate the performance of above algorithms.All of the experiments were conducted by MATLAB 2014a on a PC with 2.4GHz CPU.The corresponding parameters used in the algorithms are presented in following Table 1.Compared with existing swarm-intelligence algorithms, i.e., PIO, QPIO, PSO, we utilized path distance, cost, and altitude variations to evaluate the performance of above algorithms.All of the experiments were conducted by MATLAB 2014a on a PC with 2.4GHz CPU.The corresponding parameters used in the algorithms are presented in following Table 1.
To verify the generalization ability of the proposed algorithm in different environments, experiments, which were utilized by the mentioned algorithms, were carried out in three maps with different threat sources.The corresponding parameters of the threat sources in the three maps are provided as supplementary in the enclosure.The start coordinates of the planning mission were (18,20,2) and the target coordinates were (90,105,2).
Apparently, the planning results calculated by AOQPIO, indicated with a green line in Figures 3a, 4a and 5a, are smoother than the results of other algorithms.Convergence curves of the cost values in the three scenarios are displayed in Figures 6-8.According to the altitude value in Figures 3b, 4b and 5b, most of the waypoints generated by our algorithm are lower than those of other algorithms.Further, according to the convergence results of cost function value, our algorithm has the ability of converging more stably and quickly when compared with other three algorithms.
algorithms.Further, according to the convergence results of cost function value, our algorithm has the ability of converging more stably and quickly when compared with other three algorithms.algorithms.Further, according to the convergence results of cost function value, our algorithm has the ability of converging more stably and quickly when compared with other three algorithms.The curve of the cost iteration in map 1 shows that our algorithm convergences at the eighth iteration and the convergence value is 73.8954, which is lower than the convergence value of other algorithms.Similar experiment results obtained on the above three test scenarios verify that the performance of AOQPIO is independent of changed battlefield environments that mainly mean different enemy threats.Clearly, the modifications provided in this paper had a positive effect on AOQPIO.
To avoid the influence of random initial values, another 40 experiments in each map were carried out and the average values of path length, cost, and searching time were calculated.The statistical results are listed in Table 2. From the results in map 1, the path generated by AOQPIO is 5.9% lower than those of other algorithms, and cost decreases by at least 0.4% separately.Similar to the experiment carried out in map 1, our algorithm also has the best efficiency and can generate the shortest path in the experiments carried out in map 2 and map 3. Comprehensively comparing the statistical data of all algorithms on the three maps, AOQPIO still shows its superiority to other algorithms in terms of searching ability and efficiency.From the experiment results showed above, we can see that AOQPIO is significantly superior to other swarm-intelligence algorithms tested in this paper.The chaotic strategy indeed can improve the global search ability of the basic QPIO.Furthermore, the adaptive operator strategies help the AOQPIO to search more fully and carefully.With the adaptive compression factor, a smoother path can be generated by AOQPIO.In a word, the proposed AOQPIO has the advantage over the conventional QPIO, as well as PSO and PIO algorithms in terms of searching ability, stability, and robustness.

Conclusions
The AOQPIO algorithm for UAV path planning in 3D undulating terrain was proposed in this paper.Various improvements have been incorporated into conventional QPIO algorithms to ensure the convergence rate and avoid local optimum.By utilizing logistic mapping method, the introduced pigeon population initialization procedure improved the diversity of population and the ergodicity of

Figure 1 .
Figure 1.The image of threat computation model.

Figure 2 .
Figure 2. The image of the 3D undulating terrain environment.
generalization ability of the proposed algorithm in different environments,

Figure 2 .
Figure 2. The image of the 3D undulating terrain environment.

Figure 3 .Figure 3 .
Figure 3. (a) is the result of all algorithms and (b) is the result of optimal unmanned aerial vehicle (UAV) altitude of the experiment carried out in map 2.

Figure 3 .
Figure 3. (a) is the result of all algorithms and (b) is the result of optimal unmanned aerial vehicle (UAV) altitude of the experiment carried out in map 2.

Figure 4 .
Figure 4. (a) is the result of all algorithms and (b) is the result of optimal UAV altitude of the experiment carried out in map 2.

Figure 4 .
Figure 4. (a) is the result of all algorithms and (b) is the result of optimal UAV altitude of the experiment carried out in map 2.

Figure 4 .Figure 5 .
Figure 4. (a) is the result of all algorithms and (b) is the result of optimal UAV altitude of the experiment carried out in map 2.

Figure 5 . 16 Figure 6 .
Figure 5. (a) is the result of the experiment and (b) is the result of optimal UAV altitude of the experiment carried out in map 3. Algorithms 2019, 12, 3 13 of 16

Figure 7 .Figure 6 . 16 Figure 6 .
Figure 7.The comparative curves of iteration of the algorithms of the experiment carried out in map 2.

Figure 7 .Figure 7 .
Figure 7.The comparative curves of iteration of the algorithms of the experiment carried out in map 2.

Figure 7 .
Figure 7.The comparative curves of iteration of the algorithms of the experiment carried out in map 2.

Figure 8 .
Figure 8.The comparative curves of iteration of the algorithms of the experiment carried out in map 3.

Figure 8 .
Figure 8.The comparative curves of iteration of the algorithms of the experiment carried out in map 3.
Calculate fitness values of different pigeon individuals, SetX pi = X i , X g = argmin[calculate_ f itness(X pi )]; i , and update X pi and X g according to line (16)-(20);

Table 1 .
The parameters used in experiments of different algorithms.

Table 1 .
The parameters used in experiments of different algorithms.

Table 2 .
The statistical results of path length, cost, and running time of the algorithms in the three maps.