# A Hybrid Butterfly Optimization Algorithm for Numerical Optimization Problems

Authors: Huan Zhou; Hao-Yu Cheng; Zheng-Lei Wei; Xin Zhao; An-Di Tang; Lei Xie
Journal: Computational Intelligence and Neuroscience
 (2021)

---

## Abstract

The butterfly optimization algorithm (BOA) is a swarm-based metaheuristic algorithm inspired by the foraging behaviour and information sharing of butterflies. BOA has been applied to various fields of optimization problems due to its performance. However, BOA also suffers from drawbacks such as diminished population diversity and the tendency to get trapped in local optimum. In this paper, a hybrid butterfly optimization algorithm based on a Gaussian distribution estimation strategy, called GDEBOA, is proposed. A Gaussian distribution estimation strategy is used to sample dominant population information and thus modify the evolutionary direction of butterfly populations, improving the exploitation and exploration capabilities of the algorithm. To evaluate the superiority of the proposed algorithm, GDEBOA was compared with six state-of-the-art algorithms in CEC2017. In addition, GDEBOA was employed to solve the UAV path planning problem. The simulation results show that GDEBOA is highly competitive.


---

## Body

## 1. Introduction

Optimization problems exist in all aspects of our society, including business, engineering, and science. An optimization problem is the process of finding the best value of decision variables that satisfy the maximum or minimum objective value without violating the constraints. With the development of science and technology in these days, the optimization problems we encounter have become increasingly complex. These real-world optimization problems often involve many decision variables, complex nonlinear constraints and nonconvexity, dynamic objective functions, and expensive computational costs [1, 2]. Therefore, when we solve these problems using traditional gradient-based methods, we encounter many difficulties in achieving a satisfactory solution [3]. As the field of optimization has developed, metaheuristic algorithms have become increasingly popular. It has the property of achieving an optimal or near-optimal solution in a reasonable time and does not rely on problem-specific gradient information [4]. It is therefore widely used to solve various types of optimization problems, such as task planning [5–7], feature selection [8, 9], parameter optimization [10, 11], and image segmentation [12, 13].In the last decades, many metaheuristic algorithms have been proposed and successfully applied to different domains. These algorithms can be divided into three categories: evolutionary-based algorithms, physical-based algorithms, and swarm-based algorithms. Evolution-based algorithms are a class of algorithms that simulate the laws of evolution in nature. The genetic algorithm (GA) [14] is a widely used evolutionary-based algorithm proposed by Holland. It updates populations by simulating the natural law of survival of the fittest. With the popularity of GA and GA variants, more and more evolutionary-based algorithms are continuously being proposed, including differential evolution (DE) [15], genetic programming (GP) [16], evolutionary strategies (ES) [17], and evolutionary programming (EP) [18]. Besides these evolutionary algorithms, some novel evolutionary-based algorithms have been proposed recently, such as artificial algae algorithm (AAA) [19], biogeography-based optimization (BBO) [20], and monkey king evolutionary (MKE) [21]. Physical-based algorithms simulate the laws of physics in nature or in the universe. Simulated annealing (SA) [22] inspired by annealing phenomena in metallurgy is one of the best-known physical-based algorithms. There are other physical-based algorithms proposed, including gravitational search algorithm (GSA) [23], nuclear reaction optimization (NRO) [24], water cycle algorithm (WCA) [25], and sine cosine algorithm (SCA) [26]. Swarm-based algorithms simulate social behaviour such as self-organization and division of labour in species. The particle swarm optimization (PSO) [27]and ant colony optimization (ACO) [28] are two classic swarm-based algorithms. Inspired by these two algorithms, more and more scholars are conducting research on this subject and proposing different swarm-based algorithms. The popularity of PSO and ACO has prompted more researchers to propose new swarm-based metaheuristics. Mirjalili et al. proposed grey wolf optimizer based on the collaborative foraging of grey wolves [29]. Wang et al. proposed monarch butterfly optimization inspired by the migratory activity of monarch butterflies [30]. Inspired by the spiral foraging and parabolic foraging of tuna, Xie et al. proposed the tuna swarm optimization [31]. In addition to the above three types of algorithms, a class of human-based metaheuristics is beginning to emerge. These algorithms are inspired by the characteristics of human activity. Teaching-learning-based optimization (TLBO) [32], inspired by traditional teaching methods, is a typical example of this category of metaheuristic algorithms. Other human-based metaheuristics include social evolution and learning optimization (SELO) [33], heap-based optimizers (HBO) [34], political optimizers (PO) [35], and many others. The butterfly optimization algorithm (BOA) is a swarm-based metaheuristic algorithm proposed by Arora and Singh [36]. BOA establishes an exploitation and exploration process based on the foraging behaviour and information-sharing strategies of butterflies. Although the BOA can perform exploitation and exploration operations, the basic BOA suffers from diminished population diversity and a tendency to fall into local optimum. Meanwhile, the No free lunch theory (NFL) [37] states that no single algorithm can solve all optimization problems perfectly. These factors encourage us to further enhance and improve the performance of BOA.Metaheuristics have a common property: they find optimal solutions by exploiting and exploring the search space. Exploitation dominates and will weaken exploration. But when exploration is enhanced, exploitation is weakened. So, we need to improve algorithm performance by balancing exploitation and exploration. The improvements to the algorithm focus on three main areas. The first one is the optimization of the algorithm’s parameter settings. Fan et al. [38] adjusted the fragrance factor of BOA and proposed an adaptive fragrance factor update method to enhance the convergence of BOA. Tang et al. [39] proposed the use of chaotic mapping operators to replace the alert value of the sparrow search algorithm for the purpose of balancing exploitation and exploration. Fan et al. [32] presented a new nonlinear step-factor control parameter strategy to further enhance the global search capability of the marine predator algorithm. The second is to use some techniques from other fields to improve performance. The fractional order is an effective tool that has been used in other areas [33, 40]. Yousri et al. [34] proposed an enhanced Harris hawk optimization based on fractional-order calculus memory. Elaziz et al. [35] improved the initial population of the Harris hawk optimizer using fractional-order Gaussians and 2 × mod1 chaotic mappings. The third is to use other operators to improve the original algorithm. Wang et al. [41] proposed a hybrid metaheuristic algorithm for butterfly and flower pollination based on a reciprocal mechanism. Houssein et al. [42] proposed a variant of the slime mould algorithm with hybrid adaptive guided differential evolution in order to overcome the disadvantages of unbalanced exploitation and exploration. Inspired by these hybrid variants, this paper proposes a BOA variant with hybrid distribution estimation strategy, GDEBOA. GDEBOA uses a Gaussian probability model to describe the distribution of dominant populations and to guide the direction of evolution, improving the performance of the basic BOA. The performance of GDEBOA was evaluated on the CEC2017 test suite and compared with seven state-of-the-art algorithms. The superiority of the proposed algorithm is verified by numerical analysis, convergence analysis, stability analysis, and statistical analysis. In addition, GDEBOA is applied to the UAV route planning problem to further validate the algorithm’s ability to solve real-world optimization problems.The remainder of this paper is organized as follows. A review of the basic BOA is presented in Section2. Section 3 provides a detailed description of the proposed GDEBOA. In Section 4, the effectiveness of the proposed improvement strategy is verified using CEC 2017 test suite. Furthermore, GDEBOA is applied to solve the UAV route planning problem in Section 5. Finally, we summarize this work in Section 6 and offer directions for future research.


## 2. Butterfly Optimization Algorithm

The butterfly optimization algorithm is a swarm-based metaheuristic algorithm proposed by Arora et al. The algorithm builds a model of butterfly foraging and mating. BOA proposes three hypotheses: (1) all butterflies emit fragrance and are attracted to each other; (2) each butterfly moves randomly or towards the butterfly with the most scent; and (3) the stimulus intensity of the butterfly is determined by the landscape of the fitness function. As the butterflies move, the fragrance changes with them. All butterflies form a fragrance network, and when they do not feel the fragrance network, they fly randomly, which is called the global search phase. As the butterflies approach the butterfly with the greatest concentration of fragrance, this phase is called the local search phase. BOA solves the optimization problem through global and local search with the following mathematical model.The fragrance of the butterfly is expressed as a function of the physical intensity of the stimulus, described as follows:(1)fi=cIa,i=1,2,…,NP.where fi represents the butterfly fragrance, c represents the sensory modality, I represents the stimulus intensity, a is a power exponent with a value from 0 to 1, and NP denotes the number of butterflies. The mathematical model of the global and local search phases of BOA is represented as follows:(2)Xit+1=Xit+r2×Xbestt−Xit×fi,(3)Xit+1=Xit+r2×Xjt−Xkt×fi,where Xit denotes the position of the ith butterfly in the tth iteration, Xbestt denotes the global optimal individual, r∈0,1 is a random number, and Xjt and Xkt are the jth individual and the kth individual selected randomly. BOA constantly executes two search strategies during the search process. Therefore, a switching probability p is introduced to control the switching of the two strategies. The pseudocode for BOA is given in Algorithm 1.Algorithm 1: Butterfly optimization algorithm.
Generate initial population ofNP butterfliesInitialize parametersc, a, pWhilet < tmaxFor each butterfly doCalculate the fragrance using equation (1)End forFind the best butterflyFor each butterflyXi doGenerate a random numberr from [0,1]Ifp < rCalculate the positionXinew using equation (2)ElseCalculate the positionXinew using equation (3)End ifCalculate the fitness and select the better one from [Xinew, Xi]End forUpdate parametercEnd whileOutput the best position and fitness


## 3. Proposed GDEBOA

To overcome the shortcomings of the basic butterfly optimization algorithm, a modified butterfly optimization algorithm, called GDEBOA, is proposed in this paper. The combination of the Gaussian distribution estimation algorithm and BOA provides a solution to the problem of unbalanced exploitation and exploration capabilities that exist in BOA. Here, this paper employs a Gaussian distribution estimation strategy as an alternative to the global search strategy in BOA. The GDE is used to sample the dominant population and guide the evolutionary direction of the algorithm while increasing population diversity. The improved strategies and GDEBOA are described in detail below.


### 3.1. Gaussian Distribution Estimation

The Gaussian distribution estimation strategy represents inter-individual relationships through a probabilistic model. The strategy uses the current dominant population to calculate the probability distribution model and generates new offspring populations based on the probability distribution model sampling, eventually obtaining the optimal solution by continuous iteration. In this paper, the distribution model is estimated using a weighted maximum likelihood estimation method, and the top one-half population that performs better is taken as the dominant population. The mathematical model of this strategy is described as follows:(4)Xit+1=mean+y,y∼N0,Cov,(5)mean=Xbest+Xmeant+Xit3,(6)Covi=1NP/2∑i=1NP/2Xit+1−Xmeant×Xit−XmeantT,(7)Xmeant=∑i=1NP/2ωi×Xit,(8)ωi=ln0.5NP+0.5−lni∑i=1NP/2ln0.5NP+0.5−lni,where Xmeant denotes the weighted position of the dominant population and ω denotes the weight coefficient in the dominant population in descending order of fitness values. Cov is the weighted covariance matrix of the dominant populations. The pseudocode and flowchart of the proposed GDEBOA are shown in Algorithm 2 and Figure 1.Figure 1
Flowchart of GEDBOA.Algorithm 2: GDEBOA.
Generate initial population ofNP butterfliesInitialize parametersc, a, pWhilet < tmaxCalculate themean and Cov using equations (5) and (6)For each butterfly doCalculate the fragrance using equation (1)End forFind the best butterflyFor each butterflyXi doGenerate a random numberr from [0,1]Ifp < rCalculate the positionXinew using equation (4)ElseCalculate the positionXinew using equation (3)End ifCalculate the fitness and select the better one from [Xinew, Xi]End forUpdate parametercEnd whileOutput the best position and fitness


## 3.1. Gaussian Distribution Estimation

The Gaussian distribution estimation strategy represents inter-individual relationships through a probabilistic model. The strategy uses the current dominant population to calculate the probability distribution model and generates new offspring populations based on the probability distribution model sampling, eventually obtaining the optimal solution by continuous iteration. In this paper, the distribution model is estimated using a weighted maximum likelihood estimation method, and the top one-half population that performs better is taken as the dominant population. The mathematical model of this strategy is described as follows:(4)Xit+1=mean+y,y∼N0,Cov,(5)mean=Xbest+Xmeant+Xit3,(6)Covi=1NP/2∑i=1NP/2Xit+1−Xmeant×Xit−XmeantT,(7)Xmeant=∑i=1NP/2ωi×Xit,(8)ωi=ln0.5NP+0.5−lni∑i=1NP/2ln0.5NP+0.5−lni,where Xmeant denotes the weighted position of the dominant population and ω denotes the weight coefficient in the dominant population in descending order of fitness values. Cov is the weighted covariance matrix of the dominant populations. The pseudocode and flowchart of the proposed GDEBOA are shown in Algorithm 2 and Figure 1.Figure 1
Flowchart of GEDBOA.Algorithm 2: GDEBOA.
Generate initial population ofNP butterfliesInitialize parametersc, a, pWhilet < tmaxCalculate themean and Cov using equations (5) and (6)For each butterfly doCalculate the fragrance using equation (1)End forFind the best butterflyFor each butterflyXi doGenerate a random numberr from [0,1]Ifp < rCalculate the positionXinew using equation (4)ElseCalculate the positionXinew using equation (3)End ifCalculate the fitness and select the better one from [Xinew, Xi]End forUpdate parametercEnd whileOutput the best position and fitness


## 4. Numerical Experiment and Analysis

In order to fully validate the superior performance of the proposed GDEBOA, the algorithm was tested using the IEEE CEC2017 single-objective test function. The CEC2017 test suite consists of 28 test functions. F1 is a unimodal function with only one global optimum solution and is used to verify the local search capability of the algorithm. F2–F8 are multimodal functions and are primarily used to test the ability of an algorithm to get outside of a local optimum. F9–F17 and F18–F28 are hybrid and composite functions, respectively, and can be used to test the potential of an algorithm to solve complex optimization problems in the real world. The definition of functions and optima is given in Table1.Table 1
Descriptions of CEC 2017 test suite.
TypeNo.DescriptionFi∗Unimodal functions1Shifted and rotated bent Cigar function300Unimodal functions2Shifted and rotated Rosenbrock’s function4003Shifted and rotated Rastrigin’s function5004Shifted and rotated expanded Scaffer’s F6 function6005Shifted and rotated Lunacek bi-Rastrigin function7006Shifted and rotated noncontinuous Rastrigin’s function8007Shifted and rotated Levy function9008Shifted and rotated Schwefel’s function1000Hybrid functions9Hybrid function 1 (N = 3)110010Hybrid function 2 (N = 3)120011Hybrid function 3 (N = 3)130012Hybrid function 4 (N = 4)140013Hybrid function 5 (N = 4)150014Hybrid function 6 (N = 4)160015Hybrid function 6 (N = 5)170016Hybrid function 6 (N = 5)180017Hybrid function 6 (N = 5)190018Hybrid function 6 (N = 6)2000Composite functions19Composition function 1 (N = 3)210020Composition function 2 (N = 3)220021Composition function 3 (N = 4)230022Composition function 4 (N = 4)240023Composition function 5 (N = 5)250024Composition function 6 (N = 5)260025Composition function 7 (N = 6)270026Composition function 8 (N = 6)280027Composition function 9 (N = 3)290028Composition function 10 (N = 3)3000Seven state-of-the-art metaheuristics were used for comparison with GDEBOA, including artificial ecosystem-based optimization (AEO) [43], grey wolf optimizer (GWO) [29], Harris hawks optimization (HHO) [44], arithmetic optimization algorithm (AOA) [45], slime mould algorithm (SMA) [46], Manta ray foraging optimization (MRFO) [47], and pathfinder algorithm (PFA) [48]. In the CEC2017 test, the maximum number of iterations is 600, and the population size is 500. All the algorithm parameters were set to be the same as in the original literature, as shown in Table 2. All algorithms were run 51 times independently, and the experimental results were recorded as shown in Table 3. In this paper, the performance of GDEBOA was evaluated comprehensively by numerical analysis, convergence analysis, stability analysis, Wilcoxon test, and Friedman test. The experiments in this paper were conducted on a computer with an AMD R7 4800U processor and 16 GB RAM. Programming was performed using MATLAB R2016b.Table 2
Algorithms used for comparative analysis and their parameter settings.
AlgorithmParametersAEONo parametersGWOa=2 (linearly decreased over iterations)HHOβ=1.5,E0∈−1,1AOAMopmax=1,Mopmin=0.2,C=1,α=5,Mu=0.499SMAz=0.03MRFOS=2PFAu1=−1+2rand,u2=−1+2randTable 3
Statistical results of the nine comparison algorithms for CEC2017.
No.IndexAEOGWOHHOAOASMAMRFOPFABOAGDEBOAF1Mean9.01E − 041.74E + 041.68E + 036.91E + 044.94E − 016.61E + 004.66E + 043.82E +047.70E − 06Std.1.89E − 037.95E + 037.95E + 021.15E + 043.18E − 014.23E + 001.22E + 046.97E + 039.81E − 07Rank265934871F2Mean6.65E + 011.35E + 021.23E + 027.61E + 038.99E + 016.70E + 019.80E + 019.33E + 035.49E + 01Std.3.72E + 012.53E + 013.33E + 012.45E + 035.12E + 004.09E + 011.77E + 011.29E + 032.08E + 01Rank276843591F3Mean1.34E + 027.12E + 012.05E + 022.95E + 028.17E + 011.46E + 021.14E + 023.49E + 021.86E + 02Std.3.17E + 012.25E + 013.62E + 013.20E + 011.99E + 014.09E + 013.11E + 012.16E + 011.18E + 01Rank417825396F4Mean1.67E + 012.28E + 005.62E + 016.21E + 017.35E − 011.05E + 011.47E + 016.63E + 014.03E − 02Std.5.63E + 001.21E + 005.92E + 006.71E + 003.21E − 011.21E + 014.99E + 005.76E + 002.51E − 01Rank637824591F5Mean2.75E + 021.31E + 024.98E + 026.00E + 021.18E + 022.12E + 021.34E + 025.57E + 022.09E + 02Std.6.77E + 014.45E + 016.57E + 015.66E + 012.40E + 016.09E + 013.12E + 013.17E + 011.14E + 01Rank627915384F6Mean1.13E + 027.52E + 011.40E + 022.25E + 029.39E + 011.27E + 029.97E + 012.93E + 021.83E + 02Std.2.33E + 013.07E + 012.13E + 012.67E + 012.03E + 013.29E + 012.66E + 011.54E + 011.11E + 01Rank416825397F7Mean2.49E + 032.05E + 024.69E + 034.50E + 039.95E + 021.77E + 032.28E + 026.82E + 035.35E − 02Std.9.16E + 021.62E + 028.28E + 027.24E + 021.22E + 031.02E + 031.83E + 028.69E + 021.74E − 01Rank628745391F8Mean3.37E + 032.64E + 034.35E + 035.51E + 033.04E + 033.41E + 034.98E + 037.33E + 036.94E + 03Std.5.42E + 025.39E + 027.25E + 025.83E + 025.06E + 026.00E + 029.01E + 022.85E + 022.87E + 02Rank315724698F9Mean1.09E + 021.96E + 021.61E + 021.72E + 031.16E + 029.20E + 011.91E + 022.19E + 037.09E + 01Std.3.10E + 014.04E + 014.86E + 019.74E + 024.33E + 013.46E + 015.28E + 016.72E + 022.60E + 01Rank375842691F10Mean3.74E + 042.18E + 077.61E + 066.27E + 091.31E + 069.00E + 041.88E + 062.08E + 093.49E + 02Std.3.35E + 041.87E + 074.21E + 062.56E + 091.09E + 067.44E + 041.97E + 067.43E + 082.52E + 02Rank276943581F11Mean1.51E + 045.50E + 051.51E + 053.80E + 042.71E + 041.36E + 047.54E + 043.15E + 088.30E + 01Std.1.46E + 043.30E + 069.05E + 041.71E + 042.64E + 041.62E + 044.12E + 042.10E + 082.77E + 01Rank387542691F12Mean1.40E + 026.15E + 043.82E + 045.72E + 044.71E + 042.01E + 033.00E + 041.19E + 056.29E + 01Std.4.36E + 019.79E + 044.25E + 044.92E + 042.84E + 042.19E + 032.94E + 047.62E + 047.46E + 00Rank285763491F13Mean4.10E + 035.46E + 046.86E + 042.35E + 041.99E + 048.64E + 033.35E + 041.82E + 065.54E + 01Std.4.49E + 035.70E + 044.86E + 041.22E + 041.57E + 049.15E + 031.77E + 041.46E + 061.94E + 01Rank278543691F14Mean1.03E + 036.98E + 021.55E + 031.98E + 038.18E + 028.93E + 021.00E + 033.18E + 031.30E + 03Std.3.11E + 022.59E + 023.56E + 025.09E + 022.83E + 022.63E + 022.63E + 024.12E + 021.80E + 02Rank517823496F15Mean4.06E + 022.17E + 027.48E + 029.12E + 024.34E + 023.36E + 023.77E + 021.22E + 034.31E + 02Std.1.84E + 021.23E + 022.19E + 022.67E + 021.64E + 022.02E + 021.71E + 022.49E + 027.95E + 01Rank417862395F16Mean1.70E + 045.57E + 056.90E + 051.29E + 063.75E + 057.66E + 042.75E + 059.60E + 053.68E + 01Std.1.28E + 047.70E + 058.77E + 051.60E + 063.55E + 056.45E + 042.82E + 056.22E + 054.62E + 00Rank267953481F17Mean1.94E + 032.42E + 051.46E + 051.08E + 063.00E + 048.42E + 034.45E + 044.61E + 062.79E + 01Std.3.07E + 032.98E + 051.42E + 051.39E + 052.11E + 041.07E + 043.91E + 044.06E + 062.57E + 00Rank276843591F18Mean3.95E + 023.20E + 026.71E + 026.94E + 023.59E + 023.45E + 024.61E + 027.29E + 024.90E + 02Std.1.62E + 021.03E + 022.01E + 021.54E + 021.59E + 021.33E + 021.52E + 029.88E + 019.82E + 01Rank417832596F19Mean3.26E + 022.67E + 024.06E + 024.87E + 022.93E + 023.09E + 022.90E + 021.97E + 023.73E + 02Std.3.26E + 012.45E + 013.51E + 015.23E + 012.17E + 012.95E + 012.62E + 013.01E + 011.21E + 01Rank628945317F20Mean1.73E + 021.49E + 032.39E + 035.13E + 032.90E + 031.00E + 022.08E + 024.71E + 021.00E + 02Std.5.11E + 021.42E + 032.37E + 031.21E + 031.36E + 037.40E − 017.61E + 027.76E + 012.31E − 04Rank367982451F21Mean5.41E + 024.18E + 027.05E + 029.68E + 024.35E + 024.95E + 024.85E + 026.97E + 025.28E + 02Std.4.86E + 013.60E + 017.35E + 019.10E + 011.95E + 014.15E + 014.09E + 015.59E + 011.59E + 01Rank618924375F22Mean6.41E + 025.02E + 028.26E + 021.14E + 035.30E + 025.65E + 025.24E + 021.10E + 035.89E + 02Std.7.63E + 014.93E + 017.42E + 011.09E + 022.95E + 015.19E + 013.72E + 011.68E + 024.38E + 01Rank617934285F23Mean3.94E + 024.29E + 024.11E + 021.67E + 033.88E + 023.89E + 023.97E + 021.75E + 033.87E + 02Std.1.47E + 011.79E + 011.87E + 014.55E + 021.69E + 008.17E + 001.72E + 012.01E + 026.50E − 01Rank476823591F24Mean2.62E + 031.69E + 033.94E + 036.40E + 031.98E + 032.71E + 032.12E + 035.21E + 032.54E + 03Std.1.10E + 032.41E + 021.10E + 037.22E + 023.42E + 021.05E + 037.10E + 021.49E + 033.47E + 02Rank517926384F25Mean5.92E + 025.22E + 026.05E + 021.34E + 035.11E + 025.56E + 025.46E + 028.14E + 025.01E + 02Std.3.05E + 019.78E + 004.00E + 012.14E + 021.17E + 012.78E + 012.85E + 019.81E + 011.52E + 01Rank637925481F26Mean3.49E + 025.11E + 024.62E + 022.95E + 034.46E + 023.42E + 024.31E + 023.28E + 033.39E + 02Std.5.89E + 014.00E + 012.60E + 016.15E + 022.89E + 015.27E + 011.90E + 013.99E + 025.28E + 01Rank376852491F27Mean9.30E + 026.85E + 021.32E + 032.43E + 037.75E + 028.27E + 021.03E + 033.04E + 039.73E + 02Std.2.09E + 021.14E + 022.56E + 025.22E + 021.83E + 022.09E + 022.35E + 024.72E + 021.13E + 02Rank417823695F28Mean9.26E + 033.77E + 061.01E + 061.47E + 071.27E + 047.34E + 034.09E + 053.98E + 072.08E + 03Std.4.95E + 032.70E + 066.08E + 051.01E + 075.29E + 033.49E + 034.56E + 052.31E + 071.46E + 02Rank376842591The results in Table3 show that GDEBOA performs best on the unimodal function F1. Although GDEBOA does not consistently achieve the optimal solution, it provides the best solution among the nine algorithms. Compared to BOA, the performance of GDEBOA is significantly improved, indicating that the improvement strategy is effective in enhancing the exploitation capability. In the multimodal functions F2–F8, each algorithm performs differently. GDEBOA performs best on F2, F4, and F7. GWO provides the best solution on F3, F6, and F8. SMA achieves the optimal solution on F5. Significantly, GDEBOA outperforms BOA in all multimodal functions, which indicates that GDEBOA has a strong global search capability. The improved strategy can effectively enhance BOA’s ability to explore the solution space. The GDEBOA performs best in most hybrid and composite functions. Specifically, GDEBOA achieved satisfactory results in 7 of the 10 hybrid functions. The GDEBOA achieved optimal solutions on 5 of the 10 composite functions. Compared to BOA, GDEBOA only performed worse on F19. The analysis of the results for the hybrid and composite functions shows that GDEBOA has a good balance of exploitation and exploration capabilities and is able to solve complex optimization problems effectively.Convergence speed and convergence accuracy are important indicators of an algorithm’s performance. Convergence analysis provides information on how the algorithm has changed over a process of iterations. Figure2 shows the average fitness convergence curves for F1–F28 based on the results of all algorithms solving the test suite 51 times. The results in Figure 2 show that GDEBOA has a faster convergence speed and better convergence accuracy compared to other algorithms. Specifically, GDEBOA outperformed all comparison algorithms in terms of convergence accuracy and convergence speed on 16 functions. Notably, although GDEBOA did not perform best on all functions, GDEBOA outperformed BOA on the 27 out of 28 functions, indicating that the improvement strategy proposed in this paper better balances exploitation and exploration.Figure 2
Convergence curves of the GDEBOA and other algorithms for CEC2017.In addition, in order to analyse the distribution characteristics of GDEBOA’s solutions when solving the functions, box diagrams are shown in Figure3 based on the results of the nine algorithms solving the CEC2017 test set 51 times independently. For most of the test functions, the minimum, maximum, and median values obtained by GDEBOA are the same as the optimal solutions, which indicates that the solutions obtained by GDEBOA are more centrally distributed and more stable. Compared to BOA, GDEBOA solves the function with fewer bad values and a more concentrated distribution of solutions, indicating that GDEBOA has achieved a balance between exploitation and exploration.Figure 3
Boxplots analysis for CEC2017.Although the superiority of GDEBOA was demonstrated by comparisons of the mean and standard deviation, the literature [49, 50] demonstrates that these comparisons are not adequate. To further verify the differences between GDEBOA and other algorithms, the Wilcoxon signed-rank test was employed. Table 4 shows the results of the Wilcoxon signed-rank test with significance level P=0.05. The term R+ indicates the extent to which GDEBOA outperforms the comparison algorithm, while R− indicates the opposite. The symbol “+/=/−” indicates that GDEBOA performs better than, similar to, and worse than the comparison algorithm.Table 4
Statistical results of Wilcoxon signed-rank test for CEC2017.
AEOGWOHHOAOAP-valueR+R−WinP valueR+R−WinP-valueR+R−WinP valueR+R−Win5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+0.071907855471=5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+3.94E − 09351291−5.46E – 1011325−0.0015341001325+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.43E − 08124383+3.94E – 09351291−5.15E − 1013260+5.15E − 1013260+5.15E − 1001326−7.35E – 1061320−5.46E − 1011325−6.93E − 1013215+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1001326−5.15E – 1001326−5.15E − 1001326−5.15E − 1001326−3.96E − 071204122+5.15E – 1013260+1.40E − 09130917+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+7.35E − 1013206+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+6.52E − 052371089−1.05E – 09121314−0.0001531067259+1.42E − 08126858+0.276894547779=2.65E – 09281298−2.97E − 09129630+6.93E − 1013215+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+0.004923363963−5.53E – 09411285−2.11E − 061169157+9.14E − 09127650+1.27E − 08561270−5.15E – 1001326−8.61E − 071188138+5.15E − 1013260+0.205722798528=5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+0.264661782544=6.15E – 1031323−5.15E − 1013260+5.15E − 1013260+0.0003311046280+3.37E – 08741252−5.15 − 1013260+5.15E − 1013260+0.006378954372+5.15E – 1013260+6.93E − 1013215+5.15E − 1013260+0.035759887439+4.94E – 09391287−1.53E − 071223103+5.15E − 1013260+5.15E−1013260+3.03E – 08125472+5.15E − 1013260+5.15E − 1013260+0.285261777549=5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+0.205722528798=8.77E – 1091317−9.66E − 09127551+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+16/6/616/0/1226/0/227/0/1SMAMRFOPFABOAP-valueR+R−WinP valueR+R−WinP valueR+R−WinP valueR+R−Win5.15E − 1013260+5.15E−1013260+5.15E−1013260+5.15E−1013260+5.15E − 1013260+0.100931838488+5.15E−1013260+5.15E−1013260+5.15E − 1001326−3.25E – 071181208−5.46E − 1011325−5.15E − 1013260+7.74E − 09127947+8.77E – 1013179+5.15E − 1013260+5.15E − 1013260+5.15E − 1001326−0.880784679647=6.93E − 1051321−5.15E − 1013260+5.15E − 1001326−1.18E – 09141312−5.15E − 1001326−5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1001326−5.15E – 1001326−5.15E − 1001326−8.25E − 08123591+1.26E − 061180146+0.006562953373+5.46E − 1013251+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1013260+6.15E – 1013233+5.15E − 1013260+5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1013260+5.80E – 1013242+5.15E − 1013260+5.15E − 1013260+2.23E − 09251301−4.42E – 09371289−1.69E − 071051221−5.15E − 1013260+0.785752692634=0.004378359967−0.052342456870=5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+5.15E − 1013260+5.46E – 1013251+5.15E − 1013260+5.15E − 1013260+1.49E − 052011125−8.49E – 061881138−0.174098518808=5.15E − 1013260+5.15E − 1001326−1.05E – 09121314−5.15E − 1001326−5.15E − 1001326−5.15E − 1013260+1.31E – 051981128−7.32E − 09128046+5.15E − 1013260+5.15E − 1001326−3.87E – 052241102−1.25E − 07991227−6.93E − 1013215+1.86E − 08631263−0.002316338988−4.64E − 08801246−5.15E − 1013260+9.93E − 071185141+0.004509966360+2.94E − 071210116+5.15E − 1013260+1.08E − 08531273−0.055852867459=4.37E − 052271099−8.77E − 1013179+0.000441038288+6.53E – 1013224+9.87E − 10131511+5.15E − 1013260+5.15E − 1013260+0.425599748578=3.52E − 09129333+5.15E − 1013260+7.43E − 08891237−7.63E – 052411085−0.119709829497=5.15E − 1013260+5.15E − 1013260+5.15E – 1013260+5.15E − 1013260+5.15E − 1013260+16/1/1114/3/1116/3/927/0/1From Table4, GDEBOA outperforms the comparison algorithm on most functions. Numerically, GDEBOA is superior to MRFO on 14 functions. GDEBOA outperforms GWO and PFA on 16 functions. Compared to HHO and AOA, GDEBOA performs better on 27 and 28 functions. In particular, GDEBOA offers better solutions than BOA on all functions, except F9. In general, GDEBOA shows a superior performance.


## 5. UAV Route Planning

In this section, the application of GDEBOA to UAV route planning is discussed in detail. The UAV route planning problem aims to minimize the cost of carrying out the mission, which can be considered as a multi-constraint optimization problem. The route planning model is described in detail in the following section.


### 5.1. Cost Function

Considering the actual UAV flight scenario, we think that the UAV should reach the target as fast as possible while being free from threats during the flight. The cost function is described as the total of two functions and two constraints.(9)minJ=ω1⋅Costp+ω2⋅Costh+η⋅Penalty,Penalty=∑i=1nci,n=1,2,ci=0,Satisfing constraints,1,No satisfing constraints,where J is the total cost function, Costp denotes the flight distance cost, Costh denotes the height cost, ω1 and ω2 are the weighting factors of two costs that satisfy ω1+ω2=1, Penalty denotes the number of constraint violations, and η is the penalty function factor for converting the constrained optimization problem into an unconstrained optimization problem.The faster the UAV reaches the target, the better it is for the mission, so the path cost is represented by the sum of the route segments.(10)Costp=∑i=1NP−1Disi,where Disi is the length of the ith route segment.Besides, it is not beneficial for UAV to fly too high to avoid the threat. Therefore, the UAV needs to maintain low height flight. The corresponding height cost function is described as follows:(11)Costh=∑i=1NPzi,where zi is the height corresponding to the ith route point.Due to the limitation of UAV’s maneuverability, the turn angle and climb angle of the actual flight should be less than the theoretical maximum.(12)Constraint1=maxφij−φmax≤0,Constraint2=maxθij−θmax≤0,where φij and θij are the turn angle and climb angle between each track segment, respectively, and φmax and θmax are the corresponding theoretical maximum values.


### 5.2. Simulation Experiment

In this subsection, we will solve the route planning model using GDEBOA and compare the results with BOA, MRFO, and SMA. All the programs were written using MATLAB R2016b and run on a Windows 10 platform with AMD R7 4700U 16 GB. For fair comparison, the parameters of each algorithm are set according to the original literature. The maximum iteration is 300, the number of populations is 50, and each algorithm is executed 30 times independently. The location of the threat is shown in Table5. The best paths generated by MRFO, HHO, BOA, and GDEBOA are shown in Figure 4. The best values, mean, standard deviation, and success rate for 30 independent runs are shown in Table 6.Table 5
Threat source settings.
ThreatTypePosition (km)Radius (km)HeightThreat 1Rader(35, 20)132.8Threat 2Missile(35, 52)82.9Threat 3Artillery(52, 72)83.0Threat 4Missile(63, 45)10.72.9Threat 5Rader(78, 78)93.1Threat 6Artillery(87, 45)73.0Figure 4
The best path of each algorithm: (a) path in three-dimensional space; (b) path in two-dimensional space.
(a)(b)Table 6
Statistical results of the Friedman test.
AlgorithmMeanBsetWorstStd.Success (%)BOA9.25E + 032.11E + 024.03E + 041.03E + 0443.33MRFO6.59E + 032.20E + 023.03E + 048.51E + 0356.67SMA2.91E + 032.17E + 022.03E + 046.41E + 0383.33GDEBOA2.22E + 021.94E + 022.25E + 025.41E + 00100.00It is clear from analyzing the test results that all algorithms can give a safe flight route. This shows that the route planning model proposed in this paper is feasible. The path found by these algorithms is different as can be seen in Figure4. Obviously, GDEBOA can find a safe path at lower cost. Also, GDEBOA can consistently find safe flight paths in 30 tests, while there are failures in the other algorithms. This indicates that GDEBOA is superior to the other three algorithms. The average cost convergence curve is plotted according to the statistical results, as shown in Figure 5. GDEBOA has faster convergence speed and better convergence accuracy.Figure 5
Convergence curves for four algorithms.


## 5.1. Cost Function

Considering the actual UAV flight scenario, we think that the UAV should reach the target as fast as possible while being free from threats during the flight. The cost function is described as the total of two functions and two constraints.(9)minJ=ω1⋅Costp+ω2⋅Costh+η⋅Penalty,Penalty=∑i=1nci,n=1,2,ci=0,Satisfing constraints,1,No satisfing constraints,where J is the total cost function, Costp denotes the flight distance cost, Costh denotes the height cost, ω1 and ω2 are the weighting factors of two costs that satisfy ω1+ω2=1, Penalty denotes the number of constraint violations, and η is the penalty function factor for converting the constrained optimization problem into an unconstrained optimization problem.The faster the UAV reaches the target, the better it is for the mission, so the path cost is represented by the sum of the route segments.(10)Costp=∑i=1NP−1Disi,where Disi is the length of the ith route segment.Besides, it is not beneficial for UAV to fly too high to avoid the threat. Therefore, the UAV needs to maintain low height flight. The corresponding height cost function is described as follows:(11)Costh=∑i=1NPzi,where zi is the height corresponding to the ith route point.Due to the limitation of UAV’s maneuverability, the turn angle and climb angle of the actual flight should be less than the theoretical maximum.(12)Constraint1=maxφij−φmax≤0,Constraint2=maxθij−θmax≤0,where φij and θij are the turn angle and climb angle between each track segment, respectively, and φmax and θmax are the corresponding theoretical maximum values.


## 5.2. Simulation Experiment

In this subsection, we will solve the route planning model using GDEBOA and compare the results with BOA, MRFO, and SMA. All the programs were written using MATLAB R2016b and run on a Windows 10 platform with AMD R7 4700U 16 GB. For fair comparison, the parameters of each algorithm are set according to the original literature. The maximum iteration is 300, the number of populations is 50, and each algorithm is executed 30 times independently. The location of the threat is shown in Table5. The best paths generated by MRFO, HHO, BOA, and GDEBOA are shown in Figure 4. The best values, mean, standard deviation, and success rate for 30 independent runs are shown in Table 6.Table 5
Threat source settings.
ThreatTypePosition (km)Radius (km)HeightThreat 1Rader(35, 20)132.8Threat 2Missile(35, 52)82.9Threat 3Artillery(52, 72)83.0Threat 4Missile(63, 45)10.72.9Threat 5Rader(78, 78)93.1Threat 6Artillery(87, 45)73.0Figure 4
The best path of each algorithm: (a) path in three-dimensional space; (b) path in two-dimensional space.
(a)(b)Table 6
Statistical results of the Friedman test.
AlgorithmMeanBsetWorstStd.Success (%)BOA9.25E + 032.11E + 024.03E + 041.03E + 0443.33MRFO6.59E + 032.20E + 023.03E + 048.51E + 0356.67SMA2.91E + 032.17E + 022.03E + 046.41E + 0383.33GDEBOA2.22E + 021.94E + 022.25E + 025.41E + 00100.00It is clear from analyzing the test results that all algorithms can give a safe flight route. This shows that the route planning model proposed in this paper is feasible. The path found by these algorithms is different as can be seen in Figure4. Obviously, GDEBOA can find a safe path at lower cost. Also, GDEBOA can consistently find safe flight paths in 30 tests, while there are failures in the other algorithms. This indicates that GDEBOA is superior to the other three algorithms. The average cost convergence curve is plotted according to the statistical results, as shown in Figure 5. GDEBOA has faster convergence speed and better convergence accuracy.Figure 5
Convergence curves for four algorithms.


## 6. Conclusions

In this paper, we propose a variant of BOA using a distribution estimation strategy, called GDEBOA, to solve the global optimization problem. The performance of BOA is enhanced by using the distribution estimation strategy to sample the evolutionary information of the dominant population and to guide the direction of individual evolution. To evaluate the effectiveness of the improved strategy and the superiority of GDEBOA, it was validated using the CEC2017 test suite. It was also compared with seven advanced algorithms through numerical analysis, convergence analysis, stability analysis, and statistical tests. The simulation results show that GDEBOA balances exploitation and exploration and is competitive with other algorithms. In addition, GDEBOA is applied to solve the UAV route planning problem. Simulation results show that GDEBOA can stably find paths with better quality, demonstrating the ability of GDEBOA to solve real-world optimization problems. Of course, there is still potential for enhancement of GDEBOA. The calculation of the covariance matrix leads to an increase in the computational cost of the algorithm. Therefore, how to reduce the cost of computing while maintaining performance is something we need to investigate further.In future work, we plan to further apply GDEBOA to cooperative multi-UAV route planning and real-time route planning. Moreover, we plan to develop the multi-objective and binary versions of GDEBOA to solve optimization problems in other fields.



---

*Source: 7981670-2021-12-24.xml*