An Improved Multi-Objective Artificial Bee Colony Optimization Algorithm with Regulation Operators

To achieve effective and accurate optimization for multi-objective optimization problems, a multi-objective artificial bee colony algorithm with regulation operators (RMOABC) inspired by the intelligent foraging behavior of honey bees was proposed in this paper. The proposed algorithm utilizes the Pareto dominance theory and takes advantage of adaptive grid and regulation operator mechanisms. The adaptive grid technique is used to adaptively assess the Pareto front maintained in an external archive and the regulation operator is used to balance the weights of the local search and the global search in the evolution of the algorithm. The performance of RMOABC was evaluated in comparison with other nature inspired algorithms includes NSGA-II and MOEA/D. The experiments results demonstrated that the RMOABC approach has better accuracy and minimal execution time.


Introduction
Most of the optimization problems in the real world are multi-objective optimization problems.Multi-objective optimization is an area of multiple criteria decision making which often have multiple objectives to simultaneously optimize that are conflicting.It has been applied in many fields of science, including engineering, economics and logistics where optimal decisions need to be made in the presence of trade-offs between two or more conflicting objectives [1].If a multi-objective problem is solved by a single-objective optimization method, it can only reflect one aspect of the problem.For example, in the field of hydrological model parameters' optimization, the single-objective optimization algorithm cannot fully reflect the hydrological process and feature [2].The surveys or review articles on ABC variants and their applications also point out the limitations of single-objective optimization algorithms in solving practicable problems [3,4].Therefore, it is necessary to study efficient and stable multi-objective optimization algorithms.Researchers have carried out many interesting studies in the multi-objective optimization field to handle complex nonlinear multi-objective problems, such as the multi-objective evolutionary algorithms (MOEAs), multi-objective PSO (MPSO), and so on.Most of the multi-objective techniques have been designed based on the theories of Pareto Sort [5] and non-dominated solutions.It has been widely used in optimal control, job scheduling, mechanical design, data mining, etc. [6].
Because of the high speed of convergence in single-objective optimization, the ABC (Artificial Bee Colony) algorithm seems particularly suitable for multi-objective optimization.Therefore, in this paper we aim to propose a multi-objective artificial colony (RMOABC) algorithm with regulation operators to solve the multi-objective problems.We utilized the Adaptive Grid [7] technique to estimate the density of solutions for non-dominated solutions in the external archive and adopted Regulation operators that considerably improve the algorithm's exploratory capabilities.Six benchmark functions were selected as the objective functions to compare the performance with two highly competitive multi-objective evolution multi-objective algorithms: the NSGA-II (Non-dominated Sorting Genetic Algorithm II) [8] and the MOEA/D (Multiobjective Evolutionary Algorithm Based on Decomposition) [9].
The paper is organized as follows.Section 2 reviews the related literatures.Section 3 proposes a novel multi-objective Artificial Bee Colony (RMOABC) algorithm.Section 4 represents the optimization algorithms NSGA-II and MOEA/D and analyzes the experimental results.Finally, conclusions are presented in Section 5.

Literature Review
To solve practical problems, many efficient algorithms and methods have been used to design different classes of multi-objective optimization methods for optimizing problems with more than one objective-for example, the multi-objective particle swarm optimization (MOPSO) [5] algorithms and the multi-objective evolutionary algorithms (MOEA) [9] and so on.These algorithms have been widely used and successfully extended to deal with problems with more than one conflicting objective.The Pareto dominance is the most popularized theory for the multi-objective optimizations problems.Based on the concept of Pareto dominance, these Pareto-based approaches select non-dominated solutions to handle the trade-off between the multi-objectives in the problems.
For the multi-objective particle swarm optimization (MOPSO) [7] algorithms, Coello Coello presented a comprehensive survey of evolutionary-optimization based methods [10] and Reyes-Sierra et al. presented a useful survey on the variants of PSO for multi-objective optimization [11].A local search procedure and a flight mechanism both based on crowding distance are incorporated into the MOPSO algorithm and the computational results show that it improve the algorithm with random line search [12].Leong et al. integrated a dynamic population strategy within the multiple-swarm MOPSO and it shows results competitive with improved diversity and convergence [13].
Multi-objective evolutionary algorithms (MOEAs) are also mainstream methods.In [14], a powerful and efficient platform named MOEA Framework is provided.It is a free and open source Java library for developing and experimenting with multiobjective evolutionary algorithms (MOEAs) and other general-purpose multiobjective optimization algorithms.Guliashki et al. have presented a short review of algorithms in the fast growing area of evolutionary multiobjective optimization [15] and Zhou et al. surveyed the development of MOEAs from the algorithmic frameworks, benchmark problems, performance indicators and applications [16].Also, in [17], a new framework is presented for multi-objective Transmission Expansion Planning (TEP) to obtain the final optimal solution, and it has been applied to the real life system of the northeastern part of Iranian national 400 kV transmission grid.
In recent years, multi-objective methods based on the ABC algorithm have gradually become a hot research field.The ABC algorithm's advantages of great accuracy and satisfactory convergence speed make it suitable for the multi-objective optimization problems.Hedayatzadeh et al. designed a multi-objective artificial bee colony (MOABC) based on the Pareto theory and ε-domination notion in [18].The performance of a Pareto-based MOABC algorithm has been investigated by Akbari et al. on CEC'09 data sets in [19] and their studies showed that the algorithm could provide competitive performance.A multi-objective ABC that utilizes the Pareto-dominance concept and maintains the non-dominated solutions in an external archive was presented in [20].Omkar et al. presented a multi-objective technique named Vector Evaluated ABC (VEABC) for optimization of laminated composite components [21].Akbari also designed a multi-objective bee swarm optimization algorithm (MOBSO) and it has the ability to adaptively maintain an external archive of non-dominated solutions [22].Zhang et al. presented a hybrid multi-objective ABC (HMABC) for burdening optimization of copper strip production to solve a two-objective problem of minimizing the total cost of materials and maximizing the amount of waste material thrown into melting furnace [23].A multi-objective variant of ABC (MO-ABC) has been used to solve the frequency assignment problem in GSM networks [24].Mohammadi et al. proposed an Adaptive Multi-Objective Artificial Bee Colony (A-MOABC) Optimizer which utilizes the mechanisms of crowding distance and windowing.The experimental results indicate that the proposed approach performs well [25].Akay proposed three multi-objective Artificial Bee Colony (ABC) algorithms based on synchronous and asynchronous models using Pareto-dominance and non-dominated sorting, and the experiment results showed that synchronous multi-objective ABC using a non-dominated sorting procedure can provide good approximations of well distributed and high quality non-dominated fronts [26].Luo et al. proposed a multi-objective artificial bee colony optimization method called ε-MOABC based on performance indicators to solve multi-objective and many-objective problems [27].Kishor presented a non-dominated sorting based multi-objective artificial bee colony algorithm (NSABC) to solve multi-objective optimization problems [28].Xiang et al. proposed a dynamic, multi-colony, multi-objective artificial bee colony algorithm (DMCMOABC) by using the multi-deme model and a dynamic information exchange strategy [29].Nseef et al. put forward an adaptive multi-population artificial bee colony (ABC) algorithm for dynamic optimization problems (DOPs) [30].The experimental results show that compared with the traditional multi-objective algorithms, these variants of Multi-objective ABC can find solutions with competitive convergence and diversity within a shorter period of time.
In summary, although there are several multi-objective ABC variants, the research into artificial bee colony algorithms for multi-objective are relatively new, and the trade-off between the local search and the global search has not been discussed enough.Therefore, this paper focuses on proposing a new RMOABC algorithm based on the Pareto dominance theory to optimize multi-objective problems.The proposed method is applied on a set of well-known multi-objective problems and compared with other state-of-the-art algorithms.We intended to provide an efficient way to handle the multi-objective problems accurately and quickly.

Multi-Objective Artificial Bee Colony (RMOABC) Algorithm with Regulation Operators
The artificial bee colony (ABC) algorithm is a meta-heuristic intelligent and combinatorial optimization algorithm that was proposed by Karaboga in 2005.It has been inspired by the intelligent behavior of genuine honey bees to solve non-numerical algorithms [31,32].Each individual bee in the algorithm was taken as an agent, which could produce swarm intelligence by dividing the work and corporation among different individuals, conversion of roles, and dancing acts.Currently, ABC algorithm has become an effective method for solving complex nonlinear optimization problems.
In the ABC algorithm, the artificial bee colony is composed of three groups of bees: employed bees, onlookers and scouts.For every food source, there is just one employed bee.Each position of the nectar source represents a possible solution for the optimization problem; the nectar amount of nectar source corresponds to the quality or fitness of the corresponding solution.ABC algorithms gradually converge on optimal or near-optimal solutions through randomly searching or targeted selection of the appropriate individuals of the bee colony.Evolutionary iterations were for the most part jointly done by these three kinds of bees [33]: (1) employed bees perform a local random search in their food source annexed areas; (2) onlookers make the optimum food sources can get more opportunities based on a mechanism; (3) scouts update the stagnant food source.Overall, employed bees and onlookers make the better food sources perform random and objective evolution together.When a solution is in standstill status, scouts will start a new random exploration task.Through collaboration of these three bees, the ABC algorithm will gradually converge.

External Archive (EA)
External archive is a common mechanism of multi-objective algorithms based on Pareto theory to keep a historical record of the non-dominated solutions found along the search process [34].If a solution X 1 presents better fitness values for all the objective functions than another solution X 2 , then X 1 dominates X 2 ; If X 1 is better than X 2 in one objective function and is not worse in the others, then X 1 weakly dominates X 2 [35].
The external archive always has a fixed size, and the RMOABC algorithm has to decide whether a certain solution should be added to the external archive or not in the search process.At the beginning, the archive is empty, and the current solutions will be accepted into it.Then, the non-dominated solutions of the swarm found at each iteration will be compared with the solutions stored in the external archive one by one.If the candidate solution is dominated by an individual within the external archive, then it must be automatically discarded.Otherwise, if the candidate solutions are not dominated by the current solutions within the archive, they must be included in the external archive.If there are solutions in the archive that are dominated by the new solution, then these solutions are removed from the archive.Finally, if the external population has reached its maximum size at the end of the iteration, solutions in more crowded regions are removed from the EA using the adaptive grid mechanism.

Adaptive Grid
The measure of the spread (distribution) of solutions throughout the non-dominated vectors found so far is an important assessment metric for the multi-objective algorithm to judge how well the solutions in Pareto front are distributed.There are several methods for maintaining diversity which consists of a crowding procedure that divides objective space in a recursive manner, such as the adaptive grid [36] and the crowing distance mechanisms [7].
To produce well-distributed Pareto fronts, an adaptive grid method was utilized in this paper.The basic idea of this method is to use an external archive to store all the solutions that are non-dominated with respect to the contents of the archive.According to the non-dominated solutions in the external archive, the objective function space is mapped to a grid, and the grid is divided into many regions.Each solution is placed in a certain location in a grid.Its coordinates or the geographical location were based on the values of its objectives' values calculated by the objective functions.
This grid map can indicate and maintain the number of non-dominated solutions that reside in each grid location.If the individual inserted into the external population lies outside the current bounds of the grid, the grid will be recalculated and each individual within it has to be relocated.Thus the grid procedure is adaptive for this situation.

Regulation Operator
In the iteration of the original ABC algorithm, the current nectar source that needed to be updated is X i = (x i1 , x i2 , ..., x id ).Randomly select the nectar source X j = (x j1 , x j2 , ..., x jd ) that depended on employed bee j in bee colony, and randomly select kth in dimensional space d to calculate the formula (1): where x ij (or v ij ) denotes the jth element of X i (or V i ), and j is a random index.x k denotes another solution selected randomly from the population.And φ ij is a random number in [−1, 1].According to the calculation result, we can get a new nectar source X i = x i1 , x i2 , ..., x ik , ..., x id ) .This processing mechanism will lead to problems from insufficient search depth and the outstanding individuals cannot be fully evolved to affect the convergence performance.In addition, the algorithm also lacks protection for the optimal solution of nectar source, which may result in the abandonment of the best nectar sources, so the algorithm has difficulty converging on the global optimal solution.
To improve the exploitation and take advantage of the global best solution's information to guide the search for candidate solutions, Zhu et al. proposed a Gbest-guided Artificial Bee Colony algorithm (GABC) in [37].The algorithm adopts formula (2) to replace formula (1) in the original ABC algorithm.
where ψ ij is a random number in [0, 1.5], and y j is the optimal solution fitness value in the jth dimensional space.According to the calculation result, a new nectar source X i = x i1 , x i2 , ..., x ik , ..., x id ) could be obtained.In the GABC algorithm, comparison is taken between the new nectar source and the old nectar to retain the nectar source of the optimum quality.
The GABC algorithm looks for the high-quality nectar source through the searching of individual bees and the information sharing between the individual and their neighbors.The difference disturbance between the location of the individuals and their neighbor is affected by the φ ij and ψ ij .Also, the GABC algorithm did not fully take into account the changes in global optimization and local optimization in the evolution process.Thus, we put forward a mechanism of regulation operators into the GABC algorithm for dynamically regulating the weights of global optimization and local optimization in the evolution process to improve the depth search ability and the convergence speed.
In the searching process of our algorithm, to further control the degree of information sharing between bee individuals and the global search precision of algorithm, we introduced two dynamic regulation operators: local dynamic regulation operator k and global dynamic regulation operator r.The core formula based on dynamic regulation operator of GABC algorithm is described in formula (3).
where the local dynamic regulation operator is k = 0.5 + cos 2 ((π × i)/(2 × MFE)), the global dynamic regulation operator is r = 0.5 + sin 2 ((π × i)/(2 × MFE)), i is the current iteration number, and MFE is the maximum iteration number.The r and k should satisfy the equation r + k = 2.As the MFE is 5000, the evolution curves of r and k was shown in Figure 1.The x axis is the iteration number of the algorithm, and the y axis is the values of local regulation operator k and global regulation operator r calculated by the upper equations.
algorithm also lacks protection for the optimal solution of nectar source, which may result in the abandonment of the best nectar sources, so the algorithm has difficulty converging on the global optimal solution.To improve the exploitation and take advantage of the global best solution's information to guide the search for candidate solutions, Zhu et al. proposed a Gbest-guided Artificial Bee Colony algorithm (GABC) in [37].The algorithm adopts formula (2) to replace formula (1) in the original ABC algorithm.
where ij ψ is a random number in [0, 1.5], and yj is the optimal solution fitness value in the jth dimensional space.According to the calculation result, a new nectar source , ,..., ,..., ) could be obtained.In the GABC algorithm, comparison is taken between the new nectar source and the old nectar to retain the nectar source of the optimum quality.
The GABC algorithm looks for the high-quality nectar source through the searching of individual bees and the information sharing between the individual and their neighbors.The difference disturbance between the location of the individuals and their neighbor is affected by the ij φ and ij ψ .Also, the GABC algorithm did not fully take into account the changes in global optimization and local optimization in the evolution process.Thus, we put forward a mechanism of regulation operators into the GABC algorithm for dynamically regulating the weights of global optimization and local optimization in the evolution process to improve the depth search ability and the convergence speed.
In the searching process of our algorithm, to further control the degree of information sharing between bee individuals and the global search precision of algorithm, we introduced two dynamic regulation operators: local dynamic regulation operator k and global dynamic regulation operator r.The core formula based on dynamic regulation operator of GABC algorithm is described in formula (3).
where the local dynamic regulation operator is   It can be seen that in the early stage optimization of the algorithm, the local regulation operator k should larger than global regulation operator r, which means the algorithm is mainly focusing to improve the local search ability.But as the iteration number increases, the local search ability of algorithm should be decreased.Thus, in the latter stage of the algorithm, as the iteration number increases, the global search capacity of algorithm will be enhanced.

Main Algorithm
The proposed RMOABC algorithm is a Pareto-based algorithm.The Pareto ranking scheme could be a straightforward way to extend the original ABC approach to handle multi-objective optimization problems.RMOABC utilizes an external archive of non-dominated solutions to keep the historical record of best solutions found by an individual and incorporates the mechanism of Adaptive Grid mechanism into the ABC algorithm to maintain the solutions' diversity in the external archive.Regulation operators were also integrated into the algorithm to balance the local search and global search of the evolution process.These mechanisms would motivate convergence toward globally non-dominated solutions.The pseudocode of the RMOABC is given in Algorithm 1.This algorithm was divided into four main parts: Initialization, Update bees (SendEmployedBees, SendOnlookerBees, SendScoutBees), update the archive (UpdateExternalArchive, UpdateGrid, TrunkExternalArchive), and Return ExternalArchive.

Experimental Analyses
For the comparative analysis of the performance of RMOABC multi-objective optimization algorithm, we selected two popular and state-of-the-art multi-objective evolutionary algorithms, the NSGA-II (Non-dominated Sorting Genetic Algorithm II) and MOEA/D (Multiobjective Evolutionary Algorithm with Decomposition), for the experiments to compare with the RMOABC algorithm based on six standard benchmark functions taken from the specialized literature [38].
In this paper, all of the optimization algorithms have been implemented on the Java platform.All of our tests have been performed on an Intel (R) Core i7-4720HQ @ 2.6GHz 4 cores with 8 GB of RAM with Microsoft Windows 8 Professional Edition Version.We used the JDK (Java Development Kit) 1.7 and Eclipse 3.2 as the IDE (Integrated Development Environment).NSGA-II is one of the most widely used MOEAs and a significantly improved algorithm proposed by Deb et al. to address the deficiencies exists in the construction of a non-dominated set and the maintenance distribution strategy of the solution set of the NSGA algorithm [8].Through using a rapid non-dominated ranking scheme to reduce computational complexity, introducing elitist strategy to prevent the loss of outstanding optimal solutions, the shortcomings of NSGA have been greatly mitigated.It is shown that NSGA-II performs as well as or better than other algorithms on difficult multi-objective problems.As a non-dominated algorithm that not only has a good distribution and has a faster convergence rate, the NSGA-II algorithm has been widely used in various fields.
In the NSGA-II, the selection process at various stages of the algorithm moves toward a uniformly spread-out Pareto optimal frontier.This is guided by assigning the fitness to chromosomes based on domination and diversity.The NSGA-II algorithm was described in Algorithm 2 [8,39].[9].It decomposes a multiobjective optimization problem into a number of different single objective optimization subproblems (or simple multiobjective optimization subproblems) and then uses a population-based method to optimize these subproblems simultaneously.
Each subproblem is optimized by only using information from its several neighboring subproblems, which makes MOEA/D have lower computational complexity at each generation than MOGLS and NSGA-II.Experimental results have demonstrated that MOEA/D with simple decomposition methods outperforms or performs similarly to MOGLS and NSGA-II on multiobjective 0-1 knapsack problems and continuous multiobjective optimization problems.
The multi-objective evolutionary algorithm based on decomposition (MOEA/D) needs to decompose the MOP under consideration.The authors supposed that the Tchebycheff approach which is very trivial to modify the following MOEA/D when other decomposition methods are used.The MOEA/D algorithm was described in Algorithm 3 [9]: Compute the Euclidean distances between any two weight vectors and then work out the T closest weight vectors to each weight vector.For each i = 1, . . ., N, set B(i) = {i 1 , ..., i T }, where λ i 1 , ..., λ i T are the T closest weight vectors to λ i .
Generate an initial population x 1 , ..., x N randomly or by a problem-specific method.Set FV i = F(x i ).

Begin
For i = 1 to N, Reproduction: Randomly select two indexes k, l from B(i), and then generate a new solution y from x k and x l by using genetic operators.
Improvement: Apply a problem-specific repair/ improvement heuristic on y to produce y'.Update of z: For each j = 1, ..., m, if z j < f j (y ), then set z j = f j (y ).

Update of Neighboring Solutions:
For each index j ∈ B(i), if g te (y λ j , z) , then set x j = y and

Update of EP:
Remove from EP all the vectors dominated by F(y ).Add F(y ) to EP if no vectors in EP dominate F(y ).End For End

Parameter Settings
Appropriate parameterization of the algorithm is important for its efficiency and effectiveness.To compare the optimization algorithms and get better performance, the optimization algorithms' parameter settings in this paper essentially were adopted from previous experiments and the recommended value from the relevant literatures on the subject, NSGA-II [8], and MOEA/D [9].
The parameters' values of the RMOABC were obtained from our previous experiments and research paper [33].The parameter settings of all the optimization algorithms are shown in Table 1.The dimension of decision variables D of the problem to be resolved which is arbitrarily set to 30 to turns the problem's solving process difficult to all algorithms.It will help us to evaluate the behavior of the RMOABC algorithm regarding processing time and quality of solutions.The dimension of the objective function is 2, which means that two objectives were selected in the experiments.All of the termination conditions of the algorithms are reached at the maximum number of iterations, and the maximum number of iterations of algorithms is set to 10,000.

Benchmark Functions
The well-defined standard mathematical benchmark functions can be used as objective functions to measure and test the performance of algorithms.The multi-objective algorithms were applied to four well-known ZDT test functions (ZDT1, ZDT2, ZDT3 and ZDT6) [38] and two multi-objective optimization test instances (UF2 and UF4) for the CEC 2009 special session and competition [40] which were widely used in the literature.All these test instances are minimization of the multiple objectives.

•ZDT1
Its PF is disconnected.The two objectives are disparately scaled in the PF, f 1 is from 0 to 0.852, while f 2 is from −0.773 to 1.

•ZDT6
where g( The distribution of the Pareto solutions in the Pareto front is very non-uniform, i.e., for a set of uniformly distributed points in the Pareto set in the decision space, their images crowd in a corner of the Pareto in the objective space. •UF2 where J 1 = {j|j is odd and 2 ≤ j ≤ n}, J 2 = {j|j is even and 2 ≤ j ≤ n} and where J 1 = {j|j is odd and 2 ≤ j ≤ n}, J 2 = {j|j is even and 2 ≤ j ≤ n} and .

Evaluation Indicators of Multi-Objective Optimization Algorithms
There are many quantitative assessment methods for the performance of multi-objective optimization algorithms.We mainly evaluate the pros and cons of the multi-objective algorithms from the two aspects, convergence and distribution in this paper.Convergence refers to the approaching degree of obtained non-inferior sets with the true Pareto optimal solution sets of this problem, and distribution refers to whether the obtained non-inferior sets are evenly distributed in the optimal space.
(1) Generational Distance (GD) [41]: calculate the distance of the non-dominated solutions obtained in current iteration from the same Pareto frontier.It indicates the distance of the current searched non-dominated vector from the Pareto optimal solution set.
where n is the number of non-inferior solutions obtained by the algorithm; and g i is the minimum Euclidean distance from the i-th solution of algorithm to true Pareto front (or approximately true non-inferior frontier).The smaller the GD value indicates the closer of current solution set of the algorithm to the true Pareto front.
(2) Spacing (SP): It is a uniformity evaluation indicator of optimization algorithm to describe the distribution of solution vectors.Through compared with the true Pareto, coverage degree or distribution of non-dominated solution vectors searched by algorithm could be measured [41]. where .., m, j = i; d is the average of d i value.The smaller the SP value indicates the better distribution of the current solution set of the algorithm.
(3) Computational time (Times): times were evaluated in seconds on the same software and hardware environments for each of the algorithms to demonstrate execution efficiency of the algorithms.

Experimental Results Analysis
The multi-objective optimization algorithms all have random search procedure.Therefore, to restrict the influence of random effects on the experimental results, the simulation experiments have been independently run 30 times.In all the following examples, we report the results obtained from 30 independent runs of each algorithm.In all cases, the best average results obtained with respect to each metric are shown in boldface.(1) Generational Distance (GD) [41]: calculate the distance of the non-dominated solutions obtained in current iteration from the same Pareto frontier.It indicates the distance of the current searched non-dominated vector from the Pareto optimal solution set.
where n is the number of non-inferior solutions obtained by the algorithm; andgi is the minimum Euclidean distance from the i-th solution of algorithm to true Pareto front (or approximately true non-inferior frontier).The smaller the GD value indicates the closer of current solution set of the algorithm to the true Pareto front.
(2) Spacing (SP): It is a uniformity evaluation indicator of optimization algorithm to describe the distribution of solution vectors.Through compared with the true Pareto, coverage degree or distribution of non-dominated solution vectors searched by algorithm could be measured [41]. where ..., , ; .                It can be seen that for the uniformness of final solutions, the RMOABC algorithm is better than the NSGA-II and MOEA/D for ZDT1 and ZDT3 problems.The MOEA/D and RMOABC algorithms are better than the NSGA-II algorithms for ZDT6 problem.For the ZDT2 and UF4 test instance, our RMOABC algorithm has an overwhelmingly better performance than other algorithms.For the UF2 test instance, the three algorithms do not perform well.In particular, the NSGA-II algorithm has the biggest deviation from the Pareto Front.As shown in the figures, the NSGA-II algorithm always produced higher values than the real Pareto front.The MOEA/D algorithm only performs well for ZDT6, but it also has some points with a relatively large deviation from the real Pareto front.Thus, it can be concluded from the above results that our RMOABC algorithm with the regulation operators could converge to the true Pareto front and was able to cover the entire Pareto fronts of all the test functions.

Statistics Results
Tables 2-7 show the numerical statistical results produced by the three multi-objective algorithms in terms of the two performance metrics previously described and their computational time.The average (Mean) and standard deviation (SD) values of Pareto non-dominated sets of 30 independent algorithm runs were calculated.It can be seen that for the average computational time, compared with the other two algorithms, our RMOABC algorithm has the absolute superiority.The algorithm saves 5 to 8 times the execution time and has the smallest standard deviation value.The results also show that the RMOABC method has good performance on ZDT1, ZDT2, ZDT3, ZDT6, UF2 and UF4 test problems for the Generational Distance (GD) metric.With respect to the Spacing (SP) metric, the RMOABC performs well in ZDT1, ZDT2，UF2 and UF4 problems, and has a similar performance to the best algorithm, NSGA-II, in the ZDT6 problem.It can be seen that for the uniformness of final solutions, the RMOABC algorithm is better than the NSGA-II and MOEA/D for ZDT1 and ZDT3 problems.The MOEA/D and RMOABC algorithms are better than the NSGA-II algorithms for ZDT6 problem.For the ZDT2 and UF4 test instance, our RMOABC algorithm has an overwhelmingly better performance than other algorithms.For the UF2 test instance, the three algorithms do not perform well.In particular, the NSGA-II algorithm has the biggest deviation from the Pareto Front.As shown in the figures, the NSGA-II algorithm always produced higher values than the real Pareto front.The MOEA/D algorithm only performs well for ZDT6, but it also has some points with a relatively large deviation from the real Pareto front.Thus, it can be concluded from the above results that our RMOABC algorithm with the regulation operators could converge to the true Pareto front and was able to cover the entire Pareto fronts of all the test functions.

Statistics Results
Tables 2-7 show the numerical statistical results produced by the three multi-objective algorithms in terms of the two performance metrics previously described and their computational time.The average (Mean) and standard deviation (SD) values of Pareto non-dominated sets of 30 independent algorithm runs were calculated.It can be seen that for the average computational time, compared with the other two algorithms, our RMOABC algorithm has the absolute superiority.The algorithm saves 5 to 8 times the execution time and has the smallest standard deviation value.The results also show that the RMOABC method has good performance on ZDT1, ZDT2, ZDT3, ZDT6, UF2 and UF4 test problems for the Generational Distance (GD) metric.With respect to the Spacing (SP) metric, the RMOABC performs well in ZDT1, ZDT2, UF2 and UF4 problems, and has a similar performance to the best algorithm, NSGA-II, in the ZDT6 problem.The results shows that using the mechanisms of adaptive grid, external archive and the regulation operators not only provides a well-distributed set of non-dominated solutions, but also helps in the convergence of the algorithm to the true Pareto front.The RMOABC has a significantly lower computational time than the NSGA-II and MOEA/D algorithms, which is attributed to the adaptive grid mechanism as it can be computed faster than the crowding distance.

Conclusions
In this paper, we presented a multi-objective artificial bee colony algorithm with regulation operators (RMOABC) based on the ABC algorithm to handle multi-objective optimization problems.The proposed algorithms were implemented based on the Pareto dominance theory and incorporated an adaptive grid to control the proximity of the external archive solutions and regulation operators to maintain a trade-off between the local search and the global search.The performance of the RMOABC is evaluated based on six test functions and two metrics from the literature.The experimental study showed that the presented algorithm is highly competitive compared to the other two state-of-the-art multi-objective algorithms investigated in this work.The algorithm performs well in converging towards the Pareto front and has generated a well-distributed set of non-dominated solutions.Additionally, the low computational times required by our algorithm make it a very promising approach for problems in engineering optimization fields.

,
i is the current iteration number, and MFE is the maximum iteration number.The r and k should satisfy the equation 2 r k + = .As the MFE is 5000, the evolution curves of r and k was shown in Figure 1.The x axis is the iteration number of the algorithm, and the y axis is the values of local regulation operator k and global regulation operator r calculated by the upper equations.

Algorithm 3 .
MOEA/D Algorithm Input: • multiobjective optimization problem (MOP); • a stopping criterion; • N: The number of the subproblems considered in MOEA/D; • a uniform spread of N weight vectors: λ 1 , ..., λ N ; • T: The number of the weight vectors in the neighborhood of each weight vector.Output: an external population (EP), which is used to store nondominated solutions found during the search.Begin Initialization(); While (Stopping Criteria is not satisfied) Update(); End While Stop and output EP.End Function Initialization() Begin Set EP = Φ.

4. 3 . 1 .
Figures 2-5 show the graphical results produced by the NSGA-II, MOEA/D and RMOABC algorithms for the four ZDT test functions.Figures 6 and 7 show the graphical results produced by the NSGA-II, MOEA/D and RMOABC algorithms for the UF2 and UF4 test functions.The true Pareto front of the problem is shown as black dots, and the non-dominated fronts produced by the three multi-objective algorithms were represented by red crosses.The upper left panel of figures is for NSGA-II algorithm, the upper right panel is for the MOEA/D algorithm, the lower left panel is for the RMOABC algorithm.

Figure 3 .
Figure 3. Plots of the non-dominated front found by the three algorithms for 2-objective ZDT2 test instance.(a) NSGA-II algorithm; (b) MOEA/D algorithm; (c) RMOABC algorithm.Figure 3. Plots of the non-dominated front found by the three algorithms for 2-objective ZDT2 test instance.(a) NSGA-II algorithm; (b) MOEA/D algorithm; (c) RMOABC algorithm.

Figure 3 .
Figure 3. Plots of the non-dominated front found by the three algorithms for 2-objective ZDT2 test instance.(a) NSGA-II algorithm; (b) MOEA/D algorithm; (c) RMOABC algorithm.Figure 3. Plots of the non-dominated front found by the three algorithms for 2-objective ZDT2 test instance.(a) NSGA-II algorithm; (b) MOEA/D algorithm; (c) RMOABC algorithm.
For j = 0 to FoodNumber−1 and isDominated[i] is not true If i is not equal j and isDominated[j] is not true Verify the dominance of solution i and solution j If solution j dominates solution i Select a parameter p randomly Select neighbor k from the solution randomly for producing a mutant solution Generate a new candidate solution in the range of the parameter Calculate the fitness according to the objective functions If the new solution dominates the old one Begin Set maxValues [] = the max Fitness values of each objective of solutions in ExternArchive Set minValues [] = the min Fitness Values of each objective of solutions in ExternArchive For i = 0 to OBJNum−1 Initialize the grid's ranges to the maxValues[i]−minValues[i]

Table 1 .
The parameter settings of the three multi-objective algorithms.

Table 2 .
Results of the performance metrics for ZDT1 function of the algorithms.

Table 3 .
Results of the performance metrics for ZDT2 function of the algorithms.

Table 4 .
Results of the performance metrics for ZDT3 function of the algorithms.

Table 5 .
Results of the performance metrics for ZDT6 function of the algorithms.

Table 6 .
Results of the performance metrics for UF2 function of the algorithms.

Table 7 .
Results of the performance metrics for UF4 function of the algorithms.