Abstract

Recently, decomposition-based multiobjective evolutionary algorithms have good performances in the field of multiobjective optimization problems (MOPs) and have been paid attention by many scholars. Generally, a MOP is decomposed into a number of subproblems through a set of weight vectors with good uniformly and aggregate functions. The main role of weight vectors is to ensure the diversity and convergence of obtained solutions. However, these algorithms with uniformity of weight vectors cannot obtain a set of solutions with good diversity on some MOPs with complex Pareto optimal fronts (PFs) (i.e., PFs with a sharp peak or low tail or discontinuous PFs). To deal with this problem, an improved decomposition-based multiobjective evolutionary algorithm with adaptive weight adjustment (IMOEA/DA) is proposed. Firstly, a new method based on uniform design and crowding distance is used to generate a set of weight vectors with good uniformly. Secondly, according to the distances of obtained nondominated solutions, an adaptive weight vector adjustment strategy is proposed to redistribute the weight vectors of subobjective spaces. Thirdly, a selection strategy is used to help each subobjective space to obtain a nondominated solution (if have). Comparing with six efficient state-of-the-art algorithms, for example, NSGAII, MOEA/D, MOEA/D-AWA, EMOSA, RVEA, and KnEA on some benchmark functions, the proposed algorithm is able to find a set of solutions with better diversity and convergence.

1. Introduction

In real-world applications, there are many problems needed to simultaneously optimize multiple objectives which are typically characterized by conflicting objectives. These problems are called as multiobjective optimization problems (MOPs). A continuous optimization problem can be formulated as follows [1]: where is a -dimensional decision variable bounded in the decision space , and is the number of objective functions. is the th objective function to be minimized, defines the th inequality constraint, and defines the th equality constraint. Moreover, all the inequality and equality constraints determine a set of feasible solutions which is denoted by , and is denoted as the objective space. Because the objectives often contradict each other, the improvement of one objective may cause to the deterioration of other objectives. So, MOPs have many optimal solutions which can be called nondominated solutions [2]. Some important definitions are introduced as follows. Let , is said to be better than , if and for . If there is no other such that is better than , is called Pareto optimal solution. The set of all the Pareto optimal solutions is defined as the Pareto set (PS). The image of the PS is called the Pareto optimal front (PF) [2].

Because the number of PF may be infinite, it is impractical to obtain all the Pareto optimal solutions. Thus, the principal goal of solving MOPs is to find out a set of solutions with good diversity and convergence. Currently, multiobjective evolutionary algorithms (MOEAs) use the strategy of the population evolution to simultaneously optimize the solutions of the population in a run. MOEAs can well deal with some complex problems which are characterized with discontinuity, multimodality, and nonlinearity [3]. Nowadays, many MOEAs [430] with good performance have been proposed, such as multiobjective genetic algorithms [3], multiobjective particle swarm optimization algorithms [710], multiobjective differential evolution algorithms [10, 11], multiobjective immune clone algorithms [12], group search optimizer [13], evolutionary algorithms based on decomposition [1417], and hybrid algorithms [8, 22]. Moreover, many MOEAs are used to solve numerous applications [3134].

Recently, Zhang and Li [16, 21] introduce the decomposition approaches into MOEA and developed an outstanding MOEA, MOEA/D, which has a superior performance for many problems. MOEA/D decomposes the MOP into a number of subproblems and uses the EA to optimize these subproblems simultaneously. The two main advantages of MOEA/D are that it uses the neighbor strategy to improve the search efficiency and well maintain the diversity of obtained solutions by the given weight vectors. In the last decade, MOEA/D has attracted many research interests and many related articles [1722] have been published.

In this work, we mainly study the refinement of weight vectors in MOEA/D to enhance the diversity of obtained solutions. Zhang and Li [16] claim that the weight vectors should be selected properly to obtain the nondominated solutions evenly distributed over the true PF. The basic assumption of MOEA/D is that the set of weight vectors with good uniformity can help obtained nondominated solutions to maintain the diversity. However, recent studies have suggested that MOEA/D which uses the fixed weight vectors might not well solve MOPs with complex PFs [35].

In this paper, we develop an improved decomposition-based multiobjective evolutionary algorithm with adaptive weight vector adjustment (IMOEA/DA) to solve MOPs. The main contributions of this paper are as follows: firstly, a new method [36] based on uniform design and crowding distance [5] is used to generate a uniformity of weight vectors; secondly, some weight vectors are adaptively deleted or added according to the distances of obtained nondominated solutions to solve the problems with complex PF; thirdly, a selection strategy is used to help each subobjective space to obtain a nondominated solution (if have). The frame of decomposition-based multiobjective evolutionary algorithm with adaptive weight vector adjustment and the initialization method of weight vectors has been studied. Moreover, the research result has been presented in the conference “2017 13th International Conference on Computational Intelligence and Security (CIS)”. In this conference paper, the adaptive weight adjustment [37] is used. In this new paper, a new adaptive weight adjustment is proposed.

The rest of this paper is organized as follows: Section 2 summarizes the related works of refinements of the weight vectors. Section 3 presents the proposed algorithm IMOEA/DA in detail, while the experiment results of the proposed algorithm and the related analysis are given in Section 5; finally, Section 5 provides the conclusions and proposes the future work.

MOEA/D uses the predetermined uniformly distributed weight vectors. Recent studies have shown that the fixed weight vector used in MOEA/D might not be able to cover the whole PF very well [35]. Therefore, some researches have refined the weight vectors in MOEA/D. Gu and Liu [38] periodically create the new weight vectors according to the distribution of the current set of weight vectors. Li and Landa-Silva [35] suggest that according to the strategy, the solution of each subproblem should be a long way from the corresponding nearest neighbor to adjust each weight vector. Qi et al. [37] propose an adaptive weight adjustment which utilizes the obtained nondominated solutions to reinstall the weight vectors. In the adaptive weight adjustment, the intersection angle of the target vector of each nondominated solution and the corresponding weight vector of this nondominated solution is zero. Jiang et al. [39] develop an adaptive weight adjustment by sampling the regression curve of objective vectors of the solution in an external population.

Other MOEAs use reference points to solve MOPs. These algorithms guide solutions to converge to the reference points. The principle of algorithms based on reference points or weight vectors is the same. Jain and Deb [40] adjust the reference points in terms of the distribution of candidate solutions in the current population at each generation. Jain and Deb [40] delete reference points with an empty niche and randomly add new reference points inside each crowded reference point with a high niche count. Cheng et al. [41] design two sets of reference vectors, where one maintains uniformly distributed and the other one is adaptively adjusted. Asafuddoula et al. [42] also adopt two sets of reference vectors, where one is called active set which is adaptively adjusted and the other one is called inactive set which stores the discarded reference vectors. In this algorithm [42], the two sets of reference vectors are tuned dynamically over the course of evolution.

3. The Proposed Algorithm

In this paper, an improvement decomposition-based multiobjective evolutionary algorithm with adaptive weight vector adjustment (IMOEA/DA) is proposed to address the MOPs with complex PF. The proposed algorithm mainly consists of two parts: a new weight vector initialization method based on uniform design and crowding distance and adaptive weight vector adjustment strategy, which will be introduced in this section.

3.1. Motivation

The main goal of this paper is to use decomposition-based multiobjective evolutionary algorithm to obtain a set of nondominated solutions which evenly distribute on the true PF and have a good convergence. In this paper, we adaptively add or delete some weight vectors to achieve this goal. In decomposition-based multiobjective evolutionary algorithms, the main role of weight vectors is to improve the convergence of obtained solutions by guiding the search of subproblems. Thus, weight vectors should maintain relative stability to improve the convergence of obtained solutions. We adaptively delete or add some weight vectors by the distances of obtained nondominated solutions to maintain relative stability of weight vectors and solve the problems with complex PF. In addition, we consider that the current optimal solutions of some subproblems are dominated solution, but their optimal solutions are nondominated solution. Some of corresponding weight vectors of these subproblems are retained, and a selection strategy is used to make each subproblem obtain a nondominated solution (if have).

3.2. A New Weight Vector Initialization Method-Based Uniform Design

In this subsection, the new weight vector initialization-based [36] uniform design and crowding distance is presented. Firstly, a uniform design method is briefly shown. For a given bounded and closed set (where is the dimension of the set ), the uniform design was developed to sample some points which have a small number and are uniformly scattered on . In this paper, we only consider a specific case of and introduce the main features of uniform design. More details can be obtained by referring the literature [34].

For a given set , in general, a set of exactly uniformly scattered points on is very difficult to be found. However, there are some efficient methods that can find a set of well approximately uniformly scattered points on . The good lattice point method (GLP) [43] is one of the simple and efficient methods and can generate a set of uniformly scattered points on . The details of GLP are as follows. For given integers , , and , a integer matrix called uniform array is denoted by where , and is the remainder of . Thus, there are different integer matrices be generated by these all μ. So, for given and , they can determine a number (Table 1 lists the vales of for different values of and ) which determines an integer matrix with the smallest discrepancy among these different integer matrices. In this paper, the discrepancy is denoted as , where is the fraction of the points falling in the hyperrectangle . In practice, the greatest common divisor of and should be 1 to reduce the amount of calculation, which is because that the integer matrix with the smallest discrepancy must be determined by these [43]. Each row of matrix determiners a point of by

is given by . Then each row of matrix defines a point of by

can be considered as a set of uniformly distributed weight vectors [44]. According to (4), we can obtain that many values of the first dimension of are closed to zeros (see an illustration in Figure 1(a)), which can reduce the diversity of . We use the following method to address this problem. We use to define a set , where , , . Then we use the crowding distance to select points from as the weight vectors (see an illustration in Figure 1(b)).

3.3. Adaptive Weight Vector Adjustment

In this subsection, the adaptive weight vector adjustment is presented. The main idea of this adjustment is that, if the distance of two adjacent nondominated solutions is large, some weight vectors are added between corresponding weight vectors of these two nondominated solutions and, if the distance of two adjacent nondominated solutions is small, one or two weight vectors of these weight vectors should be deleted. This adjustment strategy uses the distances of obtained nondominated solutions to delete or add some weight vectors to solve the problems with complex PF and maintain relative stability of weight vectors. The main difference of this adaptive weight vector adjustment and the method [37] is that the adaptive weight adjustment [37] utilizes the obtained nondominated solutions to reinstall the weight vectors, and the intersection angle of the target vector of each nondominated solution and the corresponding weight vector of this nondominated solution is zero; our adjustment strategy uses the distances of obtained nondominated solutions to delete or add some weight vectors. An illustration of our adjustment strategy is shown in Figure 2.

The detail of the adaptive weight vector adjustment strategy is as follows. For the current weight vectors and current population , where is the number of solutions or weight vectors and is the current optimal solution of the corresponding subproblem of the weigh vector , we find the nondominated solutions of . For convenience, we suggest that () are the nondominated solutions of and denote as . The distances of obtained nondominated solutions of is calculated as , where and . The values of are mainly used to delete the weight vector. In addition, all and are sorted to add the weight vectors. For convenience, we use to denote the distance of obtained nondominated solutions of and , where .

The deleting strategy is as follows. If (where is the size of the initial population), weight vectors with the minimum are deleted from . Then, if , the corresponding weight vector with the minimum is deleted from . After some weight vectors whose corresponding solutions are the dominated solutions are deleted from , the adding strategy is that, if the size of the current is smaller than , find the maximum distances of obtained nondominated solutions and new weight vectors are generated as follows: where , , is the size of , and the distance of and is one of the maximum distances of obtained nondominated solutions. The condition makes the optimal solution of the new subproblem generated by the weigh vector to be nondominated solution. In other word, we do not want that the generated weight vectors locate these spaces which have no nondominated solution. The adaptive weight vector adjustment is summarized in Algorithm 1.

Require: the size of the initial population , the current weight vectors and current population
Ensure: the weight vectors
Step 1: Find the non-dominated solutions of and denote . Calculate the and .
Step 2: Deleting weight vectors:
  If , then weight vectors with the minimum are deleted from .
  While do
    The corresponding weight vector with the minimum is deleted from and recalculate the and .
Step 3: Adding weight vectors:
If then
  Find the maximum distances of obtained non-dominated solutions, and use Eq. (5) to generate the new weight vectors.
Step 4: Deleting some weight vectors of from
If >0.5 then
    Use the crowding distance to delete -0.5 weight vectors of from .

In Step 4, some weight vectors of are kept, which is to record these regions with no nondominated solution and make these subproblems to quickly find nondominated solutions (if have).

3.4. The Proposed Algorithm IMOEA/DA

IMOEA/DA uses the evolutionary framework of MOEA/D and adopts the strategy [45] of allocating different amount of computational resources to different subproblems. The major roles of this strategy are to distribute computational efforts to each subproblem and optimize these subproblems with nondominated solution to generate nondominated solutions. The major differences between MOEA/D and IMOEA/DA are that the weight vector initialization method is different and IMOEA/DA updates the weight vectors during the course of evolution. The algorithm IMOEA/DA is shown as Pseudocode 1.

Input:
 MOP (1)
 A stopping criterion
: the number of the initial weight vectors (the sub-problems)
: the number of weight vectors in the neighborhood of each weight vector,
: a set of uniformly distributed weight vectors
Output: Approximation to the PF:
Initialization: Generate an initial population randomly or by a problem-specific method; determine by a problem-specific method; determine , where are the closet weight vectors to ; set and .
Whilethe stopping criterion is not metdo
  Step 1: Update the utility function as follows:
   According to , use the 10-tournament selection method [46] to select indices which are denoted as
  Step 2: For , do
   Select mating scope: uniformly create a random from and set
   Generate offspring : set and randomly choose one index from , use the SBX operator [47] to generate a solution from and and apply the polynomial mutation operator [47] on to generate a new offspring . Set
   Update of : For , if , then set
   Update the population by the updated strategy of the literature [44].
  end for
  Step 3: If is a multiple of 50, then, use Algorithm 1 to modify the weight vectors , re-determine (where is the size of ), and randomly select solutions from to allocate the new sub-problem as their current solution.
  Step 4: Set
end while

In this work, the aggregation function is the variant of Tchebycheff approach which formulation is as follows: where is the reference point of the MOP. The optimal solution of (6) must be the Pareto optimal solution of (1). If the optimal solution of (6) is not the Pareto optimal solution of (1), there is a solution which is better than , so , . Thus, is not the optimal solution of (6), which is a contradiction.

In Step 2, we use the strategy of the literature [44] which can well maintain the diversity of obtained solutions to update the solutions.

4. Experimental Results and Discussion

In this section, some experiments are conducted to demonstrate the effectiveness of the proposed algorithm IMOEA/DA. Firstly, IMOEA/DA compares with two other algorithms: MOEA/D [16] and NSGAII [5]. Secondly, we compare IMOEA/DA with other MOEA/D with weight vector adjustment strategy: MOEA/D-AWA [37] and EMOSA [35]. Thirdly, IMOEA/DA compares with RVEA [41] and KnEA [48]. Fourthly, we study the effectiveness of IMOEA/DA on MOPs with complex PF. Fifthly, we study the effectiveness of the initialization method of weight vectors and the adaptive weight vector adjustment strategy. Moreover, the effectiveness of IMOEA/DA on many-objective test problems is studied.

4.1. Test Problems

A wide range of well-known and very challenging test problems is selected to test the performance of the proposed algorithm IMOEA/DA in the experiments. These test problems include five biobjective ZDT test problems [2], seven triobjective DTLZ problems [49], ten problems CF of CEC09 competition [50], two problems F1 and F2 [37] with a sharp peak and low tail, nine biobjective WFG test problems [51], and fifteen three-objective MaF problems [52]. To investigate the ability of IMOEA/DA on many-objective problems, two problems DTLZ5(I,m) and DTLZ4(I,m) [53] are selected as the test instances. The IMOEA/DA is implemented by using MATLAB language on a PC with Intel Xeon CPU E3-1226 (3.30 GHz for a single core and the Windows 10 operating system.

4.2. Performance Metrics

In this paper, the true Pareto optimal fronts of the selected test problems are well known. In our experiment, there are four performance metrics are used to quantificationally compare with the performances of algorithms. These three metrics are generational distance () [54], inverted generational distance () [54], and hypervolume indicator () [55]. Wilcoxon rank-sum test [56] is used in the sense of statistics to compare the mean IGD, GD, and HV of the compared algorithms. It tests whether the performance of IMOEA/DA on each test problem is better (“+”), same (“=”), or worse (“−”) than/as that of the compared algorithms at a significance level of 0.05 by a two-tailed test.

4.3. Parameter Setting

All algorithms are coded as real vectors. Distribution index is 20 in the SBX operator; distribution index is 20 and polynomial mutation probability is in the mutation operator; the initial population sizes of all algorithms are set to 105 on all test problems and 105 initial weight vectors are generated; each algorithm is run 30 times with the maximal number of function evaluations 50,000 on all test problems. MOEA/D uses the method of the literature [16] to generate the weight vectors and uses the Tchebycheff approach as the aggregate function. For IMOEA/DA and MOEA/D, the size of neighborhood list is set to , the probability of choosing mate subproblem from its neighborhood is set to 0.9.

4.4. Comparisons of IMOEA/DA with Other Multiobjective Algorithms

In this subsection, metric values obtained by IMOEA/DA and other multiobjective algorithms on test problems are shown in Tables 26, and some comparisons are made to demonstrate the performance of IMOEA/DA. The best results obtained are highlighted bold in these tables.

4.4.1. Comparisons of IMOEA/DA with NSGAII and MOEA/D

In this subsection, IMOEA/DA is compared with NSGAII and MOEA/D on these twenty-six test problems which include nine WFG problems, ten UF problems, and seven DTLZ problems, and the results obtained three algorithms which are presented in Table 2. Table 2 presents the mean and standard deviation of the HV, GD, and IGD metric values of the final solutions obtained by three algorithms on twenty-six 10-dimensional problems. We can obtain from Table 2 that the results of Wilcoxon rank-sum test show that IMOEA/DA outperforms MOEA/D and NSGAII on all twenty-six test problems and twenty-four test problems in the form of the HV and IGD metrics, respectively; these indicate that the diversity of solutions obtained by IMOEA/DA is better than NSGAII and MOEA/D, and the solutions obtained by IMEOA/DA have a good convergence; in the form of the GD metrics, IMOEA/DA outperforms NSGAII on twenty-three test problems and outperforms MOEA/D on ten test problems, and these imply that the solutions obtained by IMEOA/DA have a good convergence. Moreover, it can be seen from Table 2 that, in the form of the HV and IGD metrics, the results obtained by IMOEA/DA are better than those obtained by NSGAII and MOEA/D on all twenty-six test problems, which indicates that the final solutions obtained by IMOEA/DA have a better diversity than those obtained by NSGAII and MOEA/D and have a good convergence. We also can obtain from Table 2 that the mean values of GD metric contained by IMOEA/DA are bigger than those obtained by MOEA/D on twelve test problems, which include WFG2, WFG8, UF2, UF3, UF6, UF8, UF9, DTLZ2–DTLZ5, and DTLZ7, and are smaller than those obtained by MOEA/D on other fourteen test problem; the mean values of GD metric contained by IMOEA/DA are smaller than those obtained by NSGAII on all these problems except for problems UF6, DTLZ4, and DTLZ7, and these imply that IMOEA/DAA can obtain a set of solutions with better convergence than MOEA/D and NSGAII on most test problems. Moreover, the mean values of IGD obtained by IMOEA/DA are smaller at least 15% than those obtained by MOEA/D and NSGAII on most test problems of these twenty-six test problems. To visually show the performance of the proposed algorithm, according to the median values of IGD metric, the nondominated solutions obtained by IMEOA/DA are plotted in Figures 3 and 4. It is evident that the found approximated PF is distributed uniformly on the true PF on these twenty-six test problems except for problem UF5. These comparisons indicate that IMOEA/DA is better at maintaining the diversity of obtained solutions than MOEA/D and NSGAII, and IMOEA/DA can also obtain a set of solutions with good convergence.

MOEA/D decomposes a MOP into a number of subproblems and optimizes them simultaneously. Each subproblem is optimized by using information from its several neighboring subproblems. In MOEA/D, the current solution of each subproblem is updated by these offspring which are generated by the neighbors of this subproblem, which can improve the convergence, and which cannot be conducted to maintain the diversity. In our algorithm, each offspring is firstly classified, the current solutions of the same class as the offspring may be updated [44], and our adaptive weight vector adjustment strategy can solve MOPs with complex PFs, which help obtained solutions to maintain the diversity. Moreover, the above comparisons indicate that IMOEA/DA is better at maintaining the diversity of obtained solutions than MOEA/D. NSGAII uses the nondominated sorting and the crowding distance to rank the solutions. During selection, NSGAII uses a crowded comparison operator that takes into consideration both the nondomination rank of an individual in the population and its crowding distance (i.e., nondominated solutions are preferred over dominated solutions, but between the two solutions with the same nondomination rank, the one that resides in the less crowded region is preferred). Our algorithm uses the neighboring solutions to generate the offspring, which can generate better offspring than NSGAII. And the above comparisons declare that IMOEA/DA is better at maintaining the diversity of obtained solutions and improving the convergence than NSGAII.

4.4.2. Comparisons of IMOEA/DA with MOEA/D-AWA and EMOSA

MOEA/D-AWA [37] and EMOSA [35] are an improved MOEA/D with adaptive weight vector adjustment. It has a good performance of solving the MOPs with complex PFs. IMOEA/DA is compared with MOEA/D-AWA and EMOSA on fourteen test problems which include five ZDT problems, five DTLZ problems, two constructed problems F1-F2, and two many-objective problems DTLZ4(3,6) and DTLZ5(3,6). In this experiment, EMOSA does not test many-objective problems DTLZ4(3,6) and DTLZ5(3,6). The experimental results of MOEA/D-AWA and EMOSA are directly obtained from the original literatures to make a fair comparison, and the population sizes of IMOEA/DA are set to 100, 300, and 252 for two-objective problems, three-objective problems, and six-objective problems, respectively; the maximal number of function evaluations is set to 50,000, 75,000, and 200,000 for two-objective problems, three-objective problems, and six-objective problems, respectively, and other parameter settings are the same as Subsection 4.3.

Table 3 shows the mean and standard deviation values of IGD metric obtained by IMOEA/DA, MOEA/D-AWA, and EMOSA on these fourteen test problems. We can obtain from this table that the mean values of IGD metric obtained by IMOEA/DA are smaller than those obtained by MOEA/D-AWA on all fourteen test problems except for problem F2 and are smaller than those obtained by MOEAS on all twelve test problems except for problem F2, which indicate that the quality of the final solutions obtained by IMOEA/DA is better than those obtained by MOEA/D-AWA and EMOSA on thirteen problems and eleven problems, respectively. These comparisons imply that IMOEA/DA is better at solving MOPs with complex PFs than MOEA/D-AWA and EMOSA on most problems of these problems. According to the median values of IGD metric, the nondominated solutions obtained by IMEOA/DA on the five ZDT test problems, F1 and F2, are plotted in Figures 3 and 4, which can visually show the good performance of IMOEA/DA. These indicate that IMOEA/DA can effectively approach the true PFs.

The adaptive weight vector adjustment strategy of MOEA/D-AWA uses the obtained nondominated solutions and crowding distance to adaptively set the weight vectors; however, the crowding distance based on nondominated solutions is not a good method to maintain the uniformity of weight vectors, especially in problems with three or more objective problems. In MOEAS, each weight vector is periodically adjusted to make its solution of subproblem far from the corresponding nearest neighbor. This can well maintain the uniformity of weight vectors for two-objective problems; however, it may lose efficacy for three or more objectives problems. In our algorithm, according to the distance of neighboring nondominated solutions, the weight vectors are adaptively set, which can well maintain the uniformity of weight vectors. And the above comparisons declare that IMOEA/DA is better at maintaining the diversity of obtained solutions than MOEA/D-AWA.

4.4.3. Comparisons of IMOEA/DA with RVEA and KnEA on MaF Problems

To demonstrate the effectiveness of the proposed algorithm MaOEAIR2, the benchmark suit for CEC2018 MaOP competition [52] is chosen, which have diverse characteristics and can test the strengths and weaknesses of MOEAs. There are fifteen many-objective benchmark functions (MaF) with box constrains in the solution space in this benchmark suit. In this work, for each test problem, the number of objectives is set to 3. RVEA [41] and KnEA [48] are used to compare with our algorithm. KnEA is a knee point-driven EA to enhance the convergence performance in many-objective optimization. RVEA uses the reference vectors to decompose the original multiobjective optimization problem into a number of single-objective subproblems and elucidate user preferences to target a preferred subset of the whole Pareto front. The codes of RVEA and KnEA are obtained from PlatEMO [57]. The initial population and initial weight vectors sizes are set to 105 for these fifteen problems; each algorithm is run 30 times with the maximal number of function evaluations 100,000 on all test problems.

Table 4 shows the mean and standard deviation values of IGD metric obtained by IMOEA/DA, RVEA, and KnEA on these fifteen test problems. We can get from Table 4 that IMOEA/DA is worse than RVEA and KnEA on none test problem in the form of the HV metric; these indicate that the diversity of solutions obtained by IMOEA/DA is better than RVEA and KnEA; in the form of the GD metrics, IMOEA/DA outperforms RVEA on ten test problems and outperforms KnEA on five test problems; these imply that the solutions obtained by IMEOA/DA have a good convergence; the mean values of IGD metric obtained by IMOEA/DA are smaller than those obtained by KnEA on all fifteen test problems and are smaller than those obtained by RVEA on all fifteen test problems except for problems MaF3, MaF9, MaF11, and MaF12, which indicate that the quality of the final solutions obtained by IMOEA/DA is better than those obtained by KnEA and RVEA on fifteen problems and eleven problems, respectively. These comparisons imply that IMOEA/DA is better at solving MOPs with complex PFs than RVEA and KnEA on most problems of these problems.

4.5. Performances on MOPS with Complex PFs and Many-Objective Problems

For the problems WFG2, UF5, UF6, ZDT3, and DTLZ7 which have discontinuous PF, the above experimental results show that, though IMOEA/DA cannot well solve problem UF5, IMOEA/DA can well solve other four problems; IMOEA/DA can obtain better performances than NSGAII and MOEA/D on problems WFG2, UF5, UF6, and DTLZ7; IMOEA/D can obtain better performances than MOEA/D-AWA on problems ZDT3 and DTLZ7. These show that IMOEA/DA can obtain good performances on most problems of the MOPs with discontinuous PF.

To study the performance of IMOEA/DA on problems with sharp peak or low tail PFs, we test the problems F1 and F2. The ideal PFs of F1 and F2 are , , respectively. It can be seen from Table 3 that the performance of IMOEA/DA is worse than that of MOEA/D-AWA on problem F2 and is better than that of MOEA/D-AWA on problem F1. Figure 5 shows the distribution of the final nondominated fronts obtained by IMOEA/DA on problems F1 and F2. As shown in Figure 5, IMOEA/DA can obtain a set of nondominated solutions with good uniformity. These results suggest that the weight adjustment of IMOEA/DA does improve MOEA/D significantly in the terms of uniformity for the MOPs with complex PFs. The possible reasons for the success of IMOEA/DA are that it increases the number of subproblems of those regions with the sharp peak or low tail PFs and it treats each subproblem equally.

To test the ability of the IMOEA/DA on many-objective problems, two test problems DTLZ5(3,6) and its variation DTLZ4(3,6) are used. They are many-objective problems with low-dimensional PF in the objective space; thus, their PFs are convenient for the visual display of the distribution of solutions. The ideal PFs of DTLZ5(3,6) and DTLZ4(3,6) are described as follows: . It can be seen from Table 5 that the performance of IMOEA/DA is better than that of MOEA/D-AWA on these two problems. Figure 5 shows the distribution of the final nondominated fronts obtained by IMOEA/DA on these two problems. As shown in Figure 5, IMOEA/DA can obtain a set of nondominated solutions with good uniformity and convergence. The possible reason for the success of IMOEA/DA on many-objective problems is that IMOEA/DA with the help of the proposed weight adjustment can obtain a good diversity. The computing efforts in IMOEA/DA are evenly distributed and have no preference to the boundary solutions and nondominated solutions.

4.6. Major Contributions of the Proposed Algorithm

This paper has two major contributions; one is the initialization method of weight vectors, and the other is the adaptive weight vector adjustment strategy. In this subsection, their effects are discussed. The initialization method of weight vectors has a great effect on three or many objective problems. So, in the experiments, the test problems are three-objective problems of Subsection 4.4.1 and the parameters are the same as Subsection 4.4.1. To study the effect of the initialization method of weight vectors, we compare the performances of MOEA/D with initialization method [44], initialization method [16], and the proposed method. For convenience, these three MOEA/D algorithms are denoted as MOEA/D1, MOEA/D2, and MOEA/D3, respectively. Table 5 displays the mean and standard deviation values of HV and IGD metrics obtained by three MOEA/D algorithms on ten three-objective problems. We can obtain from Table 5 that, in the form of the HV and GD metrics, MOEA/D3 outperforms MOEA/D1 and is not worse than MOEA/D2 on all ten test problems, and the mean values of HV and IGD metrics obtained by MOEA/D3 are better than MOEA/D1 and MOEA/D2, which can indicate that the diversity of solutions obtained by MOEA/D3 is better than those obtained MOEA/D1 and MOEA/D2. These comparisons imply that the uniformity of weight vectors generated by the proposed method is preferable than the methods [16, 44].

The major role of the adaptive weight vector adjustment strategy is to maintain the diversity of obtained solutions. To identify this, IMOEA/DA is compared with IMOEA/DA without t the adaptive weight vector adjustment strategy which is denoted as IMOEA/D. In this experiment, the test problem and parameters are the same as those found in Subsection 4.4.1. Table 6 shows the mean and standard deviation values of HV, GD, and IGD metrics obtained by IMOEA/DA and IMOEA/D on these twenty-six problems. It can be seen from Table 6 that IMOEA/DA obtains the best results on all these problems in the form of the HV and IGD metrics and the results of Wilcoxon rank-sum test show that IMOEA/DA outperforms IMOEA/D on all test problems; these indicate that the diversity of solutions obtained by IMOEA/DA is better than those obtained by IMOEA/D and the proposed adaptive weight vector adjustment strategy can help obtained solutions to maintain the diversity; in the form of the GD metric, the convergence of solutions obtained by IMOEA/DA is worse than those obtained by IMOEA/D, which is because that instability in search direction of IMOEA/DA leads to a decrease in convergence speed.

5. Conclusions

In this paper, an improved decomposition-based evolutionary algorithm with adaptive weight adjustment is designed to solve multiobjective problems with complex PFs. The goal of the proposed algorithm is to enhance the diversity of obtained solutions. In this work, a new method based on uniform design and crowding distance is designed to generate a uniformity of weight vectors; an adaptive weight adjustment strategy which some weight vectors are adaptively deleted or added according to the distances of obtained nondominated solutions is proposed to adaptively change the weight vectors; a selection strategy is used to help each subobjective space to obtain a nondominated solution (if have). Moreover, the proposed algorithm tests thirty-five test problems and compares with six well-known algorithms NSGAII, MOEA/D, MOEA/D-AWA, EMOSA, RVEA, and KnEA. Simulation results show that IMOEA/DA has competitive performances on most test problems against six comparison MOEAs (i.e., NSGAII, MOEA/D, MOEA/D-AWA, EMOSA, RVEA, and KnEA. These results also imply that the proposed weight adjustment can help the proposed algorithm to obtain a good diversity for MOPs with complex PFs, and the decomposition-based multiobjective evolutionary algorithm with adaptive weight adjustment can obtain a set of solutions with good diversity on some MOPs with complex PFs (i.e., PFs with a sharp peak or low tail or discontinuous PFs). A further study of the proposed method needs to be developed for its application in optimization with complicated PFs and high-dimensional optimization problems.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 61502290, 61401263, and 61672334), China Postdoctoral Science Foundation (no. 2015M582606), Industrial Research Project of Science and Technology in Shaanxi Province (no. 2015GY016), Fundamental Research Funds for the Central Universities (nos. GK201603094 and GK201603002), Natural Science Basic Research Plan in Shaanxi Province of China (nos. 2016JQ6045 and 2015JQ6228), and Special Scientific Research Project Fund of the Education Department of Shaanxi Province (no. 16JK1381).