The paper published by Yang, X. S. (2010). “A New Metaheuristic Bat-Inspired Algorithm, in: Nature Inspired Cooperative Strategies for Optimization (NISCO 2010)”. Studies in Computational Intelligence. 284: 65–74 is implemented at https://free-thesis.com. A demo version of this is available here but we found it hard to use in any application. So developed in house full BAT code and happy to share here.
The Bats have some echolocation characteristics by using them we can develop an Algorithm for the optimization. Some rules for the Bat algorithm explained below.
- The distance can be sense by the bats using the echolocation and they have also difference between food/prey and the background barriers.
- The velocity of the Bats represented by (vi) and position of the Bat is (xi)having frequency Fmin varying wavelength and loudness A0 search for food.
- The loudness A0 also take as Amin constant value.
The Matlab coding for the Bat function done by consider various parameters. Various steps are involved for the Bat coding.
Step1. Firstly make the statement for Bat function
In above statement the input parameter is mainly a benchmark function which is represented by a ‘fobj’ and others are lb=lower bound limit and ub=upper bound limit , fmax is the maximum frequency and fmin is the minimum frequency , A= loudness of each Bat, r=pulse emission rate of each Bat, alpha and gamma are the constants for the loudness and pulse emission rate. The r0 is the initial pulse rate.
Step2. After the statement call the initialization function. The script for the initialization function written separately.
x=initialization(N, Max_iter, dim, ub, lb)
In initialization upper bound and lower bound limits are available and bats position is randomly search. Each bat has different upper and lower bound limits. There is N number of searching variable. Calculate the initial position of the search variable. We initialize the parameters of algorithm, generate and also evaluate the initial population, and then determine the best solution in the population. Here, virtual bats are moved in the search space according to updating rules of the bat algorithm.
Step3. Call the benchmark function here which is represented by the ‘fobj’and find the initial best fitness value for benchmark objective function.
The fobj function contains all the information about the benchmark function. It has 23 different benchmark function cases which have different dimension, upper bound and lower bound limits. We can randomly take any benchmark (F1 …..F23) objective function.
The initial best value is obtained
Also find the best solution for the best fitness value of the objective function.
Step4. Then start the main loop for the maximum iterations. Randomly choose the frequency as equation
Then update the velocity of the bats and position of the bats. After the update of both the terms than applied the upper and lower bound limits and update the position of bats again.
Step5. Check the condition of pulse rate emission of each bat (r). The factor 0.001 is limits the step size of the random walk.
Find the new value of the objective after check the condition. It is the new fitness value in which the loudness is too much. The loudness and emission rate updated if the new solution is improved which means the bats moving towards the optimal solution. The condition for improve the loudness and pulse emission rate is
(fitnessnew<=fitness(ii)) && (rand<A(ii)) A(ii)=alpha*A(ii) r(ii)=ro*(1-exp(-gamma*iter))
Step6. Calculate the best fitness value after optimization. Draw the convergence curve according to the best fitness value and iteration. The best position is also obtained through the optimization.
Step7. Write another script for the observation of benchmark function. Initialize the parameters like N=number of searching Bats, any benchmark function (F1…………F10).
- Call the benchmark function- in this step all the bench mark function is called.
- Call the BAT function – BAT function is calling in this step. The benchmark function initialize through the bat algorithm.
- Call the function plot for graphs and convergence curve for the benchmark function. The best position and best fitness value obtained by using the bat algorithm is shown by the curve.
In the general context of exploration and exploitation as well as genetic operators, we can analyze the roles of the main components of the standard BA. In essence, frequency tuning essentially acts as mutation because it varies the solutions mainly locally. However, if mutation is large enough, it can also leads to global search. Certain selection is carried out by applying a selection pressure that is relatively constant due to the use of the current best solution found so far. Compared with genetic algorithms, there is no explicit crossover; however, mutation varies due to the variations of loudness and pulse emission. On the other hand, the variations of loudness and pulse emission rates can also provide an auto zooming ability in the sense that exploitation becomes intensive as the search is approaching the global optimal. This essentially switches an exploitative phase to an exploitative phase automatically
The graph and curve which are obtained by using the BAT algorithm for the optimization of benchmark function.
The left hand side figure shows the parameter space for the benchmark function and right hand side figure shows the objective space of same benchmark function. The objective space curve is between the best fitness obtained value and iterations. The curve shows the minimum value of the objective function during the optimization. As shown in the curve the best fitness value near to 3.
We recommend our users to use the hybrid PSO-GWO optimization algorithm whose code is published here. https://free-thesis.com/product/hybrid-particle-swarm-and-grey-wolf-optimization/