Elite Opposition-Based Water Wave Optimization Algorithm for Global Optimization

Water wave optimization (WWO) is a novel metaheuristic method that is based on shallow water wave theory, which has simple structure, easy realization, and good performance even with a small population. To improve the convergence speed and calculation precision even further, this paper on elite opposition-based strategy water wave optimization (EOBWWO) is proposed, and it has been applied for function optimization and structure engineering design problems.There are threemajor optimization strategies in the improvement: elite opposition-based (EOB) learning strategy enhances the diversity of population, local neighborhood search strategy is introduced to enhance local search in breaking operation, and improved propagation operator provides the improved algorithm with a better balance between exploration and exploitation. EOBWWO algorithm is verified by using 20 benchmark functions and two structure engineering design problems and the performance of EOBWWO is compared against those of the state-of-the-art algorithms. Experimental results show that the proposed algorithmhas faster convergence speed, higher calculation precision, with the exact solution being even obtained on some benchmark functions, and a higher degree of stability than other comparative algorithms.


Introduction
Optimization problem is wide and varied; in many scientific and engineering computation areas, a majority of problems of people encounter can be attributed to objective optimization problem; thus, the research of optimization problem has been a very active field.In fact, the optimization method can be divided into two types of deterministic optimization and stochastic optimization; although the deterministic optimization method is relatively mature, its application condition is harsh and difficult to deal with large-scale optimization problems, which prompted the stochastic optimization method, especially development of heuristic optimization method.
In recent years, the heuristic optimization algorithm, especially the metaheuristic optimization algorithm, has been concerned by many researchers.Metaheuristic optimization algorithm originates from the simulation of various types of physical, biological, social, and other phenomena in nature to solve optimization problems.As a stochastic optimization method, the metaheuristic optimization algorithm has the advantage of simple and universal, strong robustness, suitable for parallel processing and wide application range.Due to the advantages of metaheuristic, several algorithms have been proposed recently, such as particle swarm optimization (PSO) [1], genetic algorithm (GA) [2], ant colony optimization (ACO) [3], artificial bee colony (ABC) [4], cuckoo search (CS) [5], bat algorithm (BA) [6], firefly algorithm (FA) [7], flower pollination algorithm (FPA) [8], and water wave optimization (WWO) [9].
Water wave optimization (WWO) is a relatively new metaheuristic initially proposed by Zheng in 2015 [9], inspired by the shallow water wave theory [10] for global optimization.WWO has the advantages of simple framework and thus easiness of implementation; even with a small population size it performs well [9].At present, as a new metaheuristic optimization method, WWO has been successfully applied to the optimization problems such as high speed [9] and TSP [11].

Mathematical Problems in Engineering
In order to further improve the performance of WWO, some modified approaches are introduced to strengthen its performance.Zhang et al. [12] improved on WWO; an improved version with variable population size (VC-WWO) is proposed by them, and, meanwhile, a comprehensive learning mechanism is developed in refraction operator to increase the solution diversity.Zheng and Zhang [13] developed a simplified version of WWO (Sim-WWO); in Sim-WWO, leaving out the refraction operator and in order to better balance exploration and exploitation as well as partially compensate the effect of weeding out refraction operator, a strategy of population size reducing is introduced.In order to apply WWO to combinatorial optimization problem, the traveling salesman problem (TSP), Wu et al. [11] redefined the propagation, breaking, and refraction operator based on the original WWO.In this paper, an improved water wave optimization algorithm based on elite-opposition (EOBWWO) learning strategy has been applied to function optimization and structure engineering design problems.The improvements include three parts: elite oppositionbased learning (EOBL) strategy enhances the diversity of population, local neighborhood search strategy is introduced to enhance local search in breaking operation, and improved propagation operator provides the improved algorithm with a better balance between exploration and exploitation.We tested the performance of EOBWWO on 20 benchmark functions and two structure engineering design problems.The experimental results show that proposed algorithm has signification performance advantage including a fast convergence speed and a high calculation precision; in addition, the improved algorithm is able to obtain the exact solution on some test functions.
This paper is organized into the following sections.Section 2 introduces the original WWO algorithm briefly, the detailed description of EOBWWO algorithm is presented in Section 3, simulation experiments and results discussed are described in Section 4, and finally the conclusion is given in Section 5.

Water Wave Optimization (WWO) Algorithm
Water wave optimization (WWO) algorithm is inspired by shallow water wave theory and developed by Zheng [9], where each individual in the population is analogous to the "water wave" object with a wave height ℎ and a wavelength .
Without losing generality, suppose there is a maximization problem  and its objection function is , where practical problem  can be compared with the shallow water wave model; the corresponding relation is shown as follows: Practical problem  Shallow water wave model The search space of  analogous to

󳨀 󳨀󳨀󳨀󳨀󳨀󳨀󳨀󳨀 → Seabed area
Each solution of  analogous to

󳨀 󳨀󳨀󳨀󳨀󳨀󳨀󳨀󳨀 → A water wave objection
The fitness of each solution → It is inversely proportional to vertical distance to the seabed. ( When the population is initialized, for each wave, the wave height ℎ is set to a constant ℎ max and wavelength  is generally set to 0.5.The fitness value of each water wave is inversely proportional to the vertical distance to the seabed; from this we can know that from the seabed nearer the water wave fitness value is bigger, the wave height is bigger, and the wavelength is smaller, as illustrated in Figure 1.During the process of optimization problem-solving, search globally in the solution space by simulating the propagation, breaking, and refraction operation of water waves.

Propagation.
In WWO, all water waves have to be propagated once at each generation.It is assumed that the original water wave is ,   is a new wave created by propagation operator, the dimension of the maximum value function  is , the propagation operation is shifted, and each dimension of the original water wave  is given as where  ∈ , rand(−1, 1) is used to control the propagation step which is a uniformly distributed random number fixed in [−1, 1], and () is the length of the th dimension of the search space.If the length of () is longer than the length of the th dimension of the search space, a new position will be reset randomly as where () and () are the lower bound and upper bound of th dimension of the search space and rand( ) is a random number within the range [0, 1].
After propagating, we evaluate fitness of   ; if (  ) > (),   instead of  in the population, meanwhile the wave height of   is reset to ℎ max ; otherwise,  remained, and in order to simulate energy dissipation of wave in the process of propagation, its height is decreased by one.
It is a natural phenomenon that when a wave travels from deep water to shallow water, its wave height increases and its wavelength decreases, as illustrated in Figure 1.In a bid to simulate this phenomenon, WWO uses the way in which the wavelength of each wave is updated after each generation as follows: where  is a control parameter named wavelength reduction coefficient,  max and  min are the maximum and minimum fitness values among the current population, respectively, and  is a very small positive constant to avoid division-by-zero.

Breaking.
In the water wave theory, with the energy of water wave increasing constantly, crest becomes more and more steep, and the wave breaks into a series of solitary waves when its velocity of crest exceeds the wave celerity.After propagating, WWO only performs breaking on the wave  which is a new best solution  * , which is used to improve the diversity of the population.The detailed process is as follows: first of all, we select randomly  dimensions (where  is a random number between 1 and a predefined number  max ) and perform operations on each selected dimension of original wave  to generate each dimension of solitary wave   as follows: where (0, 1) is a Gaussian random number with mean 0 and standard deviation 1 and  is breaking coefficient.If the solitary wave with the best fitness is better than  * ,   is selected instead of  * ; otherwise,  * remained.

Refraction.
In WWO, the refraction operation only performs on a wave  whose height decreases to zero to avoid search stagnation, which simulates the phenomenon that wave ray is not perpendicular to the isobath.By refraction, in the way that random number centered halfway between the original positions and  * to calculate each dimension of new wave   , the details are as follows: Followed by refraction, the wave height of   is also reset to ℎ max ; meanwhile its wavelength is updated as follows: To sum up, the role of propagation operator is to make the high fitness wave exploit small area and the low fitness wave explore large area, the breaking operator enhances the local search among the promising best waves, and the refraction operation helps avoid search stagnation and thus reduces premature convergence.The basic framework of WWO is as Algorithm 1 [9].

The EOBWWO Algorithm
In order to improve the performance including global searching and local searching abilities of WWO and obtain a better balance between exploration and exploitation even further, there are three optimization strategies applied to the original WWO; they are elite opposition-based learning (EOBL) strategy [42], local neighborhood searching (LNS) strategy [43], and improved propagation operator.

Elite Opposition-Based Learning (EOBL) Strategy.
The optimization process of WWO algorithm can be regarded as the transformation continually of its search space.When the algorithm falls into local optimum, the search space is difficult to contain the global optimal solution.Thus it is very significant to guide the current solution space approximation to the space of global optimal solution.In a bid to enhance the global search ability (i.e., exploration ability) of WWO, the elite opposition-based learning (EOBL) strategy is introduced.
Before introducing the EOBL, we should firstly explain opposition-based learning (OBL) [44].The main idea of OBL is that it generates the opposition solution of current solution, evaluates current solution and opposition solution at the same time, and chooses the better one to enter the next iteration.We assume  = ( 1 ,  2 , . . .,   ) is a point in current population ( is the dimension of search space;   ∈ [  ,   ],  = 1, 2, . . ., ), and its opposition point x = (x 1 , x2 , . . ., x ) is defined as follows: where  = 1, 2, . . ., ,  is the population size,  = 1, 2, . . ., ,  ∈ (0, 1) and  is a generalized coefficient, and [  ,   ] is the dynamic boundary of the th dimensional search space and can be obtained by the following formula: The fixed boundary is not conducive to preserve the search experience; thus we use dynamic boundary of the search space to replace the fixed boundary to preserve the search experience in order to make the opposition solution located in the search space which is narrowing.Moreover, if the operator of dynamic boundary makes x, jump out of [  ,   ], the following method can be used to reset x, : The EOBL generates the opposition population according to the elite individual and evaluates the current population and the elite population at the same time; in addition, it makes full use of the characteristics of the elite individuals to contain more useful search information than the ordinary individuals which improve the diversity of the population to certain extent.EOBL can enhance the ability of global exploration of WWO.

Local Neighborhood Search (LNS)
. WWO only performs the breaking operator on the new best solution to enhance local search around the best solution.In a bid to further enhance the local search ability to improve the convergence speed, the local neighborhood search [43] (LNS) model is added before the breaking operation.
The main idea of LNS is using the best solution found so far in a small neighborhood of the current solution rather than the entire population to update the current solution.The experience of an individual's neighborhood is considered when updating the individual's location, so that the graph of interconnections of them is called neighborhood structure.Suppose there is a WWO population  = ( 1 ,  2 , . . .,   ),   ( ∈ [1, ]) is a vector in the current population, and its dimension is .The indices of each vector are random in a bid to maintain the diversity of each neighborhood.Next, we can define the neighborhood of radius  ( is a nonzero integer and 2 + 1 < ), where  = 10, for each vector   ; that is to say, neighborhood of   consists of  − , . . .,   , . . .,  + .For analysis, we suppose that the vectors can be arranged into a ring topology according to their indices.Figure 2 illustrates the concept of local neighborhood model.In addition, the neighborhood topology is static and about the definition of the set of vectors all the time.LNS model is described in where   opt is the best vector in the neighborhood of   and ,  ∈ [ − ,  + ] ( ̸ =  ̸ = ) and  and  are the scaling factors, where ,  ∈ rand( ).In the improved version of WWO, the new best solution is updated according to (12), and the breaking operation is performed by the updated solution as where   is the best solution updated by LNS.

Improvement of WWO.
In original WWO, all water waves have to be propagated once at each generation, and Input: Define objective function (),  = ( 1 ,  2 , . . .,   ) Output: The best solution  * ; (1) Initialization: Initialize related parameters including , , , ,  max ,  min , , ,  and , initialize dynamic boundary of the search space, randomly initialize a population  of  waves (2) while stop criterion is not satisfied do (3) Update the current population with EOBL accordding to equation (9), equation (10) and equation ( 11); (4) for each  ∈  do (5) Propagate  to a new   based on equation ( 14); (6) if (  ) < () then (7) if (  ) < ( * ) then ( 8) Break  based on equation (12) and equation ( 13); (9) Update  * with   ; (10) Endif (11) Replace  with   ; ( 12) else (13) Decrease .ℎ by one; (14) if .ℎ = 0 then (15) Refract  to a new   based on equation ( 6) and equation ( 7); ( 16) endif ( 17) endif (18) Update the wavelength based on equation ( 4); ( 19) endfor ( 20) endwhile Algorithm 2: The framework of EOBWWO algorithm.the search behavior of each water wave is affected by the other waves in the group.Similar to PSO [1,45,46], an inertial weight is embedded into (2) in order to learn the past experience.Moreover, as shown in ( 2) and ( 4), propagation operator can make high fitness wave search small region and low fitness wave explore large region in global search process.In (2), the search step size is a random number fixed in the range [−1, 1], which is not very reasonable because search step size prefers fairly large step at the beginning in order to strengthen the probability of reaching the optimal regions, and, with the iteration going on, search step size should be decreased gradually to enhance the local exploitation ability.The random step size in the propagation process is improved by referencing to the method of [47]; the improved propagation operator is as follows: where  max and  min are, respectively, the maximum and minimum inertial weight, where  ∈ [0.4,1.5], iter is the current iteration number and iter max is the maximum iteration number, and  and  are two selected constants, where  = 0.02,  = 25.Improved propagation operator not only makes use of the past experience, but also makes the search step size decrease gradually with the iteration going on.The whole pseudocode of EOBWWO can be summarized as Algorithm 2.

Simulation Experiments and Results Analysis
In order to identify the effectiveness and efficiency of EOB-WWO, 20 standard test functions are applied in this section.
The detailed parameters about 20 benchmark functions [48,49] including functional form, scope, optimal solution, and the iterations are illustrated in Table 1.The 20 benchmark functions can be divided into three groups, unimodal functions ( 01 ∼ 06 ) as group 1, multimodal functions ( 07 ∼ 14 ) as group 2, and low-dimension functions ( 15 ∼ 20 ) as group 3.In the unimodal functions, the global optimum of  06 is located in a smooth, long, and narrow parabolic valley; when the traditional gradient optimization method is searched to the valley edge, it is difficult to carry out global optimization.However, it is very slow to change the value in the long and narrow area, which can be used to evaluate the performance of the algorithm.In the multimodal functions,  09 has many local minima, it is a typical nonlinear multimodal function, which has a wide range of search space, and it is generally considered to be a complex multimodal problem which is difficult to deal with.In general, unimodal functions are suitable for evaluating the exploitation; however, multimodal functions tend to be a good choice for evaluating exploration [50].The rest of this section is organized as follows: experimental setting is given in Section 4.1, experiment results of 30 dimensions and discussion are represented in Section 4.2, high-dimension test results including 100 dimensions, 1000 dimensions, and 10000 dimensions for some unimodal functions and multimodal functions are described in Section 4.3, and two design problems are shown in Section 4.4.

Experimental Setting and Comparative Methods.
The empirical analysis was conducted on a computer of Intel(R) with 3.5 GHz Xeon CPU and 8 GB of memory, the operating system is Windows 7, and the programs are written in Matlab 2012a.
The scope and dimension of variables have significant influence on the complexity of optimization.The scope of the benchmark function and the dimension of the lowdimension functions are illustrated in Table 1.The dimensions of unimodal functions and multimodal functions are, respectively, 30, 100, 1000, and 10000.
The performance of proposed EOBWWO algorithm is evaluated by comparing it to five state-of-the-art metaheuristic algorithms: ABC [4], CS [5], FPA [6], BA [8], and WWO [9]; the parameters settings of aforementioned algorithms are given in Table 2.

Experiment Results and Discussion
. In Tables 3-4, the dimension is 30, whereas the standard benchmark functions are listed in Table 1.In this paper, all function optimization experimental results of the algorithms are repeated 30 times to ensure the credibility in statistics.There are four evaluation indicators: max, min, median, and Std represent the worst fitness value, optimal fitness value, median of the test results, and standard deviation, respectively.The last column of Tables 3-5 gives the rank of the algorithms in terms of median values among the six algorithms.The minimum value, best median value, and the minimum standard deviation values among the six algorithms of each benchmark function are shown in bold.
Moreover, nonparametric Wilcoxon rank tests were conducted on the results of EOBWWO and other comparative algorithms on the 20 benchmark functions, and the test results are shown in Table 6, where the value of ℎ is 1 indicating that the performance and comparative method are statistically different with 95% confidence, and 0 implies that there is no statistical difference [9].
In Table 3, on the unimodal group, EOBWWO obtained the exact solution except function  06 and obtained the minimum standard deviation.Although ranking fifth on function  06 , standard deviation of EOBWWO is less than the other algorithms.All these mean that EOBWWO has a higher calculation precision and better stability in the optimization of the unimodal functions.
On group 2 of 8 multimodal functions, seeing from Table 4, EOBWWO can find the exact solution for  07 ,  09 ,  10 ,  11 ,  12 , and  13 , and the standard deviations of the five functions of EOBWWO are zeros.In addition, EOBWWO obtains the best median value on all functions.For functions  08 and  14 , the worst fitness value, best fitness value, median value, and standard deviations of EOBWWO are less than the other five algorithms.The multimodal functions are more complex than unimodal functions due to the local minima; thus the above analysis indicates that EOBWWO has a strong global search ability and higher calculation precision.
On group 3 of 6 low-dimension functions, results are illustrated in Table 5.For  15 , EOBWWO obtains the minimum values including the worst fitness, best value fitness, median value, and standard deviations among the comparative algorithms.For  16 , ABC, CS, WWO, and EOBWWO can find the exact solution, and the standard deviations of ABC and EOBWWO are zeros.For  17 , ABC, CS, WWO, and EOBWWO have the same median value and optimal fitness value, while the standard deviation of CS is minimal.For  18 , FPA can obtain the better fitness value and median, but the standard deviation of WWO is the best.For  19 , optimal fitness value, median value, and the standard deviation of FPA are better than those of EOBWWO.And for  20 , it is obvious that ABC, CS, WWO, and EOBWWO obtain the exact solution and the standard deviation of CS is better.Through the above analysis of Table 5, we can draw a conclusion that EOBWWO has certain advantages in dealing with lowdimension functions according to the experimental results.
In summary, as a result of introducing the three major optimization strategies in the improvement, the calculation precision of EOBWWO is better than the comparative algorithms for most benchmark functions.In addition to functions  06 , 2 )]  09 -    For unimodal function  01 , it can be seen obviously that the performance of EOBWWO outperforms the other comparative algorithms for dimensions 100, 1000, and 10000 from Table 7.With the increase of dimension, EOBWWO is still obtaining the exact solution and the standard deviation is zero on  01 .In the five comparison algorithms, as far as the median value changes in each dimension are concerned, the stability of FPA is better, but BA obtains the minimum value in each dimension.From the results of function  06 in Table 7, it is very easy to find that although the performance of EOBWWO is not very good when the dimension is 30, with the increase of dimension, not only does EOBWWO obtain the minimum value in each dimension, but also the standard deviation is the smallest and the range is not very large.In addition, for EOBWWO, the change of the order of magnitudes of median is the smallest in different dimensions.

Fitness function value
Those provide strong evidence that EOBWWO has higher performance in dealing with complex functions.Among the other comparative algorithms, the performance of BA is better, the second is CS, the third is WWO, and the fourth and fifth are BC and FPA, respectively.For multimodal functions, EOBWWO obtained the exact solution and the standard deviation is zero on  07 and  09 for the dimensions 100, 1000, and 10000.Taking into account median value, the performance of BA is better among comparative algorithms.However, considering the order of magnitudes of median, the difference between the five algorithms is not obvious.
Furthermore, some high dimensional tests of EOBWWO are also tested in functions  02 ,  04 ,  10 ,  12 , and  13 ; details of the experimental results are shown in Table 8.As shown in Tables 7 and 8 and the above analysis, EOBWWO has the ability efficiently and stably to handle high dimensional functions.

Structural Engineering Design
Examples.Many structural design problems in the real world are constrained optimization problems which are nonlinear with complex constraints and the optimal solution even does not exist in some cases.In order to evaluate the performance of EOBWWO even further, in this subsection, EOBWWO was used to solve two structural design problems: design of a compressing spring and design of a welded beam.

Test Problem 1: Design of a Tension/Compression Spring.
Design of a tension or compressing spring problem is introduced by Belegundu [14] firstly and it deals with the optimal design of tension/compression spring for a minimum weight.As shown in Figure 43, a tension/compression spring problem has three design variables: the wire diameter ( ), the mean coil diameter ( 2 ), and the number of active coils ( 3 ).The minimum weight is subject to constraints on minimum deflection, shear stress, surge frequency, and limits on outside diameter [15].A detailed description of the problem is as follows: Subject to: where the experimental parameters are set as follows: (0.05, 0.25, 2) ≤ ( 1 ,  2 ,  3 ) ≤ (1, 1.3, 15) [29].Table 9 lists the optimal solution for compression spring design obtained by EOBWWO.The results are 20 runs independently and the number of iterations of EOBWWO is 5000.
In the process of dealing with tension/compression spring  constrained optimization problem, first of all we need to determine whether the four constraints are satisfied.If these constraint conditions are all satisfied, then calculate () according to formula (15) and compare () with the original fitness values, the better result as fitness value of constrained optimization problem.Otherwise, the original fitness value remains and continues to iterate.In the process of using EOBWWO algorithm to deal with tension/compression spring constrained optimization problem, when the position of individual in current population is changed (create a new solution) and required to estimate the new solution, the steps of dealing with the tension/compression spring constrained optimization problems are as follows.

Fitness function value
Step 1. Calculate the values of the four constraint conditions (see ( 16)-( 19)) and estimate whether they are all satisfied with the constraint conditions.If these constraint conditions are all satisfied, go to Step 2; otherwise, go to Step 3.
Step 2. Calculate fitness value of new solution by formula (15), and compare new fitness value to original fitness value.According to the comparison results determine whether to update the current individual.Go to Step 4.
Step 3. Keeping the original individual that violates any constraint, go to Step 4.
Step 4. Continue performing the following operations.The reason behind keeping the individual that violates any constraint of the constraint conditions is that each iteration of the population has to implement the elite oppositionbased learning (EOBL) strategy and propagation operation; then the individual that violates any constraint may satisfy the constraints in the next iteration.As one of the most well-known design benchmark problems, many researchers have studied this problem.Belegundu [14] introduced this problem and used eight different  mathematical optimization techniques for this problem.Arora [15] solved this problem using a numerical optimization technique called a constraint correction at the constant cost.Table 10 summarized the optimal results of design of a tension/compression spring obtained by EOBWWO and other researchers.

Fitness function value
As seen from Table 10, the proposed method obtained the best design overall of 0.012665234 corresponding to  = (0.05169826, 0.356939073, 11.2760014) and the results    welded beam problem is to minimize overall cost of fabrication subject to constraints on shear stress , bending stress in the beam , buckling load on the bar   , end deflection of the beam , and side constraints.As depicted in Figure 44 [28], this problem consists of four design variables: thickness of the of using EOBWWO algorithm to deal with welded beam problem are as using EOBWWO algorithm to deal with compression spring.

Fitness function value
The results of comparing with those of other optimization algorithms reported in the literature are shown in Table 12.It can be seen from Table 12 remarkably that the proposed EOBWWO algorithm is much better than other algorithms in the design of welded beam, and the optimal solution obtained by EOBWWO is 1.69634711 corresponding to  = (0.205832588, 3.253654976, 9.0315042, 0.25962951).

Conclusions
In this paper, three strategies are added to the original WWO algorithm to improve the convergence speed and    strategy is introduced to enhance local exploitation capability via enhancing the local search around the promising optimal solution.In addition, the improved propagation operator provides the improved algorithm with a better balance between exploration and exploitation.By using the above-mentioned three strategies, EOBWWO can deal with function optimization including multimodal functions and significantly improve the convergence speed and calculation precision of the original WWO algorithm.There are various important issues for the further research topics of EOBWWO.On the one hand, structural design problems not only exist in the real world widely, but also are generally nonlinear and constrained optimization problems.Therefore,  other design problems can be resolved using EOBWWO in future research, such as multidimensional knapsack problem, permutation flow shop scheduling problem [51], and graph coloring problem.On the other hand, some improvements can be introduced to EOBWWO and WWO algorithm to enhance the ability of dealing with relevant problems.More elaborate set of parameters, such as breaking coefficient  and wavelength reduction coefficient , multiple population strategy, and combination with other optimization algorithms are some good choices.In addition, multiobjective optimization problems are also the focus in the future research.

Figure 1 :
Figure 1: Different wave shapes in deep and shallow water.

Table 1 :
18 , and  19 , the calculation precision of EOBWWO is inferior to ABC, FPA, and FPA, respectively.Moreover, EOBWWO can find the exact solution on  01 - 05 ,  07 ,  09 ,  11 ,  12 ,  13 , and  16 , and the standard deviations of these functions are zeros, which show the higher calculation precision and stronger stability of EOBWWO.EOBWWO obtains exact solution on multimodal functions  07 ,  09 ,  11 ,  12 , and  13 , which shows that EOBWWO has better global search performance.In order to show the performance of the EOBNWWO clearly, Figures 3-22 represent the convergence curves and Figures 23-42 describe the ANOVA test of global minimum of benchmark functions in Table 1.From Figures 3-22, obliviously, the convergence rate of EOBWWO is faster than other comparison algorithms including WWO on  01 - 07 , Benchmark test functions.
13, and  17 , and the exact solution of EOBWWO is obtained in some functions  01 - 05 ,  07 ,  09 - 13 , and  16 - 17 .All of these indicate that EOBWWO has a faster convergence speed and a higher calculation precision than the other comparative algorithms.Figures23-42show the ANOVA test of global minimum for  01 - 20 ; it can be easily found that the standard deviation of EOBWWO is much smaller for most functions and the standard deviation is even zero on some functions (e.g.,  01 - 05 ,  07 - 13 , and  15 - 16 ).Figures23-42imply that EOBWWO has strong stability.
01 ,  06 ) and two multimodal functions ( 07 ,  09 ) to test the 100 dimensions, 1000 dimensions, and 10000 dimensions, respectively, on the six algorithms.The results of all the algorithms about the four functions are summarized in Table7.The maximum numbers of iteration of each algorithm on each function are consistent with Table1.

Table 5 :
Experiment results of low-dimension functions for different algorithms.

Table 6 :
Statistical comparison between EOBWWO and the other five algorithms.

Table 8 :
Experiment results of high-dimension functions of EOBWWO.

Table 9 :
Statistical results of best tension/compression spring model obtained by EOBWWO.

Table 10 :
Best results of compression spring by different model.
a: mathematical optimization technique; b: numerical optimization technique; c: evolutionary programming; d: evolutionary computational technique; e: simulated annealing; f: evolution strategies; g: chaotic search, opposition-based learning, differential evolution, and quantum mechanics; h: TCA is the T-cell algorithm; i: interior search algorithm; NA: there is no relevant data.

Table 11 :
The optimal solution of the welded beam design example obtained by EOBWWO.

Table 12 :
The optimal solution of the welded beam design example using different methods.