arxiv: v3 [cs.ne] 5 Jun 2013

Transcription

1 A Hybrd Bat Algorthm Iztok Fster Jr., Dušan Fster, and Xn-She Yang Abstract Swarm ntellgence s a very powerful technque to be used for optmzaton purposes. In ths paper we present a new swarm ntellgence algorthm, based on the bat algorthm. The Bat algorthm s hybrdzed wth dfferental evoluton strateges. Besdes showng very promsng arxv: v3 [cs.ne] 5 Jun 2013 results of the standard benchmark functons, ths hybrdzaton also sgnfcantly mproves the orgnal bat algorthm. To cte paper as follows: I. Jr. Fster, D. Fster, X.-S Yang. A hybrd bat algorthm. Elektrotehnšk vestnk, 2013, n press. Unversty of Marbor, Faculty of electrcal engneerng and computer scence Smetanova 17, 2000 Marbor; Electronc address: ztok.fster@guest.arnes.s Unversty of Marbor, Faculty of mechancal engneerng Smetanova 17, 2000 Marbor; Electronc address: dusan.fster@un-mb.s Mddlesex Unversty, School of Scence and Technology London NW4 4BT, UK; Electronc address: x.yang@mdx.ac.uk 1

2 I. INTRODUCTION Nature has always been an nspraton for researchers. In the past, many new naturenspred algorthms have been developed to solve hard problems n optmzaton. In general, there are two man concepts developed n bo-nspred computaton: 1. evolutonary algorthms, 2. swarm ntellgence algorthms. The Evolutonary algorthms are optmzaton technques [7] that base on the Darwn s prncple of survval of the fttest [5]. It states that n nature the fttest ndvduals have greater chances to survve. The Evolutonary algorthms consst of the followng dscplnes: genetc algorthms, evoluton strateges, genetc programmng, evolutonary programmng, dfferental evoluton. Although all these algorthms or methods have been developed ndependently, they share smlar characterstcs (lke varaton operators, selecton operators) when solvng problems. In fact, the evolutonary algorthms are dstngushed by ther representaton of solutons. For example, the genetc algorthms [15, 16] support the bnary representaton of solutons, evoluton strateges [2, 11] and dfferental evoluton [3, 6, 29] work on real-valued solutons, genetc programmng [21] acts on programs n Lsp, whle the evolutonary programmng [13] behaves wth the fnsh state automata. The Evolutonary algorthms have been appled to a wde range of areas of optmzaton, modelng and smulaton. Essentally, the dfferental evoluton has successfully been employed n the followng areas of optmzaton: functon optmzaton [28], large-scale global optmzaton [4], graph colorng [8], chemcal-process optmzaton [1]. Swarm ntellgence s a collectve behavour of decentralzed, self-organzed systems, ether natural or artfcal. Swarm ntellgence was ntroduced by Ben n A lot of algorthms have been proposed snce then. The Swarm-ntellgence algorthms have been appled on contnuous as well as combnatoral optmzaton problems [25]. The most well-known classes of the swarm-ntellgence algorthms are: partcle swarm optmzaton, ant colony optmzaton, artfcal bee colony, frefly algorthm, cuckoo search and bat algorthm. Partcle swarm optmzaton has been successfully appled n problems of the antenna desgn [17] and electromagnetcs [27]. The Ant colony algorthms have also been used n 2

3 many areas of optmzaton [20] [26] [22]. The Artfcal bee colony shows good performance n numercal optmzaton [18] [19], n large-scale global optmzaton [10] and also n combnatoral optmzaton [24] [9] [30]. The Cuckoo search algorthm s a very strong method for functon optmzaton and also for engneerng optmzaton problems [34] [33]. The Frefly algorthm shows promsng results n functon optmzaton and provdes good results also n combnatoral optmzaton [12]. Echolocaton s an mportant feature of bat behavour. That means that bats emt a sound pulse and lsten to the echo bouncng back from obstacles whlst flyng. Ths phenomenon has nspred Yang [36] to develop the Bat Algorthm (BA). The algorthm shows good results when dealng wth lower-dmensonal optmzaton problems, but may become problematc for hgher-dmensonal problems because of ts tendng to converge very fast ntally. On the other hand, the dfferental evoluton [23] s a typcal evolutonary algorthm wth dfferental mutaton, crossover and selecton that s beng successfully appled to contnuous functon optmzaton. To mprove the bat-algorthm behavour for hgher-dmensonal problems, we hybrdzed the orgnal bat algorthm by usng dfferental-evoluton strateges. Ths Hybrdzed Bat Algorthm (HBA) was tested on a standard set of benchmark functons taken from lterature. Results of our numercal experments show that the proposed HBA can sgnfcantly mprove the performance of the orgnal BA, whch can be very useful for the future. The structure of the paper s as follows. In Secton 2, we ntroduced the orgnal BA and the dfferental evoluton algorthm and explan some bologcal foundatons of the bat behavour. In Secton 3 we descrbe our novel approach to hybrdzng the bat algorthm wth dfferental evoluton strateges. In Secton 4 we llustrate our experments and dscuss the obtaned results. The paper ends by presentng our plans for future work on HBA development. II. BAT ALGORITHM The Bat algorthm was developed by Xn-She Yang n 2010 [35]. The algorthm explots the so-called echolocaton of the bats. The bats use sonar echoes to detect and avod obstacles. It s generally known that sound pulses are transformed nto a frequency whch reflects 3

4 from obstacles. The bats navgate by usng the tme delay from emsson to reflecton. They typcally emt short, loud sound mpulses. The pulse rate s usually defned as 10 to 20 tmes per second. After httng and reflectng, the bats transform ther own pulse nto useful nformaton to gauge how far away the prey s. The bats are usng wavelengths that vary n the range from 0.7 to 17 mm or nbound frequences of khz. To mplement the algorthm, the pulse frequency and the rate have to be defned. The pulse rate can be smply determned n the range from 0 to 1, where 0 means that there s no emsson and 1 means that the bats emttng s ther maxmum [14, 31, 37]. The bat behavour can be used to formulate a new BA. Yang [35] used three generalzed rules when mplementng the bat algorthms: 1. All the bats use an echolocaton to sense the dstance and they also guess the dfference between the food/prey and background barrers n a somewhat magcal way. 2. When searchng for ther prey, the bats fly randomly wth velocty v at poston x wth fxed frequency f mn, varyng wavelength λ and loudness A 0. They can automatcally adjust the wavelength (or frequency) of ther emtted pulses and adjust the rate of pulse emsson r [0, 1], dependng on the proxmty of ther target. 3. Although the loudness can vary n many ways, we assume that t vares from a large (postve) A 0 to a mnmum constant value A mn. 4

5 Algorthm 1 Orgnal Bat Algorthm 1: Objectve functon f(x), x = (x 1,..., x d ) T 2: Intalze the bat populaton x and v for = 1... n 3: Defne pulse frequency Q [Q mn, Q max ] 4: Intalze pulse rates r and the loudness A 5: whle (t < T max ) // number of teratons 6: Generate new solutons by adjustng frequency and 7: update veloctes and locatons/solutons [Eq.(2) to (4)] 8: f(rand(0, 1) > r ) 9: Select a soluton among the best solutons 10: Generate a local soluton around the best soluton 11: end f 12: Generate a new soluton by flyng randomly 13: f(rand(0, 1) < A and f(x ) < f(x)) 14: Accept the new solutons 15: Increase r and reduce A 16: end f 17: Rank the bats and fnd the current best 18: end 19: Postprocess results and vsualzaton The orgnal BA s llustrated n Algorthm 1. In ths algorthm, the bat behavour s captured nto the ftness functon of the problem to be solved. It conssts of the followng components: ntalzaton (lnes 2-4), generaton of new solutons (lnes 6-7), local search (lnes 8-11), generaton of a new soluton by flyng randomly (lnes 12-16) and fnd the current best soluton. 5

6 Intalzaton of the bat populaton s performed randomly. Generatng new solutons s performed by movng vrtual bats accordng to the followng equatons: Q (t) = Q mn + (Q max Q mn )U(0, 1), v (t+1) x (t+1) = v t + (x t best)q (t), = x (t) + v (t), (1) where U(0, 1) s a unform dstrbuton. A random walk wth drect explotaton s used for the local search that modfes the current best soluton accordng to equaton: x (t) = best + ɛa (t) (2U(0, 1) 1), (2) where ɛ s the scalng factor, and A (t) the loudness. The local search s launched wth the proxmty dependng on pulse rate r. The term n lne 13 s smlar to the smulated annealng behavor, where the new soluton s accepted wth some proxmty dependng on parameter A. In lne wth ths, the rate of pulse emsson r ncreases and the loudness A decreases. Both characterstcs mtate the natural bats, where the rate of the pulse emsson ncreases and the loudness decreases when a bat fnds ts prey. Mathematcally, these characterstcs are captured wth the followng equatons: A (t+1) = αa (t), r (t) = r (0) [1 exp( γɛ)], (3) where α and γ are constants. Actually, α parameter plays a smlar role as the coolng factor n the smulated annealng algorthm that controls the convergence rate of ths algorthm. III. DIFFERENTIAL EVOLUTION Dfferental evoluton (DE)[29][6] s a technque for the optmzaton ntroduced by Storn and Prce n DE optmzes a problem by mantanng a populaton of canddate solutons and creates new canddate solutons by combnng the exstng ones accordng to ts smple formulae, and then keepng whchever canddate soluton has the best score or ftness on the optmzaton problem at hand. DE supports a dfferental mutaton, a 6

7 dfferental crossover and a dfferental selecton. In partcular, the dfferental mutaton randomly selects two solutons and adds a scaled dfference between these to the thrd soluton. Ths mutaton can be expressed as follows: u (t) = w (t) r0 + F (w(t) r1 w(t) r2 ), for = 1... NP, (4) where F [0.1, 1.0] denotes the scalng factor as a postve real number that scales the rate of modfcaton whle r0, r1, r2 are randomly selected vectors n the nterval 1... NP. A unform crossover s employed as a dfferental crossover by DE. The tral vector s bult out of parameter values coped from two dfferent solutons. Mathematcally, ths crossover can be expressed as follows: z,j = u (t),j rand j (0, 1) CR j = j rand, w (t),j otherwse, (5) where CR [0.0, 1.0] controls the fracton of parameters that are coped to the tral soluton. Note, the relaton j = j rand assures that the tral vector s dfferent from the orgnal soluton Y (t). Mathematcally, dfferental selecton can be expressed as follows: w (t+1) = z (t) w (t) f f(z (t) ) f(y (t) ), otherwse. (6) In the techncal sense, the crossover and mutaton can be performed n many ways n DE. Therefore, a specfc notaton was used to descrbe the varety of these methods (also strateges) n general. For example, DE/rand/1/bn denotes that the base vector s randomly selected, one vector dfference s added to t, and the number of modfed parameters n the mutaton vector follows the bnomal dstrbuton. IV. HYBRID BAT ALGORITHM As mentoned above, n ths artcle we propose a new BA, called Hybrd Bat Algorthm (HBA). It was obtaned by hybrdzng the orgnal BA usng the DE strateges. The HBA pseudo-code s llustrated n Algorthm 2. 7

8 Algorthm 2 Hybrd Bat Algorthm 1: Objectve functon f(x), x = (x 1,..., x d ) T 2: Intalze the bat populaton x and v for = 1... n 3: Defne pulse frequency Q [Q mn, Q max ] 4: Intalze pulse rates r and the loudness A 5: whle (t < T max ) // number of teratons 6: Generate new solutons by adjustng frequency and 7: updatng veloctes and locatons/solutons [Eq.(2) to (4)] 8: f(rand(0, 1) > r ) 9: Modfy the soluton usng DE/rand/1/bn 10: end f 11: Generate a new soluton by flyng randomly 12: f(rand(0, 1) < A and f(x ) < f(x)) 13: Accept the new solutons 14: Increase r and reduce A 15: end f 16: Rank the bats and fnd the current best 17: end 18: Postprocess results and vsualzaton As a result, HBA dffers from the orgnal BA n lne 9, where soluton s modfed usng DE/rand/1/bn strategy. V. EXPERIMENTS AND RESULTS The goal of our experments was to show that HBA sgnfcantly mproves the results of the orgnal BA. For ths purpose, two BAs were mplemented accordng to specfcatons gven n Algorthms 1 and 2 so that a well-selected set of test functons n the lterature are used for optmzaton benchmarks. The parameters of the two BAs were the same. The dmenson of the problem sgnfcantly affects the results optmzaton. To test the mpact of the problem dmenson on the results, three dfferent 8

9 sets of dmensons were taken nto account,.e., D = 10, D = 20, and D = 30. The functons wth dmenson D = 10 were lmted to maxmally 1,000 generatons, the functons wth dmenson D = 20 to twce as much, whle the functons wth dmenson D = 30 to 3,000. The ntal loudness was set to A 0 = 0.5 and so was also the ntal pulse rate (r 0 = 0.5). The frequency was taken from nterval Q [0.0, 2.0]. The algorthms optmzed each functon 25 tmes and results were measured accordng to the best, worst, mean, and medum values n these runs. The standard devaton of the mean values was calculated as well. A. Test sute Our test sute conssts of fve standard functons taken from lterature [32]. The functons n ths test sute are represented below. 1. The Grewangk s functon The am of the functon s overcomng falures to optmze each varable ndependently. Ths functon s multmodal, snce the number of local optma ncreases wth the dmensonalty. After dmensonaltes are suffcently hgh (n > 30), multmodalty seems to dsappear and the problem turns to be unmodal. f 1 ( x) = D =1 cos ( x ) + D =1 where 600 x 600. The global mnmum of the functon s at 0. x 2 + 1, (7) The Rosenbrock s functon As n case of the Rastrgn s functon, the Rosenbrock s functon too, has ts value 0 at the global mnmum. The global optmum s nsde the parabolc, narrow-shaped flat valley. Varables are strongly dependent on each other, snce t s dffcult to converge the global optmum. where x f 2 ( x ) = D (x +1 x 2 ) 2 + (x 1) 2, (8) =1 9

10 3. Sphere functon where x f 3 ( x ) = D =1 x 2, (9) 4. The Rastrgn s functon Based of the Sphere s functon, the Rastrgn s functon adds cosne modulaton to create many local mnma. Because of ths feature, the functon s multmodal. Its global mnmum s at 0. where x f 4 ( x ) = n 10 + n =1 (x 2 10 cos(2πx )), (10) 5. The Ackley s functon The complexty of ths functon s moderate because of the exponental term that covers ts surface wth numerous local mnma. It s based on the gradent slope. Only the algorthm that uses the gradent steepest descent wll be trapped n a local optmum. The search strategy, analyzng a wder area, wll be able to cross the valley among the optma and acheve better results. f 5 (x) = n 1 =1 [20+e 20e (x x2 ) e 0.5(cos(2πx +1)+cos(2πx )) ], where x The global mnmum of ths functon s at 0. (11) B. PC confguraton Confguraton of PC used n our experments s as follows: HP Pavlon g4, processor Intel(R) Core(TM) 2.40 GHz, 10

11 memory 8 GB, mplemented n C++. TABLE I: The results of experments Alg. D Value f 1 f 2 f 3 f 4 f 5 BA Best 3.29E E E E E+01 Worst 1.73E E E E E Mean 8.30E E E E E+01 Medan 3.91E E E E E+00 StDev 6.94E E E E E+01 Best 8.77E E E E E+02 Worst 1.43E E E E E Mean 1.46E E E E E+02 Medan 1.90E E E E E+05 StDev 1.64E E E E E+01 Best 1.58E E E E E+02 Worst 4.18E E E E E Mean 1.51E E E E E+02 Medan 4.66E E E E E+06 StDev 1.52E E E E E+01 Best 2.25E E E E E-04 Worst 3.97E E E E E Mean 3.18E E E E E+01 Medan 8.66E E E E E+00 StDev 1.14E E E E E+01 Best 1.01E E E E E-05 Worst 2.96E E E E E+01 HBA 20 Mean 8.56E E E E E-05 Medan 3.60E E E E E+02 StDev 2.17E E E E E+01 Best 6.38E E E E E-04 Worst 3.57E E E E E Mean 6.42E E E E E-03 Medan 5.99E E E E E+02 StDev 3.12E E E E E+01 11

12 C. Experment results The results of our extensve numercal experments are summarzed n Table I. The table represents results of BA and HBA algorthms (column 1) solvng the test sute of fve functons (denoted as f 1, f 2, f 3, f 4, f 5 ) wth dmensons D = 10, 20 and 30, respectvely, accordng to the best, worst, mean, medan, and standard devaton values. The results of the HBA algorthm show that ths algorthm sgnfcantly mproves the results of the orgnal BA algorthm accordng to almost all the measures safe for the standard devaton n some cases (e.g., when usng Ackley s functon). There was no statstcal analyss of the results made because they are evdently better for HBA than for BA. To show how the results of the two algorthms (.e., BA and HBA) vary wth the dmensons of the functons, the mean value of functons f 1 and f 3 wth dmensons D = 10, D = 20, and D = 30 are presented n Fgs The logarthmc scale s used to dsplay the mean value on the y-axs. The hgher the mean value, the more dffcult the functon s to solve. Mean value 1.0e e e e e e-08 BA HBA 1.0e-10 D=10 D=20 D=30 Dmensons FIG. 1: Mean value of functon f 1 wth varous dmensons. As seen from Fg. 1, the best results are obtaned by optmzng functon f 1 wth the dmenson D = 20, whle the worst by optmzng the same functon wth the hghest dmenson D = 30. It s nterestng to note that the mean value of functon f 5 wth dmenson D = 10 s most dffcult to obtan for the HBA algorthm, whle the same functon but wth dmenson D = 20 s the easest to solve. Oppostely, the results of the orgnal BA algorthm show that by ncreasng the dmensons, the results become worsen. 12

13 Mean value 1.0e e e e e e-08 BA HBA 1.0e-10 D=10 D=20 D=30 Dmensons FIG. 2: Mean value of functon f 3 wth varous dmensons. Mean value 1.0e e e e e e-08 BA HBA 1.0e-10 D=10 D=20 D=30 Dmensons FIG. 3: Mean value of functon f 5 wth varous dmensons. As seen from Fg. 2, solvng functon f 3 becomes more effcent when the dmensonalty of the problem ncreases. As a result, the most dffcult functon to solve s functon f 3 wth dmenson D = 30. Comparng the results obtaned wth the BA and the HBA algorthm shows that the results of HBA sgnfcantly outperform the results of the orgnal BA algorthm by optmzng the functons f 1, f 3, and f 5. Functons f 2 and f 4 of the HBA algorthm are also better, but the dfference s not sgnfcant. VI. CONCLUSION In ths paper we mprove the bat algorthm, BA, by developng ts new varant, the socalled hybrd bat algorthm, HBA. HBA s a hybrd of BA wth DE strateges. As shown wth our experments, HBA mproves sgnfcantly the orgnal varant of BA. In future, HBA 13

14 wll be tested on a large-scale global optmzaton. Our testng wll be extended by usng more dverse test functon sets and by deepenng our parametrc study. [1] BV Babu and R. Angra. Modfed dfferental evoluton (mde) for optmzaton of non-lnear chemcal processes. Computers & chemcal engneerng, 30(6): , [2] T. Bäck. Evolutonary algorthms n theory and practce - evoluton strateges, evolutonary programmng, genetc algorthms. Oxford Unversty Press, [3] J. Brest, S. Grener, B. Boskovc, M. Mernk, and V. Zumer. Self-adaptng control parameters n dfferental evoluton: A comparatve study on numercal benchmark problems. Evolutonary Computaton, IEEE Transactons on, 10(6): , [4] J. Brest, A. Zamuda, I. Fster, and M.S. Maucec. Large scale global optmzaton usng selfadaptve dfferental evoluton algorthm. In Evolutonary Computaton (CEC), 2010 IEEE Congress on, pages 1 8. IEEE, [5] C. Darwn. The orgn of speces. John Murray, London, UK, [6] S. Das and P.N. Suganthan. Dfferental evoluton: A survey of the state-of-the-art. Evolutonary Computaton, IEEE Transactons on, 15(1):4 31, [7] A.E. Eben and J.E. Smth. Introducton to Evolutonary Computng. Sprnger-Verlag, Berln, [8] I. Fster and J. Brest. Usng dfferental evoluton for the graph colorng. In Dfferental Evoluton (SDE), 2011 IEEE Symposum on, pages 1 7. IEEE, [9] I. Fster, I. Fster, and J. Brest. A hybrd artfcal bee colony algorthm for graph 3-colorng. Swarm and Evolutonary Computaton, pages 66 74, [10] I. Fster, I. Fster Jr, and J.B.V. Zumer. Memetc artfcal bee colony algorthm for large-scale global optmzaton. In Evolutonary Computaton (CEC), 2012 IEEE Congress on, pages 1 8. IEEE, [11] I. Fster, M. Mernk, and B. Flpč. Graph 3-colorng wth a hybrd self-adaptve evolutonary algorthm. Computatonal optmzaton and applcatons, pages 1 32, [12] I. Fster Jr, X.S. Yang, I. Fster, and J. Brest. Memetc frefly algorthm for combnatoral optmzaton. arxv preprnt arxv: , [13] L. J. Fogel, A. J. Owens, and M. J. Walsh. Artfcal Intellgence through Smulated Evoluton. 14

15 John Wley, New York, USA, [14] A.H. Gandom, X.S. Yang, A.H. Alav, and S. Talatahar. Bat algorthm for constraned optmzaton tasks. Neural Computng & Applcatons, pages 1 17, [15] D. Goldberg. Genetc Algorthms n Search, Optmzaton, and Machne Learnng. Addson- Wesley, MA, [16] John H. Holland. Adaptaton n Natural and Artfcal Systems: An Introductory Analyss wth Applcatons to Bology, Control and Artfcal Intellgence. MIT Press, Cambrdge, MA, USA, [17] N. Jn and Y. Rahmat-Sam. Advances n partcle swarm optmzaton for antenna desgns: Real-number, bnary, sngle-objectve and multobjectve mplementatons. Antennas and Propagaton, IEEE Transactons on, 55(3): , [18] D. Karaboga and B. Basturk. A powerful and effcent algorthm for numercal functon optmzaton: artfcal bee colony (abc) algorthm. Journal of Global Optmzaton, 39(3): , [19] D. Karaboga and B. Basturk. On the performance of artfcal bee colony (abc) algorthm. Appled Soft Computng, 8(1): , [20] P. Korošec, J. Šlc, and B. Flpč. The dfferental ant-stgmergy algorthm. Informaton Scences, [21] John R. Koza. Genetc programmng 2 - automatc dscovery of reusable programs. Complex adaptve systems. MIT Press, [22] D. Merkle, M. Mddendorf, and H. Schmeck. Ant colony optmzaton for resource-constraned project schedulng. Evolutonary Computaton, IEEE Transactons on, 6(4): , [23] F. Ner and V. Trronen. Recent advances n dfferental evoluton: a survey and expermental analyss. Artfcal Intellgence Revew, 33(1 2):61 106, [24] Q.K. Pan, M. Fath Tasgetren, P.N. Suganthan, and T.J. Chua. A dscrete artfcal bee colony algorthm for the lot-streamng flow shop schedulng problem. Informaton scences, 181(12): , [25] RS Parpnell and HS Lopes. New nspratons n swarm ntellgence: a survey. Internatonal Journal of Bo-Inspred Computaton, 3(1):1 16, [26] R.S. Parpnell, H.S. Lopes, and A.A. Fretas. Data mnng wth an ant colony optmzaton algorthm. Evolutonary Computaton, IEEE Transactons on, 6(4): ,

16 [27] J. Robnson and Y. Rahmat-Sam. Partcle swarm optmzaton n electromagnetcs. Antennas and Propagaton, IEEE Transactons on, 52(2): , [28] Y. Sh, H. Teng, and Z. L. Cooperatve co-evolutonary dfferental evoluton for functon optmzaton. Advances n natural computaton, pages , [29] R. Storn and K. Prce. Dfferental evoluton a smple and effcent heurstc for global optmzaton over contnuous spaces. Journal of global optmzaton, 11(4): , [30] M.F. Tasgetren, Q.K. Pan, P.N. Suganthan, and A.H.L. Chen. A dscrete artfcal bee colony algorthm for the permutaton flow shop schedulng problem wth total flowtme crteron. In Evolutonary Computaton (CEC), 2010 IEEE Congress on, pages 1 8. IEEE, [31] P.W. Tsa, J.S. Pan, B.Y. Lao, M.J. Tsa, and V. Istanda. Bat algorthm nspred algorthm for solvng numercal optmzaton problems. Appled Mechancs and Materals, 148: , [32] X.-S. Yang. Appendx a: Test problems n optmzaton. In X.-S. Yang, edtor, Engneerng Optmzaton, pages John Wley & Sons, Inc., Hoboken, NJ, USA, [33] X.-S. Yang, S. Deb. Cuckoo search va lévy flghts. In Nature & Bologcally Inspred Computng, NaBIC World Congress on, pages IEEE, [34] X.-S. Yang, S. Deb. Engneerng optmsaton by cuckoo search. Internatonal Journal of Mathematcal Modellng and Numercal Optmsaton, 1(4): , [35] X.S. Yang. A new metaheurstc bat-nspred algorthm. Nature Inspred Cooperatve Strateges for Optmzaton (NICSO 2010), pages 65 74, [36] X.S. Yang. Bat algorthm for mult-objectve optmsaton. Internatonal Journal of Bo- Inspred Computaton, 3(5): , [37] X.S. Yang. Revew of meta-heurstcs and generalsed evolutonary walk algorthm. Internatonal Journal of Bo-Inspred Computaton, 3(2):77 84, Updated 30 Aprl