Improved_Local_Search_in_Artificial_Bee_Colony_usi
Improved_Local_Search_in_Artificial_Bee_Colony_usi
net/publication/232608310
Improved Local Search in Artificial Bee Colony using Golden Section Search
CITATIONS READS
41 270
3 authors, including:
All content following this page was uploaded by Dr. Tarun Kumar Sharma on 15 October 2014.
Abstract –Artificial bee colony (ABC), an optimization algorithm is a recent addition to the family of population based
search algorithm. ABC has taken its inspiration from the collective intelligent foraging behavior of honey bees. In this
study we have incorporated golden section search mechanism in the structure of basic ABC to improve the global
convergence and prevent to stick on a local solution. The proposed variant is termed as ILS-ABC. Comparative
numerical results with the state-of-art algorithms show the performance of the proposal when applied to the set of
unconstrained engineering design problems. The simulated results show that the proposed variant can be successfully
applied to solve real life problems.
Keywords –Artificial Bee Colony; ABC; Optimization; Golden Section; Metaheuristics;
food sources randomly based on some external clue [30]. randomly to SN number of employed bees and their
This intelligent foraging mechanism can be better fitness is determined.
explained using the graphical representation shown in
Employed Bees Stage
figure 1.
The duty of employed bees is to determine new food
source say, vi with the help of the food source xi assigned
to it during the initialization phase. The equation used is:
vij x ij ij ( xij x kj ) (2)
In order to improve the performance of basic ABC, The performance of ILS-ABC is evaluated on test bed
golden section search [31] is introduced in the structure of four engineering design real life problems taken from
of basic ABC. This process helps in balancing the literature.
exploration and exploitation. The poor balance may result The performance of proposed algorithms is evaluated
in weak optimization method which may suffer from on a test bed of seven unconstrained real life problems
trapping in a local optima, stagnation and premature that are common in various fields of engineering designs.
convergence. Introduction of golden section search helps These are: (A) Optimization of Transistor Modeling (B)
in improving the global convergence and to prevent to Optimal capacity of gas production facility (C) optimal
stick on a local food source of the ABC. Hence the thermohydraulic performance of an artificially roughened
proposed variant is termed as ILS-ABC. In the proposed air heater and (D) Design of a gear train. All the problems
algorithm the onlooker bee replaces the food sources in considered in the present study are highly nonlinear in
basic ABC (Equation 2) using the following Equation: nature. Mathematical models of the problems are
vij x ij Fz ( xij x kj ) (4) discussed below:
where Fz is scaling factor depends on the value of a & b (A) Optimization of transistor modeling [32]: The
discussed below. The generation of a new food sources is objective function of this problem provides a least-
seen as an opaque procedure depending on the scale sum-of-squares approach to the solution of a set of
factor Fz. nine simultaneous nonlinear equations, which arise
Scaling factor becomes an important aspect that is to in the context of transistor modeling. The
be controlled in order to guarantee a high quality food mathematical model of the transistor design is given
source (solution) that can have an important role in the by,
4 2 2
succeeding generations. Minimize f ( x) 2 ( k k )
k 1
A local search procedure is applied to Fz in order to
Where
detect the scale factor value which guarantees a food
k (1 x1 x 2 ) x 3 {exp[ x5 ( g1k g 3k x7 10 3 g 5k x8 10 3 )] 1}g 5 k g 4 k x 2
sources with a high performance. This can be gained by
k (1 x1 x 2 ) x4 {exp[ x6 ( g 1k g 2 k g 3 k x 7 10 3 g 4 k x9 10 3 )] 1}g 5 k x1 g 4 k
minimizing the function f ( Fz ) in the decision space
x1x3 x2 x4
[−1, 1]. The meaning of the negative scale factor is the
inversion of the search direction. If the search in the Subject to: xi 0 , i = 1, 2, … , 9
negative direction succeeds, the corresponding positive And the numerical constants g ik are given by the matrix
value is assigned to the food sources for the subsequent
0.485 0.752 0.869 0.982
generations.
The Scale factor golden section search applies the 0 .369 1 . 254 0 .703 1.455
Golden section search to the scale factor in order to 5.2095 10.0677 22.9274 20.2153
generate a high quality food sources. This scheme
23.3037 101.779 111.461 191.267
processes in the interval [a = −1, b = 1] and generates two 28.5132 111.8467 134.3884 211.4823
intermediate points:
ba ba
Fz1 b ; Fz2 a (B) Optimal capacity of gas production facilities [33]:
This is the problem of determining the optimum
1 5 capacity of production facilities that combine to
where is the golden ratio. The evaluated make an oxygen producing and storing system. The
2 mathematical model of this problem is given by,
values of f ( Fz1 ) and f ( Fz2 ) are then compared and if x2 0.85
f ( x) 61.8 5.72 x1 0.2623[(40 x1 ) ln( )]
f ( Fz1 ) < f ( Fz2 ) then Fz2 replaces b and this procedure is Minimize 200
repeated in the new smaller interval [a, b]. The process is x2 0.75
0.087(40 x1 ) ln( ) 700.23x2
described as below: 200
Where R M
0.95 x 2
0.53
; GH 4.5(e ) 0.28 (0.7) 0.57 , implies the mean of the best fitness values obtained by
e x1x3 ( f / 2) 1/ 2
f ( f s f r ) / 2 f s 0.079x3 0.25 the algorithm till the stopping criteria is satisfied,
; ; ; Number of function evaluations (NFE) tells about the
0.53
f r 2(0.95 x3 2.5 * ln(1 / 2 x1 ) 2 3.75) 2 speed of an algorithm. Smaller NFE shows that the
;
algorithm converges more quickly. Numerical results
Subject to: 0.02 x1 0.8 , 10 x 2 40 , 3000 x3 20000 based on these performance measures are given in Table
(D) Design of gear train [32]: This problem is to 1 and Table 2. From Table 1, which gives the average
optimize the gear ratio for the compound gear train. fitness function value, NFE and execution time, we see
It is to be designed such that the gear ratio is as close that in terms of average fitness function value and time
as possible to 1/6.931. For each gear the number of all the algorithms gave more or less similar results
teeth must be between 12 and 60. Since the number although in some cases the proposed algorithms gave a
of teeth is to be an integer, the variables must be marginally better performance than PSO, DE and ABC.
integers. The mathematical model of gear train However if we compare the NFE in Table 2 the superior
design is given by, performance of the proposed algorithms become more
evident.
2 2
1 T T 1 xx Acceleration rate (AR) [35] is used to compare the
Minimize f ( x) d b 1 2 convergence speeds between ILS-ABC and other
6. 931 T T
a f 6. 931 x x
3 4
algorithms. It is defined as follows:
Subject to: 12 xi 60 i 1, 2,3, 4 NFE one a lg orithm NFE other a lg orithm
AR %
[ x1 , x 2 , x3 , x 4 ] [Td , Tb , Ta ,T f ] , xi’s should be integers. Ta, NFE one a lg orithm
Tb , Td, and Tf are the number of teeth on gears A, B, D From Table 2 we can see that the proposed ILS-ABC
and F respectively. gives the better results for every problem in the
comparison to the PSO, DE and ABC. Further from the
Table 2 it is clearly analysed that the proposed ILS-ABC
is faster than PSO by 38.6%, DE by 27.7% and faster
than ABC by 31.6%.
Table 1. Simulation results of four engineering design real life
problems
Figure 2. Design of Gear Train
Source
Item PSO DE ABC ILS-ABC
Result
5. Experimental settings, Result analyses and (A) Optimization of Transistor Modeling
x1 0.9010 0.9013 0.9009 0.9002 0.90
Comparison criterion x2 0.8841 0.5878 0.5224 0.4519 0.45
x3 4.0386 3.7812 1.0764 1.0352 1.0
5.1. Experimental Settings x4 4.1488 3.9021 1.9494 2.045 2.0
x5 5.2436 5.1962 7.8536 5.495 8.0
x6 9.9326 11.2697 8.8364 9.734 8.0
ABC and the proposed variant ILS-ABC are x7 0.1009 0.0979 4.7712 0.103 5.0
implemented on Dev-C++ and the experiments are x8 1.0599 1.1053 1.0074 1.012 1.0
conducted on a computer with 2.00 GHz Intel (R) core x9 0.8066 0.6799 1.8545 0.610 2.0
(TM) 2 duo CPU and 2- GB of RAM. For each problem, f(x) 0.0695 0.0618 0.0113 0.0304 NA
NFE 22195 19784 17901 14932 NA
all the algorithms independently run 25 times. Time 0.863 0.415 0.401 0.345 NA
The population size, MCN (maximum cycle numbers) (B) Optimal Capacity of Gas Production Facilities
and limit are fixed to 40, 10000 and 100 respectively. x1 17.5 17.5 17.5 17.5 17.5
And the reported results are the means and standard x2 600 593 600 600 465
deviations of the statistical experimental data. The f(x) 169.844 169.996 169.012 169.731 173.76
NFE 342 324 319 264 NA
random numbers are generated using inbuilt rand () Time 0.02 0.01 0.02 0.01 NA
function with same seed for every algorithm. (C) Optimal Thermo hydraulic Performance of an Artificially
For each algorithm, the stopping criteria is to Roughened Air Heater
terminate the search process when one of the following x1 0.05809 0.18005 0.029589 0.07111 0.052
conditions is satisfied: (i) the maximum number of x2 10 10 9.012 9.876 10
x3 10400.2 4918 5028.4 4220 10258
generations is reached (assumed 1000 generations), f(x) 4.2142 4.0923 4.2166 4.0703 4.182
(ii) | fmax - fmin | < 10-4 where f(x) is the value of objective NFE 6207 5190 4425 3026 NA
function. Time 0.3 0.3 0.3 0.2 NA
(D) Design of Gear Train
5.2. Result Analyses and Comparison Criteria x1 13 16 19 18.75 18
x2 31 19 16 13.6 22
x3 57 49 44 44 45
The performances of proposed ILS-ABC is x4 49 51 49 52 60
demonstrated on a set of four nonlinear engineering f(x) 9.98e-11 2.78e-11 2.78e-11 5.59e-07 5.7e-06
Gear ratio 0.14429 0.14428 0.14428 0.14424 0.14666
design problems and the numerical results are compared
Error (%) 0.007398 0.000497 0.000447 0.000401 1.65
with the basic PSO, DE, and ABC. Standard performance NFE 480 340 648 252 NA
measures like average fitness function values, average Time 0.1 0.01 0.02 0.01 NA
number of function evaluations and time taken to
simulate are calculated. Average fitness function value
5
[26] T.K. Sharma, M. Pant, J.C. Bansal, Artificial Bee Colony with Conference on Bio-Inspired Computing: Theories
Mean Mutation Operator for Better Exploitation, in: IEEE World
Congress on Computational Intelligence (CEC), Brisbane,
and Application (BIC-TA 2012) and SoCProS-2012
Australia, 2012, pp. 3050 - 3056. (International Conference on Soft Computing for
[27] T.K. Sharma, M. Pant, J.C. Bansal, Some Modifications to Problem Solving).
Enhance the Performance of Artificial Bee Colony, in: IEEE
World Congress on Computational Intelligence (CEC), Brisbane, V.P. Singh is Director, SCET Saharanpur, and retired
Australia, 2012, pp. 3454 - 3461.
Professor from the department of Paper Technology, IIT
[28] D. Karaboga, B. Gorkemli, C. Ozturk, N. Karaboga: A
comprehensive survey: artificial bee colony (ABC) algorithm and
Roorkee. His research area includes Mathematical
applications, Artif Intell Rev 2011, DOI 10.1007/s10462-012- Simulation and Modeling. He has number of publications
9328-0. in journal of repute.
[29] E. Bonabeau, M. Dorigo G. Theraulaz, Swarm Intelligence: From
Natural to Artificial Systems, Oxford University Press, Inc., New
York, NY, USA, 1999.
[30] T.D. Seeley, The Wisdom of the Hive: The Social Physiology of
Honey Bee Colonies, Harvard University Press, 1995.
[31] J. Kiefer, “Sequential minimax search for a maximum, in
Proceedings of the American Mathematical Society, 4, 1953, pp.
502–506.
[32] B. Babu, New optimization techniques in engineering. Springer-
Verlag, Berlin Heidelberg, 2004.
[33] C. Beightler, D. Phillips. Applied geometric programming. John
Wiley and Sons, New York, 1976.
[34] B. Prasad, J. Saini. Optimal thermo hydraulic performance of
artificially roughened solar air heater, Journal Solar Energy,
47(1991) 91–96.
[35] S. Rahnamayan, H.R. Tizhoosh, and Salama, M. M. Ali,
Opposition-Based Differential Evolution. IEEE Transactions on
Evolutionary Computation, 12:1(2008) 64– 79.