A Multiobjective Optimization Algorithm Based on Discrete Bacterial Colony Chemotaxis

Bacterial colony chemotaxis algorithm was originally developed for optimal problem with continuous space. In this paper the discrete bacterial colony chemotaxis (DBCC) algorithm is developed to solve multiobjective optimization problems. The basic DBCC algorithm has the disadvantage of being trapped into the local minimum. Therefore, some improvements are adopted in the new algorithm, such as adding chaos transfer mechanism when the bacterium choose their next locations and the crowding distance operation to maintain the population diversity in the Pareto Front. The definition of chaos transfer mechanism is used to retain the elite solution produced during the operation, and the definition of crowding distance is used to guide the bacteria for determinate variation, thus enabling the algorithm obtain well-distributed solution in the Pareto optimal set. The convergence properties of the DBCC strategy are tested on some test functions. At last, some numerical results are given to demonstrate the effectiveness of the results obtained by the new algorithm.


Introduction
In the field of optimization, many researchers have been inspired by biological processes such as evolution [1,2] or the food-searching behavior of ants [3] to develop new optimal methods such as evolutionary algorithms or ant codes.These techniques have been found to perform better than the classical heuristic or gradient-based optimization methods, especially when faced with the problem of optimizing multimodal, nondifferentiable, or discontinuous functions.Examples of applying these biologically inspired strategies in the field of engineering range from aerodynamic design [4] to job-shop scheduling problems [5].Another biologically inspired optimization method is the chemotaxis algorithm, pioneered by Bremermann [6] and his coworkers [7,8], that is proposed by analogy to the way bacteria react to chemoattractants in concentration gradients.This algorithm was analyzed for chemotaxis in a three-dimensional (3D) gradient by Bremermann [6] and employed for the training of neural networks [7,8].This strategy is evaluated on a number of test functions for local and global optimization, compared with other optimization techniques, and applied to the problem of inverse airfoil design [9].A similar algorithm is the guided accelerated random search technique [10], which was applied to optimize parameters in flight control systems, [11] and to optimize perceptrons [12].The BCC algorithm is a new colony intelligence optimization algorithm which was introduced in [13].This novel algorithm considers not only the chemotactical strategy but also the communication between the colony members, and the performance has improved greatly.
Most real-world optimization problems require making decisions involving two or more goals.When these goals are the minimization or maximization of the objective functions they are typically referred to as multiobjective optimization (MO).From the 1950s, a variety of methods known as classical have been developed for the solution of multiobjective optimal problems (MOP).These methods are based on formal logic or mathematical programming [14].Another interesting biological process that has been already implemented as an optimization technique is the BCC algorithm.This novel technique exposed the potential of implementing bacterial chemotaxis as a distributed optimization process, recognizing that in natural colonies, it is the interaction and communication between bacteria the mechanism that enables them to develop biologically advantageous patterns.Many real-world binary-discrete optimization problems require making decisions involving two or more goals that typically are in contradiction with each other.So in this paper the DBCC algorithm is developed to solve MOP.The basic DBCC algorithm has the disadvantage of being trapped into the local minimum and also cannot maintain the population diversity in the Pareto optimal set.Therefore, some improvements are adopted in the new algorithm, as adding arithmetic operator when the bacteria choose their next locations and taking into account the chaos transfer mechanism.
This paper is organized as follows.Section 2 describes the basic theorem of MOP.The model is shown in Section 3. The convergence properties of the developed DBCC strategy which is applied on test functions are given by comparing with other algorithms in Section 4. Finally, some conclusions are drawn in Section 5.

Preparative Theorem of MOP
A MOP [15] is defined as the problem of finding a vector of decision variables that satisfies some restrictions and optimizes a vector function, whose elements represent the values of the functions.A MOP may be formulated as follows [16]: where  is the vector of discrete decision variables  = ( 1 ,  2 , . . .,   )  and   () are the  objective functions.
The  inequalities   and the  equalities ℎ  are known as constraint functions.
For MOP, instead of a single optimal, there is a set of optimal solutions known as Pareto optimal front (POF).Any solution of this set represents a balance between the objective functions; therefore, it is not possible to say that there is other solution in the search space which is superior to this one when all objectives are considered.In the minimization MOP, Pareto optimality can be mathematically defined as follows [17].Definition 1. Pareto dominate: given two candidate solutions  and  from , where  is the feasible region, the vector () is said to dominate vector () (denoted by ≺y) if and only if, Denoted by ¬ ≺  if and only if, The set of solutions that satisfy all constraints is known as the Pareto optimal set and the fitness values corresponding to these solutions form the Pareto front or trade-off surface in objective space.In most cases it is not easy to find analytical expressions for the line or curve that contains the POF; thus, commonly optimal solutions points and the objective functions values in them are calculated.In order to find optimal solutions, there are two goals that any multiobjective optimization algorithm (MOA) seeks to achieve [18]: (1) guide the search towards the global Pareto optimal region, (2) maintain the population diversity in the Pareto optimal front.
Generate  bacteria randomly by using binary-discrete code and make it turn to real code.Compute the velocity V.In this model, the velocity is assumed to be a scalar constant value V = 1.Evaluate the fitness of the bacteria and store the nondominated solutions in the Pareto optimal set.
Step 2. Individual bacterial optimization and find the new position ⃗  1 .
(a) Compute the duration of the trajectory  from the distribution of a random variable with an exponential probability density function The time  is given by where  0 is minimal mean time, ⃗  pre and ⃗  cur are the position of a bacterium in the previous and current step, respectively,   pr  = 1, 2, . . .,  is the difference between the actual and the previous function value,  pr is vector connecting the previous and the actual position in the parameter space, and  is dimensionless parameter.
The angle   between the previous and the new direction obeys Gaussian distribution, for turning left or right, respectively, where   ∈ [0 0 , 180 0 ], the expectation value  = (  ), and the variance   = √VAR(  ) with where  pr is the duration of the previous step.(c) Compute the position ⃗  1 of the bacterium.The length of the path  is given by The normalized new direction vector ⃗   with | ⃗   | = 1 is multiplied by  to obtain the displacement vector ⃗  1 such that the new location of the bacterium is Step 3. Optimize by the bacterial colony and find the new position ⃗  2 .
The individual bacterial acquires information about its environment; compute the best neighbor center ⃗   .We can get the new position by where ⃗  = (0, 2) ⋅ dis( ⃗   , Center( ⃗   )), where (0, 2) is the random number within the interval (0, 2) and dis(, ) is the distance connecting the vector  and .Step 4. Optimize by the crowding distance operator and find the new position ⃗  3 .
The new position of the individual bacterial with crowding distance operator is given by where  ∈ [−1.2, 1.2] is assumed to be a real number and ⃗  pr is the biggest crowding distance position of bacterium in the previous step.The count of ⃗  3 is given no more than  * 5% + 2.
Step 5.The reference colony ⃗  ref is produced and bacteria move to the new location [19].
Step 6.If the chosen target precision  end is not reached, then go to Step 2; otherwise end.In practical optimization, the calculation precision  can be adjusted adaptively with  next =  begin /, where  > 1 is the precision update constant.
As the initial precision is reached, the parameters are adapted to another precision (defined by the number of parameter changes) and the search continues, until this new precision is reached.A given precision  end is reached if the difference between function values found by the bacterium is smaller than  end        pr      <  end  = 1, 2, . . ., for a given number  end .
When the        pr      <    = 1, 2, . . ., , (18) the difference between the actual and the previous function value is reached in a row   that has assumed before.The bacterial colony has to move to other positions and neglects information that the bacterial colony has held in previous step.This has been called bacterial migration.We modify the algorithm with chaos transfer mechanism (CTM).Therefore, one has the following: So the algorithm can escape from the local optimal region.In the other side, sometimes the bacterium jumps out of the area of the global optimal and is never able to get back.To prevent the bacteria strategy from leaving the global optimal regions, the algorithm will search each of the bacterial migration positions, compare the function values with that in the POF, and then take the better one into the POF.Flowchart of the MOADBCC algorithm is given as shown in Figure 1.

Tests Results of the MOADBCC
MOADBCC with operators, MOADBCC without operators, and NSGA-II are, respectively, applied on simulation, among which NSGA-II applies binary coding to obtain Pareto Figure 2: The result of MOADBCC with and without operators on SCH.
optimal solution set.The parameter selection of MOADBCC algorithm can refer to [6].The initial precision is set as 2.0 with a precision update constant  = 1.05.The cross probability of NSGA-II is 0.9, and the variation probability is 0.01.The colony scales of all these three algorithms are 100 and the test function iterates 100 times except for the ZDT series function which iterates 200 times.Test condition was AMD Athlon(tm) IIX2 255 Processor, 2 GB Memory, Microsoft Windows XP operation system; the procedure was tested in Matlab 7.1.This paper applies five test functions to examine the performance of the algorithm, including the research achievement SCH of Schaffer [20], DErB, and ZDT1 and ZDT2 [21][22][23] of Deb and SRN [24,25] of test function Srinivas with constraint condition,.Binary coding and discretization treatment are applied by test function while the optimization is conducted in the discrete area.Among them, ZDT1 and ZDT2 and DEB apply 10 coding and SRN and SCH apply 20 coding.Table 1 shows test functions.
In order to evaluate the comprehensive performance of the algorithm, CM, and RVI performance evaluation indexes are applied in this paper to measure the advantages of the algorithm performance.
(1) Convergence index (CM) [26] is used to measure the intimacy level between the solutions obtained by the algorithm and real solution set.The smaller the CM index value is, the more intimate the obtained solutions are with real POF (2) Effective solution proportion index (RVI) refers to the proportion of nondominant solutions in all solutions, which is used to examine the convergence degree of final solutions.The bigger the RVI index value is, the better performance the algorithm is where  means that the nondominant solutions minus the bacteria amount in same location.The algorithm operates for 30 times independently for each index to obtain their mean value Mean and standard deviation Stdev.Tables 2   MODBCC and NSGA-II is equal to each other in terms of RVI index.However, the performance of MOADBCC is better than NSGA-II in terms of CM index, indicating it has better performance than NSGA-II algorithm.It is due to the fact that the operators proposed by this paper make it easier for bacteria to jump out of local convergence and accelerate the convergence speed of the algorithm, thus enabling it to obtain more effective solutions with same iterations and increasing the value of RVI.
In order to further compare the performances among these three algorithms, this paper listed the Figures of simulation results

Conclusions
In this paper, a discrete bacterial colony chemotaxis (DBCC) algorithm is developed by chaos transfer mechanism and crowding distance strategy and used to solve multiobjective optimization problems.An additional term containing the chaos transfer mechanism draws the search away from local optima.Incorporation of crowding distance strategy also improves performance.The performance on tests is given to demonstrate the validity of the proposed algorithm.In fact, the results given in this paper not only show the effectiveness of the DBCC for the test functions, but also generally support the use of DBCC for other practical problems.Hence, there are future works including parameter settings and more applications of the proposed algorithmic framework to other optimization problems.

Figure 3 :Figure 4 :Figure 5 :Figure 6 :
Figure 3: The result of MOADBCC with and without operators on DEB.
of three algorithms corresponding to SCH, DEB, and FON test functions, shown in Figures 2, 3, 4, 5, 6, and 7. From the simulation Figures, those results show that MOADBCC with operators is better than the other two algorithms in view of test problems in the Pareto Front, as it can converge to the Pareto Front more effectively and obtain more well-distributed solutions.

Table 2 :
Comparison table of CM.

Table 3 :
Comparison table of RVI.
to 3 show the mean and Stdev of CM and RVI index for MOADBCC with operators, MOADBCC without operators, and NSGA-II respectively.From the table of comparisons, the results show that all indexes of MOADBCC with operators are better than that of MOADBCC without operators.The performance of