Abstract
Multi-objective particle swarm optimization (MOPSO) algorithms are shown to have enormous potential in solving multi-objective optimization problems (MOPs). However, most MOPSO is difficult to balance the exploration and exploitation, which may cause some problems to find true Pareto fronts when tackling some complex MOPs. A multi-objective decomposition particle swarm optimization based on completion-checking (C-DMOPSO) is improved in this paper. The updating mode of velocity is changed dynamically according to the algorithm’s evolutionary process, which balances the exploration and exploitation effectively. In addition, simulated binary crossover and opposition-based learning are adopted to improve the diversity, and the archive set strategy is added to store the optimal solutions. Furthermore, polynomial mutation is performed in archive. The effectiveness of the proposed algorithm is tested by nineteen standard functions, including ZDT, DTLZ and UF, and the experimental results show that C-DMOPSO performs better on most of test problems.
Introduction
There are many problems require the optimization of several conflicting objectives simultaneously in most practical engineering applications, which affect each other and with high computational complexity. These problems are called multi-objective optimization problems (MOPs). As a single solution is not available to minimize all objectives, the optimal solution of MOPs is a set with many trade-off solutions. Multi-objective evolutionary algorithms (MOEAs) are easy to understand and implement without much specific information, so they are very popular in solving MOPs [4] and developed in real-word application such as rolling schedules [11], sensor networking and path planning [29]. Particle swarm algorithm (PSO) [12] is a group intelligent algorithm for simulating bird flight. The main idea of PSO algorithm is closely related to evolutionary algorithm: PSO uses a swarm mode which makes it to simultaneously search large region in the solution space to optimize objective functions just like evolutionary algorithms. Particles in the swarm are moved to find the optimal value by global search and local search, and the algorithm is simple and feasible in dealing with single objective optimization problems (SOPs), especially in optimization of optimal power flow [22] and manufacturing process parameters [21]. This has motivated researchers to extended it to MOPs [8].
The methods for solving MOPs can be divided into three categories: dominance-based, decomposition-based and indicator-based methods [7]. Among them, the dominance-based and the decomposition-based method are widely used. Dominance-based approach is a fundamental method to deal with MOPs, and NSGAII algorithm [6] is one of a remarkable MOEA based on the dominance-based framework. However, it is ineffective when dealing with many-objective optimization problems or the problems with complicated Pareto front. In decomposition-based method, the MOP is transformed into a series of scalar functions so that objective function can be optimized by solving single objective problem directly. One of the famous representative MOEAs based on decomposition is MOEA/D [27]. This approach has received much attention because it can be used to simplify MOPs. MOPSO/D [19], dMOPSO [25] are two multi-objective PSO algorithms based on decomposition. Because of its low complexity and good performance, the multi-objective PSO algorithm based on decomposition is widely studied. The third type of MOEA is indicator-based method, it mainly updates offspring according to searching for information depends on the performance of the indicator such as the hypervolume [2] and R2 [20], which are used to measure the fitness of a solution by assessing its contribution to the both convergence and diversity of a MOEA.
In addition, PSO possesses different improvements in many aspects recently [23]. In the aspect of parameter adjustment, the adjustment based on Bayesian techniques [26], exponential decreasing inertia weight strategy [16] are proposed as the algorithm with the fixed inertia weight usually cannot get satisfactory results. Inspired by the idea of time-varying inertia weight, learning factors changed with time is adopted in [17]. In the population structure, Gülcü [9] proposed a new multi-population parallel collaborative particle swarm optimization algorithm, which made it easy to the effective communication between the particles. Nevertheless, the use of multiple computers to exchange information increases the difficulty on the realization of the algorithm. In the ways of searching, Lin, Q [15] and others proposed a new particle swarm update formula, which can balance the population diversity and convergence ability. DDMOPSO [1] combined the idea of decomposition and dominance, which makes the particle swarm converge to the real front without using genetic operations. However, when dealing with high-dimensional problems, the algorithm is slow to converge and inefficient.
Therefore, when using PSO to deal with MOPs, the key is to update the individual optimal location and global optimal location, and make sure to: a) converge as fast as possible, b) gain a good distribution on the Pareto front, c) avoid falling into local optimum which result in premature convergence. Considering the optimized form of PSO, it can be divided into some aspects as follows: the expansion of population topology, the improvement of particle searching, the study of hybrid strategy PSO and the adjustment of parameters.
Existing studies have shown that a single speed update mode often leads to over-exploration or over-exploitation, which result in poor convergence or premature convergence. For MOPs, using the dominance method alone may cause the high complexity and slow convergence rate, while using the decomposition method alone may generate poor distribution on the Pareto front. Therefore, combining the two methods is a key factor to the effectiveness of the algorithm.
Based on the existing research, main improvements of this paper are listed bellow: For PSO, the speed update formula is improved for the problem that exploration and exploitation is difficult to balance in the process. The completion-checking factor is proposed in order to speed up the convergence rate of the algorithm. The speed update pattern is changed and the particle neighborhood size is adjusted dynamically by checking the algorithm process. In order to solve the problem that the population is easy to fall into local optimum during the updating process, simulated binary crossover (SBX) and opposition-based learning (OBL) are adopted in this paper to maintain the diversity of particle swarm. For the problem that it is difficult to make full use of the good information after updating of each generation, the archive set is introduced to preserve the good solutions in population. What’s more, polynomial mutation (PM) is used to the archive set to avoid the premature convergence.
The remainder of this paper is organized as follows. Section 2 provides the background and related work for understanding the rest of the paper. Section 3 explains the proposed approach in detail. In Section 4, the results is shown of comparative study. Finally, Section 5 provides the conclusions and some possible paths for future work.
Background and related work
Multi-objective optimization
Without loss of generality, a continuous or unconstrained multi-objective optimization problem can be described as follows:
PSO is a metaheuristic which aims to deal with continuous and unconstrained nonlinear optimization problems. PSO simulates the movements of a flock of birds to find food. In PSO, a particle is expressed by a solution x i ∈ Ω D , i = 1, 2, . . . , N. A population can be defined by one swarm with N particles. The swarm is evolved by updating both the velocity vt+1 and the position of each particle xt+1. The velocity of each particle is updated by the following equation:
Zhang [27] proposed a multi-objective algorithm based on decomposition, which transforms a MOP into a series of SOPs. In [13], the results show that different decomposition methods produce great influence on the performance of algorithm.
After considering the characteristics of each decomposition method, C-DMOPSO adopts Tchebycheff decomposition method. The optimal solutions of the aggregate function is the Pareto optimal solutions of MOPs, and the corresponding aggregate function formula is as follows:
The proposed C-DMOPSO algorithm combines the dominance and decomposition approach. The former is used to provide selection pressure in updating and maintenance of archive, while the latter transforms MOPs into a set of scalar aggregation problems. In order to make clear of C-DMOPSO, it consist of four parts, i.e. an update mode based on completion-checking factor, population evolution strategy, updating and maintenance of archive. In the following paragraphs, the implementation details of each component in C-DMOPSO will be explained step-by-step.
Dynamic update mode based on completion-checking factor
The iterative formula is the fundamental driving force for the stepwise optimization of the swarm. In the standard PSO, the speed and position of the child particles are determined according to Equations (2) and (3) in each generation, then the individual optimum and global optimum are updated gradually. In the updating process, the particles need to make better use of the current generation and historical information in order to speed up the convergence rate of the algorithm. Many literatures make further research on the velocity formula [15], and proposes different improvement methods to help the particles choose the next flight mode. Considering the characteristics of particle’s flight, this paper proposes a completion-checking factor ξ to help the particle select the information they needed to update their velocity.
Adaptation mechanism of evolutionary operator should be used to solve searching issue during different stages in evolutionary process [10]. In evolutionary process, the number of non-dominated solutions are augmented with the iterations increasing. The completion-checking factor ξ is determined by the size of non-dominated solutions in current archive set. When dealing with complex problems, the number of non-dominated solutions is very scant at the beginning of the algorithm, so the exploration is needed to augment the number of non-dominated solutions as quickly as possible. When the number of non-dominated solutions reaches a certain extent, it will be considered that the completion degree is better and the corresponding ξ is bigger. At this point, the global excellent information is needed to guide the next flight of particles during the updating of the particles. In addition, in order to avoid falling into local optimum and losing diversity, favorable perturbation from the archive set is imposed in particles’ flight. According to above idea, the particle’s update mode of velocity is proposed in Algorithm 1:
1:
2: if ξ now < ξ
3:
4: else
5: if rand < ms,
6:
7: else
8:
9: end
10: end
11:
In Algorithm 1,
At the beginning of the iteration, the number of non-dominated solutions is small with a small value of ξ, and the search pattern of the population are supposed to focus on the exploration. Namely, the individual history optimum is used to guide the particles’ flight. With the iterations increasing, the number of non-dominated solutions are increasing quickly to make the population close to the real Pareto front. By this time, the value of ξ increases gradually. In this case, the searching pattern should be changed to local search, that is, the global optimal guide the particles’ flight. Besides, the perturbation term is added to make the particles jump out of the local optimum in order to avoid premature convergence.
In addition, ξ also determines the neighborhood size of the particles, and T (t) = ξ * T (t - 1). When ξ is small, the neighborhood size is large, so that the neighborhood is used to find the global optimum. While, the neighborhood size decreases with a large ξ, then the population change to local search to approximate the optimal front at this point.
Population evolution strategy
Crossing and opposition-based learning operator
For the complex Pareto front, the population often tends to fall into local optimum during the iterative process. In order to make better use of the information in each generation, SBX [14] is adopted in C-DMOPSO to assist in generating offspring particles.
Besides, OBL is a strategy from machine intelligence which has been proven to be an effective technique on PSO [28]. It generates a reverse solution which is far from the local optimal to guide the population out of the local optimum region and the searching range is extended to facilitate the global optimum.
A parent particle intersects with a offspring particle to produce two new particles. Between two new particles, a particle is randomly selected for executing subsequent OBL operator. The pseudocode is listed in Algorithm 2.
1: X = POP
2:
3: Generate a random integer j in [1, N];
4:
5: Generate a random integer k in [1, 2];
6: X inew = OBL (new k );
7:
Re-initialize the particles
In order to speed up convergence of the algorithm, the re-initialization strategy is adopted [25]. If the particle does not improve its individual optimum during the flight, its age is increased by one. When the particle’s age reaches the maximum age predefined, the position of the particle is re-initialized. The individual optimum and global optimum of the current generation is utilized with the Gaussian distribution, which makes the particle fly to the favorable direction. The formula is shown in Equation (5).
In PSO, a particle is an independent agent which can search the objective space according to the experience of its own and its companions. Both pbest and gbest play key roles in guiding the search of particles. Therefore, how to choose pbest and gbest is a key problem when designing a MOPSO algorithm.
On the one hand, the pbest of a particle needs to be updated to make better use of its personal knowledge. As shown in Fig. 1, the curve represents a true Pareto front. In the event that there are too many particles in swarm, A will be removed based on the dominance method because D dominates A. However, it is apparent that A is a useful particle because it locates a relatively sparse area. The diversity may be lost if the algorithm updates the pbest under the dominance mechanism. As a consequence, different from other algorithms, the decomposition method is adopted in C-DMOPSO instead of non-dominated sorting, which uniformly generate N vectors in the target space. The particle swarm is divided into N subproblems, and the particles are updated step by step under the direction vector.

Illustration of why not use dominance mechanism to update pbest.
The individual optimal value of the particle plays a guiding role in the exploration of the whole population. For the updating of pbest, the Tchebycheff approach is adopted in C-DMOPSO to construct the aggregate function. In the ideal state, the individual optimal value of any particle in the population is updated as shown in Fig. 2(a). In the figure, the dotted line represents the direction vector in the target space, the real curve represents the real front of the function, and the hollow dots are the particles in the existing population. In the process of pbest updating, A is the position of the particle at the last generation, B is the position of this particle current generation, and C is the position of this particle’s individual optimum is to run next generation. The particle is updated along the direction vector of the corresponding subproblem in each generation, the individual optimal value is updated when the aggregate function of the offspring particle is better than the original individual optimal value at this direction vector. Otherwise, the original individual optimal value unchanged. The pseudocode is listed in Algorithm 3.
1:
2:
3: pbest i = x i
4: age i = 0
5:
6: age i = age i + 1
7:
8:
On the other hand, the global best of particles is used to guide exploitation of particles. The information of neighborhood particles is used to update the gbest of particles. The updating trajectory of global optimal value of a particle within the population is shown in Fig. 2(b). In the process of gbest updating, A is the position of the gbest at the last generation, B is the position of gbest current generation, and C is the position of the gbest is to run next generation. The gbest are updated by the optimal information of the corresponding neighborhood particles. Unlike the SOP, the gbest of a MOP is not a particle, but the optimal particle with N particles in the corresponding neighborhood region. And their fitness is not a value, but a Pareto approximation front.

Update pbest and gbest.
The selection of the archive set: in order to distinguish good solutions from the whole population, an archive set is selected to record the information of the elite particles in population according to the idea of dominance relation. Most of the existing literatures [1, 15] choose the non-dominated solutions of each generation population to join the archive set. In this paper, the selection of elite particles obeys the following rules: the pbest and gbest of the population are sorted by the dominance relation after each iteration, only the first front’s particles are selected as the elite particles to add the archive set.
The maintenance of the archive set: as the increase of the iteration, the number of non-dominated solutions is increasing. Hence, the crowding distance strategy is used to maintain the archive scale. When the archive set exceed the predetermined scale, the crowding distance between the particles is calculated, and the particles with the smallest crowding distance are removed. The above steps are repeated until the archive size is maintained at the predetermined value.
Archive set variation: premature convergence phenomenon often emerged in PSO when it updates to a certain extent when dealing with complex problems. In order to jump out of premature convergence and improve the diversity, polynomial mutation(PM) [3] is imposed on the archive.
Algorithm 4 lists the pseudocode of updating and maintenance of archive.
1:
2: archive i = PM (archive i )
3:
4:
5: Execute non-dominated sorting on archive;
6: Calculate the crowding distance d i for each particle in archive;
7: Remove the archive i with d min ;
8:
1:
2: Initialize velocity and position of the population;
3: Set pbest = gbest = POP and archive = Φ;
4: Initialize the reference point
5: Initialize the N weight vectors and the particle neighborhood;
6:
7:
8:
9:
10: Update the particle according to Algorithm 1 and Eq.(3)
11:
12: Re-initialized the particles according to Eq.(5)
13:
14:
15: Execute Algorithm 2 of population evolution strategy;
16: Calculate the fitness and update the z*;
17: Update the pbest according to Algorithm 3;
18: Update the gbest according to the neighborhood information;
19: archive = archive ⋃ pbest ⋃ gbest, execute Algorithm 4;
20:
21: Return the final result in the archive.
The complete C-DMOPSO algorithm
The above subsections have described the procedures of velocity updating mode, evolution strategy, the updating and maintenance of archive, which compose the main components of C-DMOPSO. The complete C-DMOPSO algorithm is illustrated in Algorithm 5, where N is the size of population and archive, evaluation refers to the number of function evaluations, and FEs represents the maximum number of function evaluations.
Experimental results
In order to assess the performance of the proposed approach, the results of C-DMOPSO are compared with respect to those obtained by MOEA/D [27], dMOPSO [25], NSGAII [6], MOPSO [5] and a newly-proposed algorithm MIMOPSO [24].
Test instances
Nineteen test problems are adopted to estimate the performance of algorithms, whose Pareto fronts have different characteristics including convexity, concavity, connection, disconnection and multi-modality. They are two-objective test suites of Zitzler-Deb Thiele(ZDT), three-objective problems of Deb-Thiele-Laumanns Zitzler(DTLZ) and CEC’09 benchmark problems(UF1-UF10). There are 30 decision variables for ZDT1-ZDT3 and UF1-UF10. ZDT6 is tested using 10 decision variables. For DTLZ2 and DTLZ4, 12 variables are adopted. DTLZ1 and DTLZ7 are tested with 7 and 22 decision variables respectively.
Experimental settings
The parameters of C-DMOPSO are as follows: inertia weight w start = 0.9, w end = 0.4, learning factors c1 = c2 = c3 = 1.5; neighborhood size T = 20, crossover probability pc = 0.5, mutation probability pm = 0.05, distribution index η c = η m = 20; agemax = 2.
ZDT: m = 2, N = 100, FEs = 30000; DTLZ: m = 3, N = 210, FEs = 210000; UF: N = 300, FEs = 300000. For each test function, all the algorithms run 30 times independently.
Parameter settings of all compared algorithms are shown in Table 1.
Parameter settings of all compared algorithms
Parameter settings of all compared algorithms
The inverse generation distance(IGD) [18] can measure both convergence and diversity of an algorithm. The lower the value, the better the convergence and diversity of the Pareto solution set, namely the approximate Pareto front of C-DMOPSO gained are closed to the true Pareto front. Let P be a proper representation of the true Pareto front, the IGD for a set of approximated solutions P′ is calculated as:
In order to verify the validity of the proposed algorithm, five kinds of contrastive algorithms are selected: MOEA/D, dMOPSO, NSGAII, MOPSO and MIMOPSO. Furthermore, for the purpose of testing the validity of the updating mode based on completion-checking factor, the algorithm with the original velocity formula is compared with C-DMOPSO as well. Tables 2 and 3 show the mean and standard deviation(std) of the IGD performance metrics for seven algorithms, where C-D refers to C-DMOPSO and C-D1 represents the algorithm with the original velocity formula. In these tables, the best results of the mean or standard deviation for each test function are identified with bold font. In addition, the p-value(p) of Wilcoxsom rank sum test with significant level α = 0.05 are also recorded for each function, which examine that whether the IGD mean values obtained by C-DMOPSO are different from that obtained by the other algorithms statistically. Namely, when p-value is bigger than 0.05, it means that the compared IGD mean values are similar in statistics. The bold italics in tables show that there is no significant difference between C-DMOPSO and the contrastive algorithm in the corresponding test function under the Wilcoxsom rank sum test. The better/similar/worse (+/=/-) show that C-DMOPSO is superior to, similar to, and inferior to the number of the comparison algorithm on all test functions in table.
IGD values on ZDT and DTLZ
IGD values on ZDT and DTLZ
IGD values on UF
As shown in Table 2 and Fig. 3, C-DMOPSO produces a better diversity along the Pareto front on ZDT functions except ZDT6. However, the differences between C-DMOPSO and dMOPSO in the IGD values obtained are not significant on ZDT6. ZDT4 has many local Pareto fronts, while C-DMOPSO attains a better diversity compares with other five algorithms except MIMOPSO. Meanwhile, Fig. 4 indicates that C-DMOPSO also gains good distribution along the Pareto front on DTLZ problems. But on DTLZ1, C-DMOPSO performs slightly worse than MOEA/D, NSGAII and MIMOPSO. On DTLZ4, the performance of C-DMOPSO is poor and with low std value, which because there are some bad IGD value during 30 times running. It is obvious that C-DMOPSO is better than dMOPSO in 6 test functions and only worse than it on ZDT6. C-DMOPSO is significantly better than MOPSO on all 9 test problems.

Pareto front of C-DMOPSO, MOEA/D, dMOPSO, NSGAII, MIMOPSO on ZDT3 and ZDT4. (a)-(f) represent the compare algorithms on ZDT3 and (g)-(l) represent the compare algorithms on ZDT4.

Pareto front of C-DMOPSO, MOEA/D, dMOPSO, NSGAII, MIMOPSO on DTLZ1, DTLZ4 and DTLZ7, (a)-(f) represent the compare algorithms on DTLZ1, (g)-(l) represent the compare algorithms on DTLZ4, and (m)-(r) represent the compare algorithms on DTLZ7.
As can be seen in Table 3, C-DMOPSO obtains weak performance than NSGAII on UF5 and UF6 in terms of reducing IGD value. In fact, UF5 has 21 local Pareto fronts which may impede the convergence ability of C-DMOPSO, while UF6 has one isolated point and two disconnected parts which may cause difficulty in finding a good diverse set of solutions. The permissible function evaluations might not be good enough to deal with such complicated problems. In summary, the results show that the overall IGD value obtained by the proposed algorithm is superior to the other six contrast algorithms on most of the nineteen test problems. According to no free lunch theory, we can not expect the proposed algorithm is better than the contrast algorithm on each function. C-DMOPSO gains good values on most of test problems compared with C-D1, which indicates that the update mode based on completion-checking factor is effective in handling complex MOPs. The final assessment of the nineteen functions is recorded in the last column of Table 3, which verifies that the proposed algorithm outperformed among MOEA/D, dMOPSO, NSGAII, MOPSO and MIMOPSO in most of the functions with respect to the IGD.
Figures 3, 4 and 5 show the approximate Pareto fronts obtained by C-DMOPSO, MOEA/D, dMOPSO, NSGAII, MOPSO and MIMOPSO in some test problems with clearer discrimination. It can be seen that C-DMOPSO produces better approximated PF on ZDT, DTLZ, and most UF test instances in these figures which demonstrate the good distributivity of C-DMOPSO.

Pareto front of C-DMOPSO, MOEA/D, dMOPSO, NSGAII, MIMOPSO on UF1, UF7 and UF9, (a)-(f) represent the compare algorithms on UF1, (g)-(l) represent the compare algorithms on UF7, and (m)-(r) represent the compare algorithms on UF9.
To compare the convergence speed of the six algorithm, the trends in the IGD of the algorithms on some test functions with significant differences are depicted in Fig. 6, which indicates that C-DMOPSO gains a good convergence. In these figures, the x-coordinate is the number of function evaluations, while it’s noted that the y-coordinate is log (IGD) for ease of observing expediently. It is discovered from Fig. 6 that C-DMOPSO performs better than other algorithms on most of the test problems during the whole evolutionary process. For ZDT functions, C-DMOPSO gains the fastest convergence rate. DTLZ1 is a multi-modal function, which may weaken the convergence of algorithm. This is the reason for the slow convergence of the algorithms on DTLZ1. For complicated test problems, C-DMOPSO performs well except for UF2.

Variation of IGD metric value.
In conclusion, it is not difficult to find that the solutions obtained by C-DMOPSO receive a good distribution on the Pareto front and a good convergence speed during evolutionary process. Besides, C-DMOPSO does well in handling discontinuous MOPs which benefit from the searching mechanism randomly of PSO. All of these results indicate the improvement of convergence and diversity, which claim great potential of C-DMOPSO to deal with MOPs.
In this paper, an improved decomposition multi-objective particle swarm optimization algorithm based on completion-checking factor is proposed. In the updating process of particle swarm, the velocity update pattern and the size of neighborhood are determined by the completion-checking factor, which change the search pattern of the particle and improve the population convergence speed. In addition, SBX and OBL are added to enhance the global search ability in decision space. Besides, PM is performed on the non-dominated solutions in the archive to improve the diversity. After compared with the other six algorithms, it can be concluded that the C-DMOPSO can obtain better results in most test functions with respect to the IGD performance measure.
C-DMOPSO avoids over-exploration or excessive exploitation and makes a balance between convergence and diversity during evolutionary process. Besides, it combines dominance-based with decomposition-based method effectively. Actually, there is still some challenging and promising research direction on C-DMOPSO. Although the C-DMOPSO algorithm is improved in performance, the setting of parameters cannot be adjusted with the change of environment. As part of our future work, it is necessary to further study the adaptive adjustment of the parameters in C-DMOPSO. Furthermore, C-DMOPSO will be extended to many objectives optimization problems(MaOPs) in our future work.
Footnotes
Acknowledgments
The authors would like to thank the editor and reviewers for their helpful comments and suggestions to improve the quality of this paper. This work was supported by the Natural Science Foundation of Hebei Province (No. F2016203249).
