Next Article in Journal
Identifying Irregular Potatoes Using Hausdorff Distance and Intersection over Union
Previous Article in Journal
Artificial Neural Network Approach to Guarantee the Positioning Accuracy of Moving Robots by Using the Integration of IMU/UWB with Motion Capture System Data Fusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chaos-Enhanced Adaptive Hybrid Butterfly Particle Swarm Optimization Algorithm for Passive Target Localization

1
Faculty of Mechanical Engineering, University of Belgrade, 11000 Belgrade, Serbia
2
School of Electrical Engineering, University of Belgrade, 11000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(15), 5739; https://doi.org/10.3390/s22155739
Submission received: 24 June 2022 / Revised: 22 July 2022 / Accepted: 23 July 2022 / Published: 31 July 2022
(This article belongs to the Topic Wireless Sensor Networks)

Abstract

:
This paper considers the problem of finding the position of a passive target using noisy time difference of arrival (TDOA) measurements, obtained from multiple transmitters and a single receiver. The maximum likelihood (ML) estimator’s objective function is extremely nonlinear and non-convex, making it impossible to use traditional optimization techniques. In this regard, this paper proposes the chaos-enhanced adaptive hybrid butterfly particle swarm optimization algorithm, named CAHBPSO, as the hybridization of butterfly optimization (BOA) and particle swarm optimization (PSO) algorithms, to estimate passive target position. In the proposed algorithm, an adaptive strategy is employed to update the sensory fragrance of BOA algorithm, and chaos theory is incorporated into the inertia weight of PSO algorithm. Furthermore, an adaptive switch probability is employed to combine global and local search phases of BOA with the PSO algorithm. Additionally, the semidefinite programming is employed to convert the considered problem into a convex one. The statistical comparison on CEC2014 benchmark problems shows that the proposed algorithm provides a better performance compared to well-known algorithms. The CAHBPSO method surpasses the BOA, PSO and semidefinite programming (SDP) algorithms for a broad spectrum of noise, according to simulation findings, and achieves the Cramer–Rao lower bound (CRLB).

1. Introduction

Determining the location of a passive target based on time difference of arrival (TDOA) measurements from multiple transmitters and a single receiver is a key element in many technologies, such as radar or sonar, telecommunications, mobile communications [1,2], etc. In general, two groups of localization approaches, active and passive, may be distinguished. The active localization approach takes into account the scenario when in the localization the target is actively involved. However, in the second group, the target does not participate in the localization process and merely serves to reflect the transmitter’s signals [3]. The global positioning system (GPS) has been widely used to determine the position of an object in outdoor environments [4]. However, this localization system cannot provide satisfactory performance in indoor, underwater acoustics, and urban environments, where the satellite signals are unavailable [5,6]. Therefore, passive target localization has become widely used in various applications, as an effective alternative to the GPS and other conventional localization systems.
Hence, the localization of a passive target is considered in this paper, where the noisy TDOA measurements are employed. The range measurements are calculated from the difference in the time it takes for a signal coming from a transmitter via the target to a receiver and the time required for a signal coming directly from the transmitter to a receiver. Therefore, the unknown position of a target becomes difficult to estimate since the TDOA measurements are nonlinear. The noisy TDOA measurements in the line-of-sight (LOS) environment can be described as the normal-distributed Gaussian random variable. Since the probability distribution of measurement error is known, the maximum likelihood (ML) estimator can be employed to estimate a passive target’s unknown position [2]. The ML estimator’s objective function is very nonlinear and non-convex, making it challenging to find the global optimal solution. Therefore, several efficient optimization algorithms are proposed in the literature to solve this type of optimization problem. Among these algorithms, the semidefinite programming (SDP) method became widely applied as it can efficiently transform the considered problem to convex optimization problem [2]. The SDP method’s primary benefit is that it does not demand a starting solution to solve the considered optimization problem. Furthermore, it can be solved employing MATLAB toolbox CVX with the SeDuMi solver. Still, the SDP method shows some drawbacks, which reflects on the accuracy of the obtained solution, primarily when the large measurement noise is present [7]. This shows that the solving of the non-convex ML estimation problem is an important and significant challenge.
In this context, evolutionary algorithms (EAs) are proposed with the aim to successfully achieve the global optimum to the challenging ML estimation problem, and get beyond the aforementioned drawbacks [8,9]. Generally speaking, there are two steps to the optimization process for EAs: global exploration and local exploitation [10]. Here, the first stage is concerned with identifying the area of the global optimum, and the second stage is concerned with increasing the convergence speed and the solution accuracy. Therefore, to efficiently find the global optimal solution, it is important to maintain a balance between exploration and exploitation during the optimization process of EAs.
Numerous EAs are proposed to solve different optimization problems, such as genetic algorithm (GA) [11], particle swarm optimization (PSO) [12], butterfly optimization algorithm (BOA) [13], differential evolution (DE) [14], cuckoo search algorithm (CS) [15], artificial bee colony (ABC) [16] and firefly algorithm (FA) [17], etc. During recent years, BOA and PSO algorithms are successfully applied to find the global optimal solution of different complex optimization problems, due to their advantages, such as easy implementation, fast convergence, and robustness [18].
The BOA is a novel optimization algorithm, developed by Arora and Singh [13], which is inspired by the foraging and mating behavior of butterflies. During the food searching process, the butterflies emit a fragrance, and the intensity of this fragrance is proportional to the objective function value of the butterfly. Other butterflies in the swarm can sense the fragrance intensity, in order to determine the potential direction of a food source or mating partner. In general, during the optimization process, the BOA algorithm goes through two phases: the global search phase and the local search phase. In the first phase, the BOA algorithm goes through the search space and finds regions where potential global solutions exist. Conversely, in the local search phase, the BOA algorithm performs a fine search in the neighborhood of the current optimal solution. Generally, the optimization performance of the BOA algorithm is influenced by the suitable choice of two control parameters: the sensory fragrance and power exponent [13]. The BOA algorithm is successfully applied to solve a wide range of optimization problems in science and engineering, where it has demonstrated better results compared to other EAs [18,19]. It has been shown that the BOA algorithm has a strong exploitation ability; however, it suffers from drawbacks, such as premature convergence to local optima and weak global search ability [20]. Therefore, to overcome these drawbacks, a number of improved versions of the BOA algorithm have been proposed in the literature [21,22].
The PSO algorithm is another nature-inspired EA, proposed by Kennedy and Eberhart [23], which is based on the social behavior of flocks of birds searching for food. The search process of the PSO algorithm is based on the position and velocity vectors. In this regard, each particle changes its position with respect to its personal previous best position and the best position of the whole swarm. Due to the easy implementation and effectiveness, the PSO algorithm is successfully applied to solve a wide range of complex optimization problems [24,25]. However, in solving complex optimization problems the PSO algorithm shows some disadvantages such as slow convergence rate and the problem of premature convergence [26]. It is shown in the literature that the performance of the PSO algorithm strongly depends on the appropriate choice of three control parameters: the inertia weight (w), the cognitive acceleration coefficient ( c 1 ), and the social acceleration coefficient ( c 2 ) [27]. In order to overcome these drawbacks and improve the performance, a number of adaptive and self-adaptive versions of the PSO algorithm are developed [28,29]. Furthermore, in recent years, the integration of chaos theory with the PSO algorithm has proven to be an effective way to improve the optimization performance and avoid the problem of premature convergence to local optima.
During the recent years, with the development of nonlinear dynamics, chaos theory is widely employed to improve different aspects of optimization algorithms [30]. Chaos is the bounded dynamic behavior of nonlinear systems characterized by infinite unstable periodic motions [31]. In this way, chaos maps are employed as an evolution function representing the chaos behavior, which produces a bounded sequence of random numbers depending on the choice of the initial conditions. A large number of chaos maps are found in the literature, such as the sinusoidal map, Chebyshev map, tent map, sine map, logistic map, etc. [32]. Due to the pseudo-randomness, ergodicity, and irregularity of chaos maps, the integration of chaos maps into EAs has been proven to be an effective way to improve the optimization performance and solution quality [21,29].
Another way to overcome the drawbacks and improve the optimization performance is the hybridization of different EAs. In this regard, the BOA algorithm is successfully hybridized with the PSO algorithm, in HPSOBOA, with the aim to improve convergence during the evolution process [33]. Moreover, the combination of search ability of the FA and PSO algorithms is proposed in the HFPSO algorithm, to improve the convergence speed and solution accuracy [34]. Additionally, different hybridizations of the BOA algorithm with other EAs, such as ABC and BOA (BOA/ABC) [16], BOA with DE (HBODEA) [35] and BOA with PSO (HPSO) [36] can be found in the literature, which are successfully applied to solve different complex optimization problems. Therefore, the hybridization of the EAs is proven to enhance the speed of reaching optimal solution, so as to avoid the problem of prematurely converging to a solution and provide more precise solutions.
This paper proposes a hybridization of the BOA and PSO algorithms, called CAHBPSO, to efficiently tackle the presented problem of localization of a passive target. In the proposed hybrid algorithm, the global and local search phases of the BOA algorithm are incorporated into the velocity update equation of the PSO algorithm. In addition, instead of fixed-switch probability, an adaptive technique is proposed to dynamically switch between exploration and exploitation. To further improve the performance, an adaptive strategy is employed to update the sensory fragrance of the BOA algorithm. Moreover, the logistic chaos map is incorporated into the inertia weight parameter of the PSO algorithm, to maintain a trade-off between local and global abilities.
The localization performance is measured against the Cramer-Rao lower bound (CRLB), which offers a lower constraint on the unbiased estimator variance [37]. Therefore, the root mean square error (RMSE) performance of the new CAHBPSO method and the current SDP, BOA, and PSO algorithms are compared with the derived CRLB. The following is a summary of this paper’s significant contributions:
  • The problem of localization of a passive target is formulated using noisy TDOA measurements obtained from a set of transmitters and a single receiver, for the case of LOS conditions. Due to the highly nonlinear and non-convex nature of the ML estimation problem that has been formulated for the consideration localization problem, sophisticated optimization algorithms are proposed to address this complex optimization problem.
  • By converting the considered multimodal optimization problem to a problem with distinct single extremum, the SDP method, as a convex method, is employed to effectively address the ML estimation problem.
  • The enhanced CAHBPSO algorithm—a hybridization of the BOA and the PSO algorithms—is proposed, to precisely estimate the position of the passive target. To improve convergence and maintain population diversity the global and local search phases of the BOA algorithm are incorporated into the velocity update equation of the PSO algorithm. In addition, instead of fixed-switch probability, an adaptive parameter is employed to effectively maintain a trade-off between global and local search abilities throughout the iteration process. Furthermore, the sensory fragrance of the BOA algorithm is adaptively updated and logistic chaos map is incorporated into the expression for the inertia weight parameter of the PSO algorithm.
  • The Wilcoxon signed-rank test and Friedman test are employed for statistical performance comparison between CAHBPSO algorithm with several widely applied EAs on a set of CEC2014 problems. Analyzing the optimization performance, according to the statistical analysis’s findings the modifications and hybridization proposed in this paper successfully enhance the CAHBPSO algorithm.
  • The results of the numerical simulation demonstrate that the proposed CAHBPSO method outperforms SDP, BOA, and PSO algorithms in terms of localization performance and CRLB accuracy. Furthermore, according to the simulation findings, the CAHBPSO method performs the best when there is a high level of measurement noise and it is not sensitive to changes in network layout. In terms of computational complexity, the simulation results showed that the proposed algorithm provides a proper balance between localization accuracy and complexity compared to other considered algorithms.
The remainder of this paper is structured as follows. An overview of the existing literature is given in Section 2. The TDOA-based problem of passive target localization is presented in Section 3. In Section 4, the function for the problem of passive target localization based on ML estimation is derived. The formulation of the SDP method is described in Section 5. In Section 6, the conventional BOA algorithm with the improved version is described. The conventional PSO algorithm with the proposed modification of the inertia weight improved with the logistic chaos map is presented in Section 7. The proposed CAHBPSO algorithm is presented in Section 8. Section 9 gives the formulation of the CRLB for the considered passive target localization problem. Section 10 provides the numerical simulation results. The conclusions and directions for future work are drawn in Section 11. Finally, the derivation of the CRLB is given in Appendix A.

2. Background and Related Work

Accurate estimation of the location of a passive target in the presence of additive measurement noise has become an important and challenging issue [38]. The estimated target position’s accuracy is mostly determined by two factors: measurement precision and the geometric layout of receivers, transmitters, and target [39]. The two types of localization algorithms that can be roughly categorized are range-based and range-free algorithms. Range-based localization algorithms estimate a target’s position effectively by using distance or angle information between the target and receivers. These methods are based on data collected using a variety of ranging approaches, including time of arrival (TOA) [38], TDOA [17], received signal strength (RSS) [40], angle of arrival (AOA) [41] and their combinations [42]. Because of its capacity to produce high localization accuracy, the TDOA is one of the most extensively utilized techniques [43]. In contrast to range-based localization algorithms, range-free localization techniques estimate the unknown node position using connectivity and topological information [44]. Distance vector hop (DV-Hop) [45], centroid or weighted centroid technique [46], and approximation point-in-triangulation test (APIT) [47] are the most extensively used range-free algorithms. Because these methods do not require a sophisticated hardware structure to determine range measurements, they are both inexpensive and simple to build. However, when compared to range-based algorithms, range-free algorithms typically have inferior localization accuracy [48].
The nonlinear least squares (NLS) estimator [49], which is obtained by minimizing the sum of the squared measurement errors, is a widely used estimation method for obtaining the unknown position of a target. Generally, the solution in closed form for NLS estimator is not obtainable, making it a challenge to obtain the solution of NLS nonlinear and non-convex objective function. To obtain the closed-form solution of different localization problems the linear least squares (LLS) and weighted least squares (WLS) methods are often applied [50]. However, these methods do not provide the required localization accuracy, and thus can be applied to provide an initial solution to the iterative optimization methods.
In this way, a number of local search optimization algorithms are widely applied to solve different localization problems [51]. We differentiate between the two groups of optimization algorithms, whether or not they require knowledge of the gradient of objective function, both of which perform local searches around the given starting solution [52]. Hence, numerous local search algorithms, such as Gauss–Newton, Nelder–Mead, and Conjugate gradient method are employed to estimate the unknown target position by solving the NLS estimation problem. However, finding the global optimum using local search optimization algorithms highly depends on the provided initial solution [52]. This leads to the conclusion that the convergence of local search algorithms while handling multimodal optimization problems is typically not expected without suitable initialization.
Another widely used estimation method is the ML estimator, which is commonly applied when the measurement error distribution is previously known [53]. However, because of the nonlinearity and non-convexity of the ML estimator’s objective function, typical local search techniques cannot be used to address this sort of complicated optimization problem. In this regard, to formulate a convex problem, second-order cone programming (SOCP) and SDP methods are widely applied to transform ML estimation problem, overcoming the non-convexity of the ML objective function [54]. In comparison to the SOCP method, the simulation results of the comparative study of both algorithms reveal that the SDP method gives superior precision of the target location [55]. The SDP and SOCP methods, on the other hand, cannot the achieve desired precision of the estimated solution, especially in the cases when severe measurement noise is present. As a result, there is a lot of interest in improving and developing efficient EAs that can be used to identify the global best solution to the nonlinear and non-convex ML estimation problem.
Various EAs, such as CS, PSO, BOA, and FA, etc., are applied to solve different localization problems by estimating the unknown location of a target [9,18,56]. As a result, finding an effective optimization algorithm for a specific localization problem is critical in order to reduce the localization error in all situations. A significant increase on the positioning accuracy is achieved by using PSO and TDOA together [57]. Results are compared with well-known WLS and LLS techniques and the results reveal that the PSO method outperforms the well-known methods. Furthermore, the PSO algorithm based on chaos theory is proposed for the hybrid TDOA/AOA location estimation problem, where the objective function is formulated using the ML estimator [9]. The comparative analysis shows that the improved PSO algorithm has enhanced global search ability compared to the conventional optimization algorithms. For the localization of wireless sensor nodes in “concave areas”, researchers suggested a two-stage PSO technique. The approach can achieve excellent localization accuracy in wireless sensor networks while using little energy and processing resources [58]. In [18], the BOA algorithm is employed to estimate the location of nodes by minimizing the considered objective function. The results of the localization performance of BOA demonstrate more consistent and accurate location of nodes compared to widely applied methos, such as PSO and FA.
In addition, the hybridization of EAs is proven to be as an effective way to enhance the quality of the solutions and improve the optimization efficiency. In this way, using TDOA measurements a hybrid firefly algorithm (hybrid-FA), which combines the WLS and FA algorithms, is proposed in [17] for target coordinates estimation. According to the simulation findings, the hybrid-FA method outperforms other well-known localization algorithms in terms of performance and localization accuracy. In addition, various hybrid variants of EAs, such as PSO and variable neighborhood search (HPSOVNS) [59], DE and FA (HFDLA) [60] are successfully applied to enhance positioning accuracy.
Based on the preceding, an improved CAHBPSO algorithm is proposed here to increase the accuracy of the estimated target position, especially when measurement noise increases.

3. Localization Problem

This section looks at the problem of determining the unknown location of a passive target in a LOS environment using noisy TDOA measurements. As shown in Figure 1, the investigated localization system consists of one receiver at the origin of the coordinate system x r = 0 0 T ,  N transmitters with predefined known coordinates x i t = x i t y i t T R 2 , i 1 , 2 , , N , and the unknown position of the passive target at x = x y T R 2 .
In the passive target localization problem, the transmitters emit a signal, and the target reflects the signal from each of the transmitters in all directions. The TDOA measurements are then obtained by the receiver capturing the reflected signal as well as the direct signal from each of the transmitters. Furthermore, the receiver and transmitters are assumed to be perfectly synced, and it is assumed that transmitter emitted signals are reflected from the passive target across all possible directions while not engaging directly with the receiver and transmitters [37]. When modelling the measurement error it is valid to assume that its probability distribution is Gaussian, which holds in LOS conditions [2]. Therefore, the noisy TDOA measurements are
t i = 1 c x i t x 2 + x r x 2 x i t x r 2 + n ¯ i , i 1 , 2 , , N ,
where · 2 is the Euclidean distance in two dimensions, c denotes a constant which is equal to the known speed of light, and n ¯ i represents measurement noise whose underlying probability distribution, according to the assumption, is Gaussian. The range measurements r i i = 1 N , can be derived by multiplying Equation (1) with constant c, as follows
r i = c · t i = x i t x 2 + x r x 2 x i t x r 2 + n i , i 1 , 2 , , N ,
where n i = c n ¯ i follows the zero-mean Gaussian distribution. Then, introducing a new variable r ˜ i = r i + x i t x r 2 , and substituting it into the range measurements in Equation  (2), the following expression can be obtained
r ˜ i = x i t x 2 + x r x 2 + n i , i 1 , 2 , , N .
Then, the vector form of Equation (3) can be expressed as
r ˜ = d x + n ,
where
d x = x 1 t x 2 + x r x 2 x N t x 2 + x r x 2 ,
and n = n 1 , n 2 , , n N T is the vector of zero-mean Gaussian noise. It is assumed that all elements of n are independent and identically distributed. In this regard, the covariance matrix can be obtained as C = E n n T = σ 2 I N [37], where E · denotes the expectation operator and I N is N × N identity matrix. Then, the vector d x can be rewritten as
d x = H g x ,
where
g x = x r x x 1 t x x N t x T ,
H = 1 N I N = 1 1 0 0 1 0 0 0 1 0 0 1 ,
in which 1 N denotes the unit column vector of length N. Then, Equation (4) can be reformulated as
r ˜ = H g x + n .
As a result, the goal of this paper is to solve the nonlinear and non-convex ML estimation problem and to effectively estimate the unknown position of the passive target x based on noisy TDOA data.

4. Maximum Likelihood Estimator

The Maximum likelihood estimator can be successfully employed to determine the unknown coordinates of the passive target by determining the extremum of the likelihood function. Under the assumption that TDOA measurements are independent and identically distributed Gaussian zero-mean random variables, the likelihood function L r ˜ x of the obtained TDOA measurements can be expressed as
L r ˜ x = f r ˜ x = = 1 2 π N / 2 det C 1 / 2 exp 1 2 r ˜ H g x T C 1 r ˜ H g x ,
where f r ˜ x is the probability density function of the measurements. Then, taking the logarithm of the likelihood function yields
ln L r ˜ x = k 1 2 σ 2 r ˜ H g x T r ˜ H g x ,
where k = ln 1 2 π N / 2 det C 1 / 2 can be neglected as it does not depend on x . In this regard, the estimated position of the passive target x ^ is obtained as the solution of the following non-convex optimization problem
x ^ = argmin x R 2 J M L x ,
where the ML objective function corresponding to this may be written as
J M L x = r ˜ H g x T r ˜ H g x .
Figure 2 shows the corresponding contour plot of the ML objective function J M L x .
According to the plot of the the objective function J M L x in Figure 2 it is concluded that the J M L x is nonlinear function which has multiple local optima. Furthermore, it is noticed that the position of the global minimum of this function correlates to the coordinates of the unknown position of the target. As a result, advanced optimization methods, discussed in reminder of the paper, are required to obtain the global optimal solution.

5. Semidefinite Programming Method

This section presents the SDP approach to deal with non-convexity of the ML estimation problem by transforming it into a convex optimization problem in order to solve the passive target localization problem [2]. As described, the considered localization problem leads to the non-convex and multimodal ML objective function. In order to try to solve this problem, the SDP method is developed for the considered passive target TDOA-based localization problem, which converts the objective function J M L x to a convex function. The ML estimation problem Equation (12) with regard to x can be expressed as follows:
x ^ = arg min x R 2 r ˜ H g x T r ˜ H g x .
The objective function of the considered ML optimization problem can be transformed into a set of linear equations, by squaring both sides of Equation (2) and introducing an additional variable between the receiver and the target R r = x r x , as follows
2 x r x i t T x 2 r ˜ i R r + r ˜ i 2 + x r 2 x i t 2 = 2 r ˜ i R r n i , i 1 , 2 , , N ,
where the second-order term of the noise n i 2 is neglected. Introducing the variable θ , as
θ = x R r T ,
the Equation (15) can be expressed in the linear-matrix from
b A θ = m ,
where
A = 2 x r x 1 t y r y 1 r r ˜ 1 x r x 2 t y r y 2 r r ˜ 2 x r x N t y r y N r r ˜ N ,
b = r ˜ 1 2 + x r 2 + y r 2 x 1 t 2 y 1 t 2 r ˜ 2 2 + x r 2 + y r 2 x 2 t 2 y 2 t 2 r ˜ N 2 + x r 2 + y r 2 x N t 2 y N t 2 ,
and
m = 2 r ˜ 1 R r n 1 r ˜ 2 R r n 2 r ˜ N R r n N .
Then, based on Equation (16) through Equation (20), the following WLS optimization problem can be formulated as
min θ J θ = min θ b A θ T W b A θ
where W R N × N is a symmetric weighting matrix. Under the sufficiently small measurement noise, the symmetric weighting matrix can be approximated as
W = E m m T 1 = D T C D 1 ,
where
D = diag 2 r ˜ 1 R r , 2 r ˜ 2 R r , , 2 r ˜ N R r ,
It should be noted that since the measurement noise m from Equation (20) is Gaussian distributed and due to the linear relationship in Equation (17), the objective function of the ML estimator, given in Equation (13), is equivalent to that of the WLS estimator in Equation (21) [61].
Then, by introducing the range between x r and x , denoted by R r = x r x 2 as the equality constraint, the ML estimation problem in Equation (14) is expressed as
min x f o b = min θ b A θ T W b A θ s . t . R r = x r x 2 .
After corresponding algebraic manipulation, the objective function of the optimization problem f o b in Equation (24) becomes
f o b = i = 1 N b i a i θ 2 r ˜ i R r 2 = i = 1 N b i r ˜ i R r a ˜ i θ r ˜ i R r ν i R r r ˜ i R r 2 ,
where b i denotes the ith element of vector b , a i is the ith row of matrix A . Here, a i = a ˜ i ν i denotes a block matrix, where the submatrices are a ˜ i = a i 1 a i n and ν i = a i n + 1 .
To transform the objective function, the matrix property x T Ax = Tr x x T A is employed. After introduction of the matrix notation z = x T x , the Equation (25) can be rewritten as
x 1 T i = 1 N a ˜ i a ˜ i T i r ˜ i R r 2 i = 1 N a ˜ i ν i R r b i r ˜ i R r 2 i = 1 N a ˜ i T ν i R r b i r ˜ i R r 2 i = 1 N ν i R r b i 2 r ˜ i R r 2 x 1 = = Tr P x T 1 x 1 = Tr P z x x T 1 ,
in which
P = i = 1 N a ˜ i a ˜ i T i r ˜ i R r 2 i = 1 N a ˜ i ν i R r b i r ˜ i R r 2 i = 1 N a ˜ i T ν i R r b i r ˜ i R r 2 i = 1 N ν i R r b i 2 r ˜ i R r 2 ,
and Tr · denotes the trace of a square matrix.
Then, based on Equations (26) and (27) the optimization problem in Equation (24) can be equivalently written as
min x , z Tr P z x x T 1 s . t . R r = x r x 2 z = x T x .
The constraint in Equation (28) can be reformulated using the notation z = x T x , as follows
R r 2 = x r x 2 2 = x r T x r x r T x x T x r + x T x R r 2 = Tr I n x r x r T x r T x r x T 1 x 1 = Tr I n x r x r T x r T x r z x x T 1 .
Finding the global optimal solution of the optimization problem in Equation (28) is difficult due to the non-convex equality constraint z = x T x . Therefore, the equality constraint z = x T x is relaxed to a convex constraint as  
z x T x 0 ,
Then, after applying the Schur complement [62], the constraint can be equivalently rewritten as a linear matrix inequality, as follows
z x x T 1 0 .
Hence, the obtained matrix in Equation (31) is symmetric and positive semidefinite matrix. Finally, the optimization problem in Equation (28) can be phrased as follows, based on the above relaxation.
min x , z , R r Tr P z x x T 1 s . t . R r 2 = Tr I n x r x r T x r T x r z x x T 1 z x x T 1 0 .
It should be noted, that for each fixed R r , the obtained problem in Equation (32) becomes convex optimization problem. However, this optimization problem becomes non-convex, when the Equation (32) is solved with respect to variables { x , z , R r } , e.g., when the value of R r is not previously known. Therefore, it is important to find the optimal value of R r , for which the problem in Equation (32) becomes convex. Thus, in this paper, the golden section optimization algorithm (GSA) [52] is employed in order to determine the optimal value of R r . In this regard, two points within the interval R l r , R u r can be calculated using the equations
R 1 r = R l r + 1 φ R u r R l r , R 2 r = R l r + φ R u r R l r
where φ = 1 + 5 / 2 represents the golden ratio. Afterwards, for each point R j r , j = 1 , 2 obtained in Equation (33), the solutions x j and z j of the optimization problem in Equation (32) are found, and the objective function can be evaluated at these points. If f o b x 1 , z 1 , R 1 r < f o b x 2 , z 2 , R 2 r , then the optimal point belongs to the interval R l r , R 2 r , otherwise if f o b x 1 , z 1 , R 1 r > f o b x 2 , z 2 , R 2 r the solution belongs to R 1 r , R u r . Alternatively, if f o b x 1 , z 1 , R 1 r = f o b x 2 , z 2 , R 2 r the boundaries are reduced to R l r = R 1 r and R u r = R 2 r , and the optimization process is repeated. The values of R 1 r and R 2 r are calculated iteratively until the difference R 1 r R 2 r ε is less than a predefined positive number ε .
Therefore, the procedure for determining the optimal value of R r , for which the considered SDP optimization problem becomes convex problem can be stated as follows:
Step 1:
Initialize the R l r and R u r .
Step 2:
Calculate the points R 1 r and R 2 r according to Equation (33).
Step 3:
Solve Equation (32) with R r = R 1 r and R r = R 2 r .
Step 4:
If f o b x 1 , z 1 , R 1 r < f o b x 2 , z 2 , R 2 r , then the interval becomes R l r , R 2 r , otherwise interval is R 1 r , R u r . If f o b x 1 , z 1 , R 1 r = f o b x 2 , z 2 , R 2 r boundaries become R l r = R 1 r and R u r = R 2 r and go to Step 2.
Step 5:
Repeat Steps 2–4 until R 1 r R 2 r ε is satisfied.

6. Butterfly Optimization Algorithm and the Proposed Improved Version

6.1. Conventional BOA Algorithm

The butterfly optimization algorithm is a novel nature-inspired metaheuristic algorithm, where the search process is inspired by the food foraging behavior and the process of mating between butterflies [13]. The BOA is based on three assumptions:
  • All butterflies are said to release some fragrance in order to attract one another.
  • Each butterfly either moves randomly or towards the butterfly with the strongest fragrance (i.e., the best butterfly in the current generation)
  • The stimulus intensity of a butterfly is proportional to the objective function value.
In general, the optimization process of the BOA algorithm can be divided into three stages: initialization, iteration, and the final optimization stage. In the first stage, the values of control parameters of the algorithm are defined, and the set of initial solutions is randomly generated within the upper and lower bounds. In the iteration stage, the search for the global optimal solution is performed. Firstly, the algorithm calculates the objective function value of each butterfly, and then butterflies generate the fragrance at their positions. The intensity of the fragrance sensed by the butterfly, φ i , can be described as a function of the physical intensity of stimulus, as follows
φ i = c I a ,
where c denotes the sensory fragrance, a represents the power exponent, and I is the stimulus intensity, which is proportional to the objective function value. Here, coefficients c and a are assigned in the range 0 , 1 . In the case a = 1 , there is no fragrance absorption, while when a = 0 , other butterflies cannot detect the fragrance produced by any butterfly. Stimulus intensity I for ith butterfly can be determined as
I = f x i G ,
where x i G denotes the position of the ith butterfly in the Gth generation, and f is the objective function of the considered optimization problem.
During the search for the global optimum, the BOA algorithm goes through two key phases, the global search phase and the local search phase. In the first phase, the butterfly is updating its position according to the global best solution
x i G + 1 = x i G + r 2 × g * x i G × φ i ,
where g * represents the best solution found in the current iteration (i.e., butterfly with the strongest fragrance) and r represents uniform random number in the range r 0 , 1 . On the other hand, the local search phase of the BOA algorithm can be written as
x i G + 1 = x i G + r 2 × x j G x k G × φ i ,
where x j G and x k G are the jth and kth butterfly positions within the search space, respectively.
The switch probability p is used to transition between the global and local search phases of the algorithm, as shown below
x i G + 1 = x i G + r 2 × g * x i G × φ i , if r p < p x i G + 1 = x i G + r 2 × x j G x k G × φ i , otherwise
where r p denotes uniformly generated random number in the range 0 , 1 . Therefore, the global search phase of the BOA algorithm, defined in Equation (36), is applied if r p < p . Otherwise, the local search phase, given in Equation (37), is employed to search for the optimal solution. The above-mentioned process is repeated until the stopping criteria is satisfied and the optimal solution is obtained, in the final optimization stage of the algorithm.

6.2. Improved BOA Algorithm

The conventional BOA algorithm provides an excellent local search ability; however, it suffers from premature convergence to local optima due to its poor exploitation ability [20]. Therefore, there is a need to modify the trade-off between the global and local search of the BOA in order to solve multimodal and complex optimization problems.
Analyzing Equation (34), it is evident that sensory fragrance c is one of the most important parameters in BOA, which guides the movement of butterflies during the search process by enabling each butterfly to sense the fragrances emitted by other butterflies. Therefore, using a constant value for parameter c is not suitable for complex optimization problems. Using a small value of c, during the entire search process, can result in premature convergence. On the other hand, in the early stage of the search process, a large constant value of sensory fragrance c leads to a high probability of missing the area of global optimal solution, which reduces the optimization performance. In this regard, it is critical to adjust the value of the sensory fragrance c adaptively during the search process in order to enhance the balance between exploration and exploitation ability.
Therefore, to provide an appropriate balance between global and local search abilities, in this paper, an adaptive sensory fragrance c G + 1 has been introduced, which can be described as follows
c G + 1 = 1 1 + exp c G G 0.2 G max ,
where G max denotes the maximum number of generations. In this respect, the changes of proposed adaptive sensory fragrance c G + 1 , defined in Equation (39), versus the number of generations is illustrated in Figure 3.
According to Figure 3, it can be observed that the adaptive sensory fragrance c G has smaller value at the beginning of the search process, and achieves a larger value with the increase of generations. Therefore, in the early stage of the search process, the smaller value of c G can enhance the global exploration ability and prevent premature convergence. On the other hand, in the later stage of search process, a larger value of c G can improve exploitation ability and convergence speed of the algorithm.

7. Particle Swarm Optimization and the Proposed Improved Version

7.1. Conventional PSO Algorithm

In the PSO algorithm [23], the search process is performed based on velocity and position vectors. Each particle in a swarm of N P particles at Gth generation, has a position vector x i G = x i , 1 G , x i , 2 G , , x i , n G T and velocity vector v i G = v i , 1 G , v i , 2 G , , v i , n G T , i 1 , 2 , , N P , in the n-dimensional space. In the initial generation, the position and velocity vectors are randomly generated within the upper and lower bounds. During the optimization process, each particle moves through the search space with velocity v i G , which depends on the personal previous best position p i G = p i , 1 G , p i , 2 G , , p i , n G T and best position discovered by the whole population g G = g 1 G , g 2 G , , g n G T . Therefore, the position and velocity vectors of each particle at G + 1 th generation are updated as follows
v i G + 1 = v i G + c 1 r 1 p i G x i G + c 2 r 2 g G x i G ,
x i G + 1 = x i G + v i G + 1 ,
where c 1 and c 2 are the cognitive and social acceleration coefficients, respectively, which are commonly set to 2 [26]. Here, r 1 and r 2 are two distinct random numbers uniformly distributed in the range 0 , 1 . In order to ensure the balance between the exploration and exploitation abilities during the optimization process, the inertia weight ω G is introduced into the PSO algorithm. Therefore, the velocity vector v i G + 1 of each particle is updated using the linear-decreasing inertia [63], as follows
v i G + 1 = ω G v i G + c 1 r 1 p i G x i G + c 2 r 2 g G x i G ,
ω G = ω max G G max ω max ω min ,
where ω max and ω min denote the maximum and minimum values of the inertia weight, respectively. In the literature, the maximum and minimum value of inertia weight is commonly set to ω max = 0.9 and ω min = 0.4 [63].

7.2. Chaos Enhanced PSO Algorithm

The conventional PSO algorithm exhibits the issue of premature convergence to local optima, which may affect its optimization performance in solving complex optimization problems [26]. In this regard, the introduction of chaos theory into the PSO algorithm can improve the optimization performance, by modifying the PSO algorithm to escape more easily from local optima [30,64,65]. Therefore, in this paper, the logistic chaos map is implemented to dynamically adjust the value of inertia weight, in order to maintain the appropriate balance between global exploration and local exploitation abilities during the optimization process.

Chaotic Dynamic Inertia Weight

It is well known in the literature that the inertia weight has an important role in maintaining the balance between exploitation and exploration abilities during the optimization process [27]. Analyzing Equation (43), it is evident that a large value of inertia weight, in the early stage of the search process, can improve the global search abilities of the PSO algorithm and prevent the problem of premature convergence. On the other hand, a smaller value of inertia weight, in the later stage of search process, can improve the local search ability of the PSO algorithm. A number of different inertia weight strategies are proposed in the literature, among which the linear-decreasing inertia weight given in Equation (43) is widely employed. According to the analysis [27,66], choosing the appropriate inertia weight strategy for the current optimization problem depends on the properties of the objective function. In order to obtain an appropriate balance between exploration and exploitation skills for solving complex optimization problems, additional modifications to the inertia weight are necessary.
In recent years, the introduction of chaos theory is emerging as a powerful approach for improving the optimization performance of different metaheuristic algorithms [64,65]. Chaos is a bounded dynamic behavior that can be observed in certain nonlinear dynamic systems. In this regard, the chaotic behavior can be represented with chaos maps, which produce a bounded sequence of random numbers depending on initial condition. In this way, the well-known logistic chaotic map [65], which has the properties of ergodicity, non-repetition and irregularity, is employed in this paper, to adjust the value of inertia weight. The expression for logistic map is given as
x c G + 1 = a ω x c G 1 x c G , G = 1 , 2 , , G max .
where a ω denotes the control parameter. It should be noted that for certain initial conditions, e.g., x c 0 0 , 0.25 , 0.5 , 0.75 , 1 and a ω = 4 , the logistic chaos map exhibits chaotic behavior.
In Figure 4, the unbounded chaotic behavior produced by logistic chaotic map, given in Equation (44) is presented as a function of generation number G, for the x c ( 0 ) = 0.9 and a ω = 4 .
Furthermore, to better control the inertia weight parameter, the logistic chaos map is upper-bounded by introducing the term e G / G max . Therefore, the chaotic dynamic inertia weight ω c G + 1 can be calculated as follows
ω c G = e G G max · x c G
In this respect, Figure 5 illustrates the changes of the proposed chaotic dynamic inertia weight ω c G + 1 , defined in Equation (45), with the increase of generations.
As can be seen from Figure 5, the chaotic dynamic inertia weight ω c G + 1 , defined in Equation (45), is upper-bounded by the exponential factor e G / G max , and the value of ω c G + 1 gradually decreases with the increase of generations. In the early stage of the search process, ω c G + 1 has a larger value, which is suitable for enhancing the global search ability and finding the region of the global optimal solution. In the later stage of search process, a smaller value of ω c G + 1 can enhance the exploitation ability and improve the convergence towards the global optimum. This shows that the proposed chaotic dynamic inertia weight ω c G + 1 provides an effective balance between global exploration and local exploitation abilities, and thus improves the optimization performance of the PSO algorithm.

8. Chaos Enhanced Adaptive Hybrid Butterfly Particle Swarm Optimization Algorithm

To efficiently solve the considered complex passive target localization problem, in this section, CAHBPSO algorithm is introduced, as a hybridization of the BOA with PSO algorithm. In order to overcome the problem of premature convergence and enhance the exploration and exploitation abilities of PSO algorithm, in this paper, the global search and local search phases of conventional BOA algorithm are incorporated into the velocity update equation of the PSO algorithm. Instead of fixed switch probability, we propose an adaptive technique to dynamically adjust between global exploration and local exploitation abilities during the optimization process. In addition, the proposed adaptive strategy for updating the value of the sensory fragrance of the BOA algorithm is introduced into the CAHBPSO algorithm. Furthermore, the chaotic dynamic inertia weight, proposed in this paper, is incorporated into the hybrid algorithm with the aim to improve convergence and efficiently maintain the population diversity. Therefore, to improve the optimization performance and achieve a more efficient algorithm for complex optimization problems, the following two modifications are proposed.
Firstly, the exploration of the PSO algorithm is enhanced by replacing the term c 2 r 2 g G x i G in Equation (42) with r 2 × g * x i G × φ i from the global search phase of the BOA algorithm, given in Equation (36). Furthermore, v i G in Equation (42) is substituted with x i G from Equation (36), and linear inertia weight ω G is replaced with chaotic dynamic inertia weight ω c G + 1 . The equation for updating the ith butterfly’s location can thus be expressed as
x i G + 1 = ω c G + 1 x i G + c 1 r 1 p i G x i G + r 2 × g * x i G × φ i .
Next, further modification is done with the aim to enhance the local search ability PSO algorithm. In this regard, we exchange the term c 1 r 1 p i G x i G in Equation (42) with the term r 2 × x j G x k G × φ i from the local search phase of the BOA algorithm, given in Equation (37). Furthermore, v i G in Equation (42) is substituted with x i G from Equation (37) and the proposed chaotic dynamic inertia weight ω c G + 1 is introduced. Therefore, the new butterfly position is determined by
x i G + 1 = ω c G + 1 x i G + r 2 × x j G x k G × φ i + c 2 r 2 g G x i G .
In this paper, instead of using a fixed switch probability, an adaptive mechanism is proposed to dynamically adjust between global exploration and local exploitation during the optimization process. In this regard, the adaptive switch probability p G + 1 can be described as follows
p G + 1 = f m e a n G f b e s t G f w o r s t G f b e s t G ,
where f m e a n G , f b e s t G and f w o r s t G denote the mean, best, and worst values butterflies achieved in the previous generation in terms of objective function, respectively.
The trade of between global and local search abilities during the optimization process can be managed by changing the parameter p G + 1 , according to Equation (48). In the following analysis we can consider two two extreme cases. Firstly, the parameter p G + 1 is near to 1, indicating that the algorithm’s global exploration capabilities has to be improved since the diversity in the population is low. As a result, Equations (36) and (46) will be picked at random with a probability of 0.5, with the goal of improving global exploration and locating the global optimal solution region. In the second case, parameter p G + 1 is close to 0, indicating that all the solutions are near the global optimal solution, which shows that local exploitation must be improved. Therefore, Equations (37) and (47) will be randomly selected, with the probability of 0.5, to enhance exploitation ability and improve the convergence speed.
In this regard, the position of the ith butterfly can be updated based on the value of parameter p G + 1 according to the pseudocode shown in Algorithm 1.
Algorithm 1 Position update of the ith butterfly of the proposed hybrid CAHBPSO algorithm.
    if  p G + 1 > 0.5   then
        if  r a n d > 0.5  then
             x i G + 1 = x i G + r 2 × g * x i G × φ i
        else
             x i G + 1 = ω c G + 1 x i G + c 1 G r 1 p i G x i G + r 2 × g * x i G × φ i
        end if
    else if  p G + 1 0.5   then
        if  r a n d > 0.5  then
             x i G + 1 = x i G + r 2 × x j G x k G × φ i
        else
             x i G + 1 = ω c G + 1 x i G + r 2 × x j G x k G × φ i + c 2 G r 2 g G x i G
        end if
    end if
The proposed modifications, introduced in the CAHBPSO algorithm, provide an effective balance between exploration and exploitation abilities during the optimization process. Furthermore, these modifications are effective in overcoming the problem of premature convergence. In this way, the pseudocode of the proposed CAHBPSO algorithm is presented in Algorithm 2, for the considered passive target localization problem.
Algorithm 2 Pseudo-code of the proposed CAHBPSO algorithm.
    Generate initial population
    Determine stimulus intensity I i = f x i 0
    Set the value of parameters c 1 , c 2 , a, G m a x , a ω , ω ( 0 )
    Initialize values of p i ( 0 ) and g ( 0 )
    while stopping criteria not met do
        Calculate adaptive sensory modality function c G according to Equation (39)
        for each butterfly i in the population do
            Calculate fragrance φ i = c I a
        end for
        Find the best butterfly
        Calculate chaotic inertia weight according to Equation (45)
        for each butterfly i in the population do
            if  p G + 1 > 0.5  then
               if  r a n d > 0.5  then
                    x i G + 1 = x i G + r 2 × g * x i G × φ i
               else
                    x i G + 1 = ω c G + 1 x i G + c 1 G r 1 p i G x i G + r 2 × g * x i G × φ i
               end if
            else if  p G + 1 0.5  then
               if  r a n d > 0.5  then
                    x i G + 1 = x i G + r 2 × x j G x k G × φ i
               else
                    x i G + 1 = ω c G + 1 x i G + r 2 × x j G x k G × φ i + c 2 G r 2 g G x i G
               end if
            end if
        end for
        Find the global best solution g ( G )
        Determine the personal best solution according to:
         p i G + 1 = p i G , if p i G x i G x i G , otherwise
    end while

9. Cramer–Rao Lower Bound

The CRLB for the passive target localization problem provides a lower bound on the covariance matrix of any unbiased estimator [37]. Therefore, in this paper, the CRLB is used as a benchmark to evaluate the performance of the considered estimators. The derivation of the CRLB can be obtained from the inverse of the Fisher information matrix (FIM) F x , which can be defined as
F x = E ln f r ˜ x x ln f r ˜ x x T = E 2 ln f r ˜ x x x T .
Thus, the components of the FIM are given as follows
F x = F 11 F 12 F 21 F 22 ,
where the corresponding elements can be obtained as
F 11 = 1 σ 2 i = 1 N x x r x x r + x x i t x x i t 2 2 ,
F 12 = F 21 = 1 σ 2 i = 1 N x x r x x r + x x i t x x i t 2 × y y r x x r + y y i t x x i t 2 ,
F 22 = 1 σ 2 i = 1 N y y r x x r + y y i t x x i t 2 2 .
The derivations of Equations (51)–(53) are given in Appendix A. Then, the relationship between the variance and CRLB can be defined as
E x ^ x x ^ x T Tr F x 1 = CRLB x ,
where x ^ denotes the estimated value of x .

10. Experimental Study

This section present the experimental results conducted in order to evaluate the localization accuracy of the proposed CAHBPSO algorithm and compare the optimization performance with the well-known algorithms in the literature on a set of CEC2014 benchmark problems using the statistical analysis. Therefore, the obtained results are outlined in the following two subsections.

10.1. Statistical Evaluation of CAHBPSO Method against the CEC2014 Benchmark

This section outlines the results of the statistical comparison of the optimization performance between the proposed CAHBPSO algorithm and widely applied algorithms in the literature, including PSO [23], BOA [13], SHADE [67] and HPSOBOA [33] on a set of CEC2014 benchmark problems. The CEC2014 benchmark problems consist of 30 single-objective real-parameter numerical optimization problems, where D = 10 , 30 , 50 , and 100 are the considered dimensions of the search space, which are defined in [68]. The considered CEC2014 benchmark problems can be classified into four groups:
  • f 1 f 3 unimodal optimization problems;
  • f 4 f 16 simple multimodal objective functions;
  • f 17 f 22 hybrid objective functions, in which variables are subdivided and various basic functions are applied to each subset;
  • f 23 f 30 composition functions, which provide continuity around the optimal solution and merge the properties of sub-functions.
The metric of the obtained solution error f x ^ f x * was established to make a statistical comparison of the optimization performance of CAHBPSO and other widely applied algorithms, where x ^ represents the global optimum of each algorithm produced in a single run and x * denotes the solution to the CEC2014 benchmark problem, which is previously known. For each of the test functions, each algorithm was run 51 times with the termination threshold set at 10,000D and the swarm size was N P = 100 particles. The search space for each objective function is defined as [ 100 100 ] D . As a consequence, the obtained experimental findings are examined and compared using two nonparametric statistical hypothesis tests, such as the Wilcoxon signed-rank test and the Friedman test, to perform statistical assessment of optimization performance.
To perform pair-wise comparison of optimization performance between the proposed CAHBPSO and other well-known algorithms, the Wilcoxon signed-rank test is applied. In this regard, statistical comparison findings can reveal if the first algorithm outperforms the second method statistically. The significance level for this statistical test was set at 0.05. In the findings, R + means the total of rankings in which the first algorithm exceeded the second, whereas R denotes the sum of ranks in which the second algorithm outperformed the first algorithm. The Wilcoxon signed-rank test’s null hypothesis asserts that “there is no difference between the mean findings of the two samples” [69]. The alternative hypothesis, on the other hand, asserts that “there is a difference in the mean results of the two samples”. Therefore, using the p value and comparing it with the significance level α , the null hypothesis can be rejected when p α . In this regard, the obtained results are denoted with signs +, ≈, − according to the result of the statistical test. Here, plus sign (+) denotes that the first algorithm had significantly better optimization performance than second one, minus sign (−) indicates that the first algorithm performed significantly worse than the first one, while the sign (≈) denotes that there is no significant difference in optimization performance between two considered algorithms.
The Friedman test is used in this study to determine the substantial difference between the optimization performances of the studied methods. In order to determine the rankings of all examined algorithms across each CEC2014 objective function over each of the search space dimensions D, the Friedman test is utilized. As a result, the algorithm with the lowest rank performs the best in terms of optimization, while the method with the highest rank performs the worst. The Friedman test’s null hypothesis is that “there is no difference among the performance of all algorithms”, whereas the alternative hypothesis is that “there is a difference among the performance of all algorithms” [69]. When the p value is less than or equal to α = 0.05 , the null hypothesis can be rejected.
The numerical simulations results of the proposed CAHBPSO and other considered algorithms, computed over 51 independent trials on a set of CEC2014 benchmark problems are presented in Table 1. The results are presented in terms of mean (Mean) and standard deviation (STD) of the best obtained objective function value. Furthermore, the sign is added to indicate if the examined method outperforms (+), performs similarly (≈), or worse (−) than the suggested CAHBPSO method.
Based on the findings in Table 1, it is clear that the proposed CAHBPSO algorithm outperforms other studied algorithms in terms of “Mean” and “STD” values on the vast majority of benchmark functions. However, compared to the PSO algorithm, the proposed algorithm achieved a worse performance on simple multimodal functions f 4 , f 13 on dimension D = 30 , f 12 , f 16 on dimensions D = 50 , 100 , and f 11 on D = 100 . Furthermore, a worse performance is observed on hybrid function f 22 on dimension D = 100 and on composite functions f 24 , f 27 on dimensions D = 30 , 50 , 100 , as well as functions f 24 , f 25 , f 26 on dimensions D = 100 , function f 26 over dimension D = 10 , and f 28 on dimensions D = 30 . Compared to the BOA algorithm, the proposed CAHBPSO algorithm showed a worse performance on composite functions f 23 and f 27 on all dimensions, f 24 , f 25 , f 29 , f 30 on dimensions D = 30 , 50 , 100 , and function f 26 on dimensions D = 50 , 100 . When compared to the SHADE and FA algorithms, it is clear that CAHBPSO outperforms them on all functions across all dimensions. Only on the hybrid and composite functions f 18 , D = 10 and f 26 , D = 50 did the proposed algorithm perform worse than the SHADE algorithm. As a result, when compared to the SHADE and FA algorithms, the suggested CAHBPSO algorithm delivers the best results. The proposed algorithm performed similarly to or slightly worse than the BOA and PSO algorithms on problems with greater dimensions, such as D = 30 , 50 and 100 on composite functions. However, the suggested method surpassed the BOA and PSO algorithms on the majority of test functions.Compared to the HPSOBOA algorithm, it is observed that the proposed algorithm achieved worse performance over all considered dimensions on functions f 6 , f 10 , f 11 , f 13 , f 14 , f 28 and f 30 , while the CAHBPSO algorithm outperformed HPSOBOA method on functions f 15 , f 16 , f 17 , f 18 , f 19 , f 20 , f 21 , f 22 and f 29 over all dimensions. On other functions the results are mixed, depending on the dimensionality of the problem. It can be observed that when dimension D = 10 is concerned the CAHBPSO algorithm achieved the best result in the majority of the results compared to HPSOBOA. As a conclusion, the CAHPBSO method’s improvements and hybridization are confirmed, and the CAHBPSO algorithm effectively discovers potential solutions across all CEC2014 benchmark tests and dimensions.
The obtained numerical results are analyzed using Wilcoxon signed-rank test, in order to perform the pair-wise optimization performance comparison between the proposed CAHBPSO and other considered algorithms. Therefore, the results of statistical comparison using Wilcoxon’s signed-rank test on CEC2014 benchmark problems are presented in Table 2.
From the results in Table 2, it can be observed that the proposed CAHBPSO algorithm has significantly better performance than the SHADE, FA, BOA, and PSO algorithms on all considered dimensions. Compared to the BOA and PSO algorithms, on D = 30 , 50 , and 100 an increase is observed in cases where the proposed algorithm had worse performance. When comparing to the HPSOBOA algorithm, we observe that the proposed CAHBPSO algorithm achieved better performance for D = 10 , while for other dimensions the algorithms had similar performance. However, in all considered cases, the CAHBPSO algorithm achieved higher R + values than R compared to considered algorithms.
Furthermore, the Friedman test is used to assess whether there is a significant statistical difference in the optimization performances of the algorithms considered. Table 3 displays the average ranks obtained by the Friedman test of examined algorithms on various CEC2014 problems over all dimensions. The algorithm’s top rank is bolded, while the second best is underlined.
From the statistical comparison results presented in Table 3, it can be observed that the proposed CAHBPSO algorithm achieves the best performance among the considered algorithms and has the lowest rank in all considered cases. The HPSOBOA algorithm had the second-best performance. It is observed that the obtained p values for the considered Friedman test are less than the significance level α = 0.05 in all dimensions over all considered cases. As a result, the hypothesis suggests that there is a considerable discrepancy in the performance of the algorithms under consideration. Based on the statistical analysis findings, it is determined that the proposed hybridization of BOA and PSO algorithms increases optimization performance while overcoming the disadvantages of both techniques. Furthermore, the enhanced approach for updating butterfly positions, as well as the addition of a chaotic map to the inertia weight and adaptive sensory fragrance, considerably increase the CAHBPSO algorithm capabilities, demonstrating the efficacy of proposed improvements.

10.2. Localization Performance of the Proposed CAHBPSO Algorithm

In this section, simulations are performed to verify and compare the localization performance of the proposed CAHBPSO with the existing algorithms, such as SDP, PSO, BOA, and the well-known WLS method [50], for the considered passive target localization problem in LOS environment. Here, the derived CRLB is utilized as a performance benchmark in the following simulations to assess the localization capabilities in terms of RMSE metric.
The simulations are performed to evaluate the localization accuracy of all considered algorithms under the noisy TDOA measurements. Furthermore, the performance of different algorithms is assessed in relation to changes in transmitter, target, and receiver geometric configurations. The numerical simulations are carried out using a passive target localization system consisting of one receiver and four transmitters positioned at known locations, as well as a passive target located at various locations. Three simulation scenarios are considered in this regard, depending on the chosen coordinates where the passive target is placed: (i) the passive target is surrounded by the four transmitters; (ii) the target is outside the hull which are transmitters-forming; and (iii) the passive target is deployed randomly over an area 80 × 80 m 2 for each simulation run. The identical receiver and transmitter setup is used in each simulation scenario, with the receiver at x r = 0 0 T , and four transmitters creating a convex hull at x 1 t = 80 80 T m, x 2 t = 80 80 T m , x 3 t = 80 80 T m , and x 4 t = 80 80 T m. Therefore, to evaluate the localization accuracy of the considered algorithms, the RMSE is employed, which is defined as
RMSE = 1 N m n = 1 N x ^ n x 2 2 ,
where x denotes the is the actual target location, x ^ n represents the position of a target that is estimated using one of the considered algorithms and N m = 1000 is the number of Monte Carlo simulation repetitions for a given variance of measurement noise σ 2 .
Figure 6 shows the results of the first simulation scenario, in which the true position of the passive target is x = 20 30 T m. The RMSE performances of the examined algorithms are displayed in relation to measurement noise p = 10 log σ 2 and compared to the CRLB.
The RMSEs obtained using the CAHBPSO method achieve the CRLB for the whole studied ranges of p, as shown in Figure 6. Furthermore, the RMSE performances of the WLS, SDP, BOA, and PSO algorithms are several dBs higher than the CRLB. However, the obtained numerical results reveal that the SDP approach deviates significantly from the CRLB for large values of p (esp. p > 20 dB) when compared to examined algorithms. It is observed that the RMSEs of the well-known WLS method show the worst performance over all considered methods.
The findings of the second simulation scenario, in which the true location of the target is placed at x = 100 80 T m, are presented in Figure 7, where the RMSE performances of the studied algorithms are plotted and compared with the CRLB.
According to the results in Figure 7, the proposed CAHBPSO method achieves CRLB accuracy and outperforms the investigated algorithms with the increase of p. In comparison to the CAHBPSO method, the BOA and PSO algorithms’ localization performance has decreased, as shown by the simulation results. In any considered case of measurement noise, the RMSE performances of the WLS method do not reach the CRLB. Furthermore, the SDP method performs poorly, especially when the measurement noise is high (esp. p > 20 dB).
For the third simulation scenario, where the position of the passive target is randomly generated inside the examined region for each simulation run, the RMSEs of all evaluated algorithms as the function of p are shown and compared against the CRLB in Figure 8.
As shown in Figure 8, the proposed CAHBPSO method achieves CRLB accuracy throughout the whole range of p and outperforms the WLS, SDP, BOA, and PSO algorithms in terms of localization accuracy.
When comparing the numerical results of the various simulation scenarios represented in Figure 6, Figure 7, Figure 8, it can be seen that when the target is positioned inside the convex hull of the transmitters, the examined methods perform well. When the target is situated outside of the convex hull of the transmitters, the RMSE performance is higher. Furthermore, the RMSEs of the proposed CAHBPSO method achieve the CRLB throughout the whole range of p in every simulated scenario studied.
To further evaluate the performance of all considered algorithms, the cumulative distribution functions (CDFs) of the localization error are obtained for the variance of measurement noise σ 2 = 10   m 2 . The localization error (LE) is defined as L E = x ^ n x 2 , n { 1 , , N m } . Therefore, Figure 9 shows the corresponding CDFs in terms of the localization error, for the first simulation scenario, obtained for each algorithm.
From the results in Figure 9, it is observed that in 90% of the cases, the proposed CAHBPSO algorithm provides a localization error of less than 3.66 m, while the SDP, BOA, and PSO algorithms have localization errors less than 4.1m, 4.0 m, and 3.86 m, respectively. As a result, when compared to the other algorithms studied, the suggested CAHBPSO method has the lowest localization error according to CDFs.
Finally, the effect of increasing the number of transmitters on the localization accuracy is investigated for the second simulation scenario. In this respect, the ith transmitter is placed at the following coordinates
x i t = R cos φ i R sin φ i , i 1 , 2 , , N j
where R = 80 2 m is the radius of a circle, φ i = 2 π / N j denotes the angular separation between transmitters, and N j is the number of transmitters taken for simulation. Hence, the RMSE performances of the considered algorithms as a function of number of transmitters are depicted in Figure 10, for the variance of measurement noise σ 2 = 1 m 2 .
As shown in Figure 10, the RMSE performances of all considered algorithms significantly improve when the number of transmitters is increased from 4 to 10. Moreover, the proposed CAHBPSO algorithm provides 56% improvement in localization accuracy when the number of transmitters is increased from 4 to 20. Based on the simulation findings, it is determined that, when compared to other investigated methods, the performance of the proposed CAHBPSO algorithm is not susceptible to excessive measurement noise and changes in network topology.

Computational Complexity of the Considered Algorithms

The computational complexity of the proposed CAHBPSO and other algorithms that were taken into consideration, as well as the typical computation time required to arrive at the overall optimal solution, are analyzed and compared in this subsection. In the literature, it is shown that the computational complexity of SDP method is O 16 N 2 n + N n 2 + G max n + 1 3 + 4 2 n + 3 3.5 [7]. Furthermore, the appropriate computational complexity of PSO algorithm can be written as O G max N P n + f , while the complexity of BOA is O G max n × N P + n × f [70], where f denotes the computational complexity of evaluating the considered objective function. All butterflies are sorted according to the objective function value in one generation of the proposed CAHBPSO algorithm, where the average computing complexity of this operation is O N P l o g N P . After sorting, the time complexity of calculating the fragrance of each butterfly is O N P . Then, in the global and local search phase of the algorithm, each butterfly goes through the process of exploring the search space, where the average time complexity required to update the position of ith butterfly is O N P × n + O N P × O f . Therefore, the overall computing complexity of CAHBPSO algorithm is determined as O N P l o g N P + O N P + G max O N P × n + O N P × O f , which can be simplified to O N P l o g N P + G max O N P × n + O N P × O f .
Following that, the average computing time required for identifying the global optimal solution is examined, as another crucial element influencing the performance of the algorithms under consideration. An identical PC with a 3.2 GHz CPU and 16 GB of RAM is used for the analysis. In this regard, for each simulation scenario the average computational times required to obtain global optimum, for each considered algorithm, are shown in Table 4.
As shown in Table 4, The BOA technique has the fastest implementation, whereas the SDP method shows the longest computing time considering all evaluated methods. Furthermore, the results reveal that the proposed CAHBPSO method provides the optimum balance of localization accuracy and average computation time to attain the global optimal solution.

11. Conclusions

In this paper, the passive target localization problem based on the noisy TDOA measurements, obtained from multiple transmitters and a single receiver, is considered and investigated. In this regard, a hybridization of the BOA and PSO algorithms, named CAHBPSO algorithm, is proposed to solve the considered localization problem with high accuracy, even in the presence of large measurement noise. In the proposed algorithm, an adaptive parameter is introduced to combine global search and local search phases of the BOA algorithm with the PSO algorithm, with the aim to maintain an effective balance between exploration and exploitation abilities during the optimization process. Moreover, to improve convergence and maintain population diversity, chaos theory is incorporated into the inertia weight parameter of the PSO algorithm and an adaptive strategy has been employed to update the value of the sensory fragrance of the BOA. To assess the performance of the discussed techniques, the corresponding CRLB for the passive target localization issue is derived. In addition, a statistical analysis is carried out to compare the proposed CAHBPSO algorithm optimization performance with that of well-known algorithms in the literature on a set of CEC2014 benchmark test problems.
It can be shown from the statistical comparisons between CAHBPSO and other well-known algorithms in the literature that the hybrid method suggested in this paper demonstrates better optimization performance. Furthermore, it is concluded that the proposed CAHBPSO algorithm attains the CRLB accuracy and provides better localization performance compared to the WLS, SDP, BOA, and PSO algorithms. In addition, it is observed that the CAHBPSO algorithm shows lower sensitivity to the variations in network topology and higher localization accuracy under the high measurement noise. Finally, from the analysis of the execution time and computational complexity, it is concluded that the proposed CAHBPSO algorithm provides a proper balance between localization accuracy and complexity compared to other considered algorithms.
Future studies will aim to identify the optimal network architecture for the passive target localization problem in the presence of non-line-of-sight environment. In order to improve the performance of the optimization process, research may also concentrate on performing the analysis of the sensitivity to the parameter changes.

Author Contributions

Conceptualization M.R. and M.S. (Miloš Sedak); methodology and validation M.R.; formal analysis, software, investigation, data curation, M.R. and M.S. (Miloš Sedak); resources M.S. (Mirjana Simić) and P.P.; writing—original draft preparation M.R. and M.S. (Miloš Sedak); writing—review and editing M.S. (Mirjana Simić) and P.P.; visualization M.R.; supervision M.S. (Mirjana Simić) and P.P. All authors have read and agreed to the published version of the manuscript.

Funding

The research of M. Rosić was supported by the Serbian Ministry of Education and Science under Grant TR35029. The research of M. Sedak was supported by the Serbian Ministry of Education and Science under Grant No. TR35006. The work of M. Simić was supported in part by Serbian Ministry of Education and Science under Grant TR32028.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The logarithm of the likelihood function, which is defined in Equation (10), has the following form
ln f r ˜ x = c 1 2 σ 2 r ˜ H g x T r ˜ H g x ,
From Equation (A1), the log-likelihood function’s first partial derivative with respect to x is given as
ln f r ˜ x x = 1 2 σ 2 x r ˜ H g x T r ˜ H g x ,
where
x r ˜ H g x T r ˜ H g x = 2 g x T x H T r ˜ H g x .
Then, substituting Equation (A3) into Equation (A2), the following is obtained
ln f r ˜ x x = 1 σ 2 g x T x H T r ˜ H g x
As a result, taking the partial derivative with regard to x the following is obtained
ln f r ˜ x x = 1 σ 2 g x x T H T r ˜ H g x .
Using the above expressions, the first element on the diagonal of the FIM is obtained as follows
F 11 = E ln f r ˜ x x ln f r ˜ x x T = E 1 σ 4 g x x T H T r ˜ H g x r ˜ H g x T H g x x = 1 σ 4 g x x T H T E r ˜ H g x r ˜ H g x T H g x x = 1 σ 2 g x x T H T H g x x .
The rest of the element of the FIM matrix are obtained in the same way as above, e.g.,
F 12 = F 21 = E ln f r ˜ x x ln f r ˜ x y T = 1 σ 2 g x x T H T H g x y ,
F 22 = E ln f r ˜ x y ln f r ˜ x y T = 1 σ 2 g x y T H T H g x y .
The partial derivative of g x with respect to the components of x may be represented as according to Equation (7), as follows
g ( x ) x = x x r x x r 2 y y r x x r 2 x x 1 t x x 1 t 2 y y 1 t x x 1 t 2 x x N t x x N t 2 y y N t x x N t 2 ,
where g ( x ) / x is the N × 2 Jacobian matrix. The next expression may be readily produced by doing matrix multiplication
F 11 = 1 σ 2 g ( x ) x T H T H g ( x ) x = 1 σ 2 i = 1 N x x r x x r 2 + x x i t x x i t 2 2
The equations for the remaining elements of FIM are produced in the same way
F 12 = F 21 = 1 σ 2 g ( x ) x T H T H g ( x ) y = 1 σ 2 i = 1 N x x r x x r 2 + x x i t x x i t 2 y y r x x r 2 + y y i t x x i t 2
F 22 = 1 σ 2 g ( x ) y T H T H g ( x ) y = 1 σ 2 i = 1 N y y r x x r 2 + y y i t x x i t 2 2
Finally, Equations (51)–(53) of the corresponding elements of FIM are obtained.

References

  1. Sachs, J. Handbook of Ultra-Wideband Short-Range Sensing: Theory, Sensors, Applications; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar] [CrossRef]
  2. Chalise, B.K.; Zhang, Y.D.; Amin, M.G.; Himed, B. Target localization in a multi-static passive radar system through convex optimization. Signal Process. 2014, 102, 207–215. [Google Scholar] [CrossRef]
  3. Deak, G.; Curran, K.; Condell, J. A survey of active and passive indoor localisation systems. Comput. Commun. 2012, 35, 1939–1954. [Google Scholar] [CrossRef]
  4. Shen, L.; Zhang, Q.; Pang, J.; Xu, H.; Li, P.; Xue, D. ANTspin: Efficient Absolute Localization Method of RFID Tags via Spinning Antenna. Sensors 2019, 19, 2194. [Google Scholar] [CrossRef] [Green Version]
  5. Choi, K.H.; Ra, W.S.; Park, S.Y.; Park, J.B. Robust least squares approach to passive target localization using ultrasonic receiver array. IEEE Trans. Ind. Electron. 2014, 61, 1993–2002. [Google Scholar] [CrossRef]
  6. Noroozi, A.; Sebt, M.A. Target localization from bistatic range measurements in multi-transmitter multi-receiver passive radar. IEEE Signal Process. Lett. 2015, 22, 2445–2449. [Google Scholar] [CrossRef]
  7. Wang, G.; Li, Y.; Ansari, N. A semidefinite relaxation method for source localization using TDOA and FDOA measurements. IEEE Trans. Veh. Technol. 2013, 62, 853–862. [Google Scholar] [CrossRef]
  8. Zhang, L.; Liu, L.; Yang, X.S.; Dai, Y. A novel hybrid firefly algorithm for global optimization. PLoS ONE 2016, 11, e0163230. [Google Scholar] [CrossRef] [Green Version]
  9. Yue, Y.; Cao, L.; Hu, J.; Cai, S.; Hang, B.; Wu, H. A Novel Hybrid Location Algorithm Based on Chaotic Particle Swarm Optimization for Mobile Position Estimation. IEEE Access 2019, 7, 58541–58552. [Google Scholar] [CrossRef]
  10. Mohamed, A.W.; Almazyad, A.S. Differential evolution with novel mutation and adaptive crossover strategies for solving large scale global optimization problems. Appl. Comput. Intell. Soft Comput. 2017, 2017, 7974218. [Google Scholar] [CrossRef]
  11. Kuila, P.; Gupta, S.K.; Jana, P.K. A novel evolutionary approach for load balanced clustering problem for wireless sensor networks. Swarm Evol. Comput. 2013, 12, 48–56. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Wang, S.; Ji, G. A comprehensive survey on particle swarm optimization algorithm and its applications. Math. Probl. Eng. 2015, 2015, 931256. [Google Scholar] [CrossRef] [Green Version]
  13. Arora, S.; Singh, S. Butterfly optimization algorithm: A novel approach for global optimization. Soft Comput. 2019, 23, 715–734. [Google Scholar] [CrossRef]
  14. Storn, R.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  15. Cheng, J.; Xia, L. An effective Cuckoo search algorithm for node localization in wireless sensor network. Sensors 2016, 16, 1390. [Google Scholar] [CrossRef] [Green Version]
  16. Arora, S.; Singh, S. An Effective Hybrid Butterfly Optimization Algorithm with Artificial Bee Colony for Numerical Optimization. Int. J. Interact. Multimed. Artif. Intell. 2017, 4, 14–21. [Google Scholar] [CrossRef] [Green Version]
  17. Wu, P.; Su, S.; Zuo, Z.; Guo, X.; Sun, B.; Wen, X. Time Difference of Arrival (TDoA) Localization Combining Weighted Least Squares and Firefly Algorithm. Sensors 2019, 19, 2554. [Google Scholar] [CrossRef] [Green Version]
  18. Arora, S.; Singh, S. Node localization in wireless sensor networks using butterfly optimization algorithm. Arab. J. Sci. Eng. 2017, 42, 3325–3335. [Google Scholar] [CrossRef]
  19. Arora, S.; Anand, P. Binary butterfly optimization approaches for feature selection. Expert Syst. Appl. 2019, 116, 147–160. [Google Scholar] [CrossRef]
  20. Fan, Y.; Shao, J.; Sun, G.; Shao, X. A self-adaption butterfly optimization algorithm for numerical optimization problems. IEEE Access 2020, 8, 88026–88041. [Google Scholar] [CrossRef]
  21. Arora, S.; Singh, S. An improved butterfly optimization algorithm with chaos. J. Intell. Fuzzy Syst. 2017, 32, 1079–1088. [Google Scholar] [CrossRef]
  22. Li, G.; Shuang, F.; Zhao, P.; Le, C. An improved butterfly optimization algorithm for engineering design problems using the cross-entropy method. Symmetry 2019, 11, 1049. [Google Scholar] [CrossRef] [Green Version]
  23. Eberhart, R.; Kennedy, J. A new optimizer using particle swarm theory. In Proceedings of the Sixth International Symposium on Micro Machine and Human Science, Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar] [CrossRef]
  24. Tanweer, M.R.; Suresh, S.; Sundararajan, N. Dynamic mentoring and self-regulation based particle swarm optimization algorithm for solving complex real-world optimization problems. Inf. Sci. 2016, 326, 1–24. [Google Scholar] [CrossRef]
  25. Su, P.; Cai, C.; Song, Y.; Ma, J.; Tan, Q. A Hybrid Diffractive Optical Element Design Algorithm Combining Particle Swarm Optimization and a Simulated Annealing Algorithm. Appl. Sci. 2020, 10, 5485. [Google Scholar] [CrossRef]
  26. Wang, F.; Zhang, H.; Li, K.; Lin, Z.; Yang, J.; Shen, X.L. A hybrid particle swarm optimization algorithm using adaptive learning strategy. Inf. Sci. 2018, 436, 162–177. [Google Scholar] [CrossRef]
  27. Harrison, K.R.; Engelbrecht, A.P.; Ombuki-Berman, B.M. Self-adaptive particle swarm optimization: A review and analysis of convergence. Swarm Intell. 2018, 12, 187–226. [Google Scholar] [CrossRef] [Green Version]
  28. Liu, W.; Wang, Z.; Liu, X.; Zeng, N.; Bell, D. A novel particle swarm optimization approach for patient clustering from emergency departments. IEEE Trans. Evol. Comput. 2018, 23, 632–644. [Google Scholar] [CrossRef] [Green Version]
  29. Tian, D.; Zhao, X.; Shi, Z. Chaotic particle swarm optimization with sigmoid-based acceleration coefficients for numerical function optimization. Swarm Evol. Comput. 2019, 51, 100573. [Google Scholar] [CrossRef]
  30. Yang, D.; Li, G.; Cheng, G. On the efficiency of chaos optimization algorithms for global optimization. Chaos Solitons Fractals 2007, 34, 1366–1375. [Google Scholar] [CrossRef]
  31. Thompson, J.M.T.; Stewart, H.B. Nonlinear Dynamics and Chaos; John Wiley & Sons: Hoboken, NJ, USA, 2002. [Google Scholar]
  32. Tavazoei, M.S.; Haeri, M. Comparison of different one-dimensional maps as chaotic search pattern in chaos optimization algorithms. Appl. Math. Comput. 2007, 187, 1076–1085. [Google Scholar] [CrossRef]
  33. Zhang, M.; Long, D.; Qin, T.; Yang, J. A Chaotic Hybrid Butterfly Optimization Algorithm with Particle Swarm Optimization for High-Dimensional Optimization Problems. Symmetry 2020, 12, 1800. [Google Scholar] [CrossRef]
  34. Aydilek, I.B. A hybrid firefly and particle swarm optimization algorithm for computationally expensive numerical problems. Appl. Soft Comput. 2018, 66, 232–249. [Google Scholar] [CrossRef]
  35. Zhou, H.; Zhang, G.; Wang, X.; Ni, P.; Zhang, J. A hybrid identification method on butterfly optimization and differential evolution algorithm. Smart Struct. Syst. 2020, 26, 345–360. [Google Scholar]
  36. Zhang, X.; Guo, P.; Zhang, H.; Yao, J. Hybrid Particle Swarm Optimization Algorithm for Process Planning. Mathematics 2020, 8, 1745. [Google Scholar] [CrossRef]
  37. Li, W.; Wei, P.; Xiao, X. A robust TDOA-based location method and its performance analysis. Sci. China Ser. Inf. Sci. 2009, 52, 876–882. [Google Scholar] [CrossRef]
  38. Shen, J.; Molisch, A.F.; Salmi, J. Accurate passive location estimation using TOA measurements. IEEE Trans. Wirel. Commun. 2012, 11, 2182–2192. [Google Scholar] [CrossRef]
  39. Bishop, A.N.; Fidan, B.; Anderson, B.D.; Doğançay, K.; Pathirana, P.N. Optimality analysis of sensor-target localization geometries. Automatica 2010, 46, 479–492. [Google Scholar] [CrossRef]
  40. Hu, Y.; Leus, G. Robust differential received signal strength-based localization. IEEE Trans. Signal Process. 2017, 65, 3261–3276. [Google Scholar] [CrossRef] [Green Version]
  41. Xu, S.; Doğançay, K. Optimal sensor placement for 3-D angle-of-arrival target localization. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1196–1211. [Google Scholar] [CrossRef]
  42. Tomic, S.; Beko, M.; Dinis, R.; Bernardo, L. On target localization using combined RSS and AoA measurements. Sensors 2018, 18, 1266. [Google Scholar] [CrossRef] [Green Version]
  43. Jin, B.; Xu, X.; Zhang, T. Robust time-difference-of-arrival (TDOA) localization using weighted least squares with cone tangent plane constraint. Sensors 2018, 18, 778. [Google Scholar] [CrossRef] [Green Version]
  44. Han, G.; Jiang, J.; Zhang, C.; Duong, T.Q.; Guizani, M.; Karagiannidis, G.K. A survey on mobile anchor node assisted localization in wireless sensor networks. IEEE Commun. Surv. Tutor. 2016, 18, 2220–2243. [Google Scholar] [CrossRef]
  45. Xiao, H.; Zhang, H.; Wang, Z.; Gulliver, T.A. An RSSI based DV-hop algorithm for wireless sensor networks. In Proceedings of the 2017 IEEE Pacific Rim Conference on Communications, Computers and Signal Processing (PACRIM), Victoria, BC, Canada, 21–23 August 2017; pp. 1–6. [Google Scholar]
  46. Shen, J.; Wang, A.; Wang, C.; Hung, P.C.; Lai, C.F. An efficient centroid-based routing protocol for energy management in WSN-assisted IoT. IEEE Access 2017, 5, 18469–18479. [Google Scholar] [CrossRef]
  47. Liu, J.; Wang, Z.; Yao, M.; Qiu, Z. VN-APIT: Virtual nodes-based range-free APIT localization scheme for WSN. Wirel. Netw. 2016, 22, 867–878. [Google Scholar] [CrossRef]
  48. Halder, S.; Ghosal, A. A survey on mobile anchor assisted localization techniques in wireless sensor networks. Wirel. Netw. 2016, 22, 2317–2336. [Google Scholar] [CrossRef]
  49. Zekavat, R.; Buehrer, R.M. Handbook of Position Location: Theory, Practice and Advances; John Wiley & Sons: Hoboken, NJ, USA, 2011; Volume 27. [Google Scholar] [CrossRef]
  50. Noroozi, A.; Sebt, M.A. Weighted least squares target location estimation in multi-transmitter multi-receiver passive radar using bistatic range measurements. IET Radar Sonar Navig. 2016, 10, 1088–1097. [Google Scholar] [CrossRef]
  51. Kaur, R.; Arora, S. Nature inspired range based wireless sensor node localization algorithms. Int. J. Interact. Multimed. Artif. Intell. 2017, 4, 7–17. [Google Scholar] [CrossRef] [Green Version]
  52. Rao, S.S. Engineering Optimization: Theory and Practice; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  53. Destino, G.; Abreu, G. On the maximum likelihood approach for source and network localization. IEEE Trans. Signal Process. 2011, 59, 4954–4970. [Google Scholar] [CrossRef]
  54. Kocuk, B.; Dey, S.; Sun, X. Strong SOCP relaxations for the optimal power flow problem. Oper. Res. 2016, 64, 1177–1196. [Google Scholar] [CrossRef] [Green Version]
  55. Biswas, P.; Liang, T.C.; Toh, K.C.; Ye, Y.; Wang, T.C. Semidefinite programming approaches for sensor network localization with noisy distance measurements. IEEE Trans. Autom. Sci. Eng. 2006, 3, 360–371. [Google Scholar] [CrossRef] [Green Version]
  56. Jiang, Y.; Liu, M.; Chen, T.; Gao, L. TDOA Passive Location Based on Cuckoo Search Algorithm. J. Shanghai Jiaotong Univ. Sci. 2018, 23, 368–375. [Google Scholar] [CrossRef]
  57. Cakir, O.; Kaya, I.; Yazgan, A.; Cakir, O.; Tugcu, E. Emitter location finding using particle swarm optimization. Radioengineering 2014, 23, 252–258. [Google Scholar]
  58. Meng, Y.; Zhi, Q.; Zhang, Q.; Yao, N. A Two-Stage Particle Swarm Optimization Algorithm for Wireless Sensor Nodes Localization in Concave Regions. Information 2020, 11, 488. [Google Scholar] [CrossRef]
  59. Gumaida, B.F.; Luo, J. A hybrid particle swarm optimization with a variable neighborhood search for the localization enhancement in wireless sensor networks. Appl. Intell. 2019, 49, 3539–3557. [Google Scholar] [CrossRef]
  60. Harikrishnan, R.; Kumar, V.J.S.; Ponmalar, P.S. A Comparative Analysis of Intelligent Algorithms for Localization in Wireless Sensor Networks. Wirel. Pers. Commun. 2016, 87, 1057–1069. [Google Scholar] [CrossRef]
  61. Qu, X.; Xie, L. An efficient convex constrained weighted least squares source localization algorithm based on TDOA measurements. Signal Process. 2016, 119, 142–152. [Google Scholar] [CrossRef]
  62. Boyd, S.; Boyd, S.P.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  63. Shi, Y.; Eberhart, R. A modified particle swarm optimizer. In Proceedings of the 1998 IEEE International Conference on Evolutionary Computation Proceedings, IEEE World Congress on Computational Intelligence (Cat. No.98TH8360), Anchorage, AK, USA, 4–9 May 1998; pp. 69–73. [Google Scholar]
  64. Zhao, Q.; Li, C.; Zhu, D.; Xie, C. Coverage Optimization of Wireless Sensor Networks Using Combinations of PSO and Chaos Optimization. Electronics 2022, 11, 853. [Google Scholar] [CrossRef]
  65. Tian, D.; Shi, Z. MPSO: Modified particle swarm optimization and its applications. Swarm Evol. Comput. 2018, 41, 49–68. [Google Scholar] [CrossRef]
  66. Liu, W.; Wang, Z.; Yuan, Y.; Zeng, N.; Hone, K.; Liu, X. A novel sigmoid-function-based adaptive weighted particle swarm optimizer. IEEE Trans. Cybern. 2019, 51, 1085–1093. [Google Scholar] [CrossRef]
  67. Tanabe, R.; Fukunaga, A. Success-history based parameter adaptation for differential evolution. In Proceedings of the 2013 IEEE Congress on Evolutionary Computation, Cancun, Mexico, 20–23 June 2013; pp. 71–78. [Google Scholar]
  68. Liang, J.J.; Qu, B.Y.; Suganthan, P.N. Problem Definitions and Evaluation Criteria for the CEC 2014 Special Session and Competition on Single Objective Real-Parameter Numerical Optimization; Technical Report 201311, Computational Intelligence Laboratory, Zhengzhou University, Zhengzhou China and Technical Report; Nanyang Technological University: Singapore, 2013; Volume 635, p. 490. [Google Scholar]
  69. Derrac, J.; García, S.; Molina, D.; Herrera, F. A practical tutorial on the use of nonparametric statistical tests as a methodology for comparing evolutionary and swarm intelligence algorithms. Swarm Evol. Comput. 2011, 1, 3–18. [Google Scholar] [CrossRef]
  70. Assiri, A.S. On the performance improvement of Butterfly Optimization approaches for global optimization and Feature Selection. PLoS ONE 2021, 16, e0242612. [Google Scholar] [CrossRef]
Figure 1. Passive target localization using noisy TDOA measurements.
Figure 1. Passive target localization using noisy TDOA measurements.
Sensors 22 05739 g001
Figure 2. The contour plot of the considered ML objective function J M L x .
Figure 2. The contour plot of the considered ML objective function J M L x .
Sensors 22 05739 g002
Figure 3. The changes of the adaptive sensory fragrance c G during the optimization process.
Figure 3. The changes of the adaptive sensory fragrance c G during the optimization process.
Sensors 22 05739 g003
Figure 4. The the unbounded chaotic behavior produced by logistic chaotic map for the x c ( 0 ) = 0.9 and a ω = 4 .
Figure 4. The the unbounded chaotic behavior produced by logistic chaotic map for the x c ( 0 ) = 0.9 and a ω = 4 .
Sensors 22 05739 g004
Figure 5. The changes of the chaotic dynamic inertia weight ω c G + 1 during the search process.
Figure 5. The changes of the chaotic dynamic inertia weight ω c G + 1 during the search process.
Sensors 22 05739 g005
Figure 6. RMSEs in terms of p for a passive target at x = 20 30 T .
Figure 6. RMSEs in terms of p for a passive target at x = 20 30 T .
Sensors 22 05739 g006
Figure 7. RMSEs in terms of p for a passive target at x = 100 80 T .
Figure 7. RMSEs in terms of p for a passive target at x = 100 80 T .
Sensors 22 05739 g007
Figure 8. RMSEs in relation to measurement noise p for a third simulation scenario.
Figure 8. RMSEs in relation to measurement noise p for a third simulation scenario.
Sensors 22 05739 g008
Figure 9. CDFs of passive target localization error of the considered algorithms for the first simulation scenario.
Figure 9. CDFs of passive target localization error of the considered algorithms for the first simulation scenario.
Sensors 22 05739 g009
Figure 10. The RMSEs of the algorithms under consideration for σ 2 = 1 m 2 , as a function of the number of transmitters.
Figure 10. The RMSEs of the algorithms under consideration for σ 2 = 1 m 2 , as a function of the number of transmitters.
Sensors 22 05739 g010
Table 1. Table presents the objective function values in terms of mean and standard deviation obtained by the considered algorithms, as a result of the numerical simulation performed on 30 CEC2014 test functions.
Table 1. Table presents the objective function values in terms of mean and standard deviation obtained by the considered algorithms, as a result of the numerical simulation performed on 30 CEC2014 test functions.
CAHBPSOSHADEFABOAPSOHPSOBOA
Mean (STD) Sign
f 1 10 2.20 × 10 3 3.48 × 10 3 2.03 × 10 6 1.15 × 10 6 + 2.76 × 10 8 1.66 × 10 8 + 2.76 × 10 6 2.10 × 10 6 + 2.24 × 10 6 2.59 × 10 6 + 2.51 × 10 8 1.22 × 10 8 +
30 4.31 × 10 6 4.41 × 10 6 1.31 × 10 8 4.48 × 10 7 + 3.57 × 10 9 9.93 × 10 8 + 4.69 × 10 8 2.03 × 10 8 + 3.82 × 10 7 3.26 × 10 7 + 1.28 × 10 4 1.03 × 10 4
50 7.81 × 10 6 5.49 × 10 6 2.32 × 10 8 6.41 × 10 7 + 9.04 × 10 9 1.88 × 10 9 + 1.62 × 10 9 5.33 × 10 8 + 5.80 × 10 7 3.57 × 10 7 + 1.24 × 10 10 3.15 × 10 9 +
100 1.53 × 10 8 6.61 × 10 7 8.87 × 10 8 1.61 × 10 8 + 2.13 × 10 10 3.21 × 10 9 + 4.47 × 10 9 9.49 × 10 8 + 3.13 × 10 8 1.25 × 10 8 + 4.74 × 10 5 1.91 × 10 5
f 2 10 3.08 × 10 2 3.81 × 10 2 3.04 × 10 5 1.60 × 10 5 + 1.28 × 10 10 4.14 × 10 9 + 2.78 × 10 8 2.42 × 10 8 + 2.69 × 10 3 3.81 × 10 3 + 1.22 × 10 10 2.28 × 10 9 +
30 3.93 × 10 2 2.15 × 10 1 2.84 × 10 7 9.70 × 10 6 + 1.34 × 10 11 1.58 × 10 10 + 3.76 × 10 10 5.70 × 10 9 + 9.43 × 10 8 8.59 × 10 8 + 9.79 × 10 10 1.49 × 10 10 +
50 3.12 × 10 3 4.40 × 10 3 2.74 × 10 8 9.07 × 10 7 + 2.74 × 10 11 2.32 × 10 10 + 9.79 × 10 10 1.49 × 10 10 + 5.20 × 10 9 3.36 × 10 9 + 1.90 × 10 11 2.41 × 10 9 +
100 2.56 × 10 6 1.81 × 10 7 5.52 × 10 9 1.06 × 10 9 + 6.44 × 10 11 4.03 × 10 10 + 2.36 × 10 11 1.60 × 10 10 + 3.18 × 10 10 7.48 × 10 9 + 2.31 × 10 12 4.27 × 10 12
f 3 10 4.01 × 10 3 1.08 × 10 2 1.13 × 10 3 6.34 × 10 2 + 3.23 × 10 5 1.20 × 10 6 + 8.10 × 10 2 2.90 × 10 2 + 2.24 × 10 3 2.82 × 10 3 + 1.62 × 10 4 2.13 × 10 3 +
30 1.44 × 10 0 1.71 × 10 0 2.45 × 10 4 8.20 × 10 3 + 2.14 × 10 6 6.68 × 10 6 + 3.57 × 10 4 7.26 × 10 3 + 2.23 × 10 4 8.42 × 10 3 + 1.27 × 10 5 1.24 × 10 4 +
50 2.12 × 10 3 1.96 × 10 3 3.00 × 10 5 1.46 × 10 5 + 5.83 × 10 5 6.73 × 10 5 + 1.27 × 10 5 1.24 × 10 4 + 4.73 × 10 4 1.15 × 10 4 + 1.88 × 10 5 4.85 × 10 3 +
100 7.25 × 10 3 3.71 × 10 3 7.03 × 10 5 2.23 × 10 5 + 8.45 × 10 5 1.60 × 10 5 + 2.64 × 10 5 1.62 × 10 4 + 1.45 × 10 5 1.64 × 10 4 + 3.79 × 10 12 5.76 × 10 12
f 4 10 1.15 × 10 1 1.62 × 10 1 3.43 × 10 1 3.25 × 10 0 + 2.58 × 10 3 1.25 × 10 3 + 6.65 × 10 2 3.59 × 10 2 + 3.31 × 10 1 7.85 × 10 0 + 2.18 × 10 1 1.63 × 10 1 +
30 1.56 × 10 2 4.06 × 10 1 1.64 × 10 2 2.02 × 10 1 3.56 × 10 4 9.83 × 10 3 + 8.79 × 10 3 1.30 × 10 3 + 1.76 × 10 2 4.57 × 10 1 2.02 × 10 0 1.07 × 10 1
50 1.78 × 10 2 3.58 × 10 1 2.67 × 10 2 4.55 × 10 1 + 1.15 × 10 5 1.97 × 10 4 + 2.87 × 10 4 4.08 × 10 3 + 5.43 × 10 2 2.22 × 10 2 + 6.47 × 10 4 2.93 × 10 3 +
100 5.65 × 10 2 1.18 × 10 2 1.06 × 10 3 1.07 × 10 2 + 2.85 × 10 5 3.58 × 10 4 + 6.59 × 10 4 7.89 × 10 3 + 2.39 × 10 3 9.13 × 10 2 + 1.55 × 10 2 4.59 × 10 1
f 5 10 1.94 × 10 1 3.95 × 10 0 2.09 × 10 1 1.41 × 10 1 + 2.03 × 10 1 6.68 × 10 2 + 1.99 × 10 1 1.17 × 10 0 + 1.98 × 10 1 2.83 × 10 0 + 4.37 × 10 0 7.92 × 10 0
30 2.08 × 10 1 6.36 × 10 2 2.13 × 10 1 6.99 × 10 2 + 2.10 × 10 1 4.53 × 10 2 + 2.09 × 10 1 5.30 × 10 2 + 2.09 × 10 1 1.08 × 10 1 + 2.03 × 10 1 4.20 × 10 2
50 2.10 × 10 1 6.12 × 10 2 2.14 × 10 1 4.64 × 10 2 + 2.12 × 10 1 3.06 × 10 2 + 2.11 × 10 1 2.76 × 10 2 + 2.11 × 10 1 4.71 × 10 2 + 2.13 × 10 1 3.36 × 10 2 +
100 2.13 × 10 1 3.48 × 10 2 2.15 × 10 1 3.22 × 10 2 + 2.13 × 10 1 2.73 × 10 2 + 2.13 × 10 1 2.37 × 10 2 + 2.13 × 10 1 2.50 × 10 2 + 2.07 × 10 1 4.29 × 10 2
f 6 10 7.93 × 10 1 9.58 × 10 1 6.27 × 10 0 1.81 × 10 0 + 1.30 × 10 1 9.28 × 10 1 + 4.55 × 10 0 4.65 × 10 1 + 1.44 × 10 0 8.74 × 10 1 + 3.55 × 10 1 6.06 × 10 1
30 9.37 × 10 0 2.44 × 10 0 3.86 × 10 1 3.11 × 10 0 + 4.83 × 10 1 1.42 × 10 0 + 2.98 × 10 1 1.43 × 10 0 + 1.16 × 10 1 2.34 × 10 0 + 2.32 × 10 1 4.19 × 10 0
50 2.52 × 10 1 3.75 × 10 0 6.80 × 10 1 5.33 × 10 0 + 8.41 × 10 1 1.77 × 10 0 + 5.89 × 10 1 2.00 × 10 0 + 2.65 × 10 1 3.96 × 10 0 2.32 × 10 1 4.19 × 10 0
100 7.52 × 10 1 5.97 × 10 0 1.53 × 10 2 7.13 × 10 0 + 1.77 × 10 2 2.32 × 10 0 + 1.40 × 10 2 3.28 × 10 0 + 8.03 × 10 1 5.06 × 10 0 + 7.02 × 10 1 1.19 × 10 1
f 7 10 1.04 × 10 1 4.43 × 10 2 9.22 × 10 1 8.92 × 10 2 + 1.96 × 10 2 4.95 × 10 1 + 9.20 × 10 1 3.01 × 10 1 + 4.78 × 10 1 5.82 × 10 1 + 7.46 × 10 2 3.99 × 10 2
30 1.79 × 10 2 2.13 × 10 2 1.22 × 10 0 6.78 × 10 2 + 1.16 × 10 3 1.59 × 10 2 + 5.98 × 10 2 6.79 × 10 1 + 8.70 × 10 0 9.21 × 10 0 + 1.82 × 10 2 3.05 × 10 2 +
50 1.03 × 10 2 1.22 × 10 2 3.72 × 10 0 7.66 × 10 1 + 2.68 × 10 3 2.70 × 10 2 + 1.27 × 10 3 8.60 × 10 1 + 4.60 × 10 1 3.21 × 10 1 + 1.82 × 10 2 3.05 × 10 2
100 3.81 × 10 3 7.32 × 10 3 5.38 × 10 1 1.27 × 10 1 + 5.80 × 10 3 3.77 × 10 2 + 2.78 × 10 3 9.78 × 10 1 + 2.73 × 10 2 6.99 × 10 1 + 1.01 × 10 1 3.33 × 10 1 +
f 8 10 5.85 × 10 2 2.36 × 10 1 3.84 × 10 1 8.60 × 10 0 + 1.15 × 10 2 1.45 × 10 1 + 3.98 × 10 1 6.06 × 10 0 + 5.53 × 10 0 3.07 × 10 0 + 9.48 × 10 1 5.62 × 10 0 +
30 1.27 × 10 1 3.29 × 10 0 2.24 × 10 2 2.01 × 10 1 + 4.97 × 10 2 3.56 × 10 1 + 2.70 × 10 2 1.59 × 10 1 + 6.29 × 10 1 1.35 × 10 1 + 2.50 × 10 1 2.23 × 10 1
50 3.86 × 10 1 6.75 × 10 0 4.38 × 10 2 2.91 × 10 1 + 9.32 × 10 2 4.80 × 10 1 + 5.51 × 10 2 2.24 × 10 1 + 1.73 × 10 2 3.13 × 10 1 + 2.50 × 10 1 2.23 × 10 1
100 1.37 × 10 2 2.47 × 10 1 1.02 × 10 3 4.43 × 10 1 + 2.05 × 10 3 5.06 × 10 1 + 1.27 × 10 3 2.83 × 10 1 + 5.56 × 10 2 4.88 × 10 1 + 2.12 × 10 2 5.37 × 10 1 +
f 9 10 4.48 × 10 0 2.00 × 10 0 5.79 × 10 1 1.05 × 10 1 + 1.20 × 10 2 1.37 × 10 1 + 3.87 × 10 1 6.38 × 10 0 + 1.04 × 10 1 5.22 × 10 0 + 4.86 × 10 1 8.27 × 10 1
30 6.19 × 10 1 1.59 × 10 1 2.69 × 10 2 2.42 × 10 1 + 6.10 × 10 2 4.43 × 10 1 + 2.90 × 10 2 1.80 × 10 1 + 8.12 × 10 1 2.19 × 10 1 + 8.92 × 10 1 2.34 × 10 1
50 1.34 × 10 2 2.48 × 10 1 5.02 × 10 2 3.41 × 10 1 + 1.16 × 10 3 6.09 × 10 1 + 6.10 × 10 2 2.70 × 10 1 + 1.82 × 10 2 3.17 × 10 1 + 8.92 × 10 1 2.34 × 10 1
100 3.43 × 10 2 5.85 × 10 1 1.13 × 10 3 5.30 × 10 1 + 2.49 × 10 3 9.12 × 10 1 + 1.38 × 10 3 4.42 × 10 1 + 5.63 × 10 2 6.26 × 10 1 + 3.47 × 10 2 6.33 × 10 1 +
f 10 10 5.69 × 10 1 6.45 × 10 1 1.58 × 10 3 3.08 × 10 2 + 2.04 × 10 3 1.45 × 10 2 + 9.91 × 10 2 1.20 × 10 2 + 2.16 × 10 2 1.29 × 10 2 + 4.42 × 10 1 9.31 × 10 1
30 5.05 × 10 2 2.49 × 10 2 7.46 × 10 3 5.76 × 10 2 + 8.61 × 10 3 2.09 × 10 2 + 6.59 × 10 3 2.85 × 10 2 + 1.97 × 10 3 4.69 × 10 2 + 1.49 × 10 0 1.15 × 10 0
50 1.10 × 10 3 3.56 × 10 2 1.39 × 10 4 8.79 × 10 2 + 1.55 × 10 4 1.87 × 10 2 + 1.29 × 10 4 3.37 × 10 2 + 5.28 × 10 3 7.33 × 10 2 + 1.49 × 10 0 1.15 × 10 0
100 3.96 × 10 3 7.39 × 10 2 3.18 × 10 4 1.11 × 10 3 + 3.39 × 10 4 3.52 × 10 2 + 3.02 × 10 4 6.48 × 10 2 + 1.56 × 10 4 9.18 × 10 2 + 1.60 × 10 0 1.02 × 10 0
f 11 10 1.79 × 10 2 1.32 × 10 2 2.14 × 10 3 2.62 × 10 2 + 2.07 × 10 3 9.23 × 10 1 + 1.07 × 10 3 1.36 × 10 2 + 3.26 × 10 2 1.67 × 10 2 + 1.30 × 10 1 1.13 × 10 1
30 2.54 × 10 3 5.97 × 10 2 9.01 × 10 3 4.31 × 10 2 + 8.55 × 10 3 1.95 × 10 2 + 6.93 × 10 3 2.24 × 10 2 + 2.55 × 10 3 7.53 × 10 2 3.88 × 10 3 3.89 × 10 2
50 5.68 × 10 3 1.22 × 10 3 1.62 × 10 4 7.22 × 10 2 + 1.54 × 10 4 2.35 × 10 2 + 1.32 × 10 4 4.71 × 10 2 + 5.65 × 10 3 7.51 × 10 2 3.88 × 10 3 3.89 × 10 2
100 1.64 × 10 4 4.62 × 10 3 3.47 × 10 4 9.68 × 10 2 + 3.36 × 10 4 3.37 × 10 2 + 3.07 × 10 4 6.24 × 10 2 + 1.38 × 10 4 2.76 × 10 3 1.23 × 10 4 1.05 × 10 3
f 12 10 1.05 × 10 1 1.24 × 10 1 3.39 × 10 0 1.02 × 10 0 + 1.08 × 10 0 1.77 × 10 1 + 9.85 × 10 1 1.31 × 10 1 + 3.31 × 10 1 3.61 × 10 1 + 1.12 × 10 1 2.90 × 10 2 +
30 1.42 × 10 0 6.26 × 10 1 5.74 × 10 0 1.08 × 10 0 + 2.54 × 10 0 2.68 × 10 1 + 2.40 × 10 0 2.40 × 10 1 + 1.46 × 10 0 1.07 × 10 0 3.08 × 10 1 4.45 × 10 2
50 2.35 × 10 0 5.58 × 10 1 6.91 × 10 0 8.13 × 10 1 + 3.50 × 10 0 2.46 × 10 1 + 3.41 × 10 0 3.14 × 10 1 + 1.82 × 10 0 1.57 × 10 0 3.08 × 10 1 4.45 × 10 2
100 3.48 × 10 0 5.31 × 10 1 6.63 × 10 0 6.11 × 10 1 + 4.24 × 10 0 2.40 × 10 1 + 4.10 × 10 0 2.60 × 10 1 + 3.35 × 10 0 1.48 × 10 0 5.84 × 10 1 6.63 × 10 2
f 13 10 1.06 × 10 1 4.92 × 10 2 5.55 × 10 1 1.42 × 10 1 + 5.03 × 10 0 1.15 × 10 0 + 2.74 × 10 0 5.86 × 10 1 + 1.13 × 10 1 4.66 × 10 2 6.02 × 10 2 1.36 × 10 2
30 3.69 × 10 1 9.46 × 10 2 8.88 × 10 1 1.76 × 10 1 + 1.05 × 10 1 1.03 × 10 0 + 7.10 × 10 0 5.29 × 10 1 + 3.21 × 10 1 2.00 × 10 1 3.90 × 10 1 6.43 × 10 2
50 5.93 × 10 1 9.93 × 10 2 1.09 × 10 0 1.85 × 10 1 + 1.22 × 10 1 9.85 × 10 1 + 7.85 × 10 0 3.08 × 10 1 + 6.10 × 10 1 2.55 × 10 1 3.90 × 10 1 6.43 × 10 2
100 6.60 × 10 1 7.34 × 10 2 1.15 × 10 0 1.63 × 10 1 + 1.43 × 10 1 6.36 × 10 1 + 9.10 × 10 0 1.84 × 10 1 + 1.42 × 10 0 9.73 × 10 1 + 5.28 × 10 1 5.81 × 10 2
f 14 10 4.94 × 10 2 1.78 × 10 2 6.00 × 10 1 1.53 × 10 1 + 5.80 × 10 1 1.26 × 10 1 + 2.06 × 10 1 5.32 × 10 0 + 1.44 × 10 1 1.62 × 10 1 + 3.77 × 10 2 1.37 × 10 2
30 2.87 × 10 1 1.55 × 10 1 1.13 × 10 0 4.16 × 10 1 + 4.23 × 10 2 5.65 × 10 1 + 2.23 × 10 2 2.31 × 10 1 + 7.61 × 10 1 1.44 × 10 0 + 2.99 × 10 1 7.19 × 10 2
50 5.78 × 10 1 2.97 × 10 1 1.38 × 10 0 5.45 × 10 1 + 6.91 × 10 2 6.99 × 10 1 + 3.20 × 10 2 2.64 × 10 1 + 1.03 × 10 1 9.86 × 10 0 + 2.99 × 10 1 7.19 × 10 2
100 3.99 × 10 1 1.77 × 10 1 4.55 × 10 0 2.57 × 10 0 + 1.72 × 10 3 1.19 × 10 2 + 8.29 × 10 2 3.69 × 10 1 + 8.44 × 10 1 2.16 × 10 1 + 3.29 × 10 1 3.16 × 10 2
f 15 10 7.17 × 10 1 2.44 × 10 1 5.65 × 10 0 9.89 × 10 1 + 1.87 × 10 5 1.45 × 10 5 + 2.01 × 10 2 2.15 × 10 2 + 1.00 × 10 0 4.28 × 10 1 + 4.86 × 10 4 2.37 × 10 4 +
30 6.15 × 10 0 1.86 × 10 0 2.65 × 10 1 2.57 × 10 0 + 2.09 × 10 7 1.29 × 10 7 + 4.64 × 10 4 2.82 × 10 4 + 1.48 × 10 1 1.52 × 10 1 + 1.48 × 10 7 3.72 × 10 6 +
50 1.62 × 10 1 4.85 × 10 0 5.92 × 10 1 8.10 × 10 0 + 1.55 × 10 8 6.36 × 10 7 + 8.49 × 10 5 3.97 × 10 5 + 4.91 × 10 2 8.53 × 10 2 + 1.48 × 10 7 3.72 × 10 6 +
100 7.03 × 10 1 1.39 × 10 1 2.12 × 10 3 1.55 × 10 3 + 6.26 × 10 8 1.85 × 10 8 + 7.84 × 10 6 2.61 × 10 6 + 1.38 × 10 4 9.95 × 10 3 + 4.46 × 10 7 5.50 × 10 6 +
f 16 10 1.42 × 10 0 7.03 × 10 1 4.15 × 10 0 2.00 × 10 1 + 3.96 × 10 0 1.45 × 10 1 + 2.80 × 10 0 2.30 × 10 1 + 1.90 × 10 0 5.79 × 10 1 + 3.70 × 10 0 3.63 × 10 2 +
30 1.04 × 10 1 7.71 × 10 1 1.40 × 10 1 2.83 × 10 1 + 1.35 × 10 1 1.50 × 10 1 + 1.25 × 10 1 1.82 × 10 1 + 1.03 × 10 1 7.37 × 10 1 2.32 × 10 1 5.31 × 10 2 +
50 2.05 × 10 1 8.38 × 10 1 2.39 × 10 1 2.54 × 10 1 + 2.32 × 10 1 1.47 × 10 1 + 2.22 × 10 1 1.81 × 10 1 + 1.95 × 10 1 9.41 × 10 1 2.32 × 10 1 5.31 × 10 2 +
100 4.56 × 10 1 4.68 × 10 1 4.85 × 10 1 3.34 × 10 1 + 4.75 × 10 1 1.77 × 10 1 + 4.64 × 10 1 2.31 × 10 1 + 4.35 × 10 1 1.77 × 10 0 4.76 × 10 1 1.65 × 10 1 +
f 17 10 1.12 × 10 3 1.12 × 10 3 4.78 × 10 4 5.45 × 10 4 + 9.50 × 10 6 1.06 × 10 7 + 1.28 × 10 3 3.00 × 10 2 + 4.25 × 10 3 3.01 × 10 3 + 6.57 × 10 5 1.29 × 10 5 +
30 3.21 × 10 5 2.55 × 10 5 1.36 × 10 7 7.33 × 10 6 + 3.00 × 10 8 1.33 × 10 8 + 4.46 × 10 6 4.14 × 10 6 + 1.31 × 10 6 1.46 × 10 6 + 8.49 × 10 8 2.51 × 10 7 +
50 1.17 × 10 6 9.36 × 10 5 4.56 × 10 7 1.86 × 10 7 + 1.17 × 10 9 3.51 × 10 8 + 7.92 × 10 7 5.10 × 10 7 + 2.93 × 10 6 2.81 × 10 6 + 8.49 × 10 8 2.51 × 10 7 +
100 1.31 × 10 7 8.27 × 10 6 1.89 × 10 8 6.03 × 10 7 + 3.53 × 10 9 8.60 × 10 8 + 6.68 × 10 8 2.88 × 10 8 + 2.38 × 10 7 1.17 × 10 7 + 2.86 × 10 9 4.30 × 10 8 +
f 18 10 3.56 × 10 3 4.46 × 10 3 1.92 × 10 3 2.70 × 10 3 1.06 × 10 8 1.15 × 10 8 + 4.14 × 10 2 1.99 × 10 2 4.40 × 10 3 5.30 × 10 3 1.96 × 10 6 2.80 × 10 6 +
30 5.78 × 10 3 5.96 × 10 3 4.64 × 10 6 5.29 × 10 6 + 9.29 × 10 9 2.74 × 10 9 + 1.54 × 10 8 1.99 × 10 8 + 2.27 × 10 6 9.20 × 10 6 + 2.93 × 10 10 4.38 × 10 9 +
50 1.06 × 10 3 1.23 × 10 3 3.44 × 10 6 2.13 × 10 6 + 2.90 × 10 10 6.79 × 10 9 + 4.80 × 10 9 2.14 × 10 9 + 3.86 × 10 7 9.42 × 10 7 + 2.93 × 10 10 4.38 × 10 9 +
100 1.46 × 10 4 8.49 × 10 4 1.66 × 10 7 6.25 × 10 6 + 8.05 × 10 10 1.20 × 10 10 + 2.19 × 10 10 4.77 × 10 9 + 5.14 × 10 8 4.23 × 10 8 + 4.90 × 10 10 2.59 × 10 9 +
f 19 10 1.06 × 10 0 8.49 × 10 1 3.98 × 10 0 1.11 × 10 0 + 7.35 × 10 1 4.40 × 10 1 + 8.00 × 10 0 3.21 × 10 0 + 2.06 × 10 0 7.88 × 10 1 + 9.37 × 10 1 3.80 × 10 1 +
30 6.80 × 10 0 1.77 × 10 0 1.53 × 10 1 1.88 × 10 0 + 1.22 × 10 3 4.46 × 10 2 + 3.48 × 10 2 9.00 × 10 1 + 2.56 × 10 1 2.02 × 10 1 + 5.27 × 10 3 1.16 × 10 3 +
50 6.35 × 10 1 1.57 × 10 1 6.74 × 10 1 9.25 × 10 0 4.78 × 10 3 1.42 × 10 3 + 1.43 × 10 3 4.46 × 10 2 + 6.87 × 10 1 2.55 × 10 1 5.27 × 10 3 1.16 × 10 3 +
100 1.64 × 10 2 1.74 × 10 1 1.83 × 10 2 1.76 × 10 1 + 2.14 × 10 4 3.77 × 10 3 + 6.46 × 10 3 1.10 × 10 3 + 2.92 × 10 2 6.55 × 10 1 + 1.55 × 10 4 7.68 × 10 2 +
f 20 10 3.21 × 10 0 2.73 × 10 0 3.52 × 10 1 1.46 × 10 1 + 1.88 × 10 6 4.62 × 10 6 + 5.04 × 10 2 2.79 × 10 2 + 9.70 × 10 2 2.18 × 10 3 + 8.08 × 10 4 2.04 × 10 4 +
30 1.30 × 10 2 5.76 × 10 1 9.76 × 10 4 9.60 × 10 4 + 9.78 × 10 6 9.55 × 10 6 + 3.92 × 10 4 1.60 × 10 4 + 8.39 × 10 3 5.39 × 10 3 + 4.21 × 10 5 6.46 × 10 4 +
50 9.05 × 10 2 4.29 × 10 2 1.25 × 10 6 1.63 × 10 6 + 1.46 × 10 7 1.28 × 10 7 + 6.84 × 10 4 2.24 × 10 4 + 1.35 × 10 4 6.09 × 10 3 + 4.21 × 10 5 6.46 × 10 4 +
100 8.85 × 10 3 2.95 × 10 3 2.75 × 10 6 2.31 × 10 6 + 4.89 × 10 7 3.80 × 10 7 + 2.90 × 10 5 7.67 × 10 4 + 4.86 × 10 4 1.74 × 10 4 + 1.36 × 10 7 6.91 × 10 6 +
f 21 10 4.75 × 10 1 5.69 × 10 1 2.10 × 10 3 6.44 × 10 3 + 2.04 × 10 6 2.17 × 10 6 + 3.56 × 10 3 1.43 × 10 3 + 4.26 × 10 3 3.91 × 10 3 + 2.25 × 10 6 1.29 × 10 6 +
30 6.37 × 10 4 4.69 × 10 4 4.79 × 10 6 3.27 × 10 6 + 1.47 × 10 8 8.24 × 10 7 + 4.54 × 10 5 3.22 × 10 5 + 3.48 × 10 5 4.60 × 10 5 + 3.12 × 10 8 4.04 × 10 6 +
50 7.30 × 10 5 7.94 × 10 5 2.81 × 10 7 1.68 × 10 7 + 4.74 × 10 8 2.03 × 10 8 + 4.34 × 10 6 2.96 × 10 6 + 1.82 × 10 6 1.76 × 10 6 + 3.12 × 10 8 4.04 × 10 6 +
100 5.74 × 10 6 4.67 × 10 6 1.23 × 10 8 4.39 × 10 7 + 1.84 × 10 9 5.25 × 10 8 + 1.57 × 10 8 6.71 × 10 7 + 9.69 × 10 6 4.68 × 10 6 + 8.35 × 10 8 1.43 × 10 8 +
f 22 10 6.25 × 10 0 9.20 × 10 0 1.43 × 10 2 1.07 × 10 2 + 5.06 × 10 2 1.30 × 10 2 + 5.66 × 10 1 1.51 × 10 1 + 6.11 × 10 1 5.32 × 10 1 + 6.45 × 10 2 1.99 × 10 2
30 2.14 × 10 2 1.41 × 10 2 1.36 × 10 3 2.33 × 10 2 + 1.65 × 10 4 2.11 × 10 4 + 1.46 × 10 3 2.64 × 10 2 + 3.34 × 10 2 1.61 × 10 2 + 2.87 × 10 6 1.22 × 10 6 +
50 7.56 × 10 2 2.46 × 10 2 2.76 × 10 3 3.29 × 10 2 + 6.81 × 10 5 6.33 × 10 5 + 6.27 × 10 3 4.31 × 10 3 + 7.20 × 10 2 3.24 × 10 2 2.87 × 10 6 1.22 × 10 6 +
100 2.25 × 10 3 4.20 × 10 2 6.24 × 10 3 5.61 × 10 2 + 2.18 × 10 6 1.11 × 10 6 + 2.57 × 10 4 2.61 × 10 4 + 1.80 × 10 3 6.15 × 10 2 5.47 × 10 5 7.45 × 10 4 +
f 23 10 7.21 × 10 3 7.13 × 10 12 7.21 × 10 3 5.60 × 10 3 + 6.94 × 10 3 1.27 × 10 2 + 7.34 × 10 3 4.51 × 10 0 7.21 × 10 3 2.83 × 10 0 + 5.33 × 10 2 2.91 × 10 4 +
30 7.22 × 10 3 1.09 × 10 1 7.22 × 10 3 8.39 × 10 1 + 5.62 × 10 3 4.50 × 10 2 + 7.34 × 10 3 9.19 × 10 12 7.21 × 10 3 7.46 × 10 0 + 5.21 × 10 3 1.39 × 10 3
50 7.19 × 10 3 4.79 × 10 1 7.19 × 10 3 1.48 × 10 0 + 3.48 × 10 3 6.79 × 10 2 + 7.34 × 10 3 9.19 × 10 12 7.12 × 10 3 3.20 × 10 1 + 5.21 × 10 3 1.39 × 10 3
100 7.18 × 10 3 2.18 × 10 0 7.10 × 10 3 1.43 × 10 1 + 1.40 × 10 3 1.27 × 10 3 + 7.34 × 10 3 9.19 × 10 12 6.95 × 10 3 6.80 × 10 1 + 7.34 × 10 3 1.73 × 10 3
f 24 10 2.04 × 10 3 4.00 × 10 0 1.99 × 10 3 8.36 × 10 0 + 1.91 × 10 3 1.41 × 10 1 + 2.02 × 10 3 6.91 × 10 0 + 2.03 × 10 3 1.32 × 10 1 + 1.01 × 10 2 8.63 × 10 1 +
30 1.92 × 10 3 5.73 × 10 0 1.91 × 10 3 3.10 × 10 0 + 1.61 × 10 3 3.12 × 10 1 + 1.95 × 10 3 6.89 × 10 13 1.95 × 10 3 2.30 × 10 4 9.60 × 10 2 2.03 × 10 3
50 1.88 × 10 3 3.87 × 10 0 1.86 × 10 3 3.59 × 10 0 + 1.24 × 10 3 5.65 × 10 1 + 1.95 × 10 3 6.89 × 10 13 1.95 × 10 3 2.51 × 10 4 9.60 × 10 2 2.03 × 10 3
100 1.76 × 10 3 5.60 × 10 0 1.69 × 10 3 6.59 × 10 0 + 1.89 × 10 2 9.61 × 10 1 + 1.95 × 10 3 6.89 × 10 13 1.95 × 10 3 5.76 × 10 4 1.95 × 10 3 3.90 × 10 3
f 25 10 1.43 × 10 3 4.15 × 10 1 1.40 × 10 3 8.76 × 10 0 + 1.38 × 10 3 8.08 × 10 0 + 1.42 × 10 3 1.50 × 10 1 1.41 × 10 3 2.51 × 10 1 1.22 × 10 2 2.78 × 10 2 +
30 1.40 × 10 3 1.70 × 10 0 1.37 × 10 3 6.07 × 10 0 + 1.19 × 10 3 6.57 × 10 1 + 1.40 × 10 3 9.19 × 10 13 1.39 × 10 3 3.52 × 10 0 + 8.62 × 10 2 3.14 × 10 5
50 1.39 × 10 3 3.60 × 10 0 1.34 × 10 3 9.93 × 10 0 + 8.52 × 10 2 1.04 × 10 2 + 1.40 × 10 3 9.19 × 10 13 1.38 × 10 3 5.53 × 10 0 + 8.62 × 10 2 3.14 × 10 5
100 1.33 × 10 3 1.06 × 10 1 1.25 × 10 3 1.55 × 10 1 + 2.21 × 10 2 2.28 × 10 2 + 1.40 × 10 3 9.19 × 10 13 1.40 × 10 3 2.50 × 10 1 1.40 × 10 3 7.49 × 10 5
f 26 10 2.72 × 10 3 4.04 × 10 2 2.72 × 10 3 1.66 × 10 1 + 2.72 × 10 3 2.06 × 10 0 + 2.72 × 10 3 7.73 × 10 2 + 2.72 × 10 3 3.21 × 10 2 7.87 × 10 2 6.90 × 10 1
30 2.72 × 10 3 3.61 × 10 1 2.72 × 10 3 1.49 × 10 1 + 2.47 × 10 3 1.23 × 10 2 + 2.71 × 10 3 2.08 × 10 1 + 2.70 × 10 3 4.27 × 10 1 2.78 × 10 3 2.42 × 10 1 +
50 2.65 × 10 3 8.53 × 10 1 2.65 × 10 3 6.42 × 10 1 2.11 × 10 3 1.77 × 10 2 + 2.64 × 10 3 3.19 × 10 1 2.64 × 10 3 4.85 × 10 1 2.78 × 10 3 2.42 × 10 1 +
100 2.61 × 10 3 5.90 × 10 1 2.58 × 10 3 1.50 × 10 1 + 1.19 × 10 3 2.11 × 10 2 + 2.62 × 10 3 2.30 × 10 12 2.62 × 10 3 2.02 × 10 12 2.62 × 10 3 2.39 × 10 7
f 27 10 1.69 × 10 4 1.70 × 10 2 1.69 × 10 4 1.96 × 10 2 1.65 × 10 4 1.18 × 10 2 + 1.72 × 10 4 2.52 × 10 0 1.69 × 10 4 1.53 × 10 2 3.18 × 10 3 2.23 × 10 2
30 1.67 × 10 4 9.99 × 10 1 1.64 × 10 4 1.77 × 10 2 + 1.53 × 10 4 1.58 × 10 2 + 1.67 × 10 4 1.57 × 10 1 1.66 × 10 4 1.07 × 10 2 1.10 × 10 4 2.59 × 10 2
50 1.62 × 10 4 1.07 × 10 2 1.56 × 10 4 1.43 × 10 2 + 1.39 × 10 4 3.77 × 10 2 + 1.65 × 10 4 1.11 × 10 2 1.62 × 10 4 9.68 × 10 1 1.10 × 10 4 2.59 × 10 2
100 1.49 × 10 4 1.80 × 10 2 1.41 × 10 4 2.32 × 10 2 + 9.71 × 10 3 1.03 × 10 3 + 1.50 × 10 4 5.29 × 10 2 1.49 × 10 4 1.56 × 10 2 1.70 × 10 4 7.29 × 10 2
f 28 10 5.79 × 10 4 5.90 × 10 1 5.79 × 10 4 8.59 × 10 1 5.67 × 10 4 2.65 × 10 2 + 5.79 × 10 4 6.97 × 10 1 5.79 × 10 4 7.17 × 10 1 6.26 × 10 3 3.34 × 10 2
30 5.73 × 10 4 2.82 × 10 2 5.71 × 10 4 8.32 × 10 1 + 5.02 × 10 4 7.77 × 10 2 + 5.70 × 10 4 1.43 × 10 2 + 5.74 × 10 4 1.43 × 10 2 2.43 × 10 4 9.75 × 10 2
50 5.65 × 10 4 4.44 × 10 2 5.59 × 10 4 1.00 × 10 3 + 4.17 × 10 4 1.52 × 10 3 + 5.52 × 10 4 5.24 × 10 2 + 5.65 × 10 4 4.68 × 10 2 2.43 × 10 4 9.75 × 10 2
100 5.32 × 10 4 1.21 × 10 3 5.00 × 10 4 2.77 × 10 3 + 1.79 × 10 4 2.21 × 10 3 + 4.84 × 10 4 1.62 × 10 3 + 5.21 × 10 4 9.93 × 10 2 + 5.82 × 10 4 1.43 × 10 1
f 29 10 2.71 × 10 10 5.97 × 10 5 2.71 × 10 10 2.84 × 10 4 + 2.71 × 10 10 1.98 × 10 7 + 2.71 × 10 10 5.55 × 10 2 2.71 × 10 10 4.15 × 10 5 1.02 × 10 9 8.96 × 10 3 +
30 2.70 × 10 10 3.03 × 10 2 2.70 × 10 10 8.16 × 10 4 + 2.60 × 10 10 2.33 × 10 8 + 2.70 × 10 10 0.00 × 10 0 2.70 × 10 10 4.07 × 10 5 + 7.37 × 10 9 1.15 × 10 5 +
50 2.70 × 10 10 1.54 × 10 7 2.70 × 10 10 5.95 × 10 5 + 2.40 × 10 10 4.96 × 10 8 + 2.70 × 10 10 0.00 × 10 0 2.70 × 10 10 2.92 × 10 6 + 7.37 × 10 9 1.15 × 10 5 +
100 2.71 × 10 10 1.05 × 10 3 2.71 × 10 10 1.92 × 10 6 + 1.82 × 10 10 1.13 × 10 9 + 2.71 × 10 10 0.00 × 10 0 2.70 × 10 10 2.40 × 10 7 + 2.71 × 10 10 1.27 × 10 5 +
f 30 10 3.02 × 10 9 1.73 × 10 2 3.02 × 10 9 4.99 × 10 2 + 3.02 × 10 9 1.91 × 10 5 + 3.02 × 10 9 2.77 × 10 2 + 3.02 × 10 9 4.53 × 10 2 + 4.40 × 10 8 3.37 × 10 4 +
30 3.00 × 10 9 1.15 × 10 3 3.00 × 10 9 1.48 × 10 4 + 3.00 × 10 9 4.25 × 10 6 + 3.00 × 10 9 9.63 × 10 7 3.00 × 10 9 1.72 × 10 4 + 6.15 × 10 8 3.91 × 10 3
50 3.00 × 10 9 4.67 × 10 3 3.00 × 10 9 7.08 × 10 4 + 3.00 × 10 9 1.68 × 10 7 + 3.00 × 10 9 9.63 × 10 7 3.00 × 10 9 4.38 × 10 4 + 6.15 × 10 8 3.91 × 10 3
100 3.02 × 10 9 3.52 × 10 4 3.02 × 10 9 1.23 × 10 6 + 2.51 × 10 9 1.44 × 10 8 + 3.02 × 10 9 9.63 × 10 7 3.02 × 10 9 9.39 × 10 5 + 3.02 × 10 9 2.04 × 10 4
Table 2. Wilcoxon test results for D = 10 , 30 , 50 , and 100 obtained using CAHBPSO and other evaluated algorithms on a set of CEC2014 benchmark problems.
Table 2. Wilcoxon test results for D = 10 , 30 , 50 , and 100 obtained using CAHBPSO and other evaluated algorithms on a set of CEC2014 benchmark problems.
DAlgorithms R + R p Value+Dec.
10CAHBPSO versus SHADE39273 6.08 × 10 4 2721+
CAHBPSO versus FA4650 1.86 × 10 9 3000+
CAHBPSO versus BOA37689 2.37 × 10 3 2433+
CAHBPSO versus PSO41946 3.45 × 10 5 2361+
CAHBPSO versus HPSOBOA330153 4.49 × 10 2 18120+
30CAHBPSO versus SHADE4578 4.66 × 10 8 2910+
CAHBPSO versus FA4650 1.86 × 10 9 3000+
CAHBPSO versus BOA39669 4.18 × 10 4 2406+
CAHBPSO versus PSO42441 1.82 × 10 5 2253+
CAHBPSO versus HPSOBOA234231 9.84 × 10 1 13170
50CAHBPSO versus SHADE44124 1.42 × 10 6 2910+
CAHBPSO versus FA4650 1.86 × 10 9 3000+
CAHBPSO versus BOA38778 9.52 × 10 4 2406+
CAHBPSO versus PSO39174 6.67 × 10 4 2073+
CAHBPSO versus HPSOBOA184181 3.28 × 10 1 15150
100CAHBPSO versus SHADE4650 1.86 × 10 9 3000+
CAHBPSO versus FA4650 1.86 × 10 9 3000+
CAHBPSO versus BOA38976 7.98 × 10 4 2316+
CAHBPSO versus PSO39867 3.45 × 10 4 2217+
CAHBPSO versus HPSOBOA238227 9.19 × 10 1 13170
Table 3. The obtained average ranks using Friedman test for all investigated algorithms across all functions and dimensions using CEC2014, with α = 0.05 .
Table 3. The obtained average ranks using Friedman test for all investigated algorithms across all functions and dimensions using CEC2014, with α = 0.05 .
Algorithm 10 D 30 D 50 D 100 D Mean RankingRank
CAHBPSO1.772.071.972.071.971
HPSOBOA3.402.573.002.472.862
PSO2.972.972.873.002.953
BOA3.303.603.503.573.494
SHADE3.834.033.974.073.985
FA5.735.775.705.835.766
Friedman p value 3.36 × 10 14 1.44 × 10 14 1.58 × 10 13 1.84 × 10 15
Table 4. The computation times on average for the methods under consideration.
Table 4. The computation times on average for the methods under consideration.
 SDPPSOBOACAHBPSO
Scenario 1 10.03 3.02 × 10 2 1.46 × 10 2 2.06 × 10 1
Scenario 2 9.96 1.03 × 10 2 1.13 × 10 2 2.07 × 10 1
Scenario 3 14.55 1.17 × 10 2 1.46 × 10 2 2.48 × 10 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rosić, M.; Sedak, M.; Simić, M.; Pejović, P. Chaos-Enhanced Adaptive Hybrid Butterfly Particle Swarm Optimization Algorithm for Passive Target Localization. Sensors 2022, 22, 5739. https://doi.org/10.3390/s22155739

AMA Style

Rosić M, Sedak M, Simić M, Pejović P. Chaos-Enhanced Adaptive Hybrid Butterfly Particle Swarm Optimization Algorithm for Passive Target Localization. Sensors. 2022; 22(15):5739. https://doi.org/10.3390/s22155739

Chicago/Turabian Style

Rosić, Maja, Miloš Sedak, Mirjana Simić, and Predrag Pejović. 2022. "Chaos-Enhanced Adaptive Hybrid Butterfly Particle Swarm Optimization Algorithm for Passive Target Localization" Sensors 22, no. 15: 5739. https://doi.org/10.3390/s22155739

APA Style

Rosić, M., Sedak, M., Simić, M., & Pejović, P. (2022). Chaos-Enhanced Adaptive Hybrid Butterfly Particle Swarm Optimization Algorithm for Passive Target Localization. Sensors, 22(15), 5739. https://doi.org/10.3390/s22155739

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop