# An External Selection Mechanism for Differential Evolution Algorithm

Authors: Haigang Zhang; Da Wang
Journal: Computational Intelligence and Neuroscience
 (2022)

---

## Abstract

The procedures of differential evolution algorithm can be summarized as population initialization, mutation, crossover, and selection. However, successful solutions generated by each iteration have not been fully utilized to our best knowledge. In this study, an external selection mechanism (ESM) is presented to improve differential evolution (DE) algorithm performance. The proposed method stores successful solutions of each iteration into an archive. When the individual is in a state of stagnation, the parents for mutation operation are selected from the archive to restore the algorithm’s search. Most significant of all, a crowding entropy diversity measurement in fitness landscape is proposed, cooperated with fitness rank, to preserve the diversity and superiority of the archive. The ESM can be integrated into existing algorithms to improve the algorithm’s ability to escape the situation of stagnation. CEC2017 benchmark functions are used for verification of the proposed mechanism’s performance. Experimental results show that the ESM is universal, which can improve the accuracy of DE and its variant algorithms simultaneously.


---

## Body

## 1. Introduction

Differential evolution (DE) algorithm [1], introduced by Storn and Price in 1995, is one of the most popular evolutionary algorithms. DE is easy to be understood, and the principle is simple. Still, it demonstrates good optimization ability and is used in many single-objective optimization problems successfully, including continuous optimization [2], discrete optimization [3], constrained optimization [4], and unconstrained optimization [5].Different from other meta-heuristic algorithms [6], the procedures of differential evolution algorithms can be summarized as population initialization, mutation, crossover, and selection. After each iteration, promising solutions with better fitness values are selected to survive to the next iteration. These solutions can be called successful solutions. DE repeats this procedure until a predefined termination criterion is reached. However, DE may suffer stagnation during the evolution process, which means stop generating successful solutions and converging to a fixed point [7]. When the population stagnates, taking appropriate strategies to restore the search should theoretically further improve the algorithm’s performance.Successful solutions are generally located in valleys of the fitness landscape. As the iteration progresses, successful solutions containing useful information may be discarded and underutilized. We, therefore, propose an external selection mechanism (ESM). First, ESM stores successful solutions of each iteration into an archive. When the individual is in a state of stagnation, the parents for mutation operation are selected from the archive. In this way, the population can regain diversity and restore searchability. Second, the diversity and superiority of the archive are relatively significant because they can directly affect the performance of the algorithm. Therefore, a rule is proposed to update the archive based on novel crowding entropy diversity measurement in the fitness landscape and fitness rank.To verify and analyze the effectiveness of the ESM, we performed a series of experiments and comparisons on the CEC2017 benchmark set [8], incorporating three classic DEs and five state-of-the-art DE variants. Results indicate that the ESM can effectively improve the algorithm performance without increasing the computational complexity.The rest of this study is organized as follows: Section2 introduces the canonical DE algorithm, including its typical mutation, crossover, and selection operators. Section 3 reviews related work. The proposed ESM procedures are introduced in Section 4. The effectiveness of the proposed mechanism is discussed in Section 5 based on the computational results. Conclusions and future work are illustrated in Section 6.


## 2. Differential Evolution

In this section, we introduce the basic differential evolution algorithm, including the well-known mutation strategyDE/rand/1 [1] and other widely used mutation strategies.


### 2.1. Initialization

An initial random populationP consists of NP individuals, each represented by Xit=xi,1t,xi,2t,…,xi,Dt|i=1,2,…,NP, where t=0,1,2,…,Tmax is the iteration number, and D is the number of dimensions in the solution space. In DE, uniformly distributed random functions are used to generate initial solutions xi,j0=xj,min+rand0,1⋅xj,max−xj,min, where xj,max and xj,min are the maximum and minimum boundary values, respectively, on the corresponding jth dimension.


### 2.2. Mutation

At iterationt, for each target vector Xit, a mutant vector Vit is generated according to the following:(1)DE/rand/1:Vit=Xr1t+F⋅Xr2t−Xr3t.Other widely used mutations strategies include(2)DEbest/1:Vit=Xbestt+F⋅Xr1t−Xr2t,DEbest/2:Vit=Xbestt+F⋅Xr1t−Xr2t+F⋅Xr3t−Xr4t,DErand/2:Vit=Xr1t+F⋅Xr2t−Xr3t+F⋅Xr4t−Xr5t,DEcurrent−to−best/1:Vit=Xit+F⋅Xbestt−Xit+F⋅Xr1t−Xr2t,DEcurrent−to−rand/1:Vit=Xit+F⋅Xr1t−Xit+F⋅Xr2t−Xr3t,where r1,r2,r3,r4, and r5 are mutually different integers randomly generated from the range 1,2,…,NP, and they are different from i (r1≠r2≠r3≠i). Xbestt is the individual vector with the best fitness value in the population at iteration t.


### 2.3. Crossover

We illustrate the binomial crossover, in which the target vectorXit and donor vector Vit exchange elements according to the following rules:(3)ui,jt=vi,jt,ifrand0,1<CR orj=jrand,xi,jt,otherwise,where the crossover rate CR∈0,1 is a uniformly distributed random integer in 1,D that ensures at least one element of the trial vector is inherited from the donor vector.


### 2.4. Selection

The greedy selection strategy is generally used in DE. The variable,Xit+1, is assigned when the fitness value of the trial vector, Uit, is equal to or better than Xit, and otherwise Xit is reserved:(4)Xit+1=Uit,fUit≤fXit,Xit,otherwise,where fx is the fitness function.


## 2.1. Initialization

An initial random populationP consists of NP individuals, each represented by Xit=xi,1t,xi,2t,…,xi,Dt|i=1,2,…,NP, where t=0,1,2,…,Tmax is the iteration number, and D is the number of dimensions in the solution space. In DE, uniformly distributed random functions are used to generate initial solutions xi,j0=xj,min+rand0,1⋅xj,max−xj,min, where xj,max and xj,min are the maximum and minimum boundary values, respectively, on the corresponding jth dimension.


## 2.2. Mutation

At iterationt, for each target vector Xit, a mutant vector Vit is generated according to the following:(1)DE/rand/1:Vit=Xr1t+F⋅Xr2t−Xr3t.Other widely used mutations strategies include(2)DEbest/1:Vit=Xbestt+F⋅Xr1t−Xr2t,DEbest/2:Vit=Xbestt+F⋅Xr1t−Xr2t+F⋅Xr3t−Xr4t,DErand/2:Vit=Xr1t+F⋅Xr2t−Xr3t+F⋅Xr4t−Xr5t,DEcurrent−to−best/1:Vit=Xit+F⋅Xbestt−Xit+F⋅Xr1t−Xr2t,DEcurrent−to−rand/1:Vit=Xit+F⋅Xr1t−Xit+F⋅Xr2t−Xr3t,where r1,r2,r3,r4, and r5 are mutually different integers randomly generated from the range 1,2,…,NP, and they are different from i (r1≠r2≠r3≠i). Xbestt is the individual vector with the best fitness value in the population at iteration t.


## 2.3. Crossover

We illustrate the binomial crossover, in which the target vectorXit and donor vector Vit exchange elements according to the following rules:(3)ui,jt=vi,jt,ifrand0,1<CR orj=jrand,xi,jt,otherwise,where the crossover rate CR∈0,1 is a uniformly distributed random integer in 1,D that ensures at least one element of the trial vector is inherited from the donor vector.


## 2.4. Selection

The greedy selection strategy is generally used in DE. The variable,Xit+1, is assigned when the fitness value of the trial vector, Uit, is equal to or better than Xit, and otherwise Xit is reserved:(4)Xit+1=Uit,fUit≤fXit,Xit,otherwise,where fx is the fitness function.


## 3. Related Work

Most research on improving DE has focused on mutation operator [9]. In recent decades, there have been many mutation operators with distinct search performance that have been proposed. Fan et al. [10] proposed a triangular mutation strategy, which was considered a local search operator. Zhang et al. [11] introduced a new DE variant, JADE, improving optimization performance by a new mutation strategy, DE/current−to−pbest/1. DE/current−to−pbest/1 is one of the most successful mutation operators because of its relatively balanced performance between exploration and exploitation. Wang et al. [12] proposed a new mutation strategy called DE/current−to−lbest/1 based on the values near the current parameter to keep balance of exploitation and exploration capabilities during the differential evolution. A novel DE algorithm with intersect mutation operation called intersect mutation differential evolution (IMDE) was proposed [13] to further improve the performance of the standard DE algorithm. Mohamed et al. [14] proposed two mutation strategies, DE/current−to−ord_best/1 and DE/current−to−ord_pbest/1. Both of the proposed mutation strategies were based on ordering three selected vectors to achieve different search behaviors. Deng et al. [15] proposed a novel DE variant called DCDE based on a dynamic combination-based mutation operator to achieve an appropriate balance between global exploration ability and local exploitation ability.In addition, novel neighborhood-based [16], dimension-based [17], opposition-based [18], and regeneration-based [19] mechanisms also showed their effectiveness on improving the performance of DE and its variants. In this study, we focus on the related work on archive-based techniques for DE. In the multiobjective evolutionary algorithm, a subpopulation called the external archive is used to store nondominated solutions that have been found during the search [20, 21], employing an archive is almost standard in multiobjective optimization [22]. The goal of single-objective global numerical optimization is to find a global optimum in decision space of a given objective function. Recently, archive is also introduced into global numerical optimization. In JADE [11], a set of recently explored inferior solutions are archived, and their difference from the current population is considered as a promising direction toward the optimum. The archive is applied to the second difference vector of the mutation strategy in JADE. Zhou et al. [23] proposed a DE framework with guiding archive (GAR-DE), GAR-DE chose the base vector of mutation strategy from the archive to help DE escape from the situation of stagnation, and Manhattan distance metric was used to maintain the diversity of the archive. Zeng et al. [24] proposed a new selection operator (NSO) for DE, which archived not only the successful updated solutions but also the discarded trial solutions. The NSO focused on selecting appropriate solution from the archive to survive to the next generation.As can be seen from preceding explanation, the differences among archive-based mechanisms mainly consist in two aspects: what information/solutions should be stored in the archive; how to maintain and update the archive; and how to use the archive to guide the search. The most striking feature of the proposed ESM is the use of information entropy rather than distance metric [23, 25] to maintain the diversity of the archive.Shannon defined the information entropy theory [26] for the first time in 1948. In information theory, entropy is used to measure the expected value of a random variable. As the evaluation scale of information quantity of stochastic process, information entropy has been extended to a general and effective tool to solve difficult numerical problems and uncertain polynomial combinatorial optimization problems [27, 28]. An information theoretic technique was adopted to analyze the ruggedness of a continuous fitness landscape in [29, 30]. Petalas et al. [27] proposed a memetic particle swarm optimization algorithm that exploits Shannon’s information entropy for decision-making in swarm level, as well as a probabilistic decision-making scheme in particle level, for determining when and where local search is applied. As can be seen from the previous studies, entropy actually reflects the degree of chaos of system. Therefore, entropy principle can be a promising method applied to evolutionary algorithm to ensure population diversity.


## 4. The Proposed ESM

In this section, we discuss the characteristics of the successful solution archive and archive-updating procedure and describe the ESM implementation.


### 4.1. Successful Solution Archive

From Section2, we can know that DE selects one vector from the trial vector and parent vector to survive to the next generation, and the survivor is known as the successful solution. As shown in Figure 1(a), in the two-dimensional solution space of a multimodal problem, the population is evenly distributed at the initial stage. As the iteration continues, from Figure 1(b) we can conclude that superior successful solutions (red points) can be generated in possible optimal regions. However, in Figure 1(c), the random and greedy selection operation makes the population gradually converge to a local optimum. Some successful solutions may be discarded without being fully exploited because of the greedy selection operation as the iteration progresses. An appropriate mechanism to utilize the successful solutions with high diversity contribution can make the algorithm reexplore the missed potential optimal solution. To some extent, it is also a remedy for the greedy selection operation.Figure 1
Illustration of the successful solution distribution in a two-dimensional space: (a) at the initial stage; (b) at the intermediate stage; (c) at the terminal stage.
(a)(b)(c)Considering the above analysis, as depicted in Figure2, an external archive is designed to store successful solutions generated during each iteration. When an individual is in a state of stagnation, the archive is activated for parent selection of mutation operation. For instance, the classic mutation strategy DE/rand/1 can be modified to(5)Vit=Xar1t+F·Xar2t−Xar3t,iεI,where ar1,ar2,ar3∈At∧ar1≠ar2≠ar3, At is the successful solution archive at iteration t. A simple but efficient method is adopted at this point, when the number of steps for the individual’s fitness value to stop updating reaches a predetermined value Q, the individual is considered in a stagnant state [31].Figure 2
Illustration of the successful solution archive.


### 4.2. Archive Updating Procedure

As the successful solution archive has been established, the current problem is how to update the archive to keep its diversity and superiority in the whole evolution process. An archive containing diverse successful solutions can drag the individual out of stagnation without crippling the algorithm’s performance.


#### 4.2.1. Individual Diversity Measure

To our best knowledge, in single-objective optimization, diversity is usually defined by the spatial distribution of the whole population [32], but there is few diversity contribution measures for a single individual to whole population. Wang et al. [33] proposed an Euclidean distance-based diversity metric for each individual, but a huge computational cost is needed to calculate the distance between every two individuals. To utilize diversity information while reducing the computational cost, we can estimate the diversity by fitness value distribution on individual level. Fitness distance can represent the spatial distance to some extent [34].In multiobjective optimization, Deb et al. [35] proposed a crowding distance measure method to get an estimation of the density of solutions surrounding a particular solution in the population, by calculating the average distance of two points on either side of this point along each of the objectives. However, the distribution of this point is not considered, and the only average distance may not accurately reflect the crowding degree. With utilizing the characteristic of fitness distance, we can extend the crowding distance measure to single-objective optimization. For example, in Figure 3, the crowding distance of the solution A can be calculated as 3L+1L=4L (according to crowding distance calculation method in [35]), which is the same as the solution B 2L+2L=4L, whereas solution B obviously has a higher diversity contribution than solution A.Figure 3
Crowding distance presented in fitness landscape.


#### 4.2.2. Archive Update Procedure

From the above statement, this study proposes a crowding entropy diversity measurement in the fitness landscape at individual level. As shown in Figure4, we sort the archive population in the ascending order of fitness during each iteration, then, a fitness distance estimation operator can be formulated as follows:(6)Dli=fi+1−fi,(7)Dui=fi−fi−1,where fi is the fitness value of the i th individual. Dli and Dui can be further normalized, then, the crowding entropy of the i th individual can be defined as follows:(8)Pi1=DliDli+Dui,(9)Pi2=DuiDli+Dui,(10)CEi=−∑j=12Pijlog2Pij,where CEi∈0,1, the larger CE means that the individual distribution is more uniform, so it has a higher diversity contribution. According to Figure 3, we can further calculate CEA=0.8 and CEB=1, this is also consistent with their diversity contribution. Since the best and worst individuals locate in the boundary scope, there is only one neighbor around them, respectively. Therefore, we set CE1=CE2 and CENP=CENP−1 for the ESM to function normally. To keep the superiority of the archive, fitness value rank is further introduced:(11)FRi=i,i=1,2,…,NP,(12)Ri=1−FRiNP+1·CEi,where Ri represents the priority of the individual to be updated. An individual in the archive with a higher rank (i.e., far from the global optimum) or with lower entropy (i.e., located in a dense area) should have a higher probability of being updated.Figure 4
Crowding entropy presented in fitness landscape.


### 4.3. The Implementation of the ESM

Based on the ESM described above, we present the ESM updating procedure in Algorithm 1 and the complete implementation based on classic DE with mutation strategyDE/rand/1 in Algorithm 2. Different from the existing alternative for the selection of parents [33, 36], the proposed mechanism is an adaptive approach that considers the feedback of the recently successful solutions so that the robustness of the DE algorithm has the potential to be enhanced by dynamically adapting promising parents for the situation of stagnation. Moreover, the proposed ESM is computationally efficient. The time complexities of updating operations of the archive are ONP. Additionally, the selection of parents, counting of the consecutive unsuccessful updates, and stagnation detection do not increase the overall time complexity.Algorithm 1: The pseudo code of ESM updating procedure.
(1)
Input: successful solution Xi, the archive A;(2)
Output: the archive A;(3)
FORi=1 to NP(4)
CalculateRi of Ai by equations (6)–(12);(5)
END FOR(6)
Get the indexj of the smallest R;(7)
Replace the solutionAj with Xi;Algorithm 2: The pseudo code of DE/rand/1 with the proposed ESM.
(1)
Step 1: Initialization(2)
Set iterationt=0, fitness evaluations FEs=0, stagnation threshold Q=120.(3)
Initialize a population ofNP individuals P0=X10,…,XNP0 with(4)
Xi0=xi,1t,xi,2t,…,xi,Dt;(5)
Set parameterF=0.5 and CR=0.5, stagnation iterations q10=⋯=qNP0  = 0, archive(6)
A0=P0;(7)
Step 2: Evolution Iteration(8)
WHILEFEs<Max_FEs(9)
DO(10)
Step 2.1: Mutation Operation(11)
FORi=1 to NP(12)
Randomly selectr1≠r2≠r3 from the set I=1,2,…,NP;(13)
IFqit<Q(14)
Generate donor vectorsVit by (1);(15)
ELSE(16)
Generate donor vectorsVit by (5);(17)
END(18)
END FOR(19)
Step 2.2: Crossover Operation(20)
FORi=1 to NP(21)
Generate trial vectorsUit by (3);(22)
END FOR(23)
Step 2.3: Selection and Archive Update(24)
FORi=1 to NP(25)
IFfUit<fXit(26)
qit+1=0;(27)
UpdateA according to Algorithm 1 with Uit;(28)
ELSE(29)
qit+1=qit+1;(30)
END IF(31)
Generate new vectorsXit+1 by (4);(32)
END FOR(33)
Step 2.4: Increase the Count(34)
FEs=FEs+NP;(35)
t=t+1(36)
END WHILE


## 4.1. Successful Solution Archive

From Section2, we can know that DE selects one vector from the trial vector and parent vector to survive to the next generation, and the survivor is known as the successful solution. As shown in Figure 1(a), in the two-dimensional solution space of a multimodal problem, the population is evenly distributed at the initial stage. As the iteration continues, from Figure 1(b) we can conclude that superior successful solutions (red points) can be generated in possible optimal regions. However, in Figure 1(c), the random and greedy selection operation makes the population gradually converge to a local optimum. Some successful solutions may be discarded without being fully exploited because of the greedy selection operation as the iteration progresses. An appropriate mechanism to utilize the successful solutions with high diversity contribution can make the algorithm reexplore the missed potential optimal solution. To some extent, it is also a remedy for the greedy selection operation.Figure 1
Illustration of the successful solution distribution in a two-dimensional space: (a) at the initial stage; (b) at the intermediate stage; (c) at the terminal stage.
(a)(b)(c)Considering the above analysis, as depicted in Figure2, an external archive is designed to store successful solutions generated during each iteration. When an individual is in a state of stagnation, the archive is activated for parent selection of mutation operation. For instance, the classic mutation strategy DE/rand/1 can be modified to(5)Vit=Xar1t+F·Xar2t−Xar3t,iεI,where ar1,ar2,ar3∈At∧ar1≠ar2≠ar3, At is the successful solution archive at iteration t. A simple but efficient method is adopted at this point, when the number of steps for the individual’s fitness value to stop updating reaches a predetermined value Q, the individual is considered in a stagnant state [31].Figure 2
Illustration of the successful solution archive.


## 4.2. Archive Updating Procedure

As the successful solution archive has been established, the current problem is how to update the archive to keep its diversity and superiority in the whole evolution process. An archive containing diverse successful solutions can drag the individual out of stagnation without crippling the algorithm’s performance.


### 4.2.1. Individual Diversity Measure

To our best knowledge, in single-objective optimization, diversity is usually defined by the spatial distribution of the whole population [32], but there is few diversity contribution measures for a single individual to whole population. Wang et al. [33] proposed an Euclidean distance-based diversity metric for each individual, but a huge computational cost is needed to calculate the distance between every two individuals. To utilize diversity information while reducing the computational cost, we can estimate the diversity by fitness value distribution on individual level. Fitness distance can represent the spatial distance to some extent [34].In multiobjective optimization, Deb et al. [35] proposed a crowding distance measure method to get an estimation of the density of solutions surrounding a particular solution in the population, by calculating the average distance of two points on either side of this point along each of the objectives. However, the distribution of this point is not considered, and the only average distance may not accurately reflect the crowding degree. With utilizing the characteristic of fitness distance, we can extend the crowding distance measure to single-objective optimization. For example, in Figure 3, the crowding distance of the solution A can be calculated as 3L+1L=4L (according to crowding distance calculation method in [35]), which is the same as the solution B 2L+2L=4L, whereas solution B obviously has a higher diversity contribution than solution A.Figure 3
Crowding distance presented in fitness landscape.


### 4.2.2. Archive Update Procedure

From the above statement, this study proposes a crowding entropy diversity measurement in the fitness landscape at individual level. As shown in Figure4, we sort the archive population in the ascending order of fitness during each iteration, then, a fitness distance estimation operator can be formulated as follows:(6)Dli=fi+1−fi,(7)Dui=fi−fi−1,where fi is the fitness value of the i th individual. Dli and Dui can be further normalized, then, the crowding entropy of the i th individual can be defined as follows:(8)Pi1=DliDli+Dui,(9)Pi2=DuiDli+Dui,(10)CEi=−∑j=12Pijlog2Pij,where CEi∈0,1, the larger CE means that the individual distribution is more uniform, so it has a higher diversity contribution. According to Figure 3, we can further calculate CEA=0.8 and CEB=1, this is also consistent with their diversity contribution. Since the best and worst individuals locate in the boundary scope, there is only one neighbor around them, respectively. Therefore, we set CE1=CE2 and CENP=CENP−1 for the ESM to function normally. To keep the superiority of the archive, fitness value rank is further introduced:(11)FRi=i,i=1,2,…,NP,(12)Ri=1−FRiNP+1·CEi,where Ri represents the priority of the individual to be updated. An individual in the archive with a higher rank (i.e., far from the global optimum) or with lower entropy (i.e., located in a dense area) should have a higher probability of being updated.Figure 4
Crowding entropy presented in fitness landscape.


## 4.2.1. Individual Diversity Measure

To our best knowledge, in single-objective optimization, diversity is usually defined by the spatial distribution of the whole population [32], but there is few diversity contribution measures for a single individual to whole population. Wang et al. [33] proposed an Euclidean distance-based diversity metric for each individual, but a huge computational cost is needed to calculate the distance between every two individuals. To utilize diversity information while reducing the computational cost, we can estimate the diversity by fitness value distribution on individual level. Fitness distance can represent the spatial distance to some extent [34].In multiobjective optimization, Deb et al. [35] proposed a crowding distance measure method to get an estimation of the density of solutions surrounding a particular solution in the population, by calculating the average distance of two points on either side of this point along each of the objectives. However, the distribution of this point is not considered, and the only average distance may not accurately reflect the crowding degree. With utilizing the characteristic of fitness distance, we can extend the crowding distance measure to single-objective optimization. For example, in Figure 3, the crowding distance of the solution A can be calculated as 3L+1L=4L (according to crowding distance calculation method in [35]), which is the same as the solution B 2L+2L=4L, whereas solution B obviously has a higher diversity contribution than solution A.Figure 3
Crowding distance presented in fitness landscape.


## 4.2.2. Archive Update Procedure

From the above statement, this study proposes a crowding entropy diversity measurement in the fitness landscape at individual level. As shown in Figure4, we sort the archive population in the ascending order of fitness during each iteration, then, a fitness distance estimation operator can be formulated as follows:(6)Dli=fi+1−fi,(7)Dui=fi−fi−1,where fi is the fitness value of the i th individual. Dli and Dui can be further normalized, then, the crowding entropy of the i th individual can be defined as follows:(8)Pi1=DliDli+Dui,(9)Pi2=DuiDli+Dui,(10)CEi=−∑j=12Pijlog2Pij,where CEi∈0,1, the larger CE means that the individual distribution is more uniform, so it has a higher diversity contribution. According to Figure 3, we can further calculate CEA=0.8 and CEB=1, this is also consistent with their diversity contribution. Since the best and worst individuals locate in the boundary scope, there is only one neighbor around them, respectively. Therefore, we set CE1=CE2 and CENP=CENP−1 for the ESM to function normally. To keep the superiority of the archive, fitness value rank is further introduced:(11)FRi=i,i=1,2,…,NP,(12)Ri=1−FRiNP+1·CEi,where Ri represents the priority of the individual to be updated. An individual in the archive with a higher rank (i.e., far from the global optimum) or with lower entropy (i.e., located in a dense area) should have a higher probability of being updated.Figure 4
Crowding entropy presented in fitness landscape.


## 4.3. The Implementation of the ESM

Based on the ESM described above, we present the ESM updating procedure in Algorithm 1 and the complete implementation based on classic DE with mutation strategyDE/rand/1 in Algorithm 2. Different from the existing alternative for the selection of parents [33, 36], the proposed mechanism is an adaptive approach that considers the feedback of the recently successful solutions so that the robustness of the DE algorithm has the potential to be enhanced by dynamically adapting promising parents for the situation of stagnation. Moreover, the proposed ESM is computationally efficient. The time complexities of updating operations of the archive are ONP. Additionally, the selection of parents, counting of the consecutive unsuccessful updates, and stagnation detection do not increase the overall time complexity.Algorithm 1: The pseudo code of ESM updating procedure.
(1)
Input: successful solution Xi, the archive A;(2)
Output: the archive A;(3)
FORi=1 to NP(4)
CalculateRi of Ai by equations (6)–(12);(5)
END FOR(6)
Get the indexj of the smallest R;(7)
Replace the solutionAj with Xi;Algorithm 2: The pseudo code of DE/rand/1 with the proposed ESM.
(1)
Step 1: Initialization(2)
Set iterationt=0, fitness evaluations FEs=0, stagnation threshold Q=120.(3)
Initialize a population ofNP individuals P0=X10,…,XNP0 with(4)
Xi0=xi,1t,xi,2t,…,xi,Dt;(5)
Set parameterF=0.5 and CR=0.5, stagnation iterations q10=⋯=qNP0  = 0, archive(6)
A0=P0;(7)
Step 2: Evolution Iteration(8)
WHILEFEs<Max_FEs(9)
DO(10)
Step 2.1: Mutation Operation(11)
FORi=1 to NP(12)
Randomly selectr1≠r2≠r3 from the set I=1,2,…,NP;(13)
IFqit<Q(14)
Generate donor vectorsVit by (1);(15)
ELSE(16)
Generate donor vectorsVit by (5);(17)
END(18)
END FOR(19)
Step 2.2: Crossover Operation(20)
FORi=1 to NP(21)
Generate trial vectorsUit by (3);(22)
END FOR(23)
Step 2.3: Selection and Archive Update(24)
FORi=1 to NP(25)
IFfUit<fXit(26)
qit+1=0;(27)
UpdateA according to Algorithm 1 with Uit;(28)
ELSE(29)
qit+1=qit+1;(30)
END IF(31)
Generate new vectorsXit+1 by (4);(32)
END FOR(33)
Step 2.4: Increase the Count(34)
FEs=FEs+NP;(35)
t=t+1(36)
END WHILE


## 5. Experiments and Comparisons

In this section, the proposed ESM is first integrated into three representative classic DE algorithms, includingDE/rand/1, DE/best/1, and DE/current−to−best/1, to verify its efficiency. Then, five state-of-the-art DE variants, including j2020 [37], EJADE [38], BeSD [39], EBSHADE [14], and EBLSHADE [14], are introduced to incorporate with the ESM for further comparison. Finally, all experiments are conducted on the CEC2017 benchmark set [8], which contains 30 representative benchmark functions in four categories: unimodal (F1–F3), simple multimodal (F4–F10), hybrid (F11–F20), and composition (F21–F30). This benchmark set is widely used in the performance testing of evolutionary algorithms [9]. It is worth mentioning that the ESM does not change existing mechanisms in algorithms including the bound constraint handling and linear population size reduction.The parameters of all algorithms are listed in Table1. The common parameters of the incorporated algorithms are set as follows: The population size NP is 100, the stagnation tolerance Q is 120, and the maximum allowed number of function evaluations Max_FEs=10000∗D, where D is the problem dimension. Moreover, H represents the historical memory size associated with the adaptive control of scaling factors F and crossover rates CR, LP is the learning period.Table 1
Parameter settings of all compared algorithms.
No.ESMA=NP, Q=1201DE/rand/1F=0.5, CR=0.9;2DE/best/1F=0.5, CR=0.5;3DE/current−to−best/1F=0.5, CR=0.5;4j2020 [37]F=0.5, CR=0.9;5EJADE [38]uF=0.5, uCR=0.5, c=0.1, p=0.05;6BeSD [39]K=0.5;7EBSHADE [14]MF=0.5, MCR=0.5, and H=100;8EBLSHADE [14]MF=0.5, MCR=0.5, and H=100;Limited by the article space, this study presents results atD=30. For the fairness of comparison, each algorithm was implemented in MATLAB 2017b and executed over 51 independent runs on a Windows 10 64 bit desktop PC with 32 GB of RAM and a 3.0 GHz Intel Core i7-9700 processor.


### 5.1. Improving the Performance of Classic DEs and State-of-the-Art DE Variants

The following tables show the mean value (Mean) and standard deviation (Std.Dev) of the error valuefX−fX∗, where fX∗ representing the global optimum. The “+,” “−,” and “=” signs at each row of the table indicates that the ESM-based algorithm is, respectively, better than, equal to, or worse than the comparison algorithm.Tables2–4 show the performance of three classic DEs. For DE/rand/1, which is widely recognized for its diversity maintaining capability and relative low search efficiency, the ESM significantly improves its performance on 24 test functions and only cause performance degradation on two test functions, as shown in Table 2. DE/best/1 is well known as its local exploitation capability and fast convergence speed, the ESM also makes progress on 19 functions, as presented in Table 3. Relatively, the performance of ESM-DE/best/1 decreases on five functions including F13, F16, F24, F27, and F29 and keeps the same as DE/best/1 on the rest six functions. DE/current−to−best/1 is more inclined to balance global exploration and local exploitation. From Table 3, we can find that the ESM also achieves improvement on 26 functions. The results of ESM-DE/current−to−best/1 only decrease on two functions including F2 and F19. With functions F25 and F27, the two competitors keep the same.Table 2
Comparison results ofDE/rand/1 on CEC2017 test suite D=30.
DE/rand/1ESM-DE/rand/1MEanStd.DEvMEanStd.DEvF15.8903E − 012.9710E − 014.5086E − 011.7552E − 01+F25.7444E + 239.9418E + 239.2675E + 212.5025E + 22+F33.4902E + 046.1469E + 032.3752E + 043.9591E + 03+F46.5749E + 019.5494E + 006.7267E + 011.0601E + 01−F51.7859E + 029.8332E + 001.7115E + 021.1998E + 01+F61.3560E − 061.8981E − 066.2629E − 074.5960E − 07=F72.1873E + 021.0943E + 012.1104E + 028.6729E + 00+F81.8310E + 029.0762E + 001.7903E + 027.5831E + 00+F90.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F106.8115E + 032.7719E + 026.7908E + 032.7506E + 02+F117.7115E + 011.8139E + 017.6014E + 012.7738E + 01+F121.9220E + 061.3261E + 061.3328E + 051.6493E + 05+F136.9480E + 021.3966E + 024.9388E + 028.5027E + 01+F147.6534E + 016.7506E + 007.3477E + 016.1189E + 00+F157.2091E + 016.3029E + 006.0288E + 017.8284E + 00+F161.1650E + 031.6343E + 029.3936E + 021.6736E + 02+F172.8488E + 026.0664E + 017.4398E + 018.8628E + 00+F181.3888E + 045.2100E + 038.0349E + 032.2066E + 03+F193.6625E + 013.4612E + 003.6208E + 013.0908E + 00+F201.9409E + 021.4131E + 023.4660E + 016.7089E + 00+F213.7093E + 029.2962E + 003.6776E + 021.0301E + 01+F221.0000E + 021.5498E − 131.0000E + 020.0000E + 00=F235.2446E + 029.2278E + 005.1863E + 028.4044E + 00+F245.8819E + 021.2276E + 015.8615E + 028.0042E + 00+F253.8675E + 022.5661E − 023.8674E + 021.1873E − 02=F262.6453E + 039.4563E + 012.5338E + 031.1013E + 02+F274.9203E + 021.0543E + 014.9375E + 028.0730E + 00−F283.2756E + 024.6815E + 013.1464E + 023.5907E + 01+F299.0769E + 029.6800E + 017.3123E + 027.3738E + 01+F301.5880E + 043.3864E + 031.0375E + 041.4700E + 03+Total numbEr of (+/=/−): 24/4/2Table 3
Comparison results ofDE/best/1 on CEC2017 test suite D=30.
DE/best/1ESM-DE/best/1MeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F22.8258E + 101.1400E + 113.5454E + 041.6884E + 05+F36.0493E + 023.5493E + 023.1110E + 004.8214E + 00+F49.1949E + 012.7991E + 019.2586E + 012.5015E + 01=F59.2763E + 015.1334E + 014.4342E + 011.6997E + 01+F62.5865E − 026.1852E − 023.6252E − 021.1222E − 01=F71.7540E + 023.2845E + 018.2894E + 012.9691E + 01+F89.1332E + 015.5536E + 014.5522E + 011.1158E + 01+F91.3824E + 012.8355E + 011.3216E + 012.0467E + 01=F106.1379E + 035.6240E + 022.2535E + 031.6207E + 03+F116.9152E + 014.1818E + 015.9835E + 013.3934E + 01+F122.4650E + 041.7178E + 041.8722E + 041.1910E + 04+F138.0829E + 031.3281E + 041.1159E + 041.4904E + 04−F141.0157E + 022.4475E + 017.9954E + 013.5773E + 01+F151.3412E + 027.9017E + 018.4741E + 015.4255E + 01+F164.6549E + 022.8862E + 024.8784E + 022.7012E + 02−F171.6363E + 029.3907E + 011.5320E + 028.2527E + 01+F187.5687E + 045.7404E + 044.6144E + 042.9812E + 04+F196.4986E + 014.8051E + 016.0751E + 015.8037E + 01+F201.6043E + 029.2239E + 011.6611E + 029.9222E + 01=F212.8238E + 024.5823E + 012.4758E + 021.1434E + 01+F222.8353E + 032.8108E + 031.2129E + 031.6013E + 03+F234.1270E + 022.4531E + 013.9811E + 021.4860E + 01+F244.7466E + 022.8842E + 014.7812E + 021.3574E + 01−F253.9123E + 021.1628E + 013.9129E + 021.1172E + 01=F261.6222E + 032.1578E + 021.4146E + 034.2851E + 02+F275.1391E + 021.5952E + 015.1741E + 021.2805E + 01−F284.3979E + 022.9759E + 014.2865E + 022.3995E + 01+F295.7692E + 021.1599E + 025.9473E + 021.0825E + 02−F306.3732E + 033.1823E + 034.2434E + 032.5536E + 03+Total number of (+/=/−): 19/6/5Table 4
Comparison results ofDE/current−to−best/1 on CEC2017 test suiteD=30.
DE/current−to−best/1ESM-DE/current−to−best/1MeanStd.DevMeanStd.DevF14.2529E + 031.3861E + 042.4194E + 034.8764E + 03+F21.4541E + 135.2565E + 137.0969E + 133.8743E + 14−F32.5651E + 021.7770E + 021.4208E + 002.2041E + 00+F41.1469E + 021.3856E + 011.0468E + 021.8355E + 01+F51.4448E + 021.2098E + 012.9207E + 016.2955E + 00+F67.9340E − 034.4121E − 022.0523E − 031.1284E − 02+F71.7631E + 021.2016E + 017.5749E + 013.9251E + 01+F81.4187E + 021.2879E + 012.7608E + 017.8113E + 00+F94.2194E − 014.5290E − 012.8916E − 013.5412E − 01+F106.1656E + 033.9958E + 024.6688E + 031.8249E + 03+F116.8524E + 013.1083E + 015.5608E + 012.9379E + 01+F122.2269E + 041.3508E + 041.9785E + 041.3868E + 04+F136.5569E + 033.5373E + 034.7370E + 037.3061E + 03+F141.2361E + 022.9917E + 018.7387E + 012.5164E + 01+F152.1775E + 028.1693E + 011.9212E + 028.3283E + 01+F167.4089E + 022.3808E + 024.1651E + 022.2185E + 02+F171.8434E + 025.3907E + 019.5323E + 015.6810E + 01+F185.8340E + 042.6259E + 045.2731E + 043.7186E + 04+F197.9602E + 012.9727E + 011.1554E + 025.1436E + 01−F202.2812E + 028.2020E + 011.2592E + 028.0240E + 01+F213.3232E + 021.0378E + 012.3302E + 021.8782E + 01+F222.9659E + 021.0908E + 031.0053E + 021.1111E + 00+F234.6069E + 022.1824E + 013.8629E + 021.4205E + 01+F245.4998E + 021.6167E + 014.6028E + 021.5954E + 01+F254.0471E + 021.7220E + 014.0175E + 021.7238E + 01=F261.9292E + 035.1027E + 021.3090E + 032.7680E + 02+F275.1686E + 021.1760E + 015.1961E + 021.4919E + 01=F284.7064E + 023.4376E + 014.5744E + 023.0111E + 01+F297.0816E + 021.0684E + 025.2205E + 028.3692E + 01+F309.8358E + 037.4280E + 037.5623E + 033.6244E + 03+Total number of (+/=/−): 26/2/2From the above experiments, we can conclude that the ESM can significantly improve classic DE’s performance. To further verify the ESM’s universality, five representative and advanced DE variants are introduced, Tables5–9 show these variants and their competitors’ performance. For ESM-j2020, as shown in Table 5, significant improvements are achieved on 23 functions. From Tables 6 to 9, we know that the number of functions with “+” decreases compared with the results shown in Tables 2–5. Meanwhile, from Tables 6 to 9, only a few functions show performance degradation after applying the ESM. We can therefore conclude that the ESM is still competitive even in recognized excellent DE variants.Table 5
Comparison results of j2020 on CEC2017 test suiteD=30.
j2020ESM-j2020MeanStd.DevMeanStd.DevF13.9138E − 022.7831E − 024.1149E − 022.6801E − 02−F22.1525E + 227.2306E + 223.5391E + 167.4152E + 16+F33.0526E + 045.6182E + 031.5964E + 043.6765E + 03+F47.6190E + 011.0125E + 017.7537E + 011.0041E + 01−F51.5059E + 021.5882E + 013.3617E + 011.1325E + 01+F60.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F71.9215E + 021.0204E + 018.8048E + 012.4532E + 01+F81.5771E + 029.7337E + 003.5822E + 011.3602E + 01+F90.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F106.0724E + 032.8112E + 022.6253E + 039.8373E + 02+F117.1620E + 012.5123E + 011.4904E + 011.6548E + 01+F127.3213E + 056.1498E + 054.9652E + 043.2587E + 04+F133.9600E + 028.0250E + 012.6856E + 024.2577E + 01+F147.1081E + 015.4527E + 005.2870E + 011.1383E + 01+F156.0349E + 015.8258E + 003.6537E + 011.3383E + 01+F169.0601E + 021.4229E + 024.1188E + 021.7961E + 02+F179.1833E + 013.4540E + 014.3249E + 012.4777E + 01+F183.2864E + 031.5646E + 031.5206E + 038.9776E + 02+F193.3522E + 014.8067E + 001.5961E + 012.9754E + 00+F205.3644E + 012.1292E + 014.7332E + 016.8896E + 01+F213.5434E + 026.7332E + 002.3149E + 021.2889E + 01+F221.0000E + 022.2819E − 131.0000E + 028.1676E − 14=F234.9590E + 026.7701E + 003.7594E + 021.0453E + 01+F245.7509E + 028.0847E + 004.5326E + 021.1284E + 01+F253.8677E + 022.5429E − 023.8676E + 021.3848E − 02=F262.3688E + 031.1160E + 021.1769E + 031.2856E + 02+F274.9473E + 021.1847E + 014.9543E + 026.3871E + 00−F283.3192E + 024.5936E + 013.2368E + 024.8265E + 01+F297.8452E + 028.7484E + 014.8871E + 022.6830E + 01+F309.8924E + 032.2550E + 036.4595E + 031.2161E + 03+Total numbEr of (+/=/−): 23/4/3Table 6
Comparison results of EJADE on CEC2017 test suiteD=30.
EJADEESM-EJADEMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F24.5316E + 032.1794E + 041.2673E + 045.6411E + 04−F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F42.4107E + 012.9037E + 012.6740E + 012.9430E + 01−F51.8962E + 014.4630E + 001.7560E + 014.1722E + 00+F66.9845E − 072.4275E − 061.4467E − 066.4530E − 06−F74.6647E + 012.8902E + 004.5802E + 013.2092E + 00+F81.8850E + 013.1904E + 001.7218E + 015.1087E + 00+F91.9587E − 013.4505E − 011.2548E − 011.8845E − 01+F101.8437E + 033.0553E + 021.5558E + 035.3726E + 02+F114.7182E + 012.6517E + 013.9138E + 012.3463E + 01+F122.0542E + 031.8865E + 031.8759E + 031.2938E + 03+F134.9648E + 013.1191E + 013.7663E + 011.7336E + 01+F143.9862E + 011.3835E + 013.8508E + 011.1120E + 01+F154.1678E + 014.0043E + 014.2540E + 013.6838E + 01−F163.2052E + 021.8067E + 022.6554E + 021.5914E + 02+F175.7500E + 011.3465E + 015.1900E + 013.5779E + 01+F187.7806E + 015.1375E + 018.8347E + 016.3700E + 01−F191.7808E + 011.0382E + 011.9606E + 011.8146E + 01−F201.0632E + 025.7664E + 015.1910E + 015.3314E + 01+F212.1850E + 023.6330E + 002.1748E + 024.0920E + 00+F221.0008E + 024.4179E − 011.6606E + 023.6781E + 02−F233.6253E + 026.0195E + 003.6162E + 026.1887E + 00+F244.3448E + 026.1456E + 004.3460E + 025.0372E + 00=F253.8703E + 027.8427E − 013.8701E + 023.6023E − 01=F261.1198E + 037.2825E + 011.0310E + 032.7881E + 02+F275.0690E + 028.8147E + 005.0681E + 026.6749E + 00=F283.7620E + 026.0263E + 013.3258E + 025.2414E + 01+F294.6359E + 022.7079E + 014.4343E + 023.0038E + 01+F302.1586E + 031.9675E + 022.1400E + 032.0292E + 02+Total number of (+/ = /−): 18/5/7Table 7
Comparison results of BeSD on CEC2017 test suiteD=30.
BeSDESM- BeSDMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F27.2258E + 001.9085E + 013.9032E + 001.4996E + 01+F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F46.8185E + 001.6230E + 011.7131E + 012.6599E + 01−F54.5762E + 012.8543E + 012.7217E + 018.1867E + 00+F60.0000E + 000.0000E + 001.7219E − 073.7418E − 07=F71.3363E + 021.5180E + 015.1262E + 016.7117E + 00+F85.6343E + 012.9912E + 012.6864E + 017.8132E + 00+F96.1295E − 021.9377E − 016.7718E − 021.7464E − 01=F105.8972E + 033.7315E + 022.5464E + 035.6311E + 02+F112.6612E + 011.5396E + 013.7557E + 012.2784E + 01−F128.3133E + 034.3141E + 038.9538E + 037.0576E + 03−F132.2656E + 022.3629E + 026.3469E + 012.8828E + 01+F146.1777E + 011.5857E + 014.1127E + 019.6938E + 00+F156.2060E + 013.5382E + 013.8429E + 012.6456E + 01+F166.1366E + 021.7387E + 022.2941E + 021.8203E + 02+F171.0414E + 022.0838E + 013.4358E + 011.1387E + 01+F182.1008E + 029.9628E + 017.7754E + 014.7271E + 01+F194.0368E + 012.1662E + 013.0129E + 011.4941E + 01+F201.1522E + 024.3894E + 014.3747E + 015.1852E + 01+F212.3501E + 022.5203E + 012.2421E + 029.6408E + 00+F221.0000E + 021.1357E − 131.0000E + 022.2124E − 13=F233.7384E + 021.6564E + 013.7369E + 029.3329E + 00=F244.4032E + 027.1928E + 004.3973E + 028.2282E + 00+F253.8708E + 022.9529E − 013.8709E + 022.7823E − 01=F261.1447E + 031.7473E + 021.0790E + 032.7897E + 02+F275.0192E + 025.4148E + 005.0552E + 028.0174E + 00−F283.1676E + 023.8969E + 013.3033E + 024.8241E + 01−F295.9423E + 027.0224E + 014.4240E + 023.4337E + 01+F302.5624E + 033.5089E + 022.1699E + 031.3768E + 02+Total number of (+/ = /−): 18/7/5Table 8
Comparison results of EBSHADE on CEC2017 test suiteD=30.
EBSHADEESM-EBSHADEMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F21.9972E + 041.1037E + 051.6852E + 027.3798E + 02+F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F45.6852E + 011.0598E + 015.3253E + 011.7773E + 01+F51.5508E + 012.5802E + 001.6273E + 016.3554E + 00−F66.8453E − 081.0785E-074.2769E−077.4531E−07−F74.6604E + 012.7607E + 004.6489E + 014.8631E + 00+F81.6879E + 012.8304E + 001.6186E + 013.4825E + 00+F91.7544E − 028.2641E−021.4656E − 028.1599E−02+F101.6353E + 032.3269E + 021.4154E + 032.6504E + 02+F112.7724E + 012.6071E + 012.6265E + 012.5887E + 01+F121.2219E + 034.6754E + 021.2126E + 034.4003E + 02+F133.8215E + 012.0185E + 013.7691E + 011.4370E + 01+F142.8375E + 018.7055E + 003.0017E + 014.8689E + 00−F151.7726E + 011.0888E + 012.1258E + 011.2917E + 01−F162.9276E + 021.1820E + 022.0707E + 021.2958E + 02+F174.5658E + 011.0476E + 013.5719E + 011.0310E + 01=F189.0145E + 014.7859E + 019.4472E + 015.6305E + 01−F191.4666E + 018.9947E + 001.2179E + 014.1749E + 00+F206.3748E + 014.4406E + 015.2312E + 014.9270E + 01+F212.1777E + 023.4870E + 002.1847E + 024.1234E + 00−F221.0051E + 022.8400E + 001.0000E + 022.2124E−13=F233.6665E + 025.5240E + 003.6453E + 025.2131E + 00+F244.3676E + 023.7374E + 004.3642E + 026.0049E + 00=F253.8681E + 025.8187E-023.8684E + 027.5494E−02=F261.1214E + 036.7789E + 011.1160E + 037.3256E + 01+F275.0632E + 024.6020E + 005.0579E + 026.3243E + 00+F283.4497E + 025.8589E + 013.4205E + 025.8955E + 01+F294.6900E + 021.3609E + 014.4219E + 022.2017E + 01+F302.0864E + 031.3828E + 022.0213E + 031.6473E + 02+Total number of (+/ = /−): 18/6/6Table 9
Comparison results of EBLSHADE on CEC2017 test suiteD=30.
EBLSHADEESM-EBLSHADEMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F20.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F45.8562E + 013.4396E − 145.8562E + 014.0776E − 14=F56.6978E + 001.3325E + 006.9719E + 001.5418E + 00=F60.0000E + 000.0000E + 001.4398E − 084.1347E − 08=F73.7938E + 011.3728E + 003.6957E + 011.2423E + 00+F87.4627E + 001.5388E + 006.9654E + 001.6395E + 00+F90.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F101.3903E + 032.5407E + 021.2476E + 032.3487E + 02+F113.3748E + 012.8992E + 014.0898E + 012.8428E + 01−F121.1107E + 033.2096E + 021.1021E + 033.2418E + 02+F131.8932E + 017.7131E + 001.6355E + 015.4803E + 00+F142.1636E + 011.1516E + 002.1493E + 013.6187E + 00+F153.9785E + 001.8483E + 003.6417E + 001.6998E + 00+F164.1583E + 013.4857E + 012.0657E + 015.3775E + 00+F173.1880E + 015.6772E + 003.1060E + 014.8829E + 00+F182.2262E + 011.5977E + 002.2496E + 011.2879E + 00−F196.2724E + 002.2466E + 005.6099E + 001.3474E + 00+F203.2672E + 016.5190E + 002.7681E + 016.2855E + 00+F212.0835E + 021.4350E + 002.0819E + 021.4361E + 00+F221.0000E + 020.0000E + 001.0000E + 028.1676E − 14=F233.5311E + 023.0584E + 003.5407E + 022.9996E + 00−F244.2796E + 021.7733E + 004.2807E + 021.8546E + 00−F253.8672E + 021.9899E − 023.8672E + 022.4427E − 02=F269.7950E + 024.0916E + 019.7385E + 022.9174E + 01+F275.0689E + 025.5012E + 005.0580E + 024.6504E + 00+F283.5816E + 026.2765E + 013.2838E + 024.8986E + 01+F294.3621E + 027.9105E + 004.3392E + 026.5290E + 00+F301.9967E + 037.3147E + 012.0385E + 038.6201E + 01−Total number of (+/ = /−): 16/9/5We further conducted Wilcoxon signed-ranked test with a 0.05 significance [40] to statistically verify the effectiveness of the ESM. In Table 10, the p-value obtained by ESM-DE/rand/1, ESM-DE/best/1, ESM-DE/current−to−best/1, ESM-j2020, and ESM-BeSD are all less than 0.05, which indicates a significant improvement is achieved by the proposed ESM. It also can be found that ESM-EJADE, ESM-EBSHADE, and ESM-EBLSHADE still maintain their advantages in terms of the value of R+ and R–, which also reflects that the ESM achieves more promising results.Table 10
Results obtained by the Wilcoxon test.
AlgorithmsR+R−P-valueESM-DE/rand/1 vs. DE/rand/1446.518.50.039ESM-DE/best/1 vs. DE/best/1376.588.50.047ESM-DE/current−to−best/1 vs. DE/current−to−best/1423420.031ESM-j2020 vs. j2020419.545.50.038ESM-EJADE vs. EJADE326.5138.50.081ESM-BeSD vs. BeSD338.5126.50.046ESM-EBSHADE vs. EBSHADE375.589.50.075ESM-EBLSHADE vs. EBLSHADE341.5123.50.089As mentioned before, the functions of CEC2017 can be divided into four categories: unimodal functions, simple multimodal functions, hybrid functions, and composition functions. To analyze whether the performance of the ESM is related to the type of function, the performance of the ESM-based algorithms for each type of function category are presented in Table11. As shown in Table 11, the ESM performs best on hybrid functions. Following the adoption of the ESM, the algorithm has an 81.25% probability of finding a better solution but only a 16.25% probability of finding a worse solution on hybrid functions. Besides, the performance of the ESM is almost the same on simple multimodal functions and composition functions and relatively worst on unimodal functions. Table 11 demonstrates that if the optimization problem is multimodal function, hybrid function, or composition function, then the proposed selection mechanism can be considered in the DE algorithm.Table 11
The classification statistics result of proposed ESM in eight algorithms.
Functions+(%)=(%)−(%)Unimodal functions45.8341.6712.50Simple multimodal functions64.2923.2112.50Hybrid functions81.252.5016.25Composition functions62.5022.5015.00


### 5.2. Convergence Analysis

To more clearly show the influence of the ESM on algorithm convergence and avoid diversion, four representative functions from CEC2017 are extracted, including F7, F16, F20, and F24, convergence curves for these functions with the three classic DEs and their relative ESM-based versions are presented. F7 is a simple multimodal function with a characteristic of convexity. From Figure5(a), we can find that all six comparison algorithms have a similar convergence rate; this also reflects that the ESM does not cripple the convergence speed of the algorithm on F7. From Figures 5(b) to 5(d), it can be seen that when each algorithm adopts the proposed ESM, it converges to a better value than the original algorithm with faster convergence. This is mainly because the algorithm can hardly escape when it suffers stagnation. The algorithm that uses the proposed ESM has a more excellent capability to keep search efficiency.Figure 5
Convergence curves of three classic DEs and relative ESM-DEs on four 30-D CEC2017 benchmark functions: (a) F7; (b) F16; (c) F20; (d) F24.
(a)(b)(c)(d)


### 5.3. Population Diversity Analysis

As the ESM adopts the entropy-based individual diversity measure, we further calculated the diversity of the main population and archive population on F7, F16, F20, and F24 for three classic DEs. It should be pointed out that we used the concept proposed in [41] to measure population diversity:(13)xj¯=1NP∑i=1NPxj,i,Population Diversity=1NP∑i=1NP∑j=1Dxj,i−xj¯2.As illustrated in Figure6(a), the main population and archive population have an almost uniform diversity on F7. From Figures 6(b)to6(d), it can be concluded that the archive population can maintain a relatively higher diversity than the main population at the intermediate stage for DE/current−to−best/1 and DE/best/1. To sum up, the archive population can maintain relative diversity during the search process, which makes it possible for the algorithm to recover from stagnation.Figure 6
Population diversity curves of main population and archive population of three classic DEs on four 30-D CEC2017 benchmark functions: (a) F7; (b) F16; (c) F20; (d) F24.
(a)(b)(c)(d)


### 5.4. Parameter Sensitivity Analysis

As illustrated in the previous sections, the ESM contains one parameterQ required adjustment. Q is a threshold which represents the stagnation steps. Table 12 shows the B-W values that represent the variation between better (“+”) and worse (“−”) number of the mean value after applying the ESM. It can be seen that the ESM-based algorithm performance is not so sensitive to the change of Q values over a finite range. For further comparison, A wider range of Q values and its corresponding B-W values are shown in Figure 7. In general, when the Q value is too large, the overall performance of the ESM begins to deteriorate. This is mainly because too large Q value reduces the probability of the ESM intervening in the algorithm; therefore, it is difficult for the ESM to function in limited iterations. In summary, considering the overall performance of the ESM, Q is set to 120.Table 12
The B-W values of all ESM-based algorithms on CEC2017 under differentQ values.
AlgorithmsQ=100Q=110Q=120Q=130Q=140ESM-DE/rand/12221222320ESM-DE/best/11314141314ESM-DE/current−to−best/12424242220ESM-j20201821201917ESM-EJADE1012111010ESM-BeSD1110131212ESM-EBSHADE1011121011ESM-EBLSHADE910111010Total117126127119114Figure 7
The B-W values at different Q values.


### 5.5. Scalability Analysis

The dimension of the test functions governs the difficulties on finding the global optimum. Higher dimensional functions are generally more difficult to solve. To verify the relationship between the dimensionality and the performance of the proposed mechanism, we evaluate the average performance of eight ESM-based DE algorithms in CEC2017 benchmark set. Table13 shows the results (+/=/−) of the considered algorithms at D=10, 30, 50, and 100 on the test functions, for intuition, we convert the corresponding B-W values into Figure 8. We can find that the performance of the ESM-based algorithm fluctuates slightly but does not degrade significantly with the increase of problem dimension compared to the original algorithm. In Table 14, the Friedman test [40], a widely used nonparametric test in the EA community, is used to validate the performance of all algorithms based on the mean value. It is not difficult to see that the p-values from D=10 to D=100 calculated by the Friedman test are all less than 0.05. Therefore, there are significant differences in the performance of the comparison algorithms on the corresponding dimensions, and ESM-EBLSHADE gets the first rank overall.Table 13
The statistical results of all ESM-based algorithms on CEC2017 under different dimensions.
AlgorithmsD=10D=30D=50D=100ESM-DE/rand/1 vs. DE/rand/123/5/224/4/222/6/225/4/1ESM-DE/best/1 vs. DE/best/117/8/519/6/518/6/619/5/6ESM-DE/current−to−best/1 vs. DE/current−to−best/123/4/326/2/225/2/323/5/2ESM-j2020 vs. j202022/5/323/4/321/5/422/4/4ESM-EJADE vs. EJADE16/8/618/5/720/5/519/5/6ESM-BeSD vs. BeSD17/7/618/7/519/6/518/5/7ESM-EBSHADE vs. EBSHADE16/9/518/6/618/5/717/6/7ESM-EBLSHADE vs. EBLSHADE15/10/516/9/515/11/416/8/6Figure 8
The curves of B-W values for the ESM-based algorithm compared to the original algorithm.Table 14
The Friedman test results of all algorithms on CEC2017 under different dimensions.
AlgorithmsD=10D=30D=50D=100Mean rankingRankESM-EBLSHADE2.132.622.882.5710.201EBLSHADE2.983.033.213.5512.772ESM-EBSHADE4.854.785.014.8919.533EBSHADE4.985.735.555.2121.474ESM-EJADE5.896.076.245.8724.075ESM-BeSD5.316.736.846.5725.456EJADE6.516.936.956.8827.277ESM-j20207.517.627.587.4130.128BeSD7.117.938.028.5731.639j20209.8210.899.5811.2141.5010ESM-DE/rand/110.0911.0311.5911.0743.7811ESM-DE/current−to−best/110.5511.6311.2111.8545.2412ESM-DE/best/110.1111.8211.8511.9845.7613DE/best/111.5312.6512.3312.5749.0814DE/rand/112.3112.7012.5812.3449.9315DE/current−to−best/111.8913.8312.8912.5951.2016Friedman-P-value0.000.000.000.00


## 5.1. Improving the Performance of Classic DEs and State-of-the-Art DE Variants

The following tables show the mean value (Mean) and standard deviation (Std.Dev) of the error valuefX−fX∗, where fX∗ representing the global optimum. The “+,” “−,” and “=” signs at each row of the table indicates that the ESM-based algorithm is, respectively, better than, equal to, or worse than the comparison algorithm.Tables2–4 show the performance of three classic DEs. For DE/rand/1, which is widely recognized for its diversity maintaining capability and relative low search efficiency, the ESM significantly improves its performance on 24 test functions and only cause performance degradation on two test functions, as shown in Table 2. DE/best/1 is well known as its local exploitation capability and fast convergence speed, the ESM also makes progress on 19 functions, as presented in Table 3. Relatively, the performance of ESM-DE/best/1 decreases on five functions including F13, F16, F24, F27, and F29 and keeps the same as DE/best/1 on the rest six functions. DE/current−to−best/1 is more inclined to balance global exploration and local exploitation. From Table 3, we can find that the ESM also achieves improvement on 26 functions. The results of ESM-DE/current−to−best/1 only decrease on two functions including F2 and F19. With functions F25 and F27, the two competitors keep the same.Table 2
Comparison results ofDE/rand/1 on CEC2017 test suite D=30.
DE/rand/1ESM-DE/rand/1MEanStd.DEvMEanStd.DEvF15.8903E − 012.9710E − 014.5086E − 011.7552E − 01+F25.7444E + 239.9418E + 239.2675E + 212.5025E + 22+F33.4902E + 046.1469E + 032.3752E + 043.9591E + 03+F46.5749E + 019.5494E + 006.7267E + 011.0601E + 01−F51.7859E + 029.8332E + 001.7115E + 021.1998E + 01+F61.3560E − 061.8981E − 066.2629E − 074.5960E − 07=F72.1873E + 021.0943E + 012.1104E + 028.6729E + 00+F81.8310E + 029.0762E + 001.7903E + 027.5831E + 00+F90.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F106.8115E + 032.7719E + 026.7908E + 032.7506E + 02+F117.7115E + 011.8139E + 017.6014E + 012.7738E + 01+F121.9220E + 061.3261E + 061.3328E + 051.6493E + 05+F136.9480E + 021.3966E + 024.9388E + 028.5027E + 01+F147.6534E + 016.7506E + 007.3477E + 016.1189E + 00+F157.2091E + 016.3029E + 006.0288E + 017.8284E + 00+F161.1650E + 031.6343E + 029.3936E + 021.6736E + 02+F172.8488E + 026.0664E + 017.4398E + 018.8628E + 00+F181.3888E + 045.2100E + 038.0349E + 032.2066E + 03+F193.6625E + 013.4612E + 003.6208E + 013.0908E + 00+F201.9409E + 021.4131E + 023.4660E + 016.7089E + 00+F213.7093E + 029.2962E + 003.6776E + 021.0301E + 01+F221.0000E + 021.5498E − 131.0000E + 020.0000E + 00=F235.2446E + 029.2278E + 005.1863E + 028.4044E + 00+F245.8819E + 021.2276E + 015.8615E + 028.0042E + 00+F253.8675E + 022.5661E − 023.8674E + 021.1873E − 02=F262.6453E + 039.4563E + 012.5338E + 031.1013E + 02+F274.9203E + 021.0543E + 014.9375E + 028.0730E + 00−F283.2756E + 024.6815E + 013.1464E + 023.5907E + 01+F299.0769E + 029.6800E + 017.3123E + 027.3738E + 01+F301.5880E + 043.3864E + 031.0375E + 041.4700E + 03+Total numbEr of (+/=/−): 24/4/2Table 3
Comparison results ofDE/best/1 on CEC2017 test suite D=30.
DE/best/1ESM-DE/best/1MeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F22.8258E + 101.1400E + 113.5454E + 041.6884E + 05+F36.0493E + 023.5493E + 023.1110E + 004.8214E + 00+F49.1949E + 012.7991E + 019.2586E + 012.5015E + 01=F59.2763E + 015.1334E + 014.4342E + 011.6997E + 01+F62.5865E − 026.1852E − 023.6252E − 021.1222E − 01=F71.7540E + 023.2845E + 018.2894E + 012.9691E + 01+F89.1332E + 015.5536E + 014.5522E + 011.1158E + 01+F91.3824E + 012.8355E + 011.3216E + 012.0467E + 01=F106.1379E + 035.6240E + 022.2535E + 031.6207E + 03+F116.9152E + 014.1818E + 015.9835E + 013.3934E + 01+F122.4650E + 041.7178E + 041.8722E + 041.1910E + 04+F138.0829E + 031.3281E + 041.1159E + 041.4904E + 04−F141.0157E + 022.4475E + 017.9954E + 013.5773E + 01+F151.3412E + 027.9017E + 018.4741E + 015.4255E + 01+F164.6549E + 022.8862E + 024.8784E + 022.7012E + 02−F171.6363E + 029.3907E + 011.5320E + 028.2527E + 01+F187.5687E + 045.7404E + 044.6144E + 042.9812E + 04+F196.4986E + 014.8051E + 016.0751E + 015.8037E + 01+F201.6043E + 029.2239E + 011.6611E + 029.9222E + 01=F212.8238E + 024.5823E + 012.4758E + 021.1434E + 01+F222.8353E + 032.8108E + 031.2129E + 031.6013E + 03+F234.1270E + 022.4531E + 013.9811E + 021.4860E + 01+F244.7466E + 022.8842E + 014.7812E + 021.3574E + 01−F253.9123E + 021.1628E + 013.9129E + 021.1172E + 01=F261.6222E + 032.1578E + 021.4146E + 034.2851E + 02+F275.1391E + 021.5952E + 015.1741E + 021.2805E + 01−F284.3979E + 022.9759E + 014.2865E + 022.3995E + 01+F295.7692E + 021.1599E + 025.9473E + 021.0825E + 02−F306.3732E + 033.1823E + 034.2434E + 032.5536E + 03+Total number of (+/=/−): 19/6/5Table 4
Comparison results ofDE/current−to−best/1 on CEC2017 test suiteD=30.
DE/current−to−best/1ESM-DE/current−to−best/1MeanStd.DevMeanStd.DevF14.2529E + 031.3861E + 042.4194E + 034.8764E + 03+F21.4541E + 135.2565E + 137.0969E + 133.8743E + 14−F32.5651E + 021.7770E + 021.4208E + 002.2041E + 00+F41.1469E + 021.3856E + 011.0468E + 021.8355E + 01+F51.4448E + 021.2098E + 012.9207E + 016.2955E + 00+F67.9340E − 034.4121E − 022.0523E − 031.1284E − 02+F71.7631E + 021.2016E + 017.5749E + 013.9251E + 01+F81.4187E + 021.2879E + 012.7608E + 017.8113E + 00+F94.2194E − 014.5290E − 012.8916E − 013.5412E − 01+F106.1656E + 033.9958E + 024.6688E + 031.8249E + 03+F116.8524E + 013.1083E + 015.5608E + 012.9379E + 01+F122.2269E + 041.3508E + 041.9785E + 041.3868E + 04+F136.5569E + 033.5373E + 034.7370E + 037.3061E + 03+F141.2361E + 022.9917E + 018.7387E + 012.5164E + 01+F152.1775E + 028.1693E + 011.9212E + 028.3283E + 01+F167.4089E + 022.3808E + 024.1651E + 022.2185E + 02+F171.8434E + 025.3907E + 019.5323E + 015.6810E + 01+F185.8340E + 042.6259E + 045.2731E + 043.7186E + 04+F197.9602E + 012.9727E + 011.1554E + 025.1436E + 01−F202.2812E + 028.2020E + 011.2592E + 028.0240E + 01+F213.3232E + 021.0378E + 012.3302E + 021.8782E + 01+F222.9659E + 021.0908E + 031.0053E + 021.1111E + 00+F234.6069E + 022.1824E + 013.8629E + 021.4205E + 01+F245.4998E + 021.6167E + 014.6028E + 021.5954E + 01+F254.0471E + 021.7220E + 014.0175E + 021.7238E + 01=F261.9292E + 035.1027E + 021.3090E + 032.7680E + 02+F275.1686E + 021.1760E + 015.1961E + 021.4919E + 01=F284.7064E + 023.4376E + 014.5744E + 023.0111E + 01+F297.0816E + 021.0684E + 025.2205E + 028.3692E + 01+F309.8358E + 037.4280E + 037.5623E + 033.6244E + 03+Total number of (+/=/−): 26/2/2From the above experiments, we can conclude that the ESM can significantly improve classic DE’s performance. To further verify the ESM’s universality, five representative and advanced DE variants are introduced, Tables5–9 show these variants and their competitors’ performance. For ESM-j2020, as shown in Table 5, significant improvements are achieved on 23 functions. From Tables 6 to 9, we know that the number of functions with “+” decreases compared with the results shown in Tables 2–5. Meanwhile, from Tables 6 to 9, only a few functions show performance degradation after applying the ESM. We can therefore conclude that the ESM is still competitive even in recognized excellent DE variants.Table 5
Comparison results of j2020 on CEC2017 test suiteD=30.
j2020ESM-j2020MeanStd.DevMeanStd.DevF13.9138E − 022.7831E − 024.1149E − 022.6801E − 02−F22.1525E + 227.2306E + 223.5391E + 167.4152E + 16+F33.0526E + 045.6182E + 031.5964E + 043.6765E + 03+F47.6190E + 011.0125E + 017.7537E + 011.0041E + 01−F51.5059E + 021.5882E + 013.3617E + 011.1325E + 01+F60.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F71.9215E + 021.0204E + 018.8048E + 012.4532E + 01+F81.5771E + 029.7337E + 003.5822E + 011.3602E + 01+F90.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F106.0724E + 032.8112E + 022.6253E + 039.8373E + 02+F117.1620E + 012.5123E + 011.4904E + 011.6548E + 01+F127.3213E + 056.1498E + 054.9652E + 043.2587E + 04+F133.9600E + 028.0250E + 012.6856E + 024.2577E + 01+F147.1081E + 015.4527E + 005.2870E + 011.1383E + 01+F156.0349E + 015.8258E + 003.6537E + 011.3383E + 01+F169.0601E + 021.4229E + 024.1188E + 021.7961E + 02+F179.1833E + 013.4540E + 014.3249E + 012.4777E + 01+F183.2864E + 031.5646E + 031.5206E + 038.9776E + 02+F193.3522E + 014.8067E + 001.5961E + 012.9754E + 00+F205.3644E + 012.1292E + 014.7332E + 016.8896E + 01+F213.5434E + 026.7332E + 002.3149E + 021.2889E + 01+F221.0000E + 022.2819E − 131.0000E + 028.1676E − 14=F234.9590E + 026.7701E + 003.7594E + 021.0453E + 01+F245.7509E + 028.0847E + 004.5326E + 021.1284E + 01+F253.8677E + 022.5429E − 023.8676E + 021.3848E − 02=F262.3688E + 031.1160E + 021.1769E + 031.2856E + 02+F274.9473E + 021.1847E + 014.9543E + 026.3871E + 00−F283.3192E + 024.5936E + 013.2368E + 024.8265E + 01+F297.8452E + 028.7484E + 014.8871E + 022.6830E + 01+F309.8924E + 032.2550E + 036.4595E + 031.2161E + 03+Total numbEr of (+/=/−): 23/4/3Table 6
Comparison results of EJADE on CEC2017 test suiteD=30.
EJADEESM-EJADEMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F24.5316E + 032.1794E + 041.2673E + 045.6411E + 04−F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F42.4107E + 012.9037E + 012.6740E + 012.9430E + 01−F51.8962E + 014.4630E + 001.7560E + 014.1722E + 00+F66.9845E − 072.4275E − 061.4467E − 066.4530E − 06−F74.6647E + 012.8902E + 004.5802E + 013.2092E + 00+F81.8850E + 013.1904E + 001.7218E + 015.1087E + 00+F91.9587E − 013.4505E − 011.2548E − 011.8845E − 01+F101.8437E + 033.0553E + 021.5558E + 035.3726E + 02+F114.7182E + 012.6517E + 013.9138E + 012.3463E + 01+F122.0542E + 031.8865E + 031.8759E + 031.2938E + 03+F134.9648E + 013.1191E + 013.7663E + 011.7336E + 01+F143.9862E + 011.3835E + 013.8508E + 011.1120E + 01+F154.1678E + 014.0043E + 014.2540E + 013.6838E + 01−F163.2052E + 021.8067E + 022.6554E + 021.5914E + 02+F175.7500E + 011.3465E + 015.1900E + 013.5779E + 01+F187.7806E + 015.1375E + 018.8347E + 016.3700E + 01−F191.7808E + 011.0382E + 011.9606E + 011.8146E + 01−F201.0632E + 025.7664E + 015.1910E + 015.3314E + 01+F212.1850E + 023.6330E + 002.1748E + 024.0920E + 00+F221.0008E + 024.4179E − 011.6606E + 023.6781E + 02−F233.6253E + 026.0195E + 003.6162E + 026.1887E + 00+F244.3448E + 026.1456E + 004.3460E + 025.0372E + 00=F253.8703E + 027.8427E − 013.8701E + 023.6023E − 01=F261.1198E + 037.2825E + 011.0310E + 032.7881E + 02+F275.0690E + 028.8147E + 005.0681E + 026.6749E + 00=F283.7620E + 026.0263E + 013.3258E + 025.2414E + 01+F294.6359E + 022.7079E + 014.4343E + 023.0038E + 01+F302.1586E + 031.9675E + 022.1400E + 032.0292E + 02+Total number of (+/ = /−): 18/5/7Table 7
Comparison results of BeSD on CEC2017 test suiteD=30.
BeSDESM- BeSDMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F27.2258E + 001.9085E + 013.9032E + 001.4996E + 01+F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F46.8185E + 001.6230E + 011.7131E + 012.6599E + 01−F54.5762E + 012.8543E + 012.7217E + 018.1867E + 00+F60.0000E + 000.0000E + 001.7219E − 073.7418E − 07=F71.3363E + 021.5180E + 015.1262E + 016.7117E + 00+F85.6343E + 012.9912E + 012.6864E + 017.8132E + 00+F96.1295E − 021.9377E − 016.7718E − 021.7464E − 01=F105.8972E + 033.7315E + 022.5464E + 035.6311E + 02+F112.6612E + 011.5396E + 013.7557E + 012.2784E + 01−F128.3133E + 034.3141E + 038.9538E + 037.0576E + 03−F132.2656E + 022.3629E + 026.3469E + 012.8828E + 01+F146.1777E + 011.5857E + 014.1127E + 019.6938E + 00+F156.2060E + 013.5382E + 013.8429E + 012.6456E + 01+F166.1366E + 021.7387E + 022.2941E + 021.8203E + 02+F171.0414E + 022.0838E + 013.4358E + 011.1387E + 01+F182.1008E + 029.9628E + 017.7754E + 014.7271E + 01+F194.0368E + 012.1662E + 013.0129E + 011.4941E + 01+F201.1522E + 024.3894E + 014.3747E + 015.1852E + 01+F212.3501E + 022.5203E + 012.2421E + 029.6408E + 00+F221.0000E + 021.1357E − 131.0000E + 022.2124E − 13=F233.7384E + 021.6564E + 013.7369E + 029.3329E + 00=F244.4032E + 027.1928E + 004.3973E + 028.2282E + 00+F253.8708E + 022.9529E − 013.8709E + 022.7823E − 01=F261.1447E + 031.7473E + 021.0790E + 032.7897E + 02+F275.0192E + 025.4148E + 005.0552E + 028.0174E + 00−F283.1676E + 023.8969E + 013.3033E + 024.8241E + 01−F295.9423E + 027.0224E + 014.4240E + 023.4337E + 01+F302.5624E + 033.5089E + 022.1699E + 031.3768E + 02+Total number of (+/ = /−): 18/7/5Table 8
Comparison results of EBSHADE on CEC2017 test suiteD=30.
EBSHADEESM-EBSHADEMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F21.9972E + 041.1037E + 051.6852E + 027.3798E + 02+F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F45.6852E + 011.0598E + 015.3253E + 011.7773E + 01+F51.5508E + 012.5802E + 001.6273E + 016.3554E + 00−F66.8453E − 081.0785E-074.2769E−077.4531E−07−F74.6604E + 012.7607E + 004.6489E + 014.8631E + 00+F81.6879E + 012.8304E + 001.6186E + 013.4825E + 00+F91.7544E − 028.2641E−021.4656E − 028.1599E−02+F101.6353E + 032.3269E + 021.4154E + 032.6504E + 02+F112.7724E + 012.6071E + 012.6265E + 012.5887E + 01+F121.2219E + 034.6754E + 021.2126E + 034.4003E + 02+F133.8215E + 012.0185E + 013.7691E + 011.4370E + 01+F142.8375E + 018.7055E + 003.0017E + 014.8689E + 00−F151.7726E + 011.0888E + 012.1258E + 011.2917E + 01−F162.9276E + 021.1820E + 022.0707E + 021.2958E + 02+F174.5658E + 011.0476E + 013.5719E + 011.0310E + 01=F189.0145E + 014.7859E + 019.4472E + 015.6305E + 01−F191.4666E + 018.9947E + 001.2179E + 014.1749E + 00+F206.3748E + 014.4406E + 015.2312E + 014.9270E + 01+F212.1777E + 023.4870E + 002.1847E + 024.1234E + 00−F221.0051E + 022.8400E + 001.0000E + 022.2124E−13=F233.6665E + 025.5240E + 003.6453E + 025.2131E + 00+F244.3676E + 023.7374E + 004.3642E + 026.0049E + 00=F253.8681E + 025.8187E-023.8684E + 027.5494E−02=F261.1214E + 036.7789E + 011.1160E + 037.3256E + 01+F275.0632E + 024.6020E + 005.0579E + 026.3243E + 00+F283.4497E + 025.8589E + 013.4205E + 025.8955E + 01+F294.6900E + 021.3609E + 014.4219E + 022.2017E + 01+F302.0864E + 031.3828E + 022.0213E + 031.6473E + 02+Total number of (+/ = /−): 18/6/6Table 9
Comparison results of EBLSHADE on CEC2017 test suiteD=30.
EBLSHADEESM-EBLSHADEMeanStd.DevMeanStd.DevF10.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F20.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F30.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F45.8562E + 013.4396E − 145.8562E + 014.0776E − 14=F56.6978E + 001.3325E + 006.9719E + 001.5418E + 00=F60.0000E + 000.0000E + 001.4398E − 084.1347E − 08=F73.7938E + 011.3728E + 003.6957E + 011.2423E + 00+F87.4627E + 001.5388E + 006.9654E + 001.6395E + 00+F90.0000E + 000.0000E + 000.0000E + 000.0000E + 00=F101.3903E + 032.5407E + 021.2476E + 032.3487E + 02+F113.3748E + 012.8992E + 014.0898E + 012.8428E + 01−F121.1107E + 033.2096E + 021.1021E + 033.2418E + 02+F131.8932E + 017.7131E + 001.6355E + 015.4803E + 00+F142.1636E + 011.1516E + 002.1493E + 013.6187E + 00+F153.9785E + 001.8483E + 003.6417E + 001.6998E + 00+F164.1583E + 013.4857E + 012.0657E + 015.3775E + 00+F173.1880E + 015.6772E + 003.1060E + 014.8829E + 00+F182.2262E + 011.5977E + 002.2496E + 011.2879E + 00−F196.2724E + 002.2466E + 005.6099E + 001.3474E + 00+F203.2672E + 016.5190E + 002.7681E + 016.2855E + 00+F212.0835E + 021.4350E + 002.0819E + 021.4361E + 00+F221.0000E + 020.0000E + 001.0000E + 028.1676E − 14=F233.5311E + 023.0584E + 003.5407E + 022.9996E + 00−F244.2796E + 021.7733E + 004.2807E + 021.8546E + 00−F253.8672E + 021.9899E − 023.8672E + 022.4427E − 02=F269.7950E + 024.0916E + 019.7385E + 022.9174E + 01+F275.0689E + 025.5012E + 005.0580E + 024.6504E + 00+F283.5816E + 026.2765E + 013.2838E + 024.8986E + 01+F294.3621E + 027.9105E + 004.3392E + 026.5290E + 00+F301.9967E + 037.3147E + 012.0385E + 038.6201E + 01−Total number of (+/ = /−): 16/9/5We further conducted Wilcoxon signed-ranked test with a 0.05 significance [40] to statistically verify the effectiveness of the ESM. In Table 10, the p-value obtained by ESM-DE/rand/1, ESM-DE/best/1, ESM-DE/current−to−best/1, ESM-j2020, and ESM-BeSD are all less than 0.05, which indicates a significant improvement is achieved by the proposed ESM. It also can be found that ESM-EJADE, ESM-EBSHADE, and ESM-EBLSHADE still maintain their advantages in terms of the value of R+ and R–, which also reflects that the ESM achieves more promising results.Table 10
Results obtained by the Wilcoxon test.
AlgorithmsR+R−P-valueESM-DE/rand/1 vs. DE/rand/1446.518.50.039ESM-DE/best/1 vs. DE/best/1376.588.50.047ESM-DE/current−to−best/1 vs. DE/current−to−best/1423420.031ESM-j2020 vs. j2020419.545.50.038ESM-EJADE vs. EJADE326.5138.50.081ESM-BeSD vs. BeSD338.5126.50.046ESM-EBSHADE vs. EBSHADE375.589.50.075ESM-EBLSHADE vs. EBLSHADE341.5123.50.089As mentioned before, the functions of CEC2017 can be divided into four categories: unimodal functions, simple multimodal functions, hybrid functions, and composition functions. To analyze whether the performance of the ESM is related to the type of function, the performance of the ESM-based algorithms for each type of function category are presented in Table11. As shown in Table 11, the ESM performs best on hybrid functions. Following the adoption of the ESM, the algorithm has an 81.25% probability of finding a better solution but only a 16.25% probability of finding a worse solution on hybrid functions. Besides, the performance of the ESM is almost the same on simple multimodal functions and composition functions and relatively worst on unimodal functions. Table 11 demonstrates that if the optimization problem is multimodal function, hybrid function, or composition function, then the proposed selection mechanism can be considered in the DE algorithm.Table 11
The classification statistics result of proposed ESM in eight algorithms.
Functions+(%)=(%)−(%)Unimodal functions45.8341.6712.50Simple multimodal functions64.2923.2112.50Hybrid functions81.252.5016.25Composition functions62.5022.5015.00


## 5.2. Convergence Analysis

To more clearly show the influence of the ESM on algorithm convergence and avoid diversion, four representative functions from CEC2017 are extracted, including F7, F16, F20, and F24, convergence curves for these functions with the three classic DEs and their relative ESM-based versions are presented. F7 is a simple multimodal function with a characteristic of convexity. From Figure5(a), we can find that all six comparison algorithms have a similar convergence rate; this also reflects that the ESM does not cripple the convergence speed of the algorithm on F7. From Figures 5(b) to 5(d), it can be seen that when each algorithm adopts the proposed ESM, it converges to a better value than the original algorithm with faster convergence. This is mainly because the algorithm can hardly escape when it suffers stagnation. The algorithm that uses the proposed ESM has a more excellent capability to keep search efficiency.Figure 5
Convergence curves of three classic DEs and relative ESM-DEs on four 30-D CEC2017 benchmark functions: (a) F7; (b) F16; (c) F20; (d) F24.
(a)(b)(c)(d)


## 5.3. Population Diversity Analysis

As the ESM adopts the entropy-based individual diversity measure, we further calculated the diversity of the main population and archive population on F7, F16, F20, and F24 for three classic DEs. It should be pointed out that we used the concept proposed in [41] to measure population diversity:(13)xj¯=1NP∑i=1NPxj,i,Population Diversity=1NP∑i=1NP∑j=1Dxj,i−xj¯2.As illustrated in Figure6(a), the main population and archive population have an almost uniform diversity on F7. From Figures 6(b)to6(d), it can be concluded that the archive population can maintain a relatively higher diversity than the main population at the intermediate stage for DE/current−to−best/1 and DE/best/1. To sum up, the archive population can maintain relative diversity during the search process, which makes it possible for the algorithm to recover from stagnation.Figure 6
Population diversity curves of main population and archive population of three classic DEs on four 30-D CEC2017 benchmark functions: (a) F7; (b) F16; (c) F20; (d) F24.
(a)(b)(c)(d)


## 5.4. Parameter Sensitivity Analysis

As illustrated in the previous sections, the ESM contains one parameterQ required adjustment. Q is a threshold which represents the stagnation steps. Table 12 shows the B-W values that represent the variation between better (“+”) and worse (“−”) number of the mean value after applying the ESM. It can be seen that the ESM-based algorithm performance is not so sensitive to the change of Q values over a finite range. For further comparison, A wider range of Q values and its corresponding B-W values are shown in Figure 7. In general, when the Q value is too large, the overall performance of the ESM begins to deteriorate. This is mainly because too large Q value reduces the probability of the ESM intervening in the algorithm; therefore, it is difficult for the ESM to function in limited iterations. In summary, considering the overall performance of the ESM, Q is set to 120.Table 12
The B-W values of all ESM-based algorithms on CEC2017 under differentQ values.
AlgorithmsQ=100Q=110Q=120Q=130Q=140ESM-DE/rand/12221222320ESM-DE/best/11314141314ESM-DE/current−to−best/12424242220ESM-j20201821201917ESM-EJADE1012111010ESM-BeSD1110131212ESM-EBSHADE1011121011ESM-EBLSHADE910111010Total117126127119114Figure 7
The B-W values at different Q values.


## 5.5. Scalability Analysis

The dimension of the test functions governs the difficulties on finding the global optimum. Higher dimensional functions are generally more difficult to solve. To verify the relationship between the dimensionality and the performance of the proposed mechanism, we evaluate the average performance of eight ESM-based DE algorithms in CEC2017 benchmark set. Table13 shows the results (+/=/−) of the considered algorithms at D=10, 30, 50, and 100 on the test functions, for intuition, we convert the corresponding B-W values into Figure 8. We can find that the performance of the ESM-based algorithm fluctuates slightly but does not degrade significantly with the increase of problem dimension compared to the original algorithm. In Table 14, the Friedman test [40], a widely used nonparametric test in the EA community, is used to validate the performance of all algorithms based on the mean value. It is not difficult to see that the p-values from D=10 to D=100 calculated by the Friedman test are all less than 0.05. Therefore, there are significant differences in the performance of the comparison algorithms on the corresponding dimensions, and ESM-EBLSHADE gets the first rank overall.Table 13
The statistical results of all ESM-based algorithms on CEC2017 under different dimensions.
AlgorithmsD=10D=30D=50D=100ESM-DE/rand/1 vs. DE/rand/123/5/224/4/222/6/225/4/1ESM-DE/best/1 vs. DE/best/117/8/519/6/518/6/619/5/6ESM-DE/current−to−best/1 vs. DE/current−to−best/123/4/326/2/225/2/323/5/2ESM-j2020 vs. j202022/5/323/4/321/5/422/4/4ESM-EJADE vs. EJADE16/8/618/5/720/5/519/5/6ESM-BeSD vs. BeSD17/7/618/7/519/6/518/5/7ESM-EBSHADE vs. EBSHADE16/9/518/6/618/5/717/6/7ESM-EBLSHADE vs. EBLSHADE15/10/516/9/515/11/416/8/6Figure 8
The curves of B-W values for the ESM-based algorithm compared to the original algorithm.Table 14
The Friedman test results of all algorithms on CEC2017 under different dimensions.
AlgorithmsD=10D=30D=50D=100Mean rankingRankESM-EBLSHADE2.132.622.882.5710.201EBLSHADE2.983.033.213.5512.772ESM-EBSHADE4.854.785.014.8919.533EBSHADE4.985.735.555.2121.474ESM-EJADE5.896.076.245.8724.075ESM-BeSD5.316.736.846.5725.456EJADE6.516.936.956.8827.277ESM-j20207.517.627.587.4130.128BeSD7.117.938.028.5731.639j20209.8210.899.5811.2141.5010ESM-DE/rand/110.0911.0311.5911.0743.7811ESM-DE/current−to−best/110.5511.6311.2111.8545.2412ESM-DE/best/110.1111.8211.8511.9845.7613DE/best/111.5312.6512.3312.5749.0814DE/rand/112.3112.7012.5812.3449.9315DE/current−to−best/111.8913.8312.8912.5951.2016Friedman-P-value0.000.000.000.00


## 6. Conclusion and Future Work

In DE, making full use of the successful solutions generated in iterations is a meaningful work to improve algorithm performance. In this study, an external selection mechanism (ESM) is proposed to restore the searchability of the algorithm when an individual is in a state of stagnation. The ESM mainly contains a successful solution archive mechanism and a crowding entropy diversity control strategy. It can be easily integrated into the existing DE algorithms to improve its performance further.Experiments are conducted on the CEC2017 benchmark sets cooperated with three classic DEs and five state-of-the-art DE variants. From experiment results, we can see that the ESM-based algorithm can significantly improve the original algorithm’s solution accuracy; the ESM also does not increase the computational complexity of the original algorithm due to the introduction of entropy. Further, experiments also show that ESM has a fairly positive effect on multimodal function, hybrid function, and composition function, its scalability on different problem dimensions also has a certain universal.Notably, experimental results proved that the ESM can effectively improve various DE variants’ performance, more study is required. Therefore, our future work will mainly focus on (1) improving the ESM to make it suitable for various single-objective DE; (2) researching the application of the proposed ESM in other evolutionary algorithms.



---

*Source: 4544818-2022-04-04.xml*