# An Entropy-Based Multiobjective Evolutionary Algorithm with an Enhanced Elite Mechanism

Authors: Yufang Qin; Junzhong Ji; Chunnian Liu
Journal: Applied Computational Intelligence and Soft Computing
 (2012)

---

## Abstract

Multiobjective optimization problem (MOP) is an important and challenging topic in the fields of industrial design and scientific research. Multi-objective evolutionary algorithm (MOEA) has proved to be one of the most efficient algorithms solving the multi-objective optimization. In this paper, we propose an entropy-based multi-objective evolutionary algorithm with an enhanced elite mechanism (E-MOEA), which improves the convergence and diversity of solution set in MOPs effectively. In this algorithm, an enhanced elite mechanism is applied to guide the direction of the evolution of the population. Specifically, it accelerates the population to approach the true Pareto front at the early stage of the evolution process. A strategy based on entropy is used to maintain the diversity of population when the population is near to the Pareto front. The proposed algorithm is executed on widely used test problems, and the simulated results show that the algorithm has better or comparative performances in convergence and diversity of solutions compared with two state-of-the-art evolutionary algorithms: NSGA-II, SPEA2 and the MOSADE.


---

## Body

## 1. Introduction

Optimization problems exist in all kinds of engineering and scientific areas. When there is more than one objective in an optimization problem, it is called a multiobjective optimization problem (MOP). Since these objectives are usually in conflict with each other, the goal of solving a MOP is to find a set of compromise solutions regarding all objectives rather than a best one as in single-objective  optimization problems. The solutions of MOP, also called as the Pareto-optimal solutions, are optimal in the sense that there exist no other feasible solutions which would decrease some criteria without causing the increase of at least one other criterion. Evolutionary algorithm (EA) is an optimization algorithm based on the evolution of a population. As it can search for multiple solutions in parallel, it has gained great attention from researchers. In recent years, many excellent EAs [1–4] have been proposed to solve the MOPs efficiently and MOEA has been recognized as one of the best methods to solve the MOPs.Generally, there are two performance measures in evaluating the Pareto-optimal solutions obtained by MOEA. One is the convergence measurement, which evaluates the adjacent degree between the Pareto solutions and the true optimal front. Another one is the diversity measurement, which evaluates the distribution of solutions in the objective space. In order to achieve good performance, many excellent strategies and methods have been presented in MOEA [1, 2, 5–9]. For the convergence, the elite mechanism has proved to be very helpful to accelerate the evolution of population [6]. The basic idea of the elite mechanism is that the information of good solutions, which have occurred in the progress of the evolution, is used to ensure the solution set converge to the optimal front as soon as possible. Its usual practice is that a certain number of best solutions are selected from the population as the parents to produce the good offspring [1]. However, in the early stage of the algorithm applying this strategy, because there are many dominated solutions existing in the population selected as the parents, the population cannot converge at a fast speed. In order to maintain the diversity of nondominated solutions, two main methods are applied. The first is using the grid to maintain the diversity [7]. It draws grids in the objective space and controls the number of solutions in a grid. Although this way can find the better solutions quickly, sometimes it cannot accurately reflect the global distribution of solutions because the grid position is fixed. The second one is the way based on the density [2, 8, 9]. Every solution obtains a value of the density and an outstanding density calculation can help to form a good distribution of solutions.Since 1948, Shannon [10] introduced the information theoretic entropy to measure information content of a stochastic process, which led to the establishment of the field of information theory. Then, many different applications for entropy are given in various fields. In solving the multiobjective optimization problems, Farhang-Mehr and Azarm [11] and Gunawan et al. [12] have applied the entropy to maintain the diversity of the solution set well in multiobjective problems and multilevel multiobjective problems. Wang et al. proposed the MOSADE algorithm [13], which combines the self-adaptive differential evolution and the crowding entropy-based diversity measure to obtain the nondominated solution set. In this algorithm, every solution can calculate its crowding degree through the improved the information entropy formula according to solutions’ distribution. In essence, this method is similar to the crowding distance in NSGA-II. Thus, for some three objective problems, this algorithm cannot obtain the very ideal solution set.In this paper, we propose a new MOEA to solve the MOP more effectively, in which an enhanced elitism makes the nondominated solutions play the better guide role and an entropy-based strategy is applied to preserve the diversity of the population. We call it an entropy-based multiobjective evolution algorithm with an enhanced elitism, namely, E-MOEA in brief. Specifically, we employ the enhanced elitism in which only the nondominated solutions in the union population are selected as the parents to ensure that the solution set converges to the optimal front more quickly. With the algorithm going on, the number of the nondominated solutions in union population will increase gradually. In order to keep the size of the elitist population (the maximum number of the elitist population in our algorithm is set asN) and maintain the diversity of solutions, the strategy based on entropy is applied. In this strategy, a region is determined by taking a solution as its center and the most crowded regions with the most uneven distribution of solutions are found through applying the entropy; then, in these regions, the most crowded solutions are deleted one by one. Compared with the MOSADE, the enhanced elitism can make the solution set approach to the Pareto-optimal front more easily in our algorithm, and the surface of the Pareto solution set is more uniform through the combination of the region and entropy. Experimental results on the 2-objective problems and the 3-objective problems show that the novel algorithm has better performance in both convergence and diversity, compared with NSGA-II, SPEA2 and MOSADE.


## 2.  Multiobjective Optimization Problems

In this paper, we consider the following continuous multiobjective optimization problem (continuous MOP):(1)minimizeF→=(f1(x),f2(x),…,fm(x))Tsubjecttox∈X,
where X⊂Rn is the decision space and x=(x1,…,xn)T∈Rn is the decision variable vector. F:X→Rm consists of m real value objective functions fi(x)(i=1,2,…,m), and Rm is the objective space.Leta=(a1,…,am)T, b=(b1,…,bm)T∈Rm be two vectors, and a is said to dominate b, denoted by a≺b, if ai≤bi for all i=1,…,m, and a≠b.A pointx*⊂X is called Pareto-optimal if there is no x⊂X such that F(x)≺F(x*). The set of all the Pareto-optimal vectors is called Pareto set, denoted by PS. The set of all the Pareto objective vector, PF={y⊂Rm∣y=F(x),x⊂PS}, is called the Pareto front [2, 9].


## 3. The Entropy-Based MultiobjectiveEvolutionary Algorithm with an Enhanced Elite Mechanism

In this section, we first present an enhanced elitist mechanism; then an entropy-based strategy is proposed to maintain the diversity of population. Finally, the entropy-based multiobjective evolutionary algorithm with the enhanced elitism is described.


### 3.1. The Enhanced Elitist Mechanism

Recent researches have proved that the elite mechanism is an excellent method to speed up the convergence of evolutionary algorithm. The population can produce a good offspring population through the elitism’s guide role, which achieves the rapid evolution of the population. On the basis of this idea, different forms of the elite mechanism have been proposed in some EAs [1, 14]. Among them, the most popular and effective one is used in NSGA-II, which first combines parent and offspring population then chooses a certain number of solutions from the union population as the parent population of the next iteration according to the nondominated sort and crowd distance assignment. As mentioned above, in the early running of the algorithm, a part of the dominated solutions may be selected into the elitist population so that the offspring produced by them are not good enough. Thus, we make some improvements to enhance the guiding role of good solutions when producing good offspring and we call it the enhanced elite mechanism.According to the elite mechanism, it is reasonable that the better the solutions chosen as the parents are, the better the offspring solutions which are produced by these parents are. Therefore, in order to enhance the guide of the elitism in our algorithm, we just only select all the nondomination solutions as the parents of the next iteration instead of a certain number of relatively good solutions, which may include dominated solutions because the number of nondominated solutions in the union population is less thanN during the early running of the algorithm. This enhanced mechanism can avoid several complex operations, such as the hierarchical nondominated sort and the crowd distance computation. On the other hand, all offspring solutions are generated by the nondominated population, which makes use of the best solutions to improve the efficiency of the evolution of the population. When the number of the nondominated solutions exceeds N, we delete the solutions with the worst distribution one by one through the entropy-based strategy, which accurately considers the distribution information of all solutions to obtain the elitist population with better distribution.


### 3.2. The Entropy-Based Strategy for Maintaining Diversity

In the enhanced elitist mechanism, the number of the nondominated solutions in union population will gradually increase with the algorithm going on. For keeping the maximum size of the elitist populationN and maintaining the diversity of the population, we proposed the entropy-based strategy by combining the regional information and the knowledge about entropy. The strategy deletes the most crowded solutions one by one, and all the related values will be recalculated before deciding which solutions should be deleted from the population, which can reflect the distribution of solutions dynamically and accurately. It is obvious that the key operation of the strategy is how to select the most crowded solution in nondominated population. First we select the most crowded region with the most uneven distribution through applying the information entropy, and then the solution which seriously influences the distribution of the solution set in this region is picked and removed from the elitist population. The previous process is looped until the size of the nondominated population decreases to N to form the elitist population. The strategy is described in detail as follows.For ease of operations, we order the nondominated solutions by one objective, and then, the region taking each solution as the center is defined as[fi(x)-D,fi(x)+D](i=1,2,…,m,0<D<L/2), where fi(x) represents the ith objective value of a solution x, D is a parameter which controls the area of the region, and L is the length of numerical range for the objective. Obviously, for a 2-objective problem, the region is a square; for a 3-objective problem, the region is a cube. Because the areas of those regions for all solutions are the same, there is no doubt that the more the total number of solutions existing in a region is, the more crowded the region is. Thus, it is easy to find the most crowded region based on the number of solutions in a region. If there is only one region including the most solutions, the most crowded solution in this region will be deleted. However, in most case, there may be several regions which include the same and maximum number of solutions. The region with the most uneven distribution of solutions will be selected from these regions by applying the knowledge about entropy.In light of the Shannon information theory [15], the entropy can be used to measure the uniformity of the probability distribution in a normalized system. Assume a stochastic process with n possible outcomes where the probability of the ith outcome is pi. The probability distribution of this process denoted as a probability vector P can be shown as
(2)P=[p1,p2,…,pn];∑i=1n‍pi=1,pi≥0.
This probability vector has an associated Shannon’s entropy, H, of the form
(3)H(P)=-∑i=1n‍piln⁡⁡(pi),
where piln(pi) is assumed to be zero when pi=0. This function is at its maximum. Hmax=ln(n), when all probabilities have the same value, and it is at a minimum of zero when one component of the P-vector is 1 and the rest of the entries in the P-vector are zero. Inspired by this, we compute the entropy of a region through taking the region as P and viewing every distance between two adjacent solutions as a pi. If a region has a more uniform distribution where the distances between the adjacent solutions are roughly similar, its entropy will be bigger, and when a region has a more uneven distribution where the distances between the adjacent solutions vary greatly, its entropy will be smaller.The schematic diagram for computing the entropy of a region is shown in Figure1. There are n solutions x1,…,xn sorted by an objective in this region, d1,…,dn-1 are the distances between two adjacent solutions, and dl and df are, respectively, the edge distances from two extreme solutions x1, xn, to the corresponding boundaries of the region at the objective ordered. In order to obtain a normalized system and be more reasonable to employ the entropy formula, some transforms need to be performed. Let d0=2*dl, dn=2*df, then
(4)pi=di∑i=0n‍di.
In light of (3), the entropies of those regions with the most solutions are computed and the region whose entropy is the smallest is obtained. It implies that not only this region includes the most solutions, but also its distribution of solutions in this region is the most uneven. Next, we will choose the most crowded solution from this region to delete. Here, a simple and effective method is used. In the region, first the two solutions with the smallest adjacent distance are found, and then we compare distances between them and their other adjacent solutions. Finally, the solution with smaller distance will be deleted from the population. A simple example is shown in Figure 2. x1 and x2 are found having the smallest distance, and then we compare the distance d2 between x1 and x3 with the distance d4 between x2 and x4. Since d4<d2, the solution x2 will be deleted. The main difference between our method and the archive truncation method in SPEA2 is that all solutions in our algorithm have been sorted and we only perform some comparisons among the adjacent distances. This method not only can reduce the computation, but also can keep the even solution distribution reasonably.Figure 1
Computation entropy of a grid.Figure 2
Selecting the worst solution in the grid.


### 3.3. The Framework of E-MOEA

Combing the basic evolutionary algorithm and the tradition of the method producing offspring (crossover and mutation) in genetic algorithm, we proposed the entropy-based multiobjective evolutionary algorithm with an enhanced elite mechanism (E-MOEA). The main steps are shown in the following.Parameters
We have the following:N (population size), T (maximum number of generation), P (the current population), Q (the elitist population), R (the offspring population).Step 1.
Generate an initial populationP and set t=0.Step 2.
Copy all the nondominated solutions inP to the population Q.Step 3.
If|Q|>N, reduce the size of Q by the entropy-based strategy one by one until |Q|=N.Step 4.
Ift=T, then go to Step 7.Step 5.
Execute recombination and mutation operators to theQ to obtain the offspring population R.Step 6.
P
=
Q
+
R, t=t+1, and go to Step 2,Step 7.
Output the current elitist populationQ.


## 3.1. The Enhanced Elitist Mechanism

Recent researches have proved that the elite mechanism is an excellent method to speed up the convergence of evolutionary algorithm. The population can produce a good offspring population through the elitism’s guide role, which achieves the rapid evolution of the population. On the basis of this idea, different forms of the elite mechanism have been proposed in some EAs [1, 14]. Among them, the most popular and effective one is used in NSGA-II, which first combines parent and offspring population then chooses a certain number of solutions from the union population as the parent population of the next iteration according to the nondominated sort and crowd distance assignment. As mentioned above, in the early running of the algorithm, a part of the dominated solutions may be selected into the elitist population so that the offspring produced by them are not good enough. Thus, we make some improvements to enhance the guiding role of good solutions when producing good offspring and we call it the enhanced elite mechanism.According to the elite mechanism, it is reasonable that the better the solutions chosen as the parents are, the better the offspring solutions which are produced by these parents are. Therefore, in order to enhance the guide of the elitism in our algorithm, we just only select all the nondomination solutions as the parents of the next iteration instead of a certain number of relatively good solutions, which may include dominated solutions because the number of nondominated solutions in the union population is less thanN during the early running of the algorithm. This enhanced mechanism can avoid several complex operations, such as the hierarchical nondominated sort and the crowd distance computation. On the other hand, all offspring solutions are generated by the nondominated population, which makes use of the best solutions to improve the efficiency of the evolution of the population. When the number of the nondominated solutions exceeds N, we delete the solutions with the worst distribution one by one through the entropy-based strategy, which accurately considers the distribution information of all solutions to obtain the elitist population with better distribution.


## 3.2. The Entropy-Based Strategy for Maintaining Diversity

In the enhanced elitist mechanism, the number of the nondominated solutions in union population will gradually increase with the algorithm going on. For keeping the maximum size of the elitist populationN and maintaining the diversity of the population, we proposed the entropy-based strategy by combining the regional information and the knowledge about entropy. The strategy deletes the most crowded solutions one by one, and all the related values will be recalculated before deciding which solutions should be deleted from the population, which can reflect the distribution of solutions dynamically and accurately. It is obvious that the key operation of the strategy is how to select the most crowded solution in nondominated population. First we select the most crowded region with the most uneven distribution through applying the information entropy, and then the solution which seriously influences the distribution of the solution set in this region is picked and removed from the elitist population. The previous process is looped until the size of the nondominated population decreases to N to form the elitist population. The strategy is described in detail as follows.For ease of operations, we order the nondominated solutions by one objective, and then, the region taking each solution as the center is defined as[fi(x)-D,fi(x)+D](i=1,2,…,m,0<D<L/2), where fi(x) represents the ith objective value of a solution x, D is a parameter which controls the area of the region, and L is the length of numerical range for the objective. Obviously, for a 2-objective problem, the region is a square; for a 3-objective problem, the region is a cube. Because the areas of those regions for all solutions are the same, there is no doubt that the more the total number of solutions existing in a region is, the more crowded the region is. Thus, it is easy to find the most crowded region based on the number of solutions in a region. If there is only one region including the most solutions, the most crowded solution in this region will be deleted. However, in most case, there may be several regions which include the same and maximum number of solutions. The region with the most uneven distribution of solutions will be selected from these regions by applying the knowledge about entropy.In light of the Shannon information theory [15], the entropy can be used to measure the uniformity of the probability distribution in a normalized system. Assume a stochastic process with n possible outcomes where the probability of the ith outcome is pi. The probability distribution of this process denoted as a probability vector P can be shown as
(2)P=[p1,p2,…,pn];∑i=1n‍pi=1,pi≥0.
This probability vector has an associated Shannon’s entropy, H, of the form
(3)H(P)=-∑i=1n‍piln⁡⁡(pi),
where piln(pi) is assumed to be zero when pi=0. This function is at its maximum. Hmax=ln(n), when all probabilities have the same value, and it is at a minimum of zero when one component of the P-vector is 1 and the rest of the entries in the P-vector are zero. Inspired by this, we compute the entropy of a region through taking the region as P and viewing every distance between two adjacent solutions as a pi. If a region has a more uniform distribution where the distances between the adjacent solutions are roughly similar, its entropy will be bigger, and when a region has a more uneven distribution where the distances between the adjacent solutions vary greatly, its entropy will be smaller.The schematic diagram for computing the entropy of a region is shown in Figure1. There are n solutions x1,…,xn sorted by an objective in this region, d1,…,dn-1 are the distances between two adjacent solutions, and dl and df are, respectively, the edge distances from two extreme solutions x1, xn, to the corresponding boundaries of the region at the objective ordered. In order to obtain a normalized system and be more reasonable to employ the entropy formula, some transforms need to be performed. Let d0=2*dl, dn=2*df, then
(4)pi=di∑i=0n‍di.
In light of (3), the entropies of those regions with the most solutions are computed and the region whose entropy is the smallest is obtained. It implies that not only this region includes the most solutions, but also its distribution of solutions in this region is the most uneven. Next, we will choose the most crowded solution from this region to delete. Here, a simple and effective method is used. In the region, first the two solutions with the smallest adjacent distance are found, and then we compare distances between them and their other adjacent solutions. Finally, the solution with smaller distance will be deleted from the population. A simple example is shown in Figure 2. x1 and x2 are found having the smallest distance, and then we compare the distance d2 between x1 and x3 with the distance d4 between x2 and x4. Since d4<d2, the solution x2 will be deleted. The main difference between our method and the archive truncation method in SPEA2 is that all solutions in our algorithm have been sorted and we only perform some comparisons among the adjacent distances. This method not only can reduce the computation, but also can keep the even solution distribution reasonably.Figure 1
Computation entropy of a grid.Figure 2
Selecting the worst solution in the grid.


## 3.3. The Framework of E-MOEA

Combing the basic evolutionary algorithm and the tradition of the method producing offspring (crossover and mutation) in genetic algorithm, we proposed the entropy-based multiobjective evolutionary algorithm with an enhanced elite mechanism (E-MOEA). The main steps are shown in the following.Parameters
We have the following:N (population size), T (maximum number of generation), P (the current population), Q (the elitist population), R (the offspring population).Step 1.
Generate an initial populationP and set t=0.Step 2.
Copy all the nondominated solutions inP to the population Q.Step 3.
If|Q|>N, reduce the size of Q by the entropy-based strategy one by one until |Q|=N.Step 4.
Ift=T, then go to Step 7.Step 5.
Execute recombination and mutation operators to theQ to obtain the offspring population R.Step 6.
P
=
Q
+
R, t=t+1, and go to Step 2,Step 7.
Output the current elitist populationQ.


## 4. Experimental Design and Results

In this section, a large number of experiments are conducted to test the performance of E-MOEA on the biobjective and the 3-objective problems. Specifically, our algorithm is compared with other advanced MOEAs: NSGA-II and SPEA2 which have the different strategy of constructing the elitist population. And then, the comparisons of the proposed E-MOEA and the MOSADE are presented.


### 4.1. Test Instances and Performance Metrics

In our experiment, the biobjective problem is from ZDT series: ZDT1, ZDT2, ZDT3, and ZDT6. The 3-objective problems we selected is composed of the DTLZ family of scalable test problems [16].There have been several metrics proposed for measuring the performance of the Pareto-optimal obtained by MOEAs. In our work, we choose the GD metric [17] and SP metric [18]. GD can measure the distance between the nondominated solutions obtained and the real Pareto-optimal front
(5)GD=∑i=1|Q|‍ei2|Q|,
where ei represents the Euclidean distance between the solution qi∈Q and the nearest member of p* (p* is a solution set of uniform sampling from the true Pareto-optimal front). In our experiment, we use 10000 uniformly spaced Pareto-optimal solutions as the approximation of the true Pareto front.The metric SP can be used to measure the diversity of obtained solutions. Here,e¯ is the average value of all ei, and m is the number of objectives:
(6)SP=1|Q|-1∑i=1|Q|‍e¯-ei2,ei=min⁡qj∈Q∧qj≠qi⁡∑k=1m|fk(qi)-fk(qj)|.The another indicator which is used usually to evaluate the diversity of the solution set isΔ [1]. However, this indicator works only for biobjective problems and cannot be used directly to evaluate for problems of more than two objectives. Based on the metric proposed in [19], the indicator is extended to fit problems of more than two objectives by computing the distance from a given point to its nearest neighbor. The indicator is modified as
(7)Δ=∑k=1m‍d(Ei,Ω)+∑X∈Ω|d(X,Ω)-d¯|∑k=1m‍d(Ei,Ω)+(|Ω-m|)d¯,
where Ω is a set of solutions, (Ei,…,Em) are m extreme solutions in the set of Pareto-optimal solutions, m is the number of objectives, and
(8)d(X,Ω)=min⁡Y∈Ω,Y≠X‖F(X)-F(Y)‖,d¯=1|Ω|∑X∈Ω‍d(X,Ω).


### 4.2. Comparison of E-MOEA and Other MOEAs

In this part, we will compare the E-MOEA proposed and two state-of-the-art algorithms, NSGA-II and SPEA2. All three algorithms are given real-valued decision variables. Simulated binary crossover (SBX) [20] and polynomial mutation (PM) [21] are applied with distribution indexes of ηc=2 and ηm=10, respectively. A crossover probability pc=1 and a mutation probability pm=1/n (where n is the number of decision variables for test problems) are used. In E-MOEA, D should not be too big in order to decrease the computation of our algorithm; on the contrary, our experimental results also show that if D is set too small, the algorithm also do not get good result. So for 2-objective problems, we set D=0.04, and for 3-objective problems, D=0.07. In all three MOEAs, the size of the population is 100 and the maximum number of function evaluation is 25000 for 2-objective problems. And for 3-objectives problems these two numbers are 200 and 80000, respectively.Four biobjective problems, ZDT1–3 and ZDT6, and three 3-objective problems DTLZ1–3 are used. For each test problem, 30 times runs are executed. Convergence metric GD and diversity metric SP are employed to evaluate the performance. The results are given in Table1, where Average and Sdt. dev, respectively, represent the mean and the standard deviation of indicators. In terms of the GD metric, we can clearly see that E-MOEA is nearer to the Pareto-optimal front than the others for all four 2-objective test problems. The reason for this is that the parent solutions are all made up of the nondominated solution in the current population in E-MOEA. Because of not using dominated solution, the population can approach the true Pareto-optimal front more quickly and efficiently. In light of the SP metric, for ZDT1, ZDT2, and ZDT3, E-MOEA is the best, SPEA2 is the second, and NSGA-II is the last. For ZDT6, SPEA2 is a little better than E-MOEA and they are much better than NSGA-II. The possible reason is that SPEA2 used the effective fitness assignment, however, which costs too much time. E-MOEA applies the strategy based on entropy to preserve the diversity of the population, which can evaluate the uniformity of distribution of solutions more scientifically so that the elitist population with excellent distribution is constructed.Table 1
Performance comparison of the three MOEAs on the test problems.
Test problems
MOEAs
GD
SP
Average
Std. dev
Average
Std. dev
E-MOEA
0.000221511
4.88744
e
-
05
0.00424942
4.3571
e
-
04
ZDT1
NSGA-II
0.000682375
1.39413
e
-
04
0.00837093
1.20323
e
-
03
SPEA2
0.000440864
4.85921
e
-
05
0.00673195
5.59690
e
-
04
E-MOEA
0.000117987
4.14203
e
-
05
0.00416924
3.63994
e
-
04
ZDT2
NSGA-II
0.000512546
1.338
e
-
04
0.00906621
1.9978
e
-
03
SPEA2
0.000358487
6.39963
E
-
05
0.00426378
6.47733
E
-
04
E-MOEA
0.000181589
4.08439
e
-
05
0.00463638
6.10117
e
-
04
ZDT3
NSGA-II
0.000501522
1.20519
e
-
04
0.00926551
1.5051
E
-
03
SPEA2
0.000467296
1.90819
E
-
04
0.00674565
2.00025
E
-
03
E-MOEA
3.4402
e −05
3.25
e
-
06
0.00295588
2.87438
e
-
04
ZDT6
NSGA-II
3.62327
e
-
05
3.28188
E
-
06
0.00696434
7.99786
E
-
04
SPEA2
3.55231
e
-
05
2.90700
E
-
06
0.00264811
3.26533
E
-
04
E-MOEA
0.00331697
2.602846
e
-
03
0.0126931
3.5927
e
-
03
DTLZ1
NSGA-II
0.0112144
2.227298
E
-
02
0.0254555
1.90124
E
-
02
SPEA2
0.00731387
2.51490
E
-
02
0.0134803
2.77126
E
-
02
E-MOEA
0.000606691
2.04103
e
-
05
0.0214341
1.4247
e
-
03
DTLZ2
NSGA-II
0.021829197
2.56158
e
-
03
0.0281401
2.09720
e
-
03
SPEA2
0.0024392
2.82469
e
-
04
0.0258512
7.13532
e
-
03
E-MOEA
0.000927311
2.6767
e
-
04
0.0211134
1.1026
e
-
03
DTLZ3
NSGA-II
—
—
—
—
SPEA2
—
—
—
—DTLZ serial test problems are proposed by Deb K, which can be set with the different number of objectives. In here, we choose the same settings as in [22]. For DTLZ1 and DTLZ2, all three algorithms can converge the Pareto-optimal front, but E-MOEA is the nearest to the true front, SPEA2 is the second, and NSGA-II is the last. For the diversity, E-MOEA and SPEA2 are much better than NSGA-II. Since DTLZ3 is the most difficult test problem, SPEA2 and NSGA-II cannot converge to the Pareto-optimal front. The reason is that both NSGA-II and SPEA2 always generate too much wild individuals in the obtained populations. Fortunately, E-MOEA can get solutions with good convergence and diversity for DTLZ3. Thus, Table 1 for DTLZ3 only gives the results obtained by E-MOEA, and “—" denotes that corresponding algorithms cannot obtain reasonable results.Figures3, 4, 5, and 6 show the final solutions obtained by three algorithms on four biobjective test problems. Obviously, we can see that solutions obtained by E-MOEA display the broadest and the most uniform distribution and are the nearer to the true Pareto-optimal front than the NSGA-II and SPEA2. Figures 7 and 8 clearly present the final solutions by three algorithms on DTLZ1 and DTLZ2. Figure 9 shows the final solutions by E-MOEA on DTLZ3.The final solutions obtained by three MOEAs on ZDT1.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on ZDT2.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on ZDT3.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on ZDT6.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on DTLZ1.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on DTLZ2.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2Figure 9
The final solutions obtained by E-MOEA on DTLZ3.


### 4.3. Comparison of E-MOEA and MOSADE

Both of the E-MOEA and MOSADE make use of entropy to maintain the diversity of the solution sets. In E-MOEA, we select the region with the worst distribution to keep the distribution of solutions through applying the information entropy formula. In MOSADE, the improved information entropy formula is used to update the archive, which maintains the diversity of the solution set. For these two algorithms, further experiments are conducted to compare their performances. For E-MOEA, all parameters are the same as the parameters of E-MOEA as described in Section4.2 except that the population size is 100 for all the test problems. For MOSADE, the population size and the external elitist archive size are 100, the lower and upper limits of mutant constant and crossover probability Fl=0.1, Fu=0.9, CRl=0, and CRu=0.5. The biobjective problems ZDT1–3 and ZDT6, the 3-objective problems DTLZ1–7 are considered. The GD and Δ are used as the evaluation indicators in this experiment, and results are presented in Table 2.Table 2
Performance comparison of the E-MOEA and MOSADE on the test problems.
Test problems
GD
SP
E-MOEA
MOSADE
E-MOEA
MOSADE
ZDT1
2.2151
e− 4
1.2485
e
-
3
0.19677
0.13159
±
4.88744
e
-
5
±
9.7574
e
-
5
±
1.7614
e
-
2
±
5.6921
e
-
3
ZDT2
1.1789
e− 4
9.8051
e
-
4
0.18056
0.12099
±
4.14203
e
-
5
±
4.9107
e
-
5
±
2.2941
e
-
2
±
7.9444
e
-
3
ZDT3
1.8158
e− 4
2.1620
e
-
3
0.49072
0.43783
±
4.08439
e
-
5
±
1.9962
e
-
4
±
1.5302
e
-
2
±
8.0801
e
-
3
ZDT6
3.4402
e− 5
2.5156
e
-
3
0.36294
0.13319
±
3.25
e
-
6
±
1.0967
e
-
4
±
1.5114
e
-
2
±
9.8303
e
-
3
DTLZ1
1.9635
e− 3
5.2512
e
-
3
0.24111
0.63755
±
1.0045
e
-
3
±
1.5650
e
-
4
±
1.1527
e
-
2
±
3.0464
e
-
2
DTLZ2
8.7374
e− 4
5.3344
e
-
3
0.27939
0.61323
±
5.6621
e
-
5
±
1.7429
e
-
4
±
1.1527
e
-
2
±
3.6732
e
-
2
DTLZ3
2.5954
e− 3
6.4632
e
-
3
0.30944
0.58601
±
3.0169
e
-
3
±
3.3253
e
-
4
±
0.10179
±
2.9955
e
-
2
DTLZ4
5.36887
e− 4
4.7736
e
-
3
0.51411
0.93230
±
2.0186
e
-
5
±
2.7427
e
-
4
±
0.1810
e
-
2
±
1.0819
e
-
1
DTLZ5
3.6979
e − 4
3.8930
e
-
3
0.42787
0.35529
±
1.0853
e
-
4
±
2.0137
e
-
3
±
1.5032
e
-
2
±
2.2054
e
-
1
DTLZ6
1.0800
e − 5
3.7701
e
-
3
0.43264
0.30720
±
9.9346
e
-
7
±
3.1470
e
-
4
±
1.35959
±
4.7503
e
-
2
DTLZ7
4.7858
e
-
3
2.3737
e−3
0.59107
0.85454
±
2.1121
e
-
3
±
0.001303
±
9.4664
e
-
2
±
3.8634
e
-
2The results obtained from Table2 show that these two algorithms have both the good convergence according to the GD; however, E-MOEA gets better values than MOSADE in all test problems except DTLZ7. This means the resulting Pareto fronts from E-MOEA are closer to the true optimal Pareto fronts. The solution set obtained by E-MOEA can converge more quickly to the optimal front, which may be due to the enhanced elitism applied in E-MOEA. As DTLZ7 has a wider range of the optimal solutions, the effect of convergence in E-MOEA is not better than the MOSADE which applied the differential evolution. For the indicator Δ, in the 2-objective problem, MOSADE has the more uniform distribution and in the 3-objective, E-MOEA has better results. The possible reason is that the crowding entropy diversity measure tactic in MOSADE is more effective to the test problems whose optimal front is an approximate curve. However, due to considering the spacial factor, the entropy-based strategy for maintaining diversity is more suitable to solve the problem whose optimal front is a surface. Thus, E-MOEA has better diversity for DTLZ1–4 and DTLZ7.


## 4.1. Test Instances and Performance Metrics

In our experiment, the biobjective problem is from ZDT series: ZDT1, ZDT2, ZDT3, and ZDT6. The 3-objective problems we selected is composed of the DTLZ family of scalable test problems [16].There have been several metrics proposed for measuring the performance of the Pareto-optimal obtained by MOEAs. In our work, we choose the GD metric [17] and SP metric [18]. GD can measure the distance between the nondominated solutions obtained and the real Pareto-optimal front
(5)GD=∑i=1|Q|‍ei2|Q|,
where ei represents the Euclidean distance between the solution qi∈Q and the nearest member of p* (p* is a solution set of uniform sampling from the true Pareto-optimal front). In our experiment, we use 10000 uniformly spaced Pareto-optimal solutions as the approximation of the true Pareto front.The metric SP can be used to measure the diversity of obtained solutions. Here,e¯ is the average value of all ei, and m is the number of objectives:
(6)SP=1|Q|-1∑i=1|Q|‍e¯-ei2,ei=min⁡qj∈Q∧qj≠qi⁡∑k=1m|fk(qi)-fk(qj)|.The another indicator which is used usually to evaluate the diversity of the solution set isΔ [1]. However, this indicator works only for biobjective problems and cannot be used directly to evaluate for problems of more than two objectives. Based on the metric proposed in [19], the indicator is extended to fit problems of more than two objectives by computing the distance from a given point to its nearest neighbor. The indicator is modified as
(7)Δ=∑k=1m‍d(Ei,Ω)+∑X∈Ω|d(X,Ω)-d¯|∑k=1m‍d(Ei,Ω)+(|Ω-m|)d¯,
where Ω is a set of solutions, (Ei,…,Em) are m extreme solutions in the set of Pareto-optimal solutions, m is the number of objectives, and
(8)d(X,Ω)=min⁡Y∈Ω,Y≠X‖F(X)-F(Y)‖,d¯=1|Ω|∑X∈Ω‍d(X,Ω).


## 4.2. Comparison of E-MOEA and Other MOEAs

In this part, we will compare the E-MOEA proposed and two state-of-the-art algorithms, NSGA-II and SPEA2. All three algorithms are given real-valued decision variables. Simulated binary crossover (SBX) [20] and polynomial mutation (PM) [21] are applied with distribution indexes of ηc=2 and ηm=10, respectively. A crossover probability pc=1 and a mutation probability pm=1/n (where n is the number of decision variables for test problems) are used. In E-MOEA, D should not be too big in order to decrease the computation of our algorithm; on the contrary, our experimental results also show that if D is set too small, the algorithm also do not get good result. So for 2-objective problems, we set D=0.04, and for 3-objective problems, D=0.07. In all three MOEAs, the size of the population is 100 and the maximum number of function evaluation is 25000 for 2-objective problems. And for 3-objectives problems these two numbers are 200 and 80000, respectively.Four biobjective problems, ZDT1–3 and ZDT6, and three 3-objective problems DTLZ1–3 are used. For each test problem, 30 times runs are executed. Convergence metric GD and diversity metric SP are employed to evaluate the performance. The results are given in Table1, where Average and Sdt. dev, respectively, represent the mean and the standard deviation of indicators. In terms of the GD metric, we can clearly see that E-MOEA is nearer to the Pareto-optimal front than the others for all four 2-objective test problems. The reason for this is that the parent solutions are all made up of the nondominated solution in the current population in E-MOEA. Because of not using dominated solution, the population can approach the true Pareto-optimal front more quickly and efficiently. In light of the SP metric, for ZDT1, ZDT2, and ZDT3, E-MOEA is the best, SPEA2 is the second, and NSGA-II is the last. For ZDT6, SPEA2 is a little better than E-MOEA and they are much better than NSGA-II. The possible reason is that SPEA2 used the effective fitness assignment, however, which costs too much time. E-MOEA applies the strategy based on entropy to preserve the diversity of the population, which can evaluate the uniformity of distribution of solutions more scientifically so that the elitist population with excellent distribution is constructed.Table 1
Performance comparison of the three MOEAs on the test problems.
Test problems
MOEAs
GD
SP
Average
Std. dev
Average
Std. dev
E-MOEA
0.000221511
4.88744
e
-
05
0.00424942
4.3571
e
-
04
ZDT1
NSGA-II
0.000682375
1.39413
e
-
04
0.00837093
1.20323
e
-
03
SPEA2
0.000440864
4.85921
e
-
05
0.00673195
5.59690
e
-
04
E-MOEA
0.000117987
4.14203
e
-
05
0.00416924
3.63994
e
-
04
ZDT2
NSGA-II
0.000512546
1.338
e
-
04
0.00906621
1.9978
e
-
03
SPEA2
0.000358487
6.39963
E
-
05
0.00426378
6.47733
E
-
04
E-MOEA
0.000181589
4.08439
e
-
05
0.00463638
6.10117
e
-
04
ZDT3
NSGA-II
0.000501522
1.20519
e
-
04
0.00926551
1.5051
E
-
03
SPEA2
0.000467296
1.90819
E
-
04
0.00674565
2.00025
E
-
03
E-MOEA
3.4402
e −05
3.25
e
-
06
0.00295588
2.87438
e
-
04
ZDT6
NSGA-II
3.62327
e
-
05
3.28188
E
-
06
0.00696434
7.99786
E
-
04
SPEA2
3.55231
e
-
05
2.90700
E
-
06
0.00264811
3.26533
E
-
04
E-MOEA
0.00331697
2.602846
e
-
03
0.0126931
3.5927
e
-
03
DTLZ1
NSGA-II
0.0112144
2.227298
E
-
02
0.0254555
1.90124
E
-
02
SPEA2
0.00731387
2.51490
E
-
02
0.0134803
2.77126
E
-
02
E-MOEA
0.000606691
2.04103
e
-
05
0.0214341
1.4247
e
-
03
DTLZ2
NSGA-II
0.021829197
2.56158
e
-
03
0.0281401
2.09720
e
-
03
SPEA2
0.0024392
2.82469
e
-
04
0.0258512
7.13532
e
-
03
E-MOEA
0.000927311
2.6767
e
-
04
0.0211134
1.1026
e
-
03
DTLZ3
NSGA-II
—
—
—
—
SPEA2
—
—
—
—DTLZ serial test problems are proposed by Deb K, which can be set with the different number of objectives. In here, we choose the same settings as in [22]. For DTLZ1 and DTLZ2, all three algorithms can converge the Pareto-optimal front, but E-MOEA is the nearest to the true front, SPEA2 is the second, and NSGA-II is the last. For the diversity, E-MOEA and SPEA2 are much better than NSGA-II. Since DTLZ3 is the most difficult test problem, SPEA2 and NSGA-II cannot converge to the Pareto-optimal front. The reason is that both NSGA-II and SPEA2 always generate too much wild individuals in the obtained populations. Fortunately, E-MOEA can get solutions with good convergence and diversity for DTLZ3. Thus, Table 1 for DTLZ3 only gives the results obtained by E-MOEA, and “—" denotes that corresponding algorithms cannot obtain reasonable results.Figures3, 4, 5, and 6 show the final solutions obtained by three algorithms on four biobjective test problems. Obviously, we can see that solutions obtained by E-MOEA display the broadest and the most uniform distribution and are the nearer to the true Pareto-optimal front than the NSGA-II and SPEA2. Figures 7 and 8 clearly present the final solutions by three algorithms on DTLZ1 and DTLZ2. Figure 9 shows the final solutions by E-MOEA on DTLZ3.The final solutions obtained by three MOEAs on ZDT1.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on ZDT2.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on ZDT3.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on ZDT6.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on DTLZ1.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2The final solutions obtained by three MOEAs on DTLZ2.
(a)
E-MOEA
(b)
NSGA-II
(c)
SPEA2Figure 9
The final solutions obtained by E-MOEA on DTLZ3.


## 4.3. Comparison of E-MOEA and MOSADE

Both of the E-MOEA and MOSADE make use of entropy to maintain the diversity of the solution sets. In E-MOEA, we select the region with the worst distribution to keep the distribution of solutions through applying the information entropy formula. In MOSADE, the improved information entropy formula is used to update the archive, which maintains the diversity of the solution set. For these two algorithms, further experiments are conducted to compare their performances. For E-MOEA, all parameters are the same as the parameters of E-MOEA as described in Section4.2 except that the population size is 100 for all the test problems. For MOSADE, the population size and the external elitist archive size are 100, the lower and upper limits of mutant constant and crossover probability Fl=0.1, Fu=0.9, CRl=0, and CRu=0.5. The biobjective problems ZDT1–3 and ZDT6, the 3-objective problems DTLZ1–7 are considered. The GD and Δ are used as the evaluation indicators in this experiment, and results are presented in Table 2.Table 2
Performance comparison of the E-MOEA and MOSADE on the test problems.
Test problems
GD
SP
E-MOEA
MOSADE
E-MOEA
MOSADE
ZDT1
2.2151
e− 4
1.2485
e
-
3
0.19677
0.13159
±
4.88744
e
-
5
±
9.7574
e
-
5
±
1.7614
e
-
2
±
5.6921
e
-
3
ZDT2
1.1789
e− 4
9.8051
e
-
4
0.18056
0.12099
±
4.14203
e
-
5
±
4.9107
e
-
5
±
2.2941
e
-
2
±
7.9444
e
-
3
ZDT3
1.8158
e− 4
2.1620
e
-
3
0.49072
0.43783
±
4.08439
e
-
5
±
1.9962
e
-
4
±
1.5302
e
-
2
±
8.0801
e
-
3
ZDT6
3.4402
e− 5
2.5156
e
-
3
0.36294
0.13319
±
3.25
e
-
6
±
1.0967
e
-
4
±
1.5114
e
-
2
±
9.8303
e
-
3
DTLZ1
1.9635
e− 3
5.2512
e
-
3
0.24111
0.63755
±
1.0045
e
-
3
±
1.5650
e
-
4
±
1.1527
e
-
2
±
3.0464
e
-
2
DTLZ2
8.7374
e− 4
5.3344
e
-
3
0.27939
0.61323
±
5.6621
e
-
5
±
1.7429
e
-
4
±
1.1527
e
-
2
±
3.6732
e
-
2
DTLZ3
2.5954
e− 3
6.4632
e
-
3
0.30944
0.58601
±
3.0169
e
-
3
±
3.3253
e
-
4
±
0.10179
±
2.9955
e
-
2
DTLZ4
5.36887
e− 4
4.7736
e
-
3
0.51411
0.93230
±
2.0186
e
-
5
±
2.7427
e
-
4
±
0.1810
e
-
2
±
1.0819
e
-
1
DTLZ5
3.6979
e − 4
3.8930
e
-
3
0.42787
0.35529
±
1.0853
e
-
4
±
2.0137
e
-
3
±
1.5032
e
-
2
±
2.2054
e
-
1
DTLZ6
1.0800
e − 5
3.7701
e
-
3
0.43264
0.30720
±
9.9346
e
-
7
±
3.1470
e
-
4
±
1.35959
±
4.7503
e
-
2
DTLZ7
4.7858
e
-
3
2.3737
e−3
0.59107
0.85454
±
2.1121
e
-
3
±
0.001303
±
9.4664
e
-
2
±
3.8634
e
-
2The results obtained from Table2 show that these two algorithms have both the good convergence according to the GD; however, E-MOEA gets better values than MOSADE in all test problems except DTLZ7. This means the resulting Pareto fronts from E-MOEA are closer to the true optimal Pareto fronts. The solution set obtained by E-MOEA can converge more quickly to the optimal front, which may be due to the enhanced elitism applied in E-MOEA. As DTLZ7 has a wider range of the optimal solutions, the effect of convergence in E-MOEA is not better than the MOSADE which applied the differential evolution. For the indicator Δ, in the 2-objective problem, MOSADE has the more uniform distribution and in the 3-objective, E-MOEA has better results. The possible reason is that the crowding entropy diversity measure tactic in MOSADE is more effective to the test problems whose optimal front is an approximate curve. However, due to considering the spacial factor, the entropy-based strategy for maintaining diversity is more suitable to solve the problem whose optimal front is a surface. Thus, E-MOEA has better diversity for DTLZ1–4 and DTLZ7.


## 5. Conclusion

In this paper, a novel entropy-based multiobjective evolutionary algorithm with an enhanced elite mechanism (E-MOEA) is proposed. The algorithm improves the elitism and presents a new strategy based on entropy to construct the elitist population. At first we only select the nondominated solutions in the population as the elitist solutions, and when the size of the nondominated solutions exceeds the size of population, we delete worse solutions one by one to preserve the diversity of the population through the entropy-based strategy. Experimental results on seven widely used popular test functions show that E-MOEA can obtain the solutions set with better or comparative convergence and diversity performances compared with NSGA-II, SPEA2, and MOSADE.As eliminating one solution needs to recalculate the entropies of the crowded regions, the worst time complexity of the E-MOEA isO(n3), which is same with SPEA2. The future research will be how to reduce the computational expense while keeping good performance.



---

*Source: 682372-2012-08-01.xml*