Abstract

In the last two decades, the field of global optimization has become very active, and, in this regard, many deterministic and stochastic algorithms were developed for solving various optimization problems. Among them, swarm intelligence (SI) is a stochastic algorithm that is more flexible and robust and has had the ability to find an optimum solution for high-dimensional optimization and search problems. SI-based algorithms are mainly inspired by the social behavior of fish schooling or bird flocking. Among the SI-based algorithms, Bat algorithm (BA) is one of the recently developed evolutionary algorithms. It employs an echolocation behavior of microbats by varying pulse rates of emission and loudness to perform their search process. In this paper, a modified Bat algorithm (MBA) is developed. The main focus of the MBA is to further enhance the exploration and exploitation search abilities of the original Bat algorithm. The performance of the modified Bat algorithm (MBA) is examined over the benchmark functions designed for evolutionary algorithms competition in the special session of 2005 IEEE Congress on Evolutionary Computation. The used benchmark functions include the unimodal, multimodal, and hybrid benchmark functions with high dimensionality. Furthermore, the impact analysis with respect to different values of temperatures is conducted by executing the proposed algorithm twenty-five times independently by using each benchmark function with different random seeds.

1. Introduction

Global optimization has wide applications in many fields including physics, biology, engineering, economics, and business. Examples include minimizing the cost/consumption, while maximizing the profit/product is the need of the day in every time which leads human being towards the optimization process. It has been an essential part of research in applied mathematics since six or seven decades. The advancement of computers makes it possible to solve those optimization problems that were unsolvable before. The main objective of the optimization is to maximize the efficiency of the system while minimizing its cost. A scalable optimization problem can typically and generally be formulated as follows:where is the candidate solution with decision variables, comprises objective function, and are constrained functions with inequality, and are equality constrained functions. Furthermore, a solution vector is said to be globally optimal if , , and the corresponding value of is called global minimum. If and , then is single objective bound constrained function, while if , then the problem referred to as (1) is a multiobjective problem. The decision variables of the solution vector can be continuous, semicontinuous, discrete, integer, or as one of a finite set of explicit discrete categorical labels. Optimization problems can be classified based on the type of constraints, nature of design variables, physical structure of the problem, nature of the equations involved, deterministic nature of the variables, permissible value of the design variables, separability of the functions, and number of objective functions.

In the last few decades, different kinds of linear and nonlinear optimization techniques were developed to solve various linear and nonlinear optimization problems [15]. The optimization methods are further divided into two subclasses, local search methods and global search methods, as shown in Figure 1. The major difference between the two aforementioned subclasses is that local search methods work on single solution and provide a local optimum for the given problem, while global search methods operate on a set of solutions called population and usually provide a global optimum for given optimization problem.

Evolutionary computation (EC) is the subfield of soft computing [5, 7, 8] which studies a family of evolutionary algorithms (EAs) for global optimization. EAs are inspired by biological evolution and operate on a set of solutions which are generated uniformly and randomly undergo the crossover, mutation, selection, and survival of the fittest for next generation of evolutionary process. The four classical paradigms of EAs are the genetic algorithm (GA) [9], evolution strategy [10], evolution programming [11], and genetic programming [12]. Differential evolution [13, 14], particle swarm optimization [15], harmony search algorithm [16], gravitational search algorithm [17], biogeography-based optimization [18], ant lion optimizer [19], invasive weed optimization [20], shuffled frog leaping [21], ant colony optimization [22], artificial bee colony [23], firefly algorithm [24], Bat algorithm [25], cuckoo search [26], flower pollination algorithm [27], the whale optimization algorithm [28], grey wolf optimizer [29], salp swarm algorithm [30], ant colony optimization (ACO) [31, 32], firefly algorithms [33], and Bat algorithms [34] are all almost nature-inspired and swarm-intelligence-based algorithms.

Recently, in EC communities, hybrid swarm-intelligence-based algorithms and hybridization of nature-inspired algorithms [35] have got a great deal of attention to solve various test suites of single objective optimization and search problems [3647]. However, proper settings of the control parameters include population size, selection rate, and operator probabilities, and other important intrinsic parameters in EAs play important role in their success while dealing with different optimization and search problems [4850]. In this paper, a modified Bat algorithm (MBA) is developed, aiming at further improving the exploration versus exploitation issues of the baseline Bat algorithm in terms of convergence speed, convergence precision, and global searching ability. The performance of the proposed algorithm is evaluated over the benchmark functions, which are designed for the special session of evolutionary algorithms competition of the IEEE Congress on Evolutionary Computation [51]. Furthermore, experiments are carried out with different values of temperature to show their impact on the performance of the proposed algorithm for dealing with each IEEE-CEC benchmark function [51].

The rest of the paper is organized as follows: Section 2 introduces the framework of the proposed modified Bat algorithm. Section 3 demonstrates the experimental results along with characteristics of benchmark functions used in the experiments we carried out. Section 4 finally concludes this paper.

2. Modified Bat Algorithm for Global Optimization

In the last few years, swarm intelligence- (SI-) based approaches have gained much attention in both academia and industrial applications [52]. Among them, ant colony optimization (ACO) [22], artificial bee colony (ABC) [3, 23, 5256], cuckoo search (CS) [57, 58], firefly algorithm (FA) [24, 33, 34], and cuckoo search algorithm [5759] are commonly used SI-based methods. They are inspired by collective behaviors of human beings, insects, birds flocking, the behavior of honeybees, the brooding behavior of cuckoo species, the flashing behavior of tropical fireflies, and the flocking behavior of birds. SI-based algorithms have tackled various optimization and search problems including traveling salesman problem (TSP) [58, 59], controlling unmanned vehicles, data mining, and many others.

Bat algorithm (BA) was first developed by Yang for global optimization problems [6062]. It was inspired by the echolocation behavior of microbats, with varying pulse rates of emission and loudness. Bats are the only mammals with wings, and they have advanced capability of echolocation. Bats, dolphins, and other animals use the echolocation technique to determine the location of objects with reflected sound. This allows the animals to move around in pitch darkness, so they can navigate, hunt, identify friends and enemies, and avoid obstacles. The typical range of frequencies for most bat species is to . The bats emit a very loud sound pulse and listen for the echo that bounces back from objects. Bats emit to 20 sound bursts every second. The pulse emission is going to increase up to 200 pulses per second while hunting prey. The wavelength of the ultrasonic sound bursts in the range from to . Recently, an improved version of BA with novel initialization techniques was introduced in [63]. The Bats harness three rules: (i) to sense distance, background barriers, and dissimilarity between prey/food, bats use echolocation; (ii) bats fly randomly with a certain velocity to a specific position with a fixed frequency with varying wavelengths to search for prey; and (iii) the loudness changes from a large to minimum constant value [64]. Three main traits of the microbat are employed to design the basic framework of BA as explained below:(i)Most of the species of microbats exercise echolocation to search for their prey.(ii)To search for prey, the microbats search with fixed frequency,, and varying wavelength, .(iii)According to Yang’s method, the mathematical formulation of microbats is as follows:where is the minimum frequency and is the maximum frequency used by the Bat while searching for its prey, denotes the location of the Bat in the search space, represents the velocity of the Bat, and is the parameter that controls the difference between the minimum and maximum frequencies of each Bat.

2.1. Exploration and Exploitation in the Proposed Algorithm

In this paper, equation (3) is used in the modified BA (MBA) and as described in equation (4) along with as reported in equation (5), where is the distance between the sound wave source and the target in order to bound the wave back as defined by equation (5).where . In the active sonar system, the distance between the sound wave source and the target that bounds the wave back as given in equation (5), where means the time difference between the sending sound wave and receiving echo and represents the sound speed in the air, where . The bats change their position according to equations (6) and (7) as per criteria given in flowchart 2.where move the bats around the current best solution. The impact analysis of temperature with sound speed in air was settled to be constant at at [65]. In general, the speed of sound varies at different temperature levels in the air. The environmental effect of temperature with different values is also examined in the proposed MBA as outlined in Algorithm 1. The flowchart of the proposed modified Bat algorithm (MBA) is described in Figure 2. The main steps of the proposed MBA are summarized as follows:Step 1. Initialization: spread bats uniformly and randomly within the search space of the problem.Step 2. Generate new bats by using equations (4) to (7).Step 3. Evaluate bats and update their global near best position (7).Step 4. Check the termination condition to decide whether to continue or terminate the process.

(1): population size;
(2): number of decision variables
(3): maximum Iteration.
(4): random number, : maximum Iteration.
(5)Initialize bat the population ,
(6)Compute , ,
(7)Find the current best
(8)
(9)while it ≤ Mitdo
(10) Generate new solution by flying randomly.
(11)
(12)
(13)
(14)
(15)
(16) Compute the objective function value,
(17)ifthen
(18) Accept .
(19)else
(20) Retain
(21)end if
(22) Update current best solution denoted by .
(23) Move the bats around current best solution randomly.
(24)
(25) Repeating steps 11 to 15 and update the global best solution.
(26);
(27)end while
2.2. Impact of Temperature with respect to the Speed of Sound in Air

The sound waves are longitudinal, which are propagating in air with speed described as follows:

Therefore,

The speed of sound in air is directly proportional to the square root of the absolute temperature of the air. If and are the speeds of sound at temperatures and , respectively, then

We have

Applying binomial theorem and neglecting higher-order terms, we get

We have

3. Experimental Results and Discussion

Due largely to the nature of evolutionary algorithms (EAs), their behaviors and performance are mainly studied experimentally by using various test suites of optimization and search problems. For these purposes, several continuous test suites of benchmark functions have already been designed in the existing EC community over the last few years. These benchmark functions have played a crucial role in developing and scientifically studying the behavior of various existing EAs and modified enhanced versions. The benchmark functions which have recently been designed for the special session of the 2005 IEEE Congress on Evolutionary Computation are utilized in the study conducted in this paper in order to carefully assess the proposed modified Bat algorithm (MBA) and impact analysis of various temperatures used in its framework. For more details about these instances, the readers are referred to [51]. These benchmark functions that can be categorized as to are unimodal functions, to are basic multimodal functions, and are both expanded multimodal benchmark functions, and is a hybrid composition test problem. Optimization problems are classified on the basis of landscape as follows:(i)Unimodel. If the optimization problem has only one optimum that is global and local, it is called unimodal, for example, sphere function.(ii)Multimodel. If the optimization problems have more than one optimum, those may all be global, or some of them may be local, called unimodal, for example, Ackley function.

The second last column of Table 1 provides the known optimal solution of each of the benchmark functions.

3.1. PC Configuration and Programming Environment

All experiments were performed on a computer with Intel Core CPU, 2.40 GHz and 4 G RAM, under Windows 10 Pro, 64-bit OS. The proposed algorithm and all other existing algorithms used in the comparative analysis were implemented in MATLAB R2017b programming environment. The proposed algorithm was executed 25 times independently with different random seeds by using the Matlab command rand (“state,” sum ).

3.2. Parameter Settings

All experiments were performed by solving the benchmark functions with parameters settings described as follows:(i): population size.(ii): the dimension of the search space.(iii): Maximum iteration.(iv): velocity of sound at .(v): the global best.(vi): location of the Bat in the solution space.(vii): the number of independent runs.(viii): minimum frequency.(ix): maximum frequency.(x), : % different temperature in Celsius scale.(xi): distance between sound waves source and the target.(xii) and are random numbers, respectively.(xiii): time difference between the sent sound wave and echo.

Tables 2 and 3 include the numerical results provided by the proposed modified Bat algorithm in terms of minimum objective values, average objective values, standard deviation, and maximum objective function values after independent runs of twenty-five simulations with random seeds to solve each benchmark function whose characteristics are given in Table 1. All these numerical values were obtained with different values of temperatures like , , 30, and 40 in the framework of the suggested algorithms, respectively. The last columns of Tables 2 and 3 include the average CPU time in seconds elapsed by the proposed modified Bat algorithm while solving each corresponding benchmark function.

Figure 3 presents the convergence graphs displayed by the proposed algorithm to solve the benchmark functions as denoted by to . Different values , , 30, and 40 along the x-axis represent the variation in temperature values and the evolution in minimum objective function values is shown along the y-axis. Similarly, Figure 4 displays the convergence of the benchmark denoted by to with different temperature settings like , . Figure 5 displays the convergence for the test problems denoted by to in the twenty-five independent runs of executions.

4. Conclusion

Many real-life problems can be modeled as global optimization problems. Evolutionary algorithms (EAs) are population-based optimization methods that operate on a set of uniformly and randomly generated solutions in search space as specified for each optimization problem to be solved. EAs are more flexible and usually provide a high-quality set of solutions in a single simulation run, in contrast to traditional optimization methods. This paper presents a modified Bat algorithm to deal with benchmark functions that were designed for the special session of EAs competition of 2005 IEEE Congress on Evolutionary Computation (IEEE-CEC’05) [51]. The results provided by the modified Bat algorithm (MBA) seem to be more effective for dealing with most benchmark functions in terms of diversity and proximity:(i)In MBA, the basic framework of the Bat Algorithm is modified(ii)Based on the experiments that we carried out, the speed of sound varied with varying temperature(iii)The proposed algorithm has found promising results for most problems at (iv)The multimodel problems are tackled efficiently by the proposed algorithm at

In the future, we intend to conduct the following tasks:(i)To determine the efficiency, stability, and importance of the proposed MBA by using other existing unconstrained benchmark functions and some real-world problems(ii)To examine sensitivity of the different intrinsic parameters used in the framework of the proposed algorithm in order to ensure the correct use of parameters keeping in view the structure of different problems(iii)To further evaluate the robustness and efficiency of the proposed algorithm against some state-of-the-art evolutionary algorithms

Data Availability

Data are available upon request from the first author. Furthermore, the Matlab code of the benchmark functions used to conduct this study is available on https://www3.ntu.edu.sg/home/epnsugan/.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors thank and acknowledge the support of the Research Intensified Grant Scheme (RIGS) Fasa 1/2019, Vot Number 55192/7, Universiti Malaysia Terengganu.