**ABSTRACT**

Many engineering problems are characterized with multi-objectives. Differential Evolution (DE), an Evolutionary Algorithm (EA), known to be fast and robust in numerical optimization is extended to multi-objective problems in this study. The new algorithm named Multi-objective Differential Evolution Algorithm (MDEA) adjusts the selection scheme of traditional DE to solve multi-objective problems. The algorithm also modifies the domination criteria for the population. The offspring generated in subsequent generations are improved before domination check is performed on the population in the final generation. Moreover, trial solution replaces the target solution if it is better or equal in all the objectives. The proposed algorithm is coded in MATLAB 7.0 and has been successfully applied to five common test problems and an engineering cantilever design problem. Good spread of quality Pareto optimal solutions are achieved. The algorithm produces more Pareto optimal solutions than the previous algorithms and retains the fast convergence and diversity exhibited by DE in global optimization. The algorithm is a good choice for solving many practical engineering problems with ease.

PDF Abstract XML References Citation

####
**How to cite this article**

*Journal of Applied Sciences, 9: 3652-3661.*

**DOI:**10.3923/jas.2009.3652.3661

**URL:**https://scialert.net/abstract/?doi=jas.2009.3652.3661

**INTRODUCTION**

Evolutionary Algorithms (EAs) are some of the most widely used algorithm for solving Multi-objective Optimization Problems (MOPs). EAs are applied to a wide range of problems from science to engineering design problems. They are powerful optimization algorithms for solving complex real-world MOPs in a single run. The primary reason why a problem has a multi-objective formulation is that it is not possible to have a single solution that optimizes all objectives. Therefore, an algorithm that gives a large number of alternative solutions lying on or near the Pareto-optimal front is of great practical value. In multi-objective optimization, there may not exist a solution that is best with respect to all objectives, instead, there are equally good solutions that are known as Pareto optimal solutions (Deb, 2001). A Pareto optimal set of solutions is such that when we go from any one point to another in the set, at least one objective function improves and at least one other worsens (Babu and Jehan, 2003). Neither of the solutions dominates over each other and all the sets of decision variables on the Pareto are equally good. Evolutionary algorithms have the ability to find multiple Pareto optimal solutions in one single simulation run because of their population-approach. There are two goals in multi-objective optimization. The goals are to discover solutions as close to the Pareto-front as possible and to find solutions as diverse as possible in the obtained non-dominated front.

Several studies have extended evolutionary algorithm to solve multi-objective numerical optimization problems (Deb, 2001; Deb *et al*., 2002; Madavan, 2002; Babu and Jehan, 2003; Xue *et al*., 2003; Rakesh and Babu, 2005; Babu *et al*., 2005; Robic and Filipic, 2005; Santana-Quintero and Coello, 2005; Fan *et al*., 2006; Reddy and Kumar, 2007). Over the past decade, a number of Multi-objectives Evolutionary Algorithms (MOEAs) have been suggested (Deb, 2001). Recently, researchers have extended Differential Evolution (DE) which is a family of EA to solve MOPs. It has been successfully used in solving single objective problems. Differential evolution is a very simple population-based, stochastic function minimizer and very powerful at the same time. DE handles continuous, discrete and integer variables. In terms of constraints, DE can handle multiple constraints. It has the advantages of simple structure, ease of use, speed and robustness. It does not make use of some probability distribution function in order to introduce variations into the population but instead, it uses the difference between randomly selected factors (individuals) as source of random variation for a third vector (individual), referred to as the target vector. Hence, the trial solutions which will compete with the parent solution are generated by adding the weighted difference vectors to the target vector.

Pareto Differential Evolution (PDE) proposed by (Abbass and Sarker, 2002) was the first algorithm to extend DE to MOPs. The PDE was compared to SPEA (Zitzler and Thiele, 1999) on two test problems and PDE was found to outperform it. Other studies by Madavan (2002), Xue *et al*. (2003) and Robic and Filipic (2005) have proposed Pareto Differential Evolution Approach (PDEA), Multi-Objective Differential Evolution (MODE) and Differential Evolution for Multi-Objective Optimization (DEMO), respectively. All these algorithms perform well on test problems.

Multi-Objective Differential Evolution (MODE) was proposed by Babu and Jehan (2003) and Babu *et al*. (2005). In the algorithm, the dominated solutions are removed from the population in each generation. Only the non-dominated solutions are allowed to undergo DE operations of mutation, crossover and selection. So, at each generation, the population size reduces. The offspring are placed into the population if they dominate the parents. The algorithm worked well on the test problems but the number of non-dominated solutions generated was few. The dominated solutions that are removed from the population initially can still be improved using the DE operations. During mutation stage, they may be combined with other vectors to produce offspring that may replace them instead of just removing them. In this way, at the end of the generation, the number of non-dominated solutions is increased with quality members giving wider Pareto optimal fronts. In the original DE with single objective, most of the initially generated solutions are infeasible. They undergo mutation and crossover before they are improved to be feasible. In the same way, the solutions in each generation should be allowed to mutate, recombine and undergo crossover. The solutions in the final generation are the only ones to be checked for domination. Therefore, domination check is done once instead of doing it in all the generations thereby reducing the number of function evaluation.

The algorithm proposed by Fan *et al*. (2006) is implemented by modifying the selection scheme of DE. In their selection scheme, the trial population is compared with its counterpart in the current population. If the trial candidate dominates the current population member it will survive from the tournament selection to the population of the next generation and replaces the current population otherwise the current population is retained. They suggest that if the trial solution is worse than the target solution in any one of the objectives, it should be discarded.

A new DE algorithm for Multi-Objective Optimization Problems (MOOPs) is proposed in this study and it is called multi-objective differential evolution algorithm (MDEA). The DE algorithm proposed is based on the existing DE algorithm proposed by Price and Storn (1997). The only difference is its implementation of multi-objectives.

**MULTI-OBJECTIVE DIFFERENTIAL EVOLUTION ALGORITHM (MDEA)**

MDEA combines the advantages of algorithm by Fan *et al*. (2006) and Babu and Jehan (2003) while overcoming the shortcomings of the algorithms. In this way, MDEA runs faster with more quality Pareto optimal solutions. The description of the MDEA is as follows. The vectors are randomly generated to create initial vectors and solutions to the problem. The generated solutions are allowed to undergo mutation, crossover and selection for the number of generations chosen. The solutions that evolve in the final generation are checked for domination and the dominated solutions are removed. The selection procedure of MDEA is similar to (Fan *et al*., 2006). The trial solution survives to the next generation if its objective function is better or equal in all the objectives to the target solution. MDEA will produce many non-dominated solutions on the Pareto front than the algorithm by Fan *et al*. (2006) and Babu and Jehan (2003). Moreover, the algorithm by Fan *et al*. (2006) has not been modified to handle constraints except bound constraints. MDEA can handle multiple constraints. If any of the constraints is violated, a high value (10^{8}) is added to the objective function to make the solution infeasible (Deb, 2001). In this way, the solution will not be selected when compared with other solutions because of high value. MDEA uses a constant **population structure** like traditional DE. The final population is the same like the initial population. So the number of non dominated solutions are less than or equal to population size. This is an excellent advantage of MDEA over algorithms by Fan *et al*. (2006) and Babu and Jehan (2003). Though the proposed MDEA can be used on any strategy, the strategy used in this study is DE/rand/1/bin which is the most widely used of all the ten strategies of DE.

The proposed MDEA methodology can be summarized in the following pseudo code steps:

(1) | Input the required DE parameters like population size (NP), crossover constant (CR), scaling factor (F), maximum generation, number of objectives, bound constraints etc. | |

(2) | Initialize all the vectors randomly in the limit of bound constraints. | |

(3) | Set the generation counter, G = 0 | |

(4) | Perform mutation and crossover operations as explained in the DE algorithm (Price and Storn, 1997) on all the population members. | |

• | For each parent, select three distinct vectors randomly from the current population. The selected vectors must not be the parent vector. These vectors combine to produce an offspring. So in DE, there are 3 parents that mutate to produce one offspring | |

• | Calculate new mutation vector using the expression, V_{i}(g) = X_{i3}(g) + F*(X_{i1}(g) – X_{i2}(g)) | |

• | Perform crossover using any of the two crossover method. Binary crossover method is used because of the strategy DE/rand/1/bin used in MDEA in this study | |

(5) | Evaluate each member of the population and check if it is better or equal to the parents. Replace the parents with offspring in the next generation if the offspring is better or equal to the parents otherwise, the parents proceed to the next generation. | |

(6) | Increase the generation counter, G, by 1. i.e., G = G+1. If G < GMAX, then go to step 4 above and repeat mutation, crossover and selection. If G = GMAX, then goto step 7. | |

(7) | Remove the dominated solutions in the last generation. A solution is dominated if there is another solution which is better than it in all the objectives. The method used here is naïve and slow suggested by Deb (2001). | |

(8) | Output the non-dominated solutions. |

The steps above are followed, using the pseudo code below, to code the algorithm in MATLAB 7.0 (The MathWorks Inc., USA) and is executed on a 1.7GHz, 2GB RAM PC and tested on some test problems below to demonstrate its ability to solve MOPs.

The Pseudo Code for Multi-Objective Differential Evolution Algorithm (MDEA):

*Initialize the values of D, NP, Cr, F, k (number of objective functions) and maximum generation (MAXGEN)

*input the boundary constraints of the problem

*Initialize all the vectors of the population randomly

For i = 1 to NP | ||

For j = 1 to D |

*remove the dominated solutions from the last generation using naïve and slow method proposed by Deb (2001) print the results

**EVALUATION AND RESULTS**

**Test problems:** Deb (1999) suggested systematic ways of developing test problems for multi-objective optimization. Zitzler * et al*. (2000) followed these steps and suggested well known 6 test problems. Five of the test problems (ZDT1, ZDT2, ZDT3, ZDT4 and ZDT5) which are common benchmark domains used in literatures (Deb* et al*., 2002; Madavan, 2002; Xue* et al*., 2003), are chosen for evaluating the performance of our methods. These test problems have no constraints. Each of them illustrates a different class of problems. All the problems have two objectives, f_{1}(x) and f_{2}(x) which are minimized. The description of the test functions are shown in Table 1.

**Performance measures:** The common trend in the development of many successful solution methodologies, including multi-objective **evolutionary algorithms** is to compare these in terms of their performance on various test problems. In a multi-objective optimization, there are two goals unlike single-objective optimization. The two goals are; discover solutions as close to the Pareto-optimal solutions as possible and find solutions as diverse as possible in the obtained non-dominated front. A good multi-objective evolutionary algorithm will be known if both goals are satisfied (Deb, 2001). Therefore, there is a need to have at least two performance metrics for adequately evaluating both goals of multi-objective optimization. One performance metric evaluates the progress towards the Pareto-optimal front and the other evaluates the spread of solutions. Many performance metrics have been proposed in the literatures. The descriptions of three common metrics used are given below though in this study, two performance metrics will be used to evaluate MDEA. They are convergence metric Υ and diversity metric Δ.

**Convergence metric Υ:** This metric measures the distance between the obtained non-dominated front, Q and the set, P* of the Pareto-optimal solutions. It is defined as:

(1) |

where, Q is the number of non-dominated vectors found by the algorithm being analyzed and d_{i} is the Euclidean distance (in the objective space) between the obtained non dominated front Q and the nearest member in the true Pareto front P.

Table 1: | Test problems used in this study |

**Generational Distance (GD):** This metric measures the average distance of the solutions of Q obtained from the set P* of the Pareto-optimal solutions. It is defined as:

(2) |

where, d_{i} is the Euclidean distance (measure in the objective space) between non-dominated front Q and the nearest member in the true Pareto-front P. It should be noted that an algorithm having a small value of GD is better.

**Diversity metric Δ:** This metric measures the extent of spread achieved using the non-dominated solutions. It is defined as:

(3) |

where, d_{i} is the Euclidean distance (measured in the objective space) between consecutive solutions in the obtained non-dominated front Q and is the average of these distances. The parameter d_{f} and d_{l} are the Euclidean distances between the extreme solutions of the Pareto front P* and the boundary solution of the obtained front Q. Also an algorithm having a smaller value of diversity metric Δ is better.

**Cantilever design problem:** MDEA is further tested on an engineering design problem of cantilever design (Deb, 2001). Figure 1 shows a schematic representation of a cantilever beam. Figure 1 presents the end load P, the length l and the diameter d. The problem has 2 decision variables of diameter (d) and length (*l*). The beam will carry an end load P. There are 2 conflicting objectives that should be minimized; the weight f_{1} and end deflection f_{2}. Minimizing the weight, f_{1}, will result in an optimum solution that will have small dimensions of d and l. If the dimensions are small, the beam will not be adequately rigid and the end deflection of the beam will be large. If on the other hand, the beam is minimized for end deflection, the dimensions of the beam will be large, thereby making the weight of the beam to be large.

Fig. 1: | A schematic representation of a cantilever beam |

There are 2 constraints in this design problem. They are; the maximum stress, σ_{max} is less than the allowable strength S_{y} and the end deflection, δ is smaller than a specified limit of δ_{max}. The constrained optimization problem is formulated as follows:

**Objective function 1:** Minimization of weight

(4) |

**Objective function 2:** Minimization of end deflection

(5) |

The objective functions 2 are subjected to the following constraints:

**Constraint 1:** Maximum stress

The maximum stress should be less than the allowable strength, S_{y}:

(6) |

**Constraint 2:** End deflection

The end deflection, δ should be smaller than a specified limit of δ_{max}

(7) |

**Bound constraints:**

(8) |

(9) |

where, ρ, P, d and *l* are the density, force, diameter and length respectively. The maximum stress is calculated as follows:

(10) |

The following parameter values are used:

ρ = 7800 kg m^{-3}, P = 1000 N, E = 207 GPa, S_{y}= 300 MPa and δ_{max}= 5 mm

**RESULTS**

MDEA was investigated for lower values of CR and F when the recommended values of 0.8 and 0.9 failed to achieve convergence in the high dimensional functions. MDEA performs better with lower values of CR and F than the values suggested in the practical advice on DE that F is to be a little lower or higher than 0.8 and CR to be 0.9 (Price and Storn, 2008). In test problems ZDT1, ZDT2 and ZDT3, the values of F and CR of 0.6 and 0.5, respectively give faster convergence and better diversity. The value of NP used in all the test problems is 100. The values of CR and F go down to as low as 0.3 and 0.35, respectively to make the algorithm robust enough for the test problems ZDT4 and ZDT6. The number of non-dominated solutions generated for all the test problems is 100. Therefore, all the solutions in the final generation (100) are non-dominated. This is an advantage of this algorithm over other EAs.

The statistics of the test problems are presented in Table 2-6. The convergence and diversity metric for algorithms NSGA-II (real coded), NSGA-II (binary coded), SPEA, PAES, PDEA, MODE, SDE and DEMO are taken from the literatures (Deb* et al*., 2002; Madavan, 2002; Robic and Filipic, 2005; Xue* et al*., 2003). They are used to compare the effectiveness of MDEA.

In Fig. 2a-c, the dominated solutions for NSGA, DEMO and MDEA are presented. It is found that the non dominated solutions for MDEA are of more quality than those of NSGA and DEMO.

Table 2: | Statistics of the results on test problem ZDT1 |

N/A: Not available |

Table 3: | Statistics of the results on test problem ZDT2 |

N/A: Not Available |

Table 4: | Statistics of the results on test problem ZDT3 |

Though the generated Pareto optimal front is similar to NSGA and DEMO, MDEA outperforms them. The convergence metric results in Table 2 show that MDEA has lower values in the metric than all the other algorithms except NSGA-II (binary coded). Also the diversity metric for MDEA is lower than all the other algorithms. The value is close to PDEA. This shows that MDEA is better in diversity than other algorithms.

Figure 3 shows the Pareto optimal fronts for NSGA, DEMO and MDEA for test problem ZDT2.

Fig. 2: | The results of test problem ZDT1 using (a) NSGAII, (b) DEMO and (c) MDEA |

It is found from Fig. 3a-c that MDEA is better in the uniform spread of non dominated solutions than NSGA. Though they all have the same shape, MDEA is better in diversity than NSGA but comparable to DEMO. Table 3 shows the results of the convergence and diversity metrics. It is found that MDEA is better in convergence than all the other algorithms but very close to DEMO. From the results of diversity metrics, MDEA is worse in diversity than NSGA-II (real coded), NSGA-II (binary coded), PDEA, SDE and DEMO but better than SPEA and PAES. The reason MDEA’s worse diversity in this test problem is that MDEA converged faster than the other algorithms. So the diversity is sacrificed for convergence in the test problem ZDT2.

Fig. 3: | The results of test problem ZDT2 using (a) NSGAII, (b) DEMO and (c) MDEA |

It is found in the Fig. 4a-c that MDEA outperforms both NSGA and DEMO. From the convergence and diversity results in Table 4, it is found that the value of convergence metric for MDEA is 0.001139. This value is lower than the ones obtained for other algorithms. DEMO’s convergence metric value is the closest to MDEA’s. Also the diversity metric gives 0.299354 which is the lowest for all the algorithms. DEMO obtains the closest value of 0.328873. Therefore, MDEA is the best algorithm for test problem ZDT3.

From the Fig. 5a-c, MDEA outperforms NSGA but comparable to DEMO. This is also confirmed with the results of convergence and diversity metrics in Table 5.

Fig. 4: | The results of test problem ZDT3 using (a) NSGAII, (b) DEMO and (c) MDEA |

Comparison with other algorithms shows that MDEA outperforms all of them except DEMO even in hard optimization (ZDT4) in both metrics.

Test problem ZDT6 Pareto optimal front is presented in Fig. 6a-c. MDEA is comparable to DEMO but better than NSGA in the figure. Results in Table 6 shows that MDEA outperforms all other algorithms in test problem ZDT6 in closer to those of SDE than other algorithms in both metrics. Therefore, MDEA is comparable to SDE in test problem ZDT6.

In Fig. 2-6, MDEA has achieved the two goals of multi-objective optimization which are convergence to the true Pareto-front and uniform spread of solutions along the front (Deb, 2001). From the figures, MDEA’s results are comparable to NSGA II and DEMO/parent. In all the Fig. 2-6, MDEA curve is comparable to that of DEMO but better than NSGA II in terms of convergence and diversity. It can be concluded from these results that MDEA’s performance is encouraging. The results have shown that MDEA is a good optimizer as proposed in this study. It is able to solve multi-objective optimization problems. The performance is comparable to other multi-objective Eas.

Comparison of the MDEA to MODE and NSGA on engineering test problem of cantilever design is shown in Fig. 7a and b. It is found that the non-dominated solutions generated by MDEA are comparable to those of MODE and NSGA. MDEA produces quality non dominated solutions along the Pareto front. This shows that MDEA can perform well on real-world engineering problems.

Table 5: | Statistics of the results on test problem ZDT4 |

N/A: Not available |

Table 6: | Statistics of the results on test problem ZDT6 |

N/A: Not available |

Fig. 5: | The results of test problem ZDT4 using (a) NSGAII, (b) DEMO and (c) MDEA |

Fig. 6: | The results of test problem ZDT6 using (a) NSGAII, (b) DEMO and (c) MDEA |

Fig. 7: | The results of (a) MODE and (b) NSGA for cantilever problem |

**CONCLUSION**

Multi-Objective Differential Evolution Algorithm (MDEA) is illustrated in this study. The selection scheme in the original DE is improved to accommodate multi-objective optimization. Also, dominated solutions are removed in the last generation only instead of removing them in all the generations thereby reducing the number of function evaluation. This allows the dominated solutions in every generation to be improved by mutation, crossover and selection operations of DE. It increases the number of non-dominated solutions giving us wider Pareto optimal fronts with better diversity. Moreover, DE has been shown to be capable of solving multi-objective high dimensional problem with few control parameters. The advantage of traditional DE which is ease of use is also applicable to MDEA. The test problems used to validate the algorithm show that MDEA is a good optimizer. The algorithm gives good spread of solutions maintaining diversity and convergence. At higher iterations, we have uniform distribution of solutions. Therefore, MDEA is a good alternative for solving engineering problems and other multi-objective problems. MDEA as a family of DE can also handle continuous, discrete, integer variables and multiple constraints.

####
**REFERENCES**

- Abbass, H.A. and R. Sarker, 2002. The Pareto differential evolution algorithm. Int. J. Art Int. Tools, 11: 531-552.

CrossRefDirect Link - Rakesh, A. and B.V. Babu, 2005. Non-dominated sorting differential evolution (NSDE): An extension of differential evolution for multi-objective optimization. Proceedings of the 2nd Indian International Conference on Artificial Intelligence, Dec. 20-22, Pune, India, pp: 1428-1443.

Direct Link - Babu, B.V. and M.M. Jehan, 2003. Differential evolution for multi-objective optimization. Proceedings of IEEE Congress on Evolutionary Computation, Dec. 8-12, IEEE Press, Canberra, Australia, pp: 2696-2703.

CrossRefDirect Link - Babu, B.V., J.H.S. Mubeen and G.C. Pallavi, 2005. Multi-objective differential evolution (MODE) for optimization of adiabatic styrene reactor. Chem. Eng. Sci., 60: 4822-4837.

CrossRefDirect Link - Deb, K., 1999. Multi-objective genetic algorithms: Problem difficulties and construction of test problems. Evol. Comput., 7: 205-230.

CrossRefDirect Link - Deb, K., A. Pratap, S. Agarwal and T. Meyarivan, 2002. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput., 6: 182-197.

CrossRefDirect Link - Fan, H.Y., J. Lampinen and L. Yeshayahou, 2006. An easy-to-implement differential evolution approach for multi-objective optimizations. Eng. Comput. Int. J. Computer-Aided Eng. Software, 23: 124-138.

CrossRefDirect Link - Madavan, N.K., 2002. Multi-objective optimization using a pareto differential evolution approach. Proceedings of the Evolutionary Computation on 2002, May 12-17, IEEE Computer Society, Washington, DC, USA., pp: 1145-1150.

Direct Link - Storn, R. and K. Price, 1997. Differential evolution-A simple and efficient heuristic for global optimization over continuous spaces. J. Global Optim., 11: 341-359.

CrossRefDirect Link - Reddy, M.J. and D.N. Kumar, 2007. Multiobjective differential evolution with application to reservoir system optimization. J. Comput. Civil Eng., 21: 136-146.

CrossRefDirect Link - Robic, T. and B. Filipic, 2005. DEMO: Differential evolution for multiobjective optimization. Proceedings of 3rd International Conference on Evolutionary Multi-Criterion Optimization, Mar. 9-11, Guanajuato, Mexico, pp: 912-912.

Direct Link - Santana-Quintero, L.V. and C.A.C. Coello, 2005. An algorithm based on differential evolution for multi-objective problems. Int. J. Comput. Intel. Res., 1: 151-169.

Direct Link - Xue, F., A.C. Sanderson and R.J. Graves, 2003. Pareto-based multi-objective evolution. Proceedings of Congress on Evolutionary Computation (CEC), Dec. 8-12, IEEE Press, Canberra, Australia, pp: 862-869.

CrossRefDirect Link - Zitzler, E., K. Deb and L.Thiele, 2000. Comparison of multiobjective evolutionary algorithms: Empirical results. J. Evol. Comput., 8: 173-195.

CrossRef - Zitzler, E. and L. Thiele, 1999. Multiobjective evolutionary algorithms: A comparative case study and the strength pareto approach. IEEE Trans. Evol. Comput., 3: 257-271.

CrossRefDirect Link