# Passive Target Localization Problem Based on Improved Hybrid Adaptive Differential Evolution and Nelder-Mead Algorithm

Authors: Maja B. Rosić; Mirjana I. Simić; Predrag V. Pejović
Journal: Journal of Sensors
 (2020)

---

## Abstract

This paper considers a passive target localization problem in Wireless Sensor Networks (WSNs) using the noisy time of arrival (TOA) measurements, obtained from multiple receivers and a single transmitter. The objective function is formulated as a maximum likelihood (ML) estimation problem under the Gaussian noise assumption. Consequently, the objective function of the ML estimator is a highly nonlinear and nonconvex function, where conventional optimization methods are not suitable for this type of problem. Hence, an improved algorithm based on the hybridization of an adaptive differential evolution (ADE) and Nelder-Mead (NM) algorithms, named HADENM, is proposed to find the estimated position of a passive target. In this paper, the control parameters of the ADE algorithm are adaptively updated during the evolution process. In addition, an adaptive adjustment parameter is designed to provide a balance between the global exploration and the local exploitation abilities. Furthermore, the exploitation is strengthened using the NM method by improving the accuracy of the best solution obtained from the ADE algorithm. Statistical analysis has been conducted, to evaluate the benefits of the proposed modifications on the optimization performance of the HADENM algorithm. The comparison results between HADENM algorithm and its versions indicate that the modifications proposed in this paper can improve the overall optimization performance. Furthermore, the simulation shows that the proposed HADENM algorithm can attain the Cramer-Rao lower bound (CRLB) and outperforms the constrained weighted least squares (CWLS) and differential evolution (DE) algorithms. The obtained results demonstrate the high accuracy and robustness of the proposed algorithm for solving the passive target localization problem for a wide range of measurement noise levels.


---

## Body

## 1. Introduction

Determining the position of a passive target using time of arrival (TOA) measurements has become an important issue for a number of different applications, such as radar, sonar, telecommunications, mobile communications, and navigation [1, 2]. In general, localization systems can be classified into active and passive [1]. With the active localization system, the target actively participates in the localization process. Compared to the active case, in passive localization, the target is not involved in the localization process and can only reflect or scatter the signals from the transmitter to receivers [2]. Therefore, the passive localization has been widely applied to different fields, such as robots, underwater acoustics, radar, crime-prevention, surveillance, and urban environments [3, 4].Hence, this paper proposes a passive target localization problem using noisy TOA measurements. In this way, the signal propagation time from the transmitter via the target to the receiver can be measured and used to determine range measurements, e.g., transmitter-target-receiver distances. Due to the nonlinear relationship between the target position and measurements, various efficient estimation methods have been proposed, such as nonlinear least squares (NLS) and maximum likelihood (ML) [5]. To obtain a closed-form solution, the NLS problem has been linearized using the linear least squares (LLS) and the weighted least squares (WLS) algorithms [6, 7]. In order to improve the accuracy, especially in the case of high measurement noise, the constrained weighted least squares (CWLS) is introduced [8]. Hence, the CWLS problem can be formulated as a quadratically constrained quadratic program (QCQP), which is converted to an unconstrained optimization problem using the method of Lagrange multipliers [8]. Another widely used estimation method is the ML estimator, which is commonly applied when the measurement error distribution is previously known [2]. In general, the ML objective function is highly nonlinear, with multiple local optima, and thus, the solution in closed-form cannot be obtained. Therefore, finding the global optimal solution by conventional optimization methods is difficult, due to the multimodal objective function [9]. For this reason, various efficient optimization algorithms have been derived to overcome these types of difficulties and to answer the challenges of complex optimization problems [10].Motivated by these facts, this paper proposes evolutionary algorithms (EAs) and their hybrid variants to overcome these drawbacks and improve the localization performance [11]. Generally, the optimization process of EAs consists of two phases, namely, exploration and exploitation [10]. In such a context, the first phase consists of exploring the search space and locating the region of the global optimum, while the exploitation phase investigates the promising area to refine the solution around the current global best solution. Therefore, finding an appropriate balance between the exploration and exploitation abilities during the evolution process is a significant challenge for optimization algorithms.Various EAs have been successfully applied to solve different complex optimization problems, such as the particle swarm optimization (PSO) [9], the differential evolution (DE) [12], the artificial bee colony (ABC) [13], and the cuckoo search algorithm (CS) [14]. Among these algorithms, the DE has emerged as a very efficient and robust EA for global optimization with successful applications in finding a global optimum [15].The DE is a population-based EA, which has been widely used to solve numerous optimization problems in different fields of science and engineering [15, 16]. Easy implementation, compact structure, reliability, and robustness are the main advantages of the DE algorithm [15]. However, some difficulties occur, such as a weak local search and slow convergence [17]. In general, the main factors that affect the performance are mutation and crossover operators and their corresponding control parameters, such as the scale factor (F), crossover rate (CR), and population size (NP). In this direction, various empirical guidelines have been provided for selecting the most suitable values for control parameters such as F∈0.4,1 and CR∈0,1 [18]. To obtain acceptable results for a given problem, different control parameter values at different stages of the evolution process are needed. In this way, a large F is required at the early stage of the evolution process for the purpose of preventing premature convergence, while a smaller F is preferred to accelerate the convergence at the later stage of the evolution process [19]. A large value of CR in the early stage of evolution can increase the population diversity. On the contrary, a smaller CR can improve the local exploitation ability and convergence at the later stage of evolution [19]. In this regard, different adaptive mechanisms eliminate the need for manual tuning of control parameters, which are adjusted adaptively using the feedback from the evolution process [20, 21].Furthermore, choosing the appropriate mutation operator can significantly affect the optimization performance of the DE algorithm. In this way, a number of different mutation operators have been developed, such as DE/rand/1, DE/rand/2, DE/best/1, DE/best/2, DE/current-to-best/1, DE/current-to-rand/1, and different variants of them [15]. The DE/rand/1 and DE/rand/2 generally have powerful global exploration ability, while the DE/best/1 and DE/current-to-best/1 have strong local exploitation ability [20].Hybridization of DE with other algorithms is another way to overcome the drawbacks of both algorithms and further enhance the optimization performance. Depending on the type of algorithm, the DE can be hybridized with other EAs, such as ABC, CS, and PSO [13, 22] or with different local search methods such as Powell’s method, the Hook-Jeeves (HJ), and the Nelder-Mead (NM) [23–25]. Among them, the NM method has been chosen, due to its excellent local search ability. However, the convergence of the NM method is extremely sensitive to the choice of the initial point [25], and thus, it cannot be used alone to find the global optimum of a multimodal objective function.Based on the above considerations, this paper is aimed at proposing an improved HADENM algorithm, in order to efficiently find the estimated position of a passive target. The proposed algorithm is a two-stage method, where in the first stage, the adaptive differential evolution (ADE) algorithm is used as the global optimizer, to quickly locate the promising region containing the global optimum. Then, in the second stage, the NM method is employed to improve the accuracy of the best solution obtained from the ADE algorithm. Moreover, an adaptive adjustment parameter is introduced in the mutation phase of the ADE algorithm to automatically apply the appropriate mutation operator, in order to avoid the problem of stagnation and premature convergence. In addition, to further improve the optimization performance, the control parameters of the ADE algorithm are automatically and adaptively adjusted during the evolution process. Therefore, the purpose of this paper is to propose a robust optimization algorithm for the passive target localization problem for a wide range of measurement noise levels.Finally, the Cramer-Rao lower bound (CRLB), which provides a lower bound on the variance of any unbiased estimator, has been derived and used as a benchmark to evaluate the localization performance [26]. Hence, the CRLB for the passive target localization problem is compared with the root mean square error (RMSE) performance of the proposed HADENM algorithm and the existing CWLS and DE algorithms.The contributions of this paper are summarized as follows:(i)
The passive target localization problem in LOS (line-of-sight) environment is formulated using the TOA measurements obtained from a set of receivers and a single transmitter, in Wireless Sensor Networks (WSNs)
(ii)
An improved HADENM algorithm, as the hybridization of the ADE and NM algorithms, is proposed in order to efficiently solve the passive target localization problem. Furthermore, the scale factor and the crossover rate are updated using an adaptive strategy. In addition, an adaptive mutation operator has been designed in order to maintain a balance between the exploration and exploitation. To further increase the exploitation ability, the NM method is employed with the aim to enhance the accuracy of the best solution previously obtained by the ADE algorithm
(iii)
The experiments are carried out to evaluate the benefits of the proposed modifications on the optimization performance of the proposed HADENM algorithm. Based on the statistical analysis of the proposed algorithm and its versions, it can be concluded that the modifications proposed in this paper can improve the overall optimization performance. Furthermore, the simulation results show the effectiveness of the HADENM algorithm for a wide range of measurement noise levels compared to the CWLS and DE algorithms. Additionally, the proposed algorithm attains the CRLB accuracy and shows better robustness against variations in network topologyThe structure of the rest of the paper is organized as follows. Section2 presents a review of background and related work. In Section 3, the passive target localization problem using noisy TOA measurements is presented. Section 4 gives the derivation of the ML objective function for the passive target localization problem. In Section 5, the CWLS algorithm is described. Section 6 introduces the DE algorithm followed by corresponding modifications. A local search NM method is described in Section 7. In Section 8, the HADENM algorithm is introduced to solve the considered passive target localization problem. In Section 9, the derivation of the CRLB is developed. To evaluate the localization accuracy, experimental results are presented in Section 10. The conclusion and future work are given in Section 11.


## 2. Background and Related Work

Highly accurate passive target localization can be considered as one of the most significant and challenging issues in WSNs [26]. The accuracy of the estimated target position depends on the geometric configuration of sensors and the measurement accuracy [1]. The localization algorithms in WSNs can be divided into range-based and range-free [5]. Range-based algorithms use range measurements including distance or angle between the target and the receivers, such as TOA [26], time difference of arrival (TDOA) [8], received signal strength (RSS) [27], angle of arrival (AOA) [28], or a combination of them [29]. Among them, TOA is the most commonly used algorithm for solving the localization problems, which can achieve high localization accuracy [26]. On the contrary, the range-free algorithms do not measure the distance or angle information; these algorithms use connectivity information in WSNs to estimate sensor positions [30]. Compared with range-based algorithms, the range-free algorithms do not require complex hardware structure. Hence, the range-free algorithms are cost effective and easy to implement; however, they are less accurate in estimating the sensor position [30].There are many applications, in which global positioning system (GPS) or other navigation systems cannot be used, such as in indoor, underwater acoustics, and urban environments [31]. In these cases, when GPS signals are not available or do not have sufficient accuracy, the passive localization system can be employed as an efficient alternative. In this regard, the passive target localization system has become an attractive solution for determining the unknown position of the passive target under various circumstances and environments [32]. In this way, a novel two-step algorithm for TOA passive target localization has been proposed for synchronous networks [26]. Furthermore, the two new algorithms based on belief propagation on a factor graph have been developed to localize multiple passive targets, in the case where the receiver position errors exist [2]. In addition, the two-step linear algorithm has been proposed to simultaneously estimate the position of the passive target and the unknown time offset in a quasisynchronous network using TOA measurements [33].Numerous estimation methods have been proposed to find the position of a target, where the well-known NLS estimation method is commonly employed. To obtain a closed-form solution, the LLS and WLS [6, 7] methods have been widely applied and implemented. However, these linear estimation methods cannot provide high localization accuracy of the target for different measurement noise levels [34]. To further improve the localization performance, especially in the presence of high noise levels, the CWLS algorithm is developed [8]. An additional widely used estimation method in the literature is the ML estimator, which maximizes the likelihood function of the unknown target position. However, the corresponding ML objective function is a highly nonlinear function, under the Gaussian noise assumption [9]. In order to avoid difficulties related to the multimodal nature of the ML objective function, a semidefinite programing (SDP) relaxation is applied to transform the nonconvex problem into a convex problem [35]. However, in the presence of significant measurement noise, the SDP has some disadvantages in terms of accuracy; thus, it can only provide a near-global optimal solution [35]. Thus, to solve the ML estimation problem with high accuracy and improve the convergence, for a given localization problem, researchers have developed and applied various efficient EAs [10].A number of different EAs like DE, PSO, ABC, and CS have been widely applied in order to solve the localization problem, in terms of determining the unknown position of the target [36]. In this regard, the PSO algorithm and its improved variants have been applied to accurately estimate the position of the passive target using TDOA measurements [37]. The simulation results have shown that improved variants of PSO provide better localization performance compared to the conventional PSO and well-known LLS and WLS algorithms. In addition, the cuckoo search (CS) algorithm has been used to estimate the position of the target in the passive localization system based on TDOA measurements [14]. Based on the simulation results, it can be concluded that the CS algorithm has faster convergence speed and more efficient global exploration ability compared to the PSO and Newton iteration algorithms.In recent years, a number of hybrid variants of the DE algorithm have been created in order to efficiently balance the exploration and exploitation abilities during the evolution process [23–25]. With respect to this, the direct search Powell’s method has been combined with DE (DESAP) [23], where the near-global solution obtained by the DE algorithm is improved using Powell’s method. To enhance the performance, the hybridization between DE and a local search algorithm (MLSHADE-SPA), with linear population size reduction and semiparameter adaptation, has been proposed [38]. Furthermore, the DE algorithm has been hybridized with Hook-Jeeves in distributed memetic DE algorithm (DMDE), to efficiently find the global optimum and achieve a better trade-off between the exploration and the exploitation [24]. In addition, a new reflection-based mutation operation, inspired by the reflection operations in the NM method, has been incorporated into the DE algorithm, with a combination of multiple mutation strategies based on roulette wheel selection (MM-RDE) [25].Therefore, this paper proposes an improved HADENM algorithm, based on the hybridization of ADE and NM algorithms, to solve the multimodal passive target localization problem with high accuracy even in the presence of high measurement noise.


## 3. Localization Problem

This section considers passive target localization problem in two dimensions, using noisy TOA measurements. The unknown position of a passive target can be determined using one transmitterTx and a set of receivers Rii=1N, where N is the number of receivers, such that N≥3. Let x=xyT∈R2 be the unknown position of a passive target, xir=xiryirT∈R2, ∀i∈1,2,⋯,N be the known coordinates of the ith receiver, and xT=00T∈R2 be the coordinates of the transmitter, as shown in Figure 1.Figure 1
Passive target localization system using noisy TOA measurements.The passive target localization system starts with the transmitter emitting a signal, and upon arrival, the signal is reflected by the target. The considered TOA algorithm employs the information of the absolute signal travel time from the transmitter to the receivers to obtain range measurements. Therefore, it is necessary to know the signal departure time, which is achieved by the synchronization between the transmitter and receivers [26]. In addition, the Gaussian noise is widely used in the localization algorithms under LOS environment [2]. In this case, the noisy TOA measurements are obtained as follows:
(1)ti=1cx2+x−xir2+n¯i,∀i∈1,2,⋯,N,where ⋅2 denotes the Euclidean distance, c is the speed of light, and n¯i denotes the zero-mean Gaussian measurement noise with variance σi2. Then, multiplying ti by the speed of light c, the range measurements (transmitter-target-receiver distances), denoted by rii=1N, can be written as
(2)ri=c⋅ti=di+ni,∀i∈1,2,⋯,N,where di=x2+x−xir2 is the true range measurement and ni=cn¯i is the measurement noise, which follows the Gaussian distribution N0,σni2, with zero-mean and variance σni2=c2σi2.From geometric interpretation in Figure1, in the absence of measurement noise, the true range measurement di defines an ellipse Ei with focal points placed at xT and xir, respectively. Therefore, the corresponding ellipse Ei can be defined as
(3)Ei=x,y∈R2:di=x2+y2+x−xir2+y−yir2.In this regard, the true passive target positionx=xyT∈R2 is determined at the intersection point of at least three solid line ellipses Eii=1N≥3, corresponding to the TOA measurements, as shown in Figure 1.Then, in the presence of measurement noise, the vector form of Equation (2) can be written as
(4)r=dx+n,in which the vector of true range measurements is
(5)dx=x2+x−x1r2⋮x2+x−xNr2,n=n1⋯nNT is the corresponding measurement noise vector. In this case, three or more ellipses derived from measurements do not have a unique intersection point. Thus, an estimated position of the passive target can be found inside the bounded region, i.e., the region surrounded by the bold black curve, as depicted in Figure 1. In this way, the main goal is to estimate the unknown position of the passive target based on noisy TOA measurements, which involves solving a highly nonlinear and multimodal ML estimation problem explained in the following section.


## 4. Maximum Likelihood Method

The unknown position of the passive target can be estimated through maximizing the likelihood function. Under the assumptions of Gaussian noisy TOA measurements, which are independent and identically distributed, the likelihood functionLx for the passive target position [26] can be expressed as
(6)Lx=frx=12πN/2detC1/2exp−12r−dxTC−1r−dx,where frx denotes the probability density function and C=diagσ12,⋯,σN2 is a diagonal covariance matrix.In order to simplify the maximization problem, it is often convenient to take the logarithm of the likelihood function, which can be written as(7)lnLx=c−12r−dxTC−1r−dx,where c=ln1/2πN/2detC1/2 is a constant independent of x. Consequently, the ML estimator requires the maximization of the log-likelihood function, which is equivalent to minimizing the negative logarithm of the likelihood function. Then, after neglecting the constant terms, the ML estimation problem can be formulated as
(8)minx∈R2JMLx=minx∈R2r−dxTC−1r−dx,where JMLxdenotes the corresponding ML objective function.Thus, the goal is to find the optimal solutionx^ that minimizes the objective function JMLx with respect to x, such that
(9)x^=argminx∈R2JMLx.Hence, theJMLx for the passive target localization problem is depicted in Figure 2.Figure 2
Illustration of the objective functionJMLx.Figure2 provides information about the nature of the optimization problem and shows that JMLx is a highly nonlinear and nonconvex function with multiple local optima. Therefore, finding the global optimal solution for a given optimization problem becomes a significant challenge. Thus, in order to solve this kind of complex optimization problem, it is necessary to employ different optimization algorithms, which are described in the next sections.


## 5. Constrained Weighted Least Squares Algorithm

As a well-known and widely used localization algorithm in WSNs, the CWLS algorithm is considered in this paper in order to compare and evaluate the localization performance of the proposed HADENM algorithm. The passive target position is estimated using CWLS algorithm by squaring both sides of Equation (2) and introducing Rt=x2=x2+y2 [8], which yields
(10)2xirx+2yiry−2riRt=xir2+yir2−ri2−2niRt+2rini,∀i∈1,2,⋯,N,where the second-order term of noise ni2 can be neglected for small noise levels. Hence, Equation (10) is linear, which can be written as
(11)Aθ=b+m,where
(12)A=2x1ry1r−r1⋮⋮⋮xNryNr−rN,(13)θ=xyRtT,(14)b=x1r2+y1r2−r12⋮xNr2+yNr2−rN2,(15)m=2r1−Rtn1⋯rN−RtnNT.Then, based on Equation (11) through Equation (15), the following WLS optimization problem can be formulated as
(16)minθJθ=minθAθ−bTWAθ−b,where, under the small noise assumption, W∈RN×N is a symmetric weighting matrix, which can be defined as
(17)W=EmmT−1=DTCD−1,in which D=diag2r1−Rt,⋯,2rN−Rt.It should be noted that since the measurement noisem from Equation (15) is Gaussian distributed and due to the linear relationship in Equation (11), the objective function of the ML estimator JMLx, given in Equation (8), is equivalent to that of the WLS estimator [8]. Then, by incorporating the relationship between unknown target position x and an auxiliary variable Rt as a second-order equality constraint in the WLS estimator, the following CWLS target localization problem is obtained as
(18)minθJθ=minθAθ−bTWAθ−b,s.t.θTSθ=0,where S=diag1,1,−1 is a diagonal matrix.The closed-form solution is not available due to the nonlinearity of the constrained problem. In this regard, an unconstrained optimization problem is obtained by forming the auxiliary Lagrange functionℒθ,λ, as follows:
(19)minθ,λℒθ,λ=minθ,λAθ−bTWAθ−b+λθTSθ,where λ is the Lagrange multiplier associated with the equality constraint. The corresponding optimal solution θ^CWLS to the CWLS problem can be obtained as
(20)θ^CWLS=ATWA+λS−1ATWb.Hence, the Lagrange multiplierλ is found as the root of a 4th-order polynomial [8], as follows:
(21)∑i=13αiβiζi+λ2=0,where α=ΣTSATWb=α1α2α3T, β=Σ−1ATWb, and ATWAS=Σdiagζ1,ζ2,ζ3Σ−1 [8]. Then, based on the CWLS method, the algorithm for the passive target localization can be stated as follows:Step 1.
Set the symmetric weighting matrixW as unit matrix I, i.e., W=I.Step 2.
Find all roots of Equation (21) taking into consideration only real roots.Step 3.
For eachλ, obtained in Step 2, determine subestimates of θ^CWLS using Equation (20). Then, separate subestimates for which Jθ is minimal.Step 4.
Determine a weighting matrixW using Equation (17).Step 5.
Repeat Steps 2–4, while the predefined termination criterion is satisfied.


## 6. Differential Evolution Algorithm and the Proposed Modified Version

### 6.1. Differential Evolution Algorithm

The DE is a population-based EA successfully applied for global optimization, which is introduced by Storn and Price [12]. The evolution process of the DE algorithm is based on four basic steps, i.e., initialization, mutation, crossover, and selection, which are described below.


#### 6.1.1. Initialization

The evolution process of the DE algorithm starts with an initial population ofn-dimensional NP individuals, in which each individual represents a candidate solution of the problem, i.e., xiG:xi,jG∈xi,jL,xi,jU, ∀i∈1,2,⋯,NP, ∀j∈1,2,⋯,n, where G denotes the current generation. For the considered localization problem, every individual xiG in two-dimensional space (i.e., n=2) has two variables corresponding to x and y coordinates of the target position. Here, xi,jL and xi,jU are the lower and upper bounds of the jth component of the ith individual xi,j, respectively. At G=0, every individual is randomly generated as
(22)xi,j0=xi,jL+randjxi,jU−xi,jL,where randj is a uniformly distributed random number in the range 0,1.


#### 6.1.2. Mutation

The mutation operator is employed, after initialization, to generate a new mutant vectorviG=vi,1G⋯vi,jG⋯vi,nGT for each target vector xiG. The most widely used mutation operators in the literature [15] areDE/rand/1:(23)viG=xr1G+Fxr2G−xr3G,DE/rand/2:(24)viG=xr1G+Fxr2G−xr3G+Fxr4G−xr5G,DE/best/1:(25)viG=xbestG+Fxr1G−xr2G,DE/current-to-best/1(26)viG=xiG+FxbestG−xiG+Fxr1G−xr2G,where r1, r2, r3, r4, and r5 are distinct integers randomly selected from 1,2,⋯,NP\i, xbestG is the individual with the lowest objective function value, and F∈0,1 is the scale factor.


#### 6.1.3. Crossover

In order to increase the diversity of the population, a binomial crossover operator is applied on the target vectorxiG and the corresponding mutant vector viG to produce a trial vector uiG, which can be expressed as follows:
(27)ui,jG=vi,jGifrandi,j≤CR∨j=jrandxi,jGotherwise,∀j∈1,2,⋯,n,∀i∈1,2,⋯,NP,where randi,j is an uniform number generated within 0,1, CR∈0,1 is the crossover rate, and jrand∈1,2,⋯,n is an integer selected randomly.


#### 6.1.4. Selection

After crossover, the selection operator is employed to compare the objective function values of the corresponding trial vectoruiG with the objective function value of the target vector xiG. The objective function value of the ith individual in the Gth generation has been calculated using the ML objective function in (8), as follows:
(28)fxiG=JMLxiG,(29)fuiG=JMLuiG.In this way, the vector with a better objective function value is selected, according to(30)xiG+1=uiGiffuiG≤fxiG,xiGotherwise.


### 6.2. A Modified Differential Evolution Algorithm

In order to improve the performance of the conventional DE algorithm, this subsection introduces modifications to the scale factor and crossover rate. Furthermore, an automatically adapted mutation operator is proposed, to select an appropriate mutation strategy based on the current optimization state.


#### 6.2.1. Adaptive Scale Factor

The performance of the DE algorithm depends largely on the proper choice of the scale factor, which can affect the convergence. In this context, a higher value of the scale factor, in the early stage, improves the exploration ability, which has a positive effect on the population diversity and, thus, avoids premature convergence at local optima [18, 39]. In contrast, a smaller value of the scale factor, in the later stage, improves the exploitation ability, which can enhance the convergence speed [18, 39].Based on the above analysis, an adaptive scale factorFiG, which is dynamically adapted in each generation for each individual, is introduced as follows:
(31)FiG=Fmax−Fminexp−10GfiG−fbestG/fworstG−fbestG,∀i∈1,2,⋯,NP,where Fmin and Fmax are minimum and maximum values of FiG, and fbestG, fworstG, and fiG are the objective function values of the best, worst, and ith individual, respectively.In this respect, Figure3 illustrates the changes of the proposed adaptive scale factor FiG defined in Equation (31) during the evolution process, where Fmax=0.9 and Fmin=0.5.Figure 3
The changes of the proposed adaptive scale factorFiG during the evolution process.According to Figure3, it is evident that in the early stage, the adaptive scale factor FiG has a larger value, which is suitable for exploration. During the later stage, with the increase of generations, the adaptive scale factor FiG starts to decrease, which improves exploitation and accelerates the local convergence of the algorithm.


#### 6.2.2. Adaptive Crossover Rate

According to the analysis in [18, 40], adapting the suitable value of the crossover rate can maintain the diversity of the population and improve the quality of the solution. From Equation (30), it is evident that for the large value of the crossover rate, the mutant vector viG has a greater contribution to the trial vector uiG. In this case, a higher crossover rate increases the population diversity and enhances the global search. In contrast, for a smaller crossover rate, any trial vector uiG keeps the previous state xiG with a large probability. This can further refine the trial vector, which is beneficial for improving the quality of the solution.In this regard, an adaptive crossover rateCRG is proposed here, to address the above mentioned issues, which is defined as
(32)CRG=CRmax−CRmin⋅2−exp3−Gmax/G+1+CRmin,where Gmax is the maximum number of generations and CRmax and CRmin are the maximum and minimum values of CRG, respectively.In this way, Figure4 shows the proposed adaptive crossover rate CRG versus the number of generations, where CRmax=0.9 and CRmin=0.1.Figure 4
The changes of the proposed adaptive crossover rateCRG during the evolution process.As can be seen from Figure4, the CRG has a large value in the initial stage of the evolution process and gradually decreases with the increase of generations. In this way, larger CRG can advance the population diversity and strongly enhance the exploration ability. During the later stage of the evolution process, the CRG takes a small value, which is beneficial for improving the quality of the solution.


#### 6.2.3. The Adaptive Mutation Operator

The performance of the DE algorithm, such as convergence, population diversity, and exploration ability, is greatly affected by mutation operators. The appropriate mutation operators for different evolutionary stages have been proposed here, with the aim to avoid premature convergence and prevent stagnation. In this manner, the DE/rand/1 and DE/rand/2 have strong global exploration ability, while the DE/best/1 operator has a good local exploitation ability [20]. In order to further improve the local search ability, the DE/current-to-best/1 mutation operator can be employed [20].Based on the above considerations, to dynamically adjust the global exploration and local exploitation abilities, this paper proposes an adaptive adjustment parameterδG, as follows:
(33)δG=fbestG−fmeanGfbestG−fworstG,where fmeanG is mean value of the objective function. Thus, the pseudocode of the proposed adaptive mutation operator is shown in Algorithm 1.Algorithm 1: Pseudocode of the adaptive mutation operator.
if δG>0.5 then
if rand>0.5 then
viG=xr1+FiGxr2−xr3
else
viG=xr1+FiGxr2−xr3+FiGxr4−xr5
end if
else if δG>0.5 then
if rand>0.5 then
viG=xbest+FiGxr1−xr2
else
viG=xi+FiGxbest−xi+FiGxr1−xr2
end if
end ifAs a matter of fact, the adaptive parameterδG∈0,1 has an important influence on the identification of evolutionary stages in the search process. According to Equation (33), in the early stage, the value of δG is close to 1, which indicates that the population is far from the region of the global optimum, and this corresponds to the global exploration in the search space. Hence, the DE/rand/1 and DE/rand/2 will be chosen randomly with the probability of 0.5, in order to improve the exploration and find the region of the global optimum. In the later stage, the value of δG is close to 0, which shows that the population is near the region of the global optimum, and this corresponds to the local exploitation. Therefore, the DE/best/1 and DE/current-to-best/1 will be randomly selected, with the probability of 0.5, to strengthen the exploitation ability and further improve the solution quality.


## 6.1. Differential Evolution Algorithm

The DE is a population-based EA successfully applied for global optimization, which is introduced by Storn and Price [12]. The evolution process of the DE algorithm is based on four basic steps, i.e., initialization, mutation, crossover, and selection, which are described below.


### 6.1.1. Initialization

The evolution process of the DE algorithm starts with an initial population ofn-dimensional NP individuals, in which each individual represents a candidate solution of the problem, i.e., xiG:xi,jG∈xi,jL,xi,jU, ∀i∈1,2,⋯,NP, ∀j∈1,2,⋯,n, where G denotes the current generation. For the considered localization problem, every individual xiG in two-dimensional space (i.e., n=2) has two variables corresponding to x and y coordinates of the target position. Here, xi,jL and xi,jU are the lower and upper bounds of the jth component of the ith individual xi,j, respectively. At G=0, every individual is randomly generated as
(22)xi,j0=xi,jL+randjxi,jU−xi,jL,where randj is a uniformly distributed random number in the range 0,1.


### 6.1.2. Mutation

The mutation operator is employed, after initialization, to generate a new mutant vectorviG=vi,1G⋯vi,jG⋯vi,nGT for each target vector xiG. The most widely used mutation operators in the literature [15] areDE/rand/1:(23)viG=xr1G+Fxr2G−xr3G,DE/rand/2:(24)viG=xr1G+Fxr2G−xr3G+Fxr4G−xr5G,DE/best/1:(25)viG=xbestG+Fxr1G−xr2G,DE/current-to-best/1(26)viG=xiG+FxbestG−xiG+Fxr1G−xr2G,where r1, r2, r3, r4, and r5 are distinct integers randomly selected from 1,2,⋯,NP\i, xbestG is the individual with the lowest objective function value, and F∈0,1 is the scale factor.


### 6.1.3. Crossover

In order to increase the diversity of the population, a binomial crossover operator is applied on the target vectorxiG and the corresponding mutant vector viG to produce a trial vector uiG, which can be expressed as follows:
(27)ui,jG=vi,jGifrandi,j≤CR∨j=jrandxi,jGotherwise,∀j∈1,2,⋯,n,∀i∈1,2,⋯,NP,where randi,j is an uniform number generated within 0,1, CR∈0,1 is the crossover rate, and jrand∈1,2,⋯,n is an integer selected randomly.


### 6.1.4. Selection

After crossover, the selection operator is employed to compare the objective function values of the corresponding trial vectoruiG with the objective function value of the target vector xiG. The objective function value of the ith individual in the Gth generation has been calculated using the ML objective function in (8), as follows:
(28)fxiG=JMLxiG,(29)fuiG=JMLuiG.In this way, the vector with a better objective function value is selected, according to(30)xiG+1=uiGiffuiG≤fxiG,xiGotherwise.


## 6.1.1. Initialization

The evolution process of the DE algorithm starts with an initial population ofn-dimensional NP individuals, in which each individual represents a candidate solution of the problem, i.e., xiG:xi,jG∈xi,jL,xi,jU, ∀i∈1,2,⋯,NP, ∀j∈1,2,⋯,n, where G denotes the current generation. For the considered localization problem, every individual xiG in two-dimensional space (i.e., n=2) has two variables corresponding to x and y coordinates of the target position. Here, xi,jL and xi,jU are the lower and upper bounds of the jth component of the ith individual xi,j, respectively. At G=0, every individual is randomly generated as
(22)xi,j0=xi,jL+randjxi,jU−xi,jL,where randj is a uniformly distributed random number in the range 0,1.


## 6.1.2. Mutation

The mutation operator is employed, after initialization, to generate a new mutant vectorviG=vi,1G⋯vi,jG⋯vi,nGT for each target vector xiG. The most widely used mutation operators in the literature [15] areDE/rand/1:(23)viG=xr1G+Fxr2G−xr3G,DE/rand/2:(24)viG=xr1G+Fxr2G−xr3G+Fxr4G−xr5G,DE/best/1:(25)viG=xbestG+Fxr1G−xr2G,DE/current-to-best/1(26)viG=xiG+FxbestG−xiG+Fxr1G−xr2G,where r1, r2, r3, r4, and r5 are distinct integers randomly selected from 1,2,⋯,NP\i, xbestG is the individual with the lowest objective function value, and F∈0,1 is the scale factor.


## 6.1.3. Crossover

In order to increase the diversity of the population, a binomial crossover operator is applied on the target vectorxiG and the corresponding mutant vector viG to produce a trial vector uiG, which can be expressed as follows:
(27)ui,jG=vi,jGifrandi,j≤CR∨j=jrandxi,jGotherwise,∀j∈1,2,⋯,n,∀i∈1,2,⋯,NP,where randi,j is an uniform number generated within 0,1, CR∈0,1 is the crossover rate, and jrand∈1,2,⋯,n is an integer selected randomly.


## 6.1.4. Selection

After crossover, the selection operator is employed to compare the objective function values of the corresponding trial vectoruiG with the objective function value of the target vector xiG. The objective function value of the ith individual in the Gth generation has been calculated using the ML objective function in (8), as follows:
(28)fxiG=JMLxiG,(29)fuiG=JMLuiG.In this way, the vector with a better objective function value is selected, according to(30)xiG+1=uiGiffuiG≤fxiG,xiGotherwise.


## 6.2. A Modified Differential Evolution Algorithm

In order to improve the performance of the conventional DE algorithm, this subsection introduces modifications to the scale factor and crossover rate. Furthermore, an automatically adapted mutation operator is proposed, to select an appropriate mutation strategy based on the current optimization state.


### 6.2.1. Adaptive Scale Factor

The performance of the DE algorithm depends largely on the proper choice of the scale factor, which can affect the convergence. In this context, a higher value of the scale factor, in the early stage, improves the exploration ability, which has a positive effect on the population diversity and, thus, avoids premature convergence at local optima [18, 39]. In contrast, a smaller value of the scale factor, in the later stage, improves the exploitation ability, which can enhance the convergence speed [18, 39].Based on the above analysis, an adaptive scale factorFiG, which is dynamically adapted in each generation for each individual, is introduced as follows:
(31)FiG=Fmax−Fminexp−10GfiG−fbestG/fworstG−fbestG,∀i∈1,2,⋯,NP,where Fmin and Fmax are minimum and maximum values of FiG, and fbestG, fworstG, and fiG are the objective function values of the best, worst, and ith individual, respectively.In this respect, Figure3 illustrates the changes of the proposed adaptive scale factor FiG defined in Equation (31) during the evolution process, where Fmax=0.9 and Fmin=0.5.Figure 3
The changes of the proposed adaptive scale factorFiG during the evolution process.According to Figure3, it is evident that in the early stage, the adaptive scale factor FiG has a larger value, which is suitable for exploration. During the later stage, with the increase of generations, the adaptive scale factor FiG starts to decrease, which improves exploitation and accelerates the local convergence of the algorithm.


### 6.2.2. Adaptive Crossover Rate

According to the analysis in [18, 40], adapting the suitable value of the crossover rate can maintain the diversity of the population and improve the quality of the solution. From Equation (30), it is evident that for the large value of the crossover rate, the mutant vector viG has a greater contribution to the trial vector uiG. In this case, a higher crossover rate increases the population diversity and enhances the global search. In contrast, for a smaller crossover rate, any trial vector uiG keeps the previous state xiG with a large probability. This can further refine the trial vector, which is beneficial for improving the quality of the solution.In this regard, an adaptive crossover rateCRG is proposed here, to address the above mentioned issues, which is defined as
(32)CRG=CRmax−CRmin⋅2−exp3−Gmax/G+1+CRmin,where Gmax is the maximum number of generations and CRmax and CRmin are the maximum and minimum values of CRG, respectively.In this way, Figure4 shows the proposed adaptive crossover rate CRG versus the number of generations, where CRmax=0.9 and CRmin=0.1.Figure 4
The changes of the proposed adaptive crossover rateCRG during the evolution process.As can be seen from Figure4, the CRG has a large value in the initial stage of the evolution process and gradually decreases with the increase of generations. In this way, larger CRG can advance the population diversity and strongly enhance the exploration ability. During the later stage of the evolution process, the CRG takes a small value, which is beneficial for improving the quality of the solution.


### 6.2.3. The Adaptive Mutation Operator

The performance of the DE algorithm, such as convergence, population diversity, and exploration ability, is greatly affected by mutation operators. The appropriate mutation operators for different evolutionary stages have been proposed here, with the aim to avoid premature convergence and prevent stagnation. In this manner, the DE/rand/1 and DE/rand/2 have strong global exploration ability, while the DE/best/1 operator has a good local exploitation ability [20]. In order to further improve the local search ability, the DE/current-to-best/1 mutation operator can be employed [20].Based on the above considerations, to dynamically adjust the global exploration and local exploitation abilities, this paper proposes an adaptive adjustment parameterδG, as follows:
(33)δG=fbestG−fmeanGfbestG−fworstG,where fmeanG is mean value of the objective function. Thus, the pseudocode of the proposed adaptive mutation operator is shown in Algorithm 1.Algorithm 1: Pseudocode of the adaptive mutation operator.
if δG>0.5 then
if rand>0.5 then
viG=xr1+FiGxr2−xr3
else
viG=xr1+FiGxr2−xr3+FiGxr4−xr5
end if
else if δG>0.5 then
if rand>0.5 then
viG=xbest+FiGxr1−xr2
else
viG=xi+FiGxbest−xi+FiGxr1−xr2
end if
end ifAs a matter of fact, the adaptive parameterδG∈0,1 has an important influence on the identification of evolutionary stages in the search process. According to Equation (33), in the early stage, the value of δG is close to 1, which indicates that the population is far from the region of the global optimum, and this corresponds to the global exploration in the search space. Hence, the DE/rand/1 and DE/rand/2 will be chosen randomly with the probability of 0.5, in order to improve the exploration and find the region of the global optimum. In the later stage, the value of δG is close to 0, which shows that the population is near the region of the global optimum, and this corresponds to the local exploitation. Therefore, the DE/best/1 and DE/current-to-best/1 will be randomly selected, with the probability of 0.5, to strengthen the exploitation ability and further improve the solution quality.


## 6.2.1. Adaptive Scale Factor

The performance of the DE algorithm depends largely on the proper choice of the scale factor, which can affect the convergence. In this context, a higher value of the scale factor, in the early stage, improves the exploration ability, which has a positive effect on the population diversity and, thus, avoids premature convergence at local optima [18, 39]. In contrast, a smaller value of the scale factor, in the later stage, improves the exploitation ability, which can enhance the convergence speed [18, 39].Based on the above analysis, an adaptive scale factorFiG, which is dynamically adapted in each generation for each individual, is introduced as follows:
(31)FiG=Fmax−Fminexp−10GfiG−fbestG/fworstG−fbestG,∀i∈1,2,⋯,NP,where Fmin and Fmax are minimum and maximum values of FiG, and fbestG, fworstG, and fiG are the objective function values of the best, worst, and ith individual, respectively.In this respect, Figure3 illustrates the changes of the proposed adaptive scale factor FiG defined in Equation (31) during the evolution process, where Fmax=0.9 and Fmin=0.5.Figure 3
The changes of the proposed adaptive scale factorFiG during the evolution process.According to Figure3, it is evident that in the early stage, the adaptive scale factor FiG has a larger value, which is suitable for exploration. During the later stage, with the increase of generations, the adaptive scale factor FiG starts to decrease, which improves exploitation and accelerates the local convergence of the algorithm.


## 6.2.2. Adaptive Crossover Rate

According to the analysis in [18, 40], adapting the suitable value of the crossover rate can maintain the diversity of the population and improve the quality of the solution. From Equation (30), it is evident that for the large value of the crossover rate, the mutant vector viG has a greater contribution to the trial vector uiG. In this case, a higher crossover rate increases the population diversity and enhances the global search. In contrast, for a smaller crossover rate, any trial vector uiG keeps the previous state xiG with a large probability. This can further refine the trial vector, which is beneficial for improving the quality of the solution.In this regard, an adaptive crossover rateCRG is proposed here, to address the above mentioned issues, which is defined as
(32)CRG=CRmax−CRmin⋅2−exp3−Gmax/G+1+CRmin,where Gmax is the maximum number of generations and CRmax and CRmin are the maximum and minimum values of CRG, respectively.In this way, Figure4 shows the proposed adaptive crossover rate CRG versus the number of generations, where CRmax=0.9 and CRmin=0.1.Figure 4
The changes of the proposed adaptive crossover rateCRG during the evolution process.As can be seen from Figure4, the CRG has a large value in the initial stage of the evolution process and gradually decreases with the increase of generations. In this way, larger CRG can advance the population diversity and strongly enhance the exploration ability. During the later stage of the evolution process, the CRG takes a small value, which is beneficial for improving the quality of the solution.


## 6.2.3. The Adaptive Mutation Operator

The performance of the DE algorithm, such as convergence, population diversity, and exploration ability, is greatly affected by mutation operators. The appropriate mutation operators for different evolutionary stages have been proposed here, with the aim to avoid premature convergence and prevent stagnation. In this manner, the DE/rand/1 and DE/rand/2 have strong global exploration ability, while the DE/best/1 operator has a good local exploitation ability [20]. In order to further improve the local search ability, the DE/current-to-best/1 mutation operator can be employed [20].Based on the above considerations, to dynamically adjust the global exploration and local exploitation abilities, this paper proposes an adaptive adjustment parameterδG, as follows:
(33)δG=fbestG−fmeanGfbestG−fworstG,where fmeanG is mean value of the objective function. Thus, the pseudocode of the proposed adaptive mutation operator is shown in Algorithm 1.Algorithm 1: Pseudocode of the adaptive mutation operator.
if δG>0.5 then
if rand>0.5 then
viG=xr1+FiGxr2−xr3
else
viG=xr1+FiGxr2−xr3+FiGxr4−xr5
end if
else if δG>0.5 then
if rand>0.5 then
viG=xbest+FiGxr1−xr2
else
viG=xi+FiGxbest−xi+FiGxr1−xr2
end if
end ifAs a matter of fact, the adaptive parameterδG∈0,1 has an important influence on the identification of evolutionary stages in the search process. According to Equation (33), in the early stage, the value of δG is close to 1, which indicates that the population is far from the region of the global optimum, and this corresponds to the global exploration in the search space. Hence, the DE/rand/1 and DE/rand/2 will be chosen randomly with the probability of 0.5, in order to improve the exploration and find the region of the global optimum. In the later stage, the value of δG is close to 0, which shows that the population is near the region of the global optimum, and this corresponds to the local exploitation. Therefore, the DE/best/1 and DE/current-to-best/1 will be randomly selected, with the probability of 0.5, to strengthen the exploitation ability and further improve the solution quality.


## 7. Nelder-Mead Method

The Nelder-Mead is a local search method that does not require the derivative information of the objective function [25]. In general, the NM method is described for the minimization of an objective function in n-dimensional space. The evolution process of the NM method starts with an initial simplex in the search space, which is a polyhedron with n+1 vertices, i.e., xi:1≤i≤n. The objective function is evaluated at each vertex, and then all vertices are ranked based on their objective function values, i.e., the vertex corresponding to the best objective function value is denoted as x1 and the vertex corresponding to the worst objective function value is denoted as xn+1. The NM method uses four elementary geometric transformations, called reflection, expansion, contraction, and shrinkage. Using these geometric transformations, the simplex moves through the search space towards the optimal solution. The following steps are executed repeatedly until reaching stopping criteria:Step 1.
Initialization. Given an initial vertex x10, obtained by the ADE algorithm, generate the remaining n vertices xi+10:1≤i≤n according to
(34)xi+10=x10+λei,∀i∈1,2,⋯,n,where ei is the unit vector of the ith axis and λ is the initial step size, set as λ=1.Step 2.
Sorting. Rank the vertices as follows
(35)fx1k≤⋯≤fxn+1k,where k is the current iteration. Here, the objective function value is obtained from the ML objective function in Equation (8), as follows: fxik=JMLxik.Step 3.
Reflection. Generate the vertex xrk by reflecting the vertex xn+1k, as follows:
(36)xrk=x¯k+αx¯k−xn+1k,where α>0 is the reflection coefficient, usually suggested as α=1 and
(37)x¯k=1n∑i=1nxik,is the centroid of the n best vertices. If fx1k≤fxrk<fxnk, replace xn+1k+1 with xrk and terminate the iteration. If fxrk<fx1k, continue with the expansion; otherwise, perform contraction.Step 4.
Expansion. Calculate the vertex xek in the same direction as xrk, as follows:
(38)xek=x¯k+βxrk−x¯k,where β>1 is the expansion coefficient, usually suggested as β=2. If fxek<fxrk, replace the vertex xn+1k+1 with xek and terminate the iteration; otherwise, replace xn+1k+1 with xrk and terminate the iteration.Step 5.
Contraction. After reflection, there are two possible contraction cases:Step 5.1.
Outside contraction. If fxnk≤fxrk<fxn+1k, generate the vertex xck as
(39)xck=x¯k+γxrk−x¯k,where 0<γ<1 is the contraction coefficient, usually suggested as γ=0.5. If fxck≤fxrk, replace xn+1k+1 with xck and terminate the iteration; otherwise, perform shrinkage.Step 5.2.
Inside contraction. If fxrk≥fxn+1k, generate the vertex xcck as follows:
(40)xcck=x¯k−γxrk−x¯k.
Iffxcck<fxn+1k, replace xn+1k+1 with xcck and terminate the iteration; otherwise, do shrinkage.Step 6.
Shrinkage. Perform the shrinkage on the vertices ∀i∈2,⋯,n+1 as follows:
(41)xik+1=x1k+δxik−x1k,where 0<δ<1 is the shrinkage coefficient, usually suggested as δ=0.5. The optimization process described in Steps 1–6 is repeated until the termination criteria are reached. The obtained optimal solution represents the estimated position x^ of the unknown target position x.


## 8. HADENM Algorithm

In this section, the proposed HADENM algorithm is introduced that hybridizes the ADE algorithm with the local search NM method, in order to efficiently solve the passive target localization problem given in Equation (8). In this way, the pseudocode of the HADENM algorithm is presented in Algorithm 2, for the passive target localization problem.Algorithm 2: Pseudocode of the HADENM algorithm.
Initialize the parameters: NP, Gmax, Fmax, Fmin, CRmax, CRmin, α, β, γ
Generate a uniformly distributed random population
while ADE termination condition is not satisfied do
for i=1 to NP do
Compute the adaptive parameters FiG and CRG
Compute δG
if δG>0.5 then
if rand>0.5 then
viG⟵xr1+FiGxr2−xr3
else
viG⟵xr1+FiGxr2−xr3+FiGxr4−xr5
end if
else if δG≤0.5 then
if rand>0.5 then
viG⟵xbest+FiGxr1−xr2
else
viG⟵xi+FiGxbest−xi+FiGxr1−xr2
end if
end if
if randj,i⩽CRG∨j=jrand then
ui,jG⟵vi,jG
else
ui,jG⟵xi,jG
end if
if fui,jG⩽fxi,jG then
xiG+1⟵uiG
else
xiG+1⟵xiG
end if
end for
end while
Initialize the NM with the current best solution found by ADE
while NM termination condition is not satisfied do
Sort the vertices according to fx1k≤⋯≤fxn+1k
Perform reflection
if fx1k≤fxrk<fxnk then
xn+1k+1⟵xrk
end if
Perform expansion
if fxek<fxrk then
xn+1k+1⟵xek
else
xn+1k+1⟵xrk
end if
if fxnk≤fxrk<fxn+1k then
Perform outside contraction
if fxck≤fxrk then
xn+1k+1⟵xck
else
Perform Shrinkage
end if
else if fxrk≥fxn+1k then
Perform inside contraction
if fxcck<fxn+1k then
xn+1k+1⟵xcck
else
Perform Shrinkage
end if
end if
end while
Display the results


## 9. Cramer-Rao Lower Bound

In the passive target localization problem, the CRLB can be used as a benchmark for evaluating the performance of different algorithms. The CRLB is obtained from the diagonal elements of the inverse of the Fisher information matrix (FIM) [26], denoted by Ix, which is based on the probability density function frx, defined in Equation (6), and is represented as
(42)Ix=E∂lnfrx∂x∂lnfrx∂xT=−E∂2lnfrx∂x∂xT.Based on this, the FIM can be defined as(43)Ix=IxxIxyIyxIyy,where the corresponding elements are given as follows:
(44)Ixx=∑i=1N1σni2x−xirx−xir2+xx22,(45)Ixy=Iyx=∑i=1N1σni2x−xirx−xir2+xx2y−yirx−xir2+yx2,(46)Iyy=∑i=1N1σni2y−yirx−xir2+yx22.Hence, the relationship between the variance and CRLB can be determined as(47)Ex^−xx^−xT≥trIx−1=CRLBx,where x^ is the estimated value of x.


## 10. Experimental Study

In this section, experiments are conducted to evaluate the localization performance and to perform the analysis of the benefits of the proposed modifications on the optimization performance of the HADENM algorithm. In this regard, the presentation of the experimental results is divided into two subsections, described below.


### 10.1. A Parametric Study on HADENM Algorithm

In this subsection, the experiments are carried out to evaluate the performance of the HADENM algorithm by studying the effects of the proposed adaptive scale factor and crossover rate, as well as the proposed adaptive mutation operator on the optimization performance. Furthermore, the effects of the hybridization between the ADE and NM algorithms have been analysed. Then, the experiments are performed in order to verify that the performance is enhanced after combining the previously described improvements.To evaluate the performance of the proposed HADENM algorithm, the solution error measurefx^−fx∗ has been employed, where x^ denotes the best solution of the algorithm obtained in one run and x∗ represents the well-known global optimal solution of the corresponding ML objective function. All experiments for each algorithm have been independently run 30 times and statistical results are provided. From the statistical point of view, the quality of the obtained solutions has been analysed and compared using two nonparametric statistical hypothesis tests such as the Wilcoxon signed-rank test and Friedman test.Firstly, the Wilcoxon signed-rank test can be used to determine the significant differences between two samples obtained by the algorithms. This statistic test has been applied with a significance level ofα=0.05. Here, R+ denotes the sum of ranks for the problem in which the first algorithm outperformed the second and R− is the sum of ranks for the problem where the first algorithm performed worse than the second. According to [41], in the Wilcoxon signed-rank test for the null hypothesis, it is assumed that, “there is no difference between the mean results of the two samples.” On the contrary, the alternative hypothesis is that, “there is a difference in the mean results of the two samples.” In the statistical analysis, the p value is used and compared to the significance level. Thus, the null hypothesis can be rejected when the p value is less than or equal to α=0.05. Based on the obtained results of the statistical test, one of the following three signs (+, −, and ≈) has been assigned for the comparison between any two algorithms. Thus, the plus (+) sign denotes that the first algorithm is significantly better than the second, the minus (−) sign means that the first algorithm is significantly worse than the second, and the approximation (≈) sign denotes that there is no significant difference between two algorithms.The second test is the Friedman test, which obtains the ranks of all considered algorithms over every tested function, with the aim to find the significant differences in performance between two or more algorithms. In this statistical test, the algorithm with a minimum rank value is considered as the best performing algorithm, while the one with the highest rank value is considered as the worst. According to [41], the null hypothesis for Friedman test states that, “there is no difference among the performance of all algorithms,” whereas the alternative hypothesis states that, “there is a difference among the performance of all algorithms”.In this regard, the presentation of the obtained results of the considered evaluations is divided into four sub-subsections. In the first sub-subsection, the effectiveness of the proposed adaptive scale factor and crossover rate has been evaluated. In the second sub-subsection, the effectiveness of the proposed adaptive mutation operator has been considered. The third sub-subsection considers the hybridization of ADE and NM algorithms and the overall performance improvements. Finally, in the fourth sub-subsection, the statistical results of applying the Friedman test between the considered algorithms have been analysed.


#### 10.1.1. Effectiveness of the Adaptive Scale Factor and Crossover Rate

In this sub-subsection, the experimental studies have been performed to evaluate the effectiveness and benefits of the proposed adaptive scale factorFiG (given in Equation (31)) and adaptive crossover rate CRG (given in Equation (32)) on the optimization performance of the HADENM algorithm. Firstly, the effectiveness of the adaptive scale factor is considered, where the performance of the HADENM algorithm is evaluated with two different fixed values of F, as follows:
(1)
HADENMF=0.5, which has the same operators as HADENM, except that the scale factor is set to a fixed value of F=0.5
(2)
HADENMF=0.9, which has the same operators as HADENM, except that the scale factor is set to a fixed value of F=0.9The summary of statistical results of applying the Wilcoxon test between the proposed HADENM and the above two algorithms is presented in Table1.Table 1
Results of the Wilcoxon test for HADENM, HADENMF=0.5, and HADENMF=0.9 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENMF=0.5
346
119
0.0177
+
HADENM versus HADENMF=0.9
282
183
0.3207
≈From the results in Table1, it can be seen that the proposed HADENM algorithm is significantly better than HADENMF=0.5. On the other hand, in the case of HADENM versus HADENMF=0.9, the proposed algorithm has a higher R+ value than R−. This shows the effectiveness of the proposed adaptive scale factor FiG given in Equation (31) and the idea of increasing the scale factor during the evolution process. Hence, the main purpose of this comparison is to show that the HADENM algorithm with the adaptive scale factor FiG can achieve better performance than the same algorithm with a fixed value of F.Secondly, in order to show the efficiency of the proposed adaptive crossover rateCRG in Equation (32), the performance of the HADENM algorithm is compared with three different fixed values of CR, as follows:
(1)
HADENMCR=0.1, which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.1
(2)
HADENMCR=0.5,which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.5
(3)
HADENMCR=0.9, which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.9Table2 shows the summary of statistical results of applying the Wilcoxon test between the proposed HADENM and the above three algorithms.Table 2
Results of the Wilcoxon test for HADENM, HADENMCR=0.1, HADENMCR=0.5, and HADENMCR=0.9 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENMCR=0.1
264.5
150
0.7761
≈
HADENM versus HADENMCR=0.5
253
72
0.0135
+
HADENM versus HADENMCR=0.9
259
66
0.0079
+From the results in Table2, it can be seen that the proposed HADENM algorithm has significantly better performance than HADENMCR=0.5 and HADENMCR=0.9. Furthermore, in the case of HADENM versus HADENMCR=0.1, it can be observed that the proposed algorithm provides higher R+ values than R− in all considered cases. This shows that the proposed adaptive crossover rate CRG plays a vital role in determining the optimal crossover rate value for the considered optimization problem.


#### 10.1.2. Effectiveness of the Proposed Adaptive Mutation Operator

In this sub-subsection, to evaluate the effectiveness and improvement of the proposed adaptive mutation operator on the optimization performance of the HADENM algorithm, the experimental studies have been performed. In this way, two different versions of the HADENM algorithm have been tested and compared against the proposed one, such as(1)
HADENM-1, which has the same operators as HADENM, except that only the explorative mutations of the proposed adaptive mutation operator are applied. Hence, only the mutation operators DE/rand/1 and DE/rand/2 will be chosen randomly with the probability of 0.5
(2)
HADENM-2, which has the same operators as HADENM, except that only the exploitative mutations of the proposed adaptive mutation operator are applied. In this case, the mutation operators DE/best/1 and DE/current-to-best/1 will be randomly selected, with the probability of 0.5In Table3, the summary of statistical results of applying the Wilcoxon test between the proposed HADENM, HADENM-1, and HADENM-2 algorithms is shown.Table 3
Results of the Wilcoxon test for HADENM, HADENM-1, and HADENM-2 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENM-1
314.5
105.5
0.1019
≈
HADENM versus HADENM-2
283
182
0.2915
≈According to the Wilcoxon test, given in Table3, it can be observed that the proposed HADENM algorithm provides higher R+ values than R− in both cases. The obtained results indicate that the proposed adaptive mutation operator can effectively keep the balance between exploration and exploitation abilities, thus improving the overall optimization performance of the HADENM algorithm.


#### 10.1.3. Effectiveness of the Proposed Hybridization

To study the effects of the proposed hybridization between the ADE and NM algorithms, in this sub-subsection, the following experiment has been performed. In this regard, the performance of the HADENM algorithm without using the NM method has been compared to the performance of the proposed HADENM algorithm. Thus, the summary of statistical results of applying the Wilcoxon test between the two abovementioned algorithms is shown in Table4.Table 4
Results of the Wilcoxon test for the HADENM algorithm without using the NM method and the proposed HADENM algorithm at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus ADE
465
0
0.000
+According to the obtained statistical results in Table4, it can be concluded that the proposed HADENM algorithm has better performance compared to the algorithm without using the NM method. This shows that the NM method can further enhance the quality of the obtained solution, and in this way, it can improve the optimization performance of the HADENM algorithm.


#### 10.1.4. Comparison of the Proposed Improvements

In this sub-subsection, the experimental studies have been performed in order to demonstrate the effectiveness of the proposed HADENM algorithm. As there are more than two algorithms for comparison, the overall performance of the considered algorithms has been compared using the Friedman test. In this regard, Table5 shows the average ranks according to the Friedman test for the considered algorithms, using different values of the variance of the noise 10logσni2, where the proposed HADENM algorithm is selected as the base algorithm. The best ranks are shown in bold, and the second ranks are in italics.Table 5
Average ranks computed through the Friedman test for all algorithms across different values of the variance of the noise10logσni2, at the significance level of 0.05.
Algorithm
-40
-20
20
40
Mean ranking
Rank
HADENM
3.07
3.30
3.37
2.98
3.27
1
HADENM-2
3.87
3.70
3.47
3.02
3.63
2
HADENM-1
4.00
3.93
4.10
4.67
4.22
3
HADENMCR=0.1
4.73
4.53
3.80
4.28
4.27
4
HADENMF=0.9
4.03
5.07
4.80
4.68
4.64
5
HADENMCR=0.5
4.80
4.87
5.33
5.08
4.94
6
HADENMF=0.5
5.57
5.00
5.47
5.35
5.34
7
HADENMCR=0.9
5.93
5.60
5.67
5.93
5.69
8
Friedmanp value
0.000
0.003
0.000
0.000From Table5, it can be noted that the p values computed through the Friedman test for different values of the variance of the noise 10logσni2 are less than 0.05. Therefore, this means that there is a significant difference in performance between the considered algorithms. Furthermore, it can be observed that the proposed HADENM algorithm outperforms other considered algorithms for all values of 10logσni2, which demonstrates the effectiveness of the modifications proposed in this paper.


### 10.2. Localization Performance

The simulations are performed for the passive target localization system in the LOS environment, which consists of one transmitter, four receivers located at known positions, and the passive target, which is located at a different position for each simulation scenario. The transmitter is located atxT=00T m, while the receiver positions are x1r=6060T m, x2r=60−60T m, x3r=−60−60T m, and x4r=−6060T m, respectively. To evaluate the effectiveness of the HADENM algorithm to the measurement error, three scenarios are considered: (i) the position of the target is inside the convex hull formed by the set of receivers at the position x=525T m, (ii) the target is located outside of the convex hull at the position x=8050T m, and (iii) the position of the target is deployed randomly over an area 80×80 m2, and its position is different for each simulation run.Here, the RMSE is employed to evaluate the localization accuracy, which is expressed as(48)RMSE=1Nm∑n=1Nmx^n−x22,where the estimated value of the true target position is x^n and Nm=1000 is the number of simulation runs for a given variance of the noise σni2.Figure5 shows the results of the first scenario, where the RMSEs of the considered algorithms versus p=10logσni2 are plotted and compared against the CRLB.Figure 5
RMSEs versusp for a passive target at x=525T m.From Figure5, it can be observed that the RMSE performance of the HADENM algorithm is better compared to the DE and the CWLS, and the proposed algorithm attains the CRLB for all the considered ranges of p. Moreover, the DE and the CWLS algorithms can achieve the RMSE performance several dBs above the CRLB, for small measurement noise. However, as the noise level increases, i.e., p≥20 dB, the CWLS algorithm separates faster from the CRLB accuracy.The simulation results of the second scenario are presented in Figure6, where the RMSEs of all the considered algorithms versus measurement noise p are plotted.Figure 6
RMSEs versusp for a passive target at x=8050T m.According to the simulation results in Figure6, it can be observed that the HADENM algorithm attains the CRLB accuracy and has a superior localization performance compared to the other considered algorithms. Furthermore, it is concluded that there is a significant degradation in the localization accuracy of the CWLS and DE algorithms compared to the first simulation scenario. It is also noted that, the RMSE of the CWLS algorithm diverges rapidly from the CRLB for large values of the measurement noise, i.e., when p≥5 dB.Moreover, the RMSEs of the considered algorithms versusp are plotted in Figure 7 for the third simulation scenario, where the true position of the passive target is randomly generated within the considered area for each simulation run.Figure 7
RMSEs versusp with a random passive target position.As expected, the results of the third simulation scenario, presented in Figure7, show that the HADENM algorithm attains the CRLB accuracy and provides superior performance over the CWLS and DE algorithms. Therefore, summarizing the results of the three simulation scenarios, presented in Figures 5–7, it can be concluded that the proposed HADENM algorithm has a better localization performance compared to the considered algorithms and successfully attains the CRLB for every considered scenario. Furthermore, it is observed that the proposed algorithm is less sensitive to the changes in the geometric configuration of the transmitter, receivers, and target.For a better comparison of the simulation results, the cumulative distribution functions (CDFs) of the passive target localization error of the considered algorithms are investigated, with different variances of the measurement noiseσni2. The localization error (LE) is defined as the Euclidean distance between the estimated and the true position of the passive target, i.e.,
(49)LE=x^n−x2,∀n∈1,⋯,Nm.Figure8 represents the simulation results for the second scenario with the corresponding CDFs in terms of the localization error obtained for each algorithm for different variances of the measurement noise σni2.Figure 8
CDFs of passive target localization error of the considered algorithms.From the CDFs of the localization error, depicted in Figure8, it can be observed that for a fixed CDF percentage, e.g., 90%, under the measurement noise σni2=1 m2, the HADENM algorithm has the lowest localization error of 1.2 m, while the DE and CWLS algorithms have 1.3 m and 1.5 m, respectively. When the value of the measurement noise is increased to σni2=10 m2, the HADENM algorithm has the localization error of 3.8 m, while the DE and the CWLS show localization errors of 3.9 m and 4.2 m, respectively. For the higher value of measurement noise, e.g., σni2=20 m2, the HADENM algorithm has a localization error of 5.6 m, while the DE and CWLS algorithms have 5.8 m and 6.1 m, respectively. Based on the above, it is evident that the HADENM algorithm has a smaller localization error compared to the other considered algorithms.Finally, the influence of increasing the number of receivers on the localization accuracy of different algorithms is investigated for the first simulation scenario. In this regard, an arrangement ofNj=21 receivers uniformly distributed on a circle of radius R=602 m is considered, where the coordinates of the ith receiver are obtained as follows: xir=RcosφiRsinφiT, ∀i∈1,2,⋯,Nj. In this way, Figure 9 shows the RMSE performances of the considered algorithms versus the number of receivers, when the variance of the measurement noise is σni2=1 m2.Figure 9
RMSE performances of the considered algorithms versus the number of receivers.From Figure9, it is observed that with the increase of the number of receivers from 4 to 15, the RMSE performances of the CWLS, DE, and HADENM algorithms improve significantly. Further increase in the number of receivers does not explicitly enhance the localization accuracy. Thus, the difference between the RMSE values of the considered algorithms is smaller. Therefore, based on the obtained results, it can be concluded that the proposed HADENM algorithm is robust against variations in the topology and provides a superior performance among all the considered algorithms.


#### 10.2.1. Computational Complexity of the Considered Algorithms

In this sub-subsection, the computational complexity of the HADENM algorithm and other considered algorithms is presented, and the analysis of the average computation time is also given for the comparison. It should be noted that in this paper, the computational complexity of the considered algorithms is only analysed in a single generation. It has been previously shown that the computational complexity of the CWLS algorithm isOGmaxn+1 [8], while the conventional DE algorithm has the complexity of OGmaxNPn [42]. During one generation of the proposed HADENM algorithm, all the individuals in the population are sorted according to the objective function value, where the average computational complexity of this process is ONPlogNP. Then, all individuals in the population go through mutation, crossover, and selection operations, which have the computational complexity approximately equal to ONPn+Coc+Cos [43]. Here, Coc and Cos denote the costs of crossover and selection operations in the DE algorithm, respectively. Afterwards, the accuracy of the best solution is further improved by the NM algorithm, which has the computational complexity Onlogn [44]. Based on the above considerations, the computational complexity of the HADENM algorithm can be simplified to OGmaxNPn+nlogn.Next, the average computation times in searching for the global optimum have been determined on the same computer with 3.2 GHz CPU and 16 GB of RAM. In this way, based on the considered simulation scenarios, a comparison between the average computation times of the CWLS, DE and HADENM algorithms is shown in Table6.Table 6
Average computation times of the considered localization algorithms.
CWLS (ms)
DE (ms)
HADENM (ms)
Scenario 1
0.567
7.1
8.3
Scenario 2
0.466
6.4
7.8
Scenario 3
0.926
8.6
8.7Table6 shows that the CWLS algorithm has the fastest implementation among the considered algorithms, while there is no significant difference between the HADENM and DE algorithms. Based on the results of the analysis of localization accuracy, the results in Table 6 indicate that the proposed HADENM algorithm achieves the best compromise between the localization accuracy and the average computation time in searching for the global optimum.


## 10.1. A Parametric Study on HADENM Algorithm

In this subsection, the experiments are carried out to evaluate the performance of the HADENM algorithm by studying the effects of the proposed adaptive scale factor and crossover rate, as well as the proposed adaptive mutation operator on the optimization performance. Furthermore, the effects of the hybridization between the ADE and NM algorithms have been analysed. Then, the experiments are performed in order to verify that the performance is enhanced after combining the previously described improvements.To evaluate the performance of the proposed HADENM algorithm, the solution error measurefx^−fx∗ has been employed, where x^ denotes the best solution of the algorithm obtained in one run and x∗ represents the well-known global optimal solution of the corresponding ML objective function. All experiments for each algorithm have been independently run 30 times and statistical results are provided. From the statistical point of view, the quality of the obtained solutions has been analysed and compared using two nonparametric statistical hypothesis tests such as the Wilcoxon signed-rank test and Friedman test.Firstly, the Wilcoxon signed-rank test can be used to determine the significant differences between two samples obtained by the algorithms. This statistic test has been applied with a significance level ofα=0.05. Here, R+ denotes the sum of ranks for the problem in which the first algorithm outperformed the second and R− is the sum of ranks for the problem where the first algorithm performed worse than the second. According to [41], in the Wilcoxon signed-rank test for the null hypothesis, it is assumed that, “there is no difference between the mean results of the two samples.” On the contrary, the alternative hypothesis is that, “there is a difference in the mean results of the two samples.” In the statistical analysis, the p value is used and compared to the significance level. Thus, the null hypothesis can be rejected when the p value is less than or equal to α=0.05. Based on the obtained results of the statistical test, one of the following three signs (+, −, and ≈) has been assigned for the comparison between any two algorithms. Thus, the plus (+) sign denotes that the first algorithm is significantly better than the second, the minus (−) sign means that the first algorithm is significantly worse than the second, and the approximation (≈) sign denotes that there is no significant difference between two algorithms.The second test is the Friedman test, which obtains the ranks of all considered algorithms over every tested function, with the aim to find the significant differences in performance between two or more algorithms. In this statistical test, the algorithm with a minimum rank value is considered as the best performing algorithm, while the one with the highest rank value is considered as the worst. According to [41], the null hypothesis for Friedman test states that, “there is no difference among the performance of all algorithms,” whereas the alternative hypothesis states that, “there is a difference among the performance of all algorithms”.In this regard, the presentation of the obtained results of the considered evaluations is divided into four sub-subsections. In the first sub-subsection, the effectiveness of the proposed adaptive scale factor and crossover rate has been evaluated. In the second sub-subsection, the effectiveness of the proposed adaptive mutation operator has been considered. The third sub-subsection considers the hybridization of ADE and NM algorithms and the overall performance improvements. Finally, in the fourth sub-subsection, the statistical results of applying the Friedman test between the considered algorithms have been analysed.


### 10.1.1. Effectiveness of the Adaptive Scale Factor and Crossover Rate

In this sub-subsection, the experimental studies have been performed to evaluate the effectiveness and benefits of the proposed adaptive scale factorFiG (given in Equation (31)) and adaptive crossover rate CRG (given in Equation (32)) on the optimization performance of the HADENM algorithm. Firstly, the effectiveness of the adaptive scale factor is considered, where the performance of the HADENM algorithm is evaluated with two different fixed values of F, as follows:
(1)
HADENMF=0.5, which has the same operators as HADENM, except that the scale factor is set to a fixed value of F=0.5
(2)
HADENMF=0.9, which has the same operators as HADENM, except that the scale factor is set to a fixed value of F=0.9The summary of statistical results of applying the Wilcoxon test between the proposed HADENM and the above two algorithms is presented in Table1.Table 1
Results of the Wilcoxon test for HADENM, HADENMF=0.5, and HADENMF=0.9 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENMF=0.5
346
119
0.0177
+
HADENM versus HADENMF=0.9
282
183
0.3207
≈From the results in Table1, it can be seen that the proposed HADENM algorithm is significantly better than HADENMF=0.5. On the other hand, in the case of HADENM versus HADENMF=0.9, the proposed algorithm has a higher R+ value than R−. This shows the effectiveness of the proposed adaptive scale factor FiG given in Equation (31) and the idea of increasing the scale factor during the evolution process. Hence, the main purpose of this comparison is to show that the HADENM algorithm with the adaptive scale factor FiG can achieve better performance than the same algorithm with a fixed value of F.Secondly, in order to show the efficiency of the proposed adaptive crossover rateCRG in Equation (32), the performance of the HADENM algorithm is compared with three different fixed values of CR, as follows:
(1)
HADENMCR=0.1, which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.1
(2)
HADENMCR=0.5,which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.5
(3)
HADENMCR=0.9, which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.9Table2 shows the summary of statistical results of applying the Wilcoxon test between the proposed HADENM and the above three algorithms.Table 2
Results of the Wilcoxon test for HADENM, HADENMCR=0.1, HADENMCR=0.5, and HADENMCR=0.9 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENMCR=0.1
264.5
150
0.7761
≈
HADENM versus HADENMCR=0.5
253
72
0.0135
+
HADENM versus HADENMCR=0.9
259
66
0.0079
+From the results in Table2, it can be seen that the proposed HADENM algorithm has significantly better performance than HADENMCR=0.5 and HADENMCR=0.9. Furthermore, in the case of HADENM versus HADENMCR=0.1, it can be observed that the proposed algorithm provides higher R+ values than R− in all considered cases. This shows that the proposed adaptive crossover rate CRG plays a vital role in determining the optimal crossover rate value for the considered optimization problem.


### 10.1.2. Effectiveness of the Proposed Adaptive Mutation Operator

In this sub-subsection, to evaluate the effectiveness and improvement of the proposed adaptive mutation operator on the optimization performance of the HADENM algorithm, the experimental studies have been performed. In this way, two different versions of the HADENM algorithm have been tested and compared against the proposed one, such as(1)
HADENM-1, which has the same operators as HADENM, except that only the explorative mutations of the proposed adaptive mutation operator are applied. Hence, only the mutation operators DE/rand/1 and DE/rand/2 will be chosen randomly with the probability of 0.5
(2)
HADENM-2, which has the same operators as HADENM, except that only the exploitative mutations of the proposed adaptive mutation operator are applied. In this case, the mutation operators DE/best/1 and DE/current-to-best/1 will be randomly selected, with the probability of 0.5In Table3, the summary of statistical results of applying the Wilcoxon test between the proposed HADENM, HADENM-1, and HADENM-2 algorithms is shown.Table 3
Results of the Wilcoxon test for HADENM, HADENM-1, and HADENM-2 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENM-1
314.5
105.5
0.1019
≈
HADENM versus HADENM-2
283
182
0.2915
≈According to the Wilcoxon test, given in Table3, it can be observed that the proposed HADENM algorithm provides higher R+ values than R− in both cases. The obtained results indicate that the proposed adaptive mutation operator can effectively keep the balance between exploration and exploitation abilities, thus improving the overall optimization performance of the HADENM algorithm.


### 10.1.3. Effectiveness of the Proposed Hybridization

To study the effects of the proposed hybridization between the ADE and NM algorithms, in this sub-subsection, the following experiment has been performed. In this regard, the performance of the HADENM algorithm without using the NM method has been compared to the performance of the proposed HADENM algorithm. Thus, the summary of statistical results of applying the Wilcoxon test between the two abovementioned algorithms is shown in Table4.Table 4
Results of the Wilcoxon test for the HADENM algorithm without using the NM method and the proposed HADENM algorithm at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus ADE
465
0
0.000
+According to the obtained statistical results in Table4, it can be concluded that the proposed HADENM algorithm has better performance compared to the algorithm without using the NM method. This shows that the NM method can further enhance the quality of the obtained solution, and in this way, it can improve the optimization performance of the HADENM algorithm.


### 10.1.4. Comparison of the Proposed Improvements

In this sub-subsection, the experimental studies have been performed in order to demonstrate the effectiveness of the proposed HADENM algorithm. As there are more than two algorithms for comparison, the overall performance of the considered algorithms has been compared using the Friedman test. In this regard, Table5 shows the average ranks according to the Friedman test for the considered algorithms, using different values of the variance of the noise 10logσni2, where the proposed HADENM algorithm is selected as the base algorithm. The best ranks are shown in bold, and the second ranks are in italics.Table 5
Average ranks computed through the Friedman test for all algorithms across different values of the variance of the noise10logσni2, at the significance level of 0.05.
Algorithm
-40
-20
20
40
Mean ranking
Rank
HADENM
3.07
3.30
3.37
2.98
3.27
1
HADENM-2
3.87
3.70
3.47
3.02
3.63
2
HADENM-1
4.00
3.93
4.10
4.67
4.22
3
HADENMCR=0.1
4.73
4.53
3.80
4.28
4.27
4
HADENMF=0.9
4.03
5.07
4.80
4.68
4.64
5
HADENMCR=0.5
4.80
4.87
5.33
5.08
4.94
6
HADENMF=0.5
5.57
5.00
5.47
5.35
5.34
7
HADENMCR=0.9
5.93
5.60
5.67
5.93
5.69
8
Friedmanp value
0.000
0.003
0.000
0.000From Table5, it can be noted that the p values computed through the Friedman test for different values of the variance of the noise 10logσni2 are less than 0.05. Therefore, this means that there is a significant difference in performance between the considered algorithms. Furthermore, it can be observed that the proposed HADENM algorithm outperforms other considered algorithms for all values of 10logσni2, which demonstrates the effectiveness of the modifications proposed in this paper.


## 10.1.1. Effectiveness of the Adaptive Scale Factor and Crossover Rate

In this sub-subsection, the experimental studies have been performed to evaluate the effectiveness and benefits of the proposed adaptive scale factorFiG (given in Equation (31)) and adaptive crossover rate CRG (given in Equation (32)) on the optimization performance of the HADENM algorithm. Firstly, the effectiveness of the adaptive scale factor is considered, where the performance of the HADENM algorithm is evaluated with two different fixed values of F, as follows:
(1)
HADENMF=0.5, which has the same operators as HADENM, except that the scale factor is set to a fixed value of F=0.5
(2)
HADENMF=0.9, which has the same operators as HADENM, except that the scale factor is set to a fixed value of F=0.9The summary of statistical results of applying the Wilcoxon test between the proposed HADENM and the above two algorithms is presented in Table1.Table 1
Results of the Wilcoxon test for HADENM, HADENMF=0.5, and HADENMF=0.9 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENMF=0.5
346
119
0.0177
+
HADENM versus HADENMF=0.9
282
183
0.3207
≈From the results in Table1, it can be seen that the proposed HADENM algorithm is significantly better than HADENMF=0.5. On the other hand, in the case of HADENM versus HADENMF=0.9, the proposed algorithm has a higher R+ value than R−. This shows the effectiveness of the proposed adaptive scale factor FiG given in Equation (31) and the idea of increasing the scale factor during the evolution process. Hence, the main purpose of this comparison is to show that the HADENM algorithm with the adaptive scale factor FiG can achieve better performance than the same algorithm with a fixed value of F.Secondly, in order to show the efficiency of the proposed adaptive crossover rateCRG in Equation (32), the performance of the HADENM algorithm is compared with three different fixed values of CR, as follows:
(1)
HADENMCR=0.1, which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.1
(2)
HADENMCR=0.5,which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.5
(3)
HADENMCR=0.9, which has the same operators as HADENM, except that the crossover rate is set to a fixed value of CR=0.9Table2 shows the summary of statistical results of applying the Wilcoxon test between the proposed HADENM and the above three algorithms.Table 2
Results of the Wilcoxon test for HADENM, HADENMCR=0.1, HADENMCR=0.5, and HADENMCR=0.9 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENMCR=0.1
264.5
150
0.7761
≈
HADENM versus HADENMCR=0.5
253
72
0.0135
+
HADENM versus HADENMCR=0.9
259
66
0.0079
+From the results in Table2, it can be seen that the proposed HADENM algorithm has significantly better performance than HADENMCR=0.5 and HADENMCR=0.9. Furthermore, in the case of HADENM versus HADENMCR=0.1, it can be observed that the proposed algorithm provides higher R+ values than R− in all considered cases. This shows that the proposed adaptive crossover rate CRG plays a vital role in determining the optimal crossover rate value for the considered optimization problem.


## 10.1.2. Effectiveness of the Proposed Adaptive Mutation Operator

In this sub-subsection, to evaluate the effectiveness and improvement of the proposed adaptive mutation operator on the optimization performance of the HADENM algorithm, the experimental studies have been performed. In this way, two different versions of the HADENM algorithm have been tested and compared against the proposed one, such as(1)
HADENM-1, which has the same operators as HADENM, except that only the explorative mutations of the proposed adaptive mutation operator are applied. Hence, only the mutation operators DE/rand/1 and DE/rand/2 will be chosen randomly with the probability of 0.5
(2)
HADENM-2, which has the same operators as HADENM, except that only the exploitative mutations of the proposed adaptive mutation operator are applied. In this case, the mutation operators DE/best/1 and DE/current-to-best/1 will be randomly selected, with the probability of 0.5In Table3, the summary of statistical results of applying the Wilcoxon test between the proposed HADENM, HADENM-1, and HADENM-2 algorithms is shown.Table 3
Results of the Wilcoxon test for HADENM, HADENM-1, and HADENM-2 at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus HADENM-1
314.5
105.5
0.1019
≈
HADENM versus HADENM-2
283
182
0.2915
≈According to the Wilcoxon test, given in Table3, it can be observed that the proposed HADENM algorithm provides higher R+ values than R− in both cases. The obtained results indicate that the proposed adaptive mutation operator can effectively keep the balance between exploration and exploitation abilities, thus improving the overall optimization performance of the HADENM algorithm.


## 10.1.3. Effectiveness of the Proposed Hybridization

To study the effects of the proposed hybridization between the ADE and NM algorithms, in this sub-subsection, the following experiment has been performed. In this regard, the performance of the HADENM algorithm without using the NM method has been compared to the performance of the proposed HADENM algorithm. Thus, the summary of statistical results of applying the Wilcoxon test between the two abovementioned algorithms is shown in Table4.Table 4
Results of the Wilcoxon test for the HADENM algorithm without using the NM method and the proposed HADENM algorithm at a significance level of 0.05.
Algorithm
R+
R−
p value
Dec.
HADENM versus ADE
465
0
0.000
+According to the obtained statistical results in Table4, it can be concluded that the proposed HADENM algorithm has better performance compared to the algorithm without using the NM method. This shows that the NM method can further enhance the quality of the obtained solution, and in this way, it can improve the optimization performance of the HADENM algorithm.


## 10.1.4. Comparison of the Proposed Improvements

In this sub-subsection, the experimental studies have been performed in order to demonstrate the effectiveness of the proposed HADENM algorithm. As there are more than two algorithms for comparison, the overall performance of the considered algorithms has been compared using the Friedman test. In this regard, Table5 shows the average ranks according to the Friedman test for the considered algorithms, using different values of the variance of the noise 10logσni2, where the proposed HADENM algorithm is selected as the base algorithm. The best ranks are shown in bold, and the second ranks are in italics.Table 5
Average ranks computed through the Friedman test for all algorithms across different values of the variance of the noise10logσni2, at the significance level of 0.05.
Algorithm
-40
-20
20
40
Mean ranking
Rank
HADENM
3.07
3.30
3.37
2.98
3.27
1
HADENM-2
3.87
3.70
3.47
3.02
3.63
2
HADENM-1
4.00
3.93
4.10
4.67
4.22
3
HADENMCR=0.1
4.73
4.53
3.80
4.28
4.27
4
HADENMF=0.9
4.03
5.07
4.80
4.68
4.64
5
HADENMCR=0.5
4.80
4.87
5.33
5.08
4.94
6
HADENMF=0.5
5.57
5.00
5.47
5.35
5.34
7
HADENMCR=0.9
5.93
5.60
5.67
5.93
5.69
8
Friedmanp value
0.000
0.003
0.000
0.000From Table5, it can be noted that the p values computed through the Friedman test for different values of the variance of the noise 10logσni2 are less than 0.05. Therefore, this means that there is a significant difference in performance between the considered algorithms. Furthermore, it can be observed that the proposed HADENM algorithm outperforms other considered algorithms for all values of 10logσni2, which demonstrates the effectiveness of the modifications proposed in this paper.


## 10.2. Localization Performance

The simulations are performed for the passive target localization system in the LOS environment, which consists of one transmitter, four receivers located at known positions, and the passive target, which is located at a different position for each simulation scenario. The transmitter is located atxT=00T m, while the receiver positions are x1r=6060T m, x2r=60−60T m, x3r=−60−60T m, and x4r=−6060T m, respectively. To evaluate the effectiveness of the HADENM algorithm to the measurement error, three scenarios are considered: (i) the position of the target is inside the convex hull formed by the set of receivers at the position x=525T m, (ii) the target is located outside of the convex hull at the position x=8050T m, and (iii) the position of the target is deployed randomly over an area 80×80 m2, and its position is different for each simulation run.Here, the RMSE is employed to evaluate the localization accuracy, which is expressed as(48)RMSE=1Nm∑n=1Nmx^n−x22,where the estimated value of the true target position is x^n and Nm=1000 is the number of simulation runs for a given variance of the noise σni2.Figure5 shows the results of the first scenario, where the RMSEs of the considered algorithms versus p=10logσni2 are plotted and compared against the CRLB.Figure 5
RMSEs versusp for a passive target at x=525T m.From Figure5, it can be observed that the RMSE performance of the HADENM algorithm is better compared to the DE and the CWLS, and the proposed algorithm attains the CRLB for all the considered ranges of p. Moreover, the DE and the CWLS algorithms can achieve the RMSE performance several dBs above the CRLB, for small measurement noise. However, as the noise level increases, i.e., p≥20 dB, the CWLS algorithm separates faster from the CRLB accuracy.The simulation results of the second scenario are presented in Figure6, where the RMSEs of all the considered algorithms versus measurement noise p are plotted.Figure 6
RMSEs versusp for a passive target at x=8050T m.According to the simulation results in Figure6, it can be observed that the HADENM algorithm attains the CRLB accuracy and has a superior localization performance compared to the other considered algorithms. Furthermore, it is concluded that there is a significant degradation in the localization accuracy of the CWLS and DE algorithms compared to the first simulation scenario. It is also noted that, the RMSE of the CWLS algorithm diverges rapidly from the CRLB for large values of the measurement noise, i.e., when p≥5 dB.Moreover, the RMSEs of the considered algorithms versusp are plotted in Figure 7 for the third simulation scenario, where the true position of the passive target is randomly generated within the considered area for each simulation run.Figure 7
RMSEs versusp with a random passive target position.As expected, the results of the third simulation scenario, presented in Figure7, show that the HADENM algorithm attains the CRLB accuracy and provides superior performance over the CWLS and DE algorithms. Therefore, summarizing the results of the three simulation scenarios, presented in Figures 5–7, it can be concluded that the proposed HADENM algorithm has a better localization performance compared to the considered algorithms and successfully attains the CRLB for every considered scenario. Furthermore, it is observed that the proposed algorithm is less sensitive to the changes in the geometric configuration of the transmitter, receivers, and target.For a better comparison of the simulation results, the cumulative distribution functions (CDFs) of the passive target localization error of the considered algorithms are investigated, with different variances of the measurement noiseσni2. The localization error (LE) is defined as the Euclidean distance between the estimated and the true position of the passive target, i.e.,
(49)LE=x^n−x2,∀n∈1,⋯,Nm.Figure8 represents the simulation results for the second scenario with the corresponding CDFs in terms of the localization error obtained for each algorithm for different variances of the measurement noise σni2.Figure 8
CDFs of passive target localization error of the considered algorithms.From the CDFs of the localization error, depicted in Figure8, it can be observed that for a fixed CDF percentage, e.g., 90%, under the measurement noise σni2=1 m2, the HADENM algorithm has the lowest localization error of 1.2 m, while the DE and CWLS algorithms have 1.3 m and 1.5 m, respectively. When the value of the measurement noise is increased to σni2=10 m2, the HADENM algorithm has the localization error of 3.8 m, while the DE and the CWLS show localization errors of 3.9 m and 4.2 m, respectively. For the higher value of measurement noise, e.g., σni2=20 m2, the HADENM algorithm has a localization error of 5.6 m, while the DE and CWLS algorithms have 5.8 m and 6.1 m, respectively. Based on the above, it is evident that the HADENM algorithm has a smaller localization error compared to the other considered algorithms.Finally, the influence of increasing the number of receivers on the localization accuracy of different algorithms is investigated for the first simulation scenario. In this regard, an arrangement ofNj=21 receivers uniformly distributed on a circle of radius R=602 m is considered, where the coordinates of the ith receiver are obtained as follows: xir=RcosφiRsinφiT, ∀i∈1,2,⋯,Nj. In this way, Figure 9 shows the RMSE performances of the considered algorithms versus the number of receivers, when the variance of the measurement noise is σni2=1 m2.Figure 9
RMSE performances of the considered algorithms versus the number of receivers.From Figure9, it is observed that with the increase of the number of receivers from 4 to 15, the RMSE performances of the CWLS, DE, and HADENM algorithms improve significantly. Further increase in the number of receivers does not explicitly enhance the localization accuracy. Thus, the difference between the RMSE values of the considered algorithms is smaller. Therefore, based on the obtained results, it can be concluded that the proposed HADENM algorithm is robust against variations in the topology and provides a superior performance among all the considered algorithms.


### 10.2.1. Computational Complexity of the Considered Algorithms

In this sub-subsection, the computational complexity of the HADENM algorithm and other considered algorithms is presented, and the analysis of the average computation time is also given for the comparison. It should be noted that in this paper, the computational complexity of the considered algorithms is only analysed in a single generation. It has been previously shown that the computational complexity of the CWLS algorithm isOGmaxn+1 [8], while the conventional DE algorithm has the complexity of OGmaxNPn [42]. During one generation of the proposed HADENM algorithm, all the individuals in the population are sorted according to the objective function value, where the average computational complexity of this process is ONPlogNP. Then, all individuals in the population go through mutation, crossover, and selection operations, which have the computational complexity approximately equal to ONPn+Coc+Cos [43]. Here, Coc and Cos denote the costs of crossover and selection operations in the DE algorithm, respectively. Afterwards, the accuracy of the best solution is further improved by the NM algorithm, which has the computational complexity Onlogn [44]. Based on the above considerations, the computational complexity of the HADENM algorithm can be simplified to OGmaxNPn+nlogn.Next, the average computation times in searching for the global optimum have been determined on the same computer with 3.2 GHz CPU and 16 GB of RAM. In this way, based on the considered simulation scenarios, a comparison between the average computation times of the CWLS, DE and HADENM algorithms is shown in Table6.Table 6
Average computation times of the considered localization algorithms.
CWLS (ms)
DE (ms)
HADENM (ms)
Scenario 1
0.567
7.1
8.3
Scenario 2
0.466
6.4
7.8
Scenario 3
0.926
8.6
8.7Table6 shows that the CWLS algorithm has the fastest implementation among the considered algorithms, while there is no significant difference between the HADENM and DE algorithms. Based on the results of the analysis of localization accuracy, the results in Table 6 indicate that the proposed HADENM algorithm achieves the best compromise between the localization accuracy and the average computation time in searching for the global optimum.


## 10.2.1. Computational Complexity of the Considered Algorithms

In this sub-subsection, the computational complexity of the HADENM algorithm and other considered algorithms is presented, and the analysis of the average computation time is also given for the comparison. It should be noted that in this paper, the computational complexity of the considered algorithms is only analysed in a single generation. It has been previously shown that the computational complexity of the CWLS algorithm isOGmaxn+1 [8], while the conventional DE algorithm has the complexity of OGmaxNPn [42]. During one generation of the proposed HADENM algorithm, all the individuals in the population are sorted according to the objective function value, where the average computational complexity of this process is ONPlogNP. Then, all individuals in the population go through mutation, crossover, and selection operations, which have the computational complexity approximately equal to ONPn+Coc+Cos [43]. Here, Coc and Cos denote the costs of crossover and selection operations in the DE algorithm, respectively. Afterwards, the accuracy of the best solution is further improved by the NM algorithm, which has the computational complexity Onlogn [44]. Based on the above considerations, the computational complexity of the HADENM algorithm can be simplified to OGmaxNPn+nlogn.Next, the average computation times in searching for the global optimum have been determined on the same computer with 3.2 GHz CPU and 16 GB of RAM. In this way, based on the considered simulation scenarios, a comparison between the average computation times of the CWLS, DE and HADENM algorithms is shown in Table6.Table 6
Average computation times of the considered localization algorithms.
CWLS (ms)
DE (ms)
HADENM (ms)
Scenario 1
0.567
7.1
8.3
Scenario 2
0.466
6.4
7.8
Scenario 3
0.926
8.6
8.7Table6 shows that the CWLS algorithm has the fastest implementation among the considered algorithms, while there is no significant difference between the HADENM and DE algorithms. Based on the results of the analysis of localization accuracy, the results in Table 6 indicate that the proposed HADENM algorithm achieves the best compromise between the localization accuracy and the average computation time in searching for the global optimum.


## 11. Conclusion

In this paper, the passive target localization problem has been considered, based on TOA measurements, for the system with one transmitter and a set of receivers. In order to solve this nonlinear and nonconvex localization problem even in highly noisy environments, an improved HADENM algorithm, based on the hybridization of the ADE and NM algorithms, has been proposed. An adaptive mutation operator has been developed to provide a good balance between the global exploration and local exploitation abilities. In addition, to further enhance the optimization performance, the control parameters of the ADE are adaptively updated during the evolution process. Furthermore, the exploitation ability is enhanced by the NM method, which improves the accuracy of the best solution previously obtained by the ADE algorithm. To evaluate the benefits of the proposed modifications on the optimization performance, statistical analysis has been conducted.Based on the comparison results between the HADENM algorithm and its versions, it can be concluded that the modifications proposed in this paper can improve the overall optimization performance. Furthermore, the simulation results show that the proposed HADENM algorithm provides superior localization performance, under different measurement noise conditions, compared to the DE and CWLS algorithms. Moreover, the HADENM algorithm attains the CRLB accuracy and exhibits better robustness against variations in network topology and under high measurement noise. In this way, the HADENM algorithm has the ability to provide both accuracy and robustness compared to the other considered algorithms.Future research can be focused on further improving the performance of the proposed algorithm and applying it to other complex optimization problems in WSNs.



---

*Source: 3482463-2020-02-27.xml*