Next Article in Journal / Special Issue
An Efficient Framework for Multi-Objective Risk-Informed Decision Support Systems for Drainage Rehabilitation
Previous Article in Journal
Convolutional Neural Network Based Ensemble Approach for Homoglyph Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Evolution in Robust Optimization Over Time Using a Survival Time Approach

1
Artificial Intelligence Research Center, University of Veracruz, Sebastián Camacho 5, Col. Centro, Xalapa 91000, Veracruz, Mexico
2
National Laboratory on Advanced Informatics, Rebsamen 80, Xalapa 91000, Veracruz, Mexico
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2020, 25(4), 72; https://0-doi-org.brum.beds.ac.uk/10.3390/mca25040072
Submission received: 2 October 2020 / Revised: 22 October 2020 / Accepted: 23 October 2020 / Published: 26 October 2020
(This article belongs to the Special Issue Numerical and Evolutionary Optimization 2020)

Abstract

:
This study presents an empirical comparison of the standard differential evolution (DE) against three random sampling methods to solve robust optimization over time problems with a survival time approach to analyze its viability and performance capacity of solving problems in dynamic environments. A set of instances with four different dynamics, generated by two different configurations of two well-known benchmarks, are solved. This work also introduces a comparison criterion that allows the algorithm to discriminate among solutions with similar survival times to benefit the selection process. The results show that the standard DE holds a good performance to find ROOT solutions, improving the results reported by state-of-the-art approaches in the studied environments. Finally, it was found that the chaotic dynamic, disregarding the type of peak movement in the search space, is a source of difficulty for the proposed DE algorithm.

1. Introduction

Optimization is an inherent process in various areas of study and everyday life. The search to improve processes, services, and performances has originated in different solution techniques. However, there are problems in which uncertainty is present over time, given that the solution’s environment can change at a specific time. These types of problems are named Dynamic Optimization Problems (DOPs) [1]. This study deals with dynamic problems where the environment of the problem changes over time. Various studies have been carried out to resolve DOPs through tracking moving optima (TMO), which is characterized by the search for and implementation of the global optimal-solution every time the environment changes [2,3,4].
Evolutionary algorithms, such as Differential Evolution (DE), have shown good performance to solve tracking problems [5,6,7]. However, the search and implementation of the optimum each time the environment changes may not be feasible due to different circumstances, such as time or cost.
The approach introduced in [8] tries to solve DOPs through a procedure known as robust optimization over time (ROOT). ROOT seeks to solve DOPs by looking for a good solution for multiple environments and preserve it for as long as possible, while its quality does not decrease from a pre-established threshold. The solution found is called Robust Solution Over Time (RSOT).
In this regard, Fu et al. [9] introduced different measures to characterize environmental changes. After that, the authors developed two definitions for robustness [10]. The first one was based on “Survival Time”—when a solution is considered acceptable (an aptitude threshold must be previously defined). The second definition was based on the “Average Fitness”—the solution’s average fitness is maintained during a previously defined time window. The measurements incorporate information on the concepts of robustness and consider the values of error estimators. An algorithm performance measure was suggested to find ROOT solutions. The study was carried out using a modified version of the moving peaks benchmark (mMPB) [10].
Jin et al. [11] proposed a ROOT framework that includes three variants of the Particle Swarm Optimization algorithm (PSO). PSO with a simple restart strategy (sPSO), PSO with memory scheme (memPSO), and a variant that implements the species technique SPSO. The authors applied a radial basis function as an approximation and an autoregressive model as a predictor.
On the other hand, Wang introduced the concept of robustness in a multi-objective environment, where a framework is created to find robust Pareto fronts [12]. The author adopted the dynamic multi-objective evolutionary optimization algorithm in the experiments. At the same time, Huang et al. considered the cost derived from implementing new solutions, thus addressing the ROOT problem by a multi-objective PSO (MOPSO); The Fu metric was applied in that study [13].
Yazdani et al. introduced a new semi-ROOT algorithm that looks for a new solution when the current one is not acceptable, or if the current one is acceptable but the algorithm finds a better solution, whose implementation is preferable even with the cost of change [14].
Novoa-Hernández, Pelta, and Corona analyzed the ROOT behavior using some approximation models [15]. The authors suggested that the radial basis network model with radial basis function works better for problems with a low number of peaks. However, considering all the scenarios, the SVM model with Laplace Kernel shows notably better performance to those compared in the tests carried out.
Novoa-Hernández and Amilkar in [16] reviewed different relevant contributions to ROOT. The authors analyzed papers hosted in the SCOPUS database. Concerning new methods to solve ROOT problems, Yazdani, Nguyen, and Branke proposed a new framework using a multi-population approach where sub-populations track peaks and collect information from them [17]. Adam and Yao introduced three methods to address ROOT (Mesh, Optimal in time, and Robust). The authors mentioned that they significantly improves the results obtained for ROOT in the state-of-the-art [18]. Fox, Yang and Caraffini studied different prediction methods in ROOT, including the Linear and Quadratic Regression methods, an Autoregressive model, and Support Vector Regression [19]. Finally, Liu and Liang mapped a ROOT approach to minimize the electric bus transit center’s total cost in the first stage [20].
In different studies, DE has been used to solve ROOT problems using the Average Fitness approach, achieving competitive results [21,22]. However, to the best of the authors’ knowledge, there are no studies that determine DE’s performance in solving ROOT problems with the Survival Time approach, and this is where this work precisely focuses. This research aims to present an empirical comparison of the standard DE against three random sampling methods to solve robust optimization over time problems with a survival time approach to analyze its viability and performance capacity of solving problems in dynamic environments.
The paper is organized as follows: Section 2 includes ROOT’s definition under a survival time approach, while in Section 3 the implemented methods based on random sampling are detailed. Section 4 details the standard differential evolution and the objective function used by the algorithm in the present study. Section 5 specifies the benchmark problems to be solved. After that, Section 6 specifies the experimental settings and Section 7 shows the results obtained. Finally, Section 8 summarizes the conclusions and future work.

2. Survival Time Approach

Under this approach, a threshold is predefined to specify the quality that a solution must have to be considered good or suitable to survive. Once the threshold is defined, the search begins for a solution whose fitness can remain above the threshold in as many environments as possible. In this sense, the solution is maintained until its quality does not meet the predefined expectations, and then new robust solution over time must be sought.
In Equation (1), the function F s ( x , t , V ) to calculate the survival time fitness of a solution x at time t is detailed. It measures the number of environments that a solution remains in above the threshold V.
F s ( x , t , V ) = 0 if f t ( x ) < V 1 + max { l | i { t , t + 1 , , t + l } : f i ( x ) V } , in other case

3. Random Sampling Methods

In the study presented in [18] the authors proposed three random sampling methods to solve ROOT problems, with a better performance against the state-of-the-art algorithms. The methods are described below.
In all three methods, the best solution should be searched in the current time’s solutions space, modifying the solution space when the “Robust method” is used and then using that solution next time according to the approach used (Survival Time or Average Fitness).

3.1. Mesh

This method performs random sampling in the current search space, then uses the sample with the best fitness in the current environment as a robust solution over time, using the solution found in the following times (Figure 1).

3.2. Time-Optimal

This method performs a search similar to the “Mesh” method, with the difference that the best solution found being improved using a local search (Figure 2).

3.3. Robust

This method performs a search similar to the “Mesh” method, differing in that a smoothing preprocessing of the solution space is performed before the search process. As seen in Figure 3, the solution obtained (green dot) can vary concerning the solution with better suitability in the raw environment (yellow dot).

4. Differential Evolution

In 1995, Storm and Price proposed an evolutionary algorithm to solve optimization problems in continuous spaces. DE is based on a population of solutions that, through simple recombination and mutation, evolves, thus improving individuals’ fitness [23].
Considering the fact that this work, to the best of the authors’ knowledge, is the first attempt to study DE in this type of ROOT problems (i.e., survival time), and also taking into account that the most popular DE variant (DE/rand/1/bin) has provided competitive results in ROOT problems under an average fitness approach [21,22], the algorithm used in this study is precisely the most popular variant known as DE/rand/1/bin, where ’‘rand” (random) refers to the base vector used in the mutation, ‘’bin” (binomial) refers to the crossover type used, and 1 means one vector difference computed.
The algorithm starts by randomly generating a uniformly distributed population x i , G i = 1 , , N P , where NP is the number of individuals for each generation “G”.
After that, the algorithm enters an evolution cycle until the stop condition is reached. We applied the maximum number of evaluations allowed “ M A X E v a l ” as the stop condition.
Subsequently, to adapt individuals, the algorithm performs recombination, mutation, and the replacement of each one in the current generation. One of the most popular mutations is DE/rand/1 in Equation (2), where r 0 r 1 r 2 i are the indices of individuals randomly chosen from the population, 1 is the number of differences used in the mutation and F > 0 is the scale factor.
v i , G = x r 0 , G + F ( x r 1 , G x r 2 , G )
The vector obtained v i , G is known as mutant vector, which is recombined with the target (parent) vector by binomial crossover, as detailed in Equation (3).
u i , j , G = v j , i , G , i f r a n d j C R o r j = j r a n d x j , i , G , o t h e r w i s e
In this study, the elements of the child vector (also called trial) u i , j , G are limited according to the pre-established maximum and minimum limits, also known as boundary constraints. Based on the study in [24], we use the boundary method (see Algorithm 1, line 14). In the selection process, the algorithm determines the vector that will prevail for the next generation between parent (target) and child (trial), as expressed in Equation (4).
x i , G + 1 = u i , G , i f f u i , G f x i , G x i , G , o t h e r w i s e
Algorithm 1: “DE/rand/1/bin” Algorithm. NP, M A X E v a l , CR and F are parameters defined by the user. D is the dimension of the problem.
Mca 25 00072 i001
In Equation (1), the function to obtain an individual’s fitness through the survival time approach is shown. However, the fitness obtained is not enough to differentiate similar individuals, i.e., individuals who have survived the same amount of environments. That is why, in the implemented algorithm, we consider an additional calculation to help identify better solutions.
We propose to obtain the average of the solution’s quality in the environments that have survived. This average value helps to differentiate solutions with similar survival times. Therefore, the objective function now considers both, the number of surviving environments and the performance achieved by this solution in those environments that it has survived.
Considering the fact that the maximum height of the peaks is defined at 70 (see Table 3), the objective function for the implemented algorithm is given by the result obtained in Equation (1) multiplied by 100 plus the average fitness of the solution throughout the environments it has survived.

5. Benchmark Problems

The problems tackled in this study are based on Moving Peaks Benchmark (MPB) [25] and are configured in a similar way to that used in various specialized literature publications on ROOT, and specifically as used in [18].
Two modified MPBs can be highlighted, which are described in the following subsections. The dynamics used are presented in Table 1, where Δ ϕ is the increment from time t to time t + 1 of the ϕ parameter.

5.1. Moving Peaks Benchmark 1 (MPB1)

In this benchmark, environments with conical peaks of height h ( t ) [ h m i n , h m a x ] , width w ( t ) [ w m i n , w m a x ] and center c ( t ) [ x m i n , x m a x ] are generated, where the design variable x is bounded in [ x m i n , x m a x ] . The function to generate the environment is expressed in Equation (5), where the dynamic function for height and width is given as in Table 1, while the center moves according to Equation (6). r i follows an uniform distribution of a D-dimensional sphere of radius s i , and λ [ 0 , 1 ] is a fixed parameter.
f x , a ( t ) = m a x i = 1 i = m { h i ( t ) w i ( t ) x c i ( t ) l 2 }
c i ( t + 1 ) = c i ( t ) + v i ( t + 1 ) v i ( t + 1 ) = s i ( 1 λ ) r i ( t + 1 ) + λ v i ( t ) ( 1 λ ) r i ( t + 1 ) + λ v i ( t )
In the present study, two problems generated by this benchmark are solved, with λ = 0 it implies that the movement of the peaks is random, while with λ = 1 it implies that the movement is constant in the direction v i ( t ) .

5.2. Moving Peaks Benchmark 2 (MPB2)

The set of test functions in this benchmark is described in Equation (7), where a ( t ) is the environment at time step t, h i ( t ) , w i ( t ) , c i ( t ) is the height, width and center of the i-th peak function at time t, respectively; x is the decision variable and m is the total number of peaks. h i ( t + 1 ) and w i ( t + 1 ) vary according to Table 1. An additional technique that uses a rotation matrix is used to rotate the centers [25].
f x , a ( t ) = 1 d j = 1 d m a x i = 1 i = m { h i ( t ) w i ( t ) x c i ( t ) }

6. Experimental Settings

Based on the information in Section 5, different environments are generated as test problems and they are summarized in Table 2.
The parameter settings of the problems are detailed in Table 3.
The height and width of the peaks were randomly initialized in the predefined ranges. The centers were randomly initialized within the solution space.
A survival threshold V = 50 is selected, representing the most difficult cases that have been resolved in the literature under the survival approach. The higher the survival threshold, the more difficult it is to find solutions that satisfy multiple scenarios.
The DE parameters were fine-tuned using the iRace tool [26] and they are summarized in Table 4, where NP is the population size, CR is the crossover parameter, and F is the scale factor.
For each problem, a solution is sought at each time t ( 2 , , 100 ) .
In order to evaluate an RSOT in a specific time, approximate and predictive methods have been used in the literature so that the performance of an algorithm depends on their accuracy. However, in this study, we want to know the DE behavior when solving the ROOT environments considering they had ideal predictors to evaluate the solutions. In this regard, the process to study the algorithm’s ability to find RSOT using DE at each instant of time is as follows:
  • A solution is sought according to the algorithm described in Section 4, and the measured solution value by Equation (1) is recorded.
  • Subsequently, to obtain the algorithm’s performance in the following environment, the search process is performed again using the real-environments; the best solution found is newly measured by Equation (1) and is also recorded.
  • The described procedure is carried out at each instant of time that is being recorded. Therefore, in the present study, it is not necessary to detect environmental changes to know at what point in time a solution is no longer considered good. Each time a solution is sought, the algorithm initializes its population randomly, avoiding diversity problems.

7. Results and Discussion

The results for the problems generated with dynamics 1–4 are detailed in Table 7 and graphically shown in Figures 6–9, for each one of the four dynamics. In all four figures, those labeled with (a) and (b) present the results obtained in the MPB1 problems, while those labeled with (c) and (d) refer to the MPB2 problems. In all cases, the average survival values obtained by Mesh, time-optimal and robust approaches are compared against DE.
Non-parametric statistical tests [27] were applied to the corresponding numerical results presented in Table 7. The 95 % -confidence Kruskal–Wallis and 95 % -confidence Friedman tests were applied and their obtained p-values are reported in Table 5.
To further determine differences among the compared algorithms, the 95 % -confidence Wilcoxon test was applied to pair-wise comparisons for each problem instance. The obtained p-values are reported in Table 6, where the significant improvement with a significance level α = 0.5 is shown in boldface. We can observe that the Wilcoxon test confirmed significant differences obtained in Kruskal–Wallis and Friedman Tests, most of them comparing DE/rand/1/bin versus random sampling methods, with the exception of four corresponding to problems generated with dynamic 4 (B1D4-1 and B1D4-2 both, in the comparison of DE/rand/1/bin versus Mesh and DE/rand/1/bin versus Time-optimal).
Table 7 summarizes the mean and standard deviation statistical results obtained by the compared algorithms. It can be seen that DE/rand/1/bin obtains the highest average values in all the problems that were solved. Nevertheless, the higher standard deviation values obtained by DE/rand/1/bin in problems B1D4-1 and B1D4-2 confirm those expressed by the non-parametric tests—the differences are not significant with respect to the random sampling methods. Figure 4 and Figure 5 have the box-plots for B1D4-1 and B1D4-2, where all the compared algorithms reach survival times between 1 and 3, with the exception of the Robust approach in B1D4-1, but such a difference was not significant.
With respect to the graphical results, when MPB1 (items (a) and (b)) is compared against MPB2 (items (c) and (d)) in Figure 6, Figure 7, Figure 8 and Figure 9, it is clear that MPB1 is more difficult to solve by all four approaches. However, in all cases (MPB1 and MPB2 in the four dynamics) DE is able to provide better results against the three other algorithms. Such a performance is more evident in all MPB2 instances.
Regarding MPB1 (items (a) and (b)), it is important to note that it is more difficult to find higher survival times when the peak movement is random (items (a), where λ = 0).
Another interesting behavior found is that all four compared methods are affected mainly by the random and chaotic dynamics in those MPB1 instances, the latter one being the most complex (chaotic dynamic). However, even in such a case DE was able to match and in some cases improve the survival values of the compared approaches. This source of difficulty now found motivates part of our future research.

8. Conclusions

A performance analysis of the differential evolution algorithm, with one of its original variants, called DE/rand/1/bin, when solving robust optimization over time problems with a survival time approach, was presented in this paper. Three state-of-the-art random sampling methods to solve ROOT problems were used for comparison purposes. Sixteen generated problems by two benchmarks with two configurations and four different dynamics were solved. The solutions generated by the DE were obtained using the real environments without prediction mechanisms with the aim to analyze its behavior in ideal conditions. The findings supported by the obtained results indicate that DE is a suitable algorithm to deal with this type of dynamic search space when a survival time approach is considered. Moreover, the additional criterion that was added to the DE objective function allowed the algorithm to better discriminate between similar solutions in terms of survival time. Furthermore, it was found that the combination of a chaotic dynamic with both, random and constant peak movements, is a source of difficulty that requires further analysis.
This last finding is the starting point of our future research, where more recent DE variants, such as DE/current-to-p-best, will be tested in those complex ROOT instances. Moreover, the effect of predictors in DE-based approaches will be studied.

Author Contributions

Conceptualization, J.-Y.G.-G.; methodology, E.M.-M. and J.-Y.G.-G.; software, J.-Y.G.-G.; data curation, J.-Y.G.-G.; investigation, J.-Y.G.-G.; formal analysis, J.-Y.G.-G. and E.M.-M.; validation, E.M.-M. and S.D.-I.; writting—original draft preparation, J.-Y.G.-G., E.M.-M. and S.D.-I.; writing—review and editing, S.D.-I., E.M.-M. and J.-Y.G.-G. All authors have read and agreed to the published version of the manuscript.

Funding

The first author acknowledges support from the Mexican National Council of Science and Technology (CONACyT) through a scholarship to pursue graduate studies at the University of Veracruz.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DEDifferential Evolution
MPB1Moving Peaks Benchmark 1
MPB2Moving Peaks Benchmark 2
ROOTRobust Optimization over Time
RSOTRobust Solution over Time
S.D.Standard deviation

References

  1. Nguyen, T.T.; Yang, S.; Branke, J. Evolutionary dynamic optimization: A survey of the state of the art. Swarm Evolut. Comput. 2012, 6, 1–24. [Google Scholar] [CrossRef]
  2. Dang, D.C.; Jansen, T.; Lehre, P.K. Populations Can Be Essential in Tracking Dynamic Optima. Algorithmica 2017, 78, 660–680. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, S.; Li, C. A clustering particle swarm optimizer for locating and tracking multiple optima in dynamic environments. IEEE Trans. Evolut. Comput. 2010, 14, 959–974. [Google Scholar] [CrossRef] [Green Version]
  4. Yang, S.; Yao, X. Evolutionary Computation for Dynamic Optimization Problems; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar] [CrossRef] [Green Version]
  5. Das, S.; Mullick, S.S.; Suganthan, P. Recent advances in differential evolution—An updated survey. Swarm Evolut. Comput. 2016, 27, 1–30. [Google Scholar] [CrossRef]
  6. Lin, L.; Zhu, M. Efficient Tracking of Moving Target Based on an Improved Fast Differential Evolution Algorithm. IEEE Access 2018, 6, 6820–6828. [Google Scholar] [CrossRef]
  7. Zhu, Z.; Chen, L.; Yuan, C.; Xia, C. Global replacement-based differential evolution with neighbor-based memory for dynamic optimization. Appl. Intell. 2018, 48, 3280–3294. [Google Scholar] [CrossRef]
  8. Yu, X.; Jin, Y.; Tang, K.; Yao, X. Robust optimization over time; A new perspective on dynamic optimization problems. In Proceedings of the IEEE Congress on Evolutionary Computation, Barcelona, Spain, 18–23 July 2010; pp. 1–6. [Google Scholar] [CrossRef] [Green Version]
  9. Fu, H.; Sendhoff, B.; Tang, K.; Yao, X. Characterizing environmental changes in Robust Optimization Over Time. In Proceedings of the IEEE Congress on Evolutionary Computation, Brisbane, QLD, Australia, 10–15 June 2012; pp. 1–8. [Google Scholar] [CrossRef]
  10. Fu, H.; Sendhoff, B.; Tang, K.; Yao, X. Finding Robust Solutions to Dynamic Optimization Problems. In Applications of Evolutionary Computation; Esparcia-Alcázar, A.I., Ed.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 616–625. [Google Scholar]
  11. Jin, Y.; Tang, K.; Yu, X.; Sendhoff, B.; Yao, X. A framework for finding robust optimal solutions over time. Memetic Comput. 2013, 5, 3–18. [Google Scholar] [CrossRef]
  12. Wang, H.L.G.C. The Evolutionary Algorithm to Find Robust Pareto-Optimal Solutions over Time. Math. Probl. Eng. 2014, 2014, 814210. [Google Scholar]
  13. Huang, Y.; Ding, Y.; Hao, K.; Jin, Y. A multi-objective approach to robust optimization over time considering switching cost. Inf. Sci. 2017, 394–395, 183–197. [Google Scholar] [CrossRef]
  14. Yazdani, D.; Branke, J.; Omidvar, M.N.; Nguyen, T.T.; Yao, X. Changing or Keeping Solutions in Dynamic Optimization Problems with Switching Costs. In Proceedings of the Genetic and Evolutionary Computation Conference, Kyoto, Japan, 15–19 July 2018; pp. 1095–1102. [Google Scholar] [CrossRef]
  15. Novoa-Hernández, P.; Pelta, D.A.; Corona, C.C. Approximation Models in Robust Optimization Over Time—An Experimental Study. In Proceedings of the 2018 IEEE Congress on Evolutionary Computation (CEC), Rio de Janeiro, Brazil, 8–13 July 2018; pp. 1–6. [Google Scholar] [CrossRef]
  16. Novoa-Hernández, P.; Puris, A. Robust optimization over time: A review of most relevant contributions [Optimización robusta en el tiempo: Una revisión de las contribuciones más relevantes]. Rev. Iber. Sist. Tecnol. Inf. 2019, 2019, 156–164. [Google Scholar]
  17. Yazdani, D.; Nguyen, T.T.; Branke, J. Robust Optimization Over Time by Learning Problem Space Characteristics. IEEE Trans. Evolut. Comput. 2019, 23, 143–155. [Google Scholar] [CrossRef] [Green Version]
  18. Adam, L.; Yao, X. A Simple Yet Effective Approach to Robust Optimization Over Time. In Proceedings of the 2019 IEEE Symposium Series on Computational Intelligence (SSCI), Xiamen, China, 6–9 December 2019; pp. 680–688. [Google Scholar]
  19. Fox, M.; Yang, S.; Caraffini, F. An Experimental Study of Prediction Methods in Robust optimization Over Time. In Proceedings of the 2020 IEEE Congress on Evolutionary Computation (CEC), Glasgow, UK, 19–24 July 2020; pp. 1–7. [Google Scholar]
  20. Liu, Y.; Liang, H. A ROOT Approach for Stochastic Energy Management in Electric Bus Transit Center with PV and ESS. In Proceedings of the 2019 IEEE Global Communications Conference (GLOBECOM), Waikoloa, HI, USA, 9–13 December 2019; pp. 1–6. [Google Scholar]
  21. Guzmán-Gaspar, J.; Mezura-Montes, E. Differential Evolution Variants in Robust Optimization Over Time. In Proceedings of the 2019 International Conference on Electronics, Communications and Computers (CONIELECOMP), Cholula, Mexico, 27 February–1 March 2019; pp. 164–169. [Google Scholar]
  22. Guzmán-Gaspar, J.; Mezura-Montes, E. Robust Optimization Over Time with Differential Evolution using an Average Time Approach. In Proceedings of the 2019 IEEE Congress on Evolutionary Computation (CEC), Wellington, New Zealand, 10–13 June 2019; pp. 1548–1555. [Google Scholar]
  23. Price, K.; Storn, R.M.; Lampinen, J.A. Differential Evolution a Practical Approach to Global Optimization, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  24. Juárez-Castillo, E.; Acosta-Mesa, H.G.; Mezura-Montes, E. Adaptive boundary constraint-handling scheme for constrained optimization. Soft Comput. 2019, 23, 8247–8280. [Google Scholar] [CrossRef]
  25. Li, C.; Yang, S.; Nguyen, T.T.; Yu, E.L.; Yao, X.; Jin, Y.; Beyer, H.G.; Suganthan, P.N. Benchmark Generator for CEC 2009 Competition on Dynamic Optimization. Available online: https://bura.brunel.ac.uk/bitstream/2438/5897/2/Fulltext.pdf (accessed on 24 October 2020).
  26. López-Ibáñez, M.; Dubois-Lacoste, J.; Pérez Cáceres, L.; Birattari, M.; Stützle, T. The irace package: Iterated racing for automatic algorithm configuration. Oper. Res. Perspect. 2016, 3, 43–58. [Google Scholar] [CrossRef]
  27. García, S.; Molina, D.; Lozano, M.; Herrera, F. A Study on the Use of Non-Parametric Tests for Analyzing the Evolutionary Algorithms’ Behaviour: A Case Study on the CEC’2005 Special Session on Real Parameter Optimization. J. Heuristics 2009, 15, 617–644. [Google Scholar] [CrossRef]
Figure 1. Mesh method. The yellow point is the robust solution found by the method and it is used in the next times.
Figure 1. Mesh method. The yellow point is the robust solution found by the method and it is used in the next times.
Mca 25 00072 g001
Figure 2. Time-optimal method. The green point is the robust solution found by the method, which is used in the next times.
Figure 2. Time-optimal method. The green point is the robust solution found by the method, which is used in the next times.
Mca 25 00072 g002
Figure 3. Robust method. The green point is the robust solution found by the method which is used in the next times.
Figure 3. Robust method. The green point is the robust solution found by the method which is used in the next times.
Mca 25 00072 g003
Figure 4. Boxplot of the results obtained in B1D4-1.
Figure 4. Boxplot of the results obtained in B1D4-1.
Mca 25 00072 g004
Figure 5. Boxplot of the results obtained in B1D4-2.
Figure 5. Boxplot of the results obtained in B1D4-2.
Mca 25 00072 g005
Figure 6. Results obtained using dynamic 1.
Figure 6. Results obtained using dynamic 1.
Mca 25 00072 g006
Figure 7. Results obtained using dynamic 2.
Figure 7. Results obtained using dynamic 2.
Mca 25 00072 g007
Figure 8. Results obtained using dynamic 3.
Figure 8. Results obtained using dynamic 3.
Mca 25 00072 g008
Figure 9. Results obtained using dynamic 4.
Figure 9. Results obtained using dynamic 4.
Mca 25 00072 g009
Table 1. Dynamic Functions.
Table 1. Dynamic Functions.
1. Small Step Δ ϕ = γ · ϕ · r · ϕ s e v e r i t y
2. Large Step Δ ϕ = ϕ · ( γ * s i g n ( r ) + ( γ m a x γ ) · r ) · ϕ s e v e r i t y
3. Random Δ ϕ = N ( 0 , 1 ) · ϕ s e v e r i t y
4. Chaotic ϕ t + 1 = ϕ m i n + A · ( ϕ t ϕ m i n ) · ( 1 ( ϕ t ϕ m i n ) / ϕ )
Table 2. Summary of problems.
Table 2. Summary of problems.
BenchmarkAbbreviationConfigurationDynamic ( δ )
MPB1 B 1 D δ 1 λ = 0 {1,2,3,4}
MPB1 B 1 D δ 2 λ = 1 {1,2,3,4}
MPB2 B 2 D δ a uniform start of peak distribution{1,2,3,4}
MPB2 B 2 D δ b random start of peak distribution{1,2,3,4}
Table 3. Parameters settings.
Table 3. Parameters settings.
ParameterMPB1MBP2
Number of peaks m525
Number of dimensions d22
Search range [ x m i n , x m a x ] [ 0 , 50 ] [ 25 , 25 ]
Height range [ h m i n , h m a x ] [ 30 , 70 ] [ 30 , 70 ]
Width range [ w m i n , w m a x ] [ 1 , 12 ] [ 1 , 13 ]
Angle range [ θ m i n , θ m a x ] - [ π , π ]
h e i g h t s e v e r i t y U ( 1 , 10 ) 5.0
w i d t h s e v e r i t y U ( 0.1 , 1 ) 0.5
a n g l e s e v e r i t y -1.0
Initial h50 U ( h m i n , h m a x )
Initial w6 U ( w m i n , w m a x )
Initial Angle-0
λ { 0 , 1 } -
Number of dimensions for rotation l r -2
Computational budget at each step Δ e 25002500
Table 4. Parameter settings of DE.
Table 4. Parameter settings of DE.
NPCRF
540.530.73
Table 5. Results of the 95 % -confidence Kruskal–Wallis (KW) and Friedman (F) tests. The symbol (*) after letter ‘’D” in the Problem column refers to the type of dynamic used according to columns Dynamic. A p-value less than 0.05 means that there are significant differences among the compared algorithms in such problems.
Table 5. Results of the 95 % -confidence Kruskal–Wallis (KW) and Friedman (F) tests. The symbol (*) after letter ‘’D” in the Problem column refers to the type of dynamic used according to columns Dynamic. A p-value less than 0.05 means that there are significant differences among the compared algorithms in such problems.
Problem Instancep-Value
Dynamic 1Dynamic 2Dynamic 3Dynamic 4
KWFKWFKWFKWF
B1D*-1<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
B1D*-2<0.0001<0.0001<0.0001<0.0001<0.0001<0.00010.0082<0.0001
B2D*-a<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
B2D*-b<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
Table 6. Results of the 95 % -confidence Wilcoxon signed-rank test. A p-value less than 0.05 means that exists significant differences.
Table 6. Results of the 95 % -confidence Wilcoxon signed-rank test. A p-value less than 0.05 means that exists significant differences.
ProblemAlgorithmp-Value
Dynamic 1Dynamic 2Dynamic 3Dynamic 4
B1D*-1DE/rand/1/bin versus Mesh<0.0001<0.0001<0.00010.2545
DE/rand/1/bin versus Time-optimal<0.0001<0.0001<0.00010.2881
DE/rand/1/bin versus Robust<0.0001<0.0001<0.0001<0.0001
Mesh versus Time-optimal0.11670.00190.04830.8758
Mesh versus Robust<0.00010.00060.1019<0.0001
Time-optimal versus Robust<0.0001<0.00010.0007<0.0001
B1D*-2DE/rand/1/bin versus Mesh<0.0001<0.0001<0.00010.7724
DE/rand/1/bin versus Time-optimal<0.0001<0.0001<0.00010.7531
DE/rand/1/bin versus Robust<0.0001<0.0001<0.00010.0037
Mesh versus Time-optimal0.96020.93580.88520.9984
Mesh versus Robust0.09870.90780.00870.0067
Time-optimal versus Robust0.1040.85580.01260.0069
B2D*-aDE/rand/1/bin versus Mesh<0.0001<0.0001<0.0001<0.0001
DE/rand/1/bin versus Time-optimal<0.0001<0.0001<0.0001<0.0001
DE/rand/1/bin versus Robust<0.0001<0.0001<0.0001<0.0001
Mesh versus Time-optimal0.87960.60440.84790.9221
Mesh versus Robust0.55020.9639<0.00010.761
Time-optimal versus Robust0.42220.5658<0.00010.7259
B2D*-bDE/rand/1/bin versus Mesh<0.0001<0.0001<0.0001<0.0001
DE/rand/1/bin versus Time-optimal<0.0001<0.0001<0.0001<0.0001
DE/rand/1/bin versus Robust<0.0001<0.0001<0.0001<0.0001
Mesh versus Time-optimal0.79370.98050.62420.9028
Mesh versus Robust0.72310.69220.00010.3679
Time-optimal versus Robust0.87860.67610.00140.3369
Table 7. Statistical Results obtained by each Algorithm in each one of the problem instances. Best statistical results are marked with boldface.
Table 7. Statistical Results obtained by each Algorithm in each one of the problem instances. Best statistical results are marked with boldface.
ProblemAlgorithmDynamic 1Dynamic 2Dynamic 3Dynamic 4
MeanS.D.MeanS.D.MeanS.D.MeanS.D.
B1D*-1DE/rand/1/bin10.14381.50566.13060.681610.65811.49731.56781.7595
Mesh4.93370.77722.97720.34485.68440.88940.94290.8177
Time-optimal5.05270.77563.03650.35035.79790.88570.95690.8728
Robust4.48190.78532.89470.35775.53710.98370.31760.2515
B1D*-2DE/rand/1/bin13.56333.8668.6232.138813.83943.69421.82191.9531
Mesh10.3294.05755.61621.861910.58873.94881.72691.8693
Time-optimal10.33094.03285.61041.846110.57313.9231.72351.8584
Robust10.14373.92415.65331.816610.22463.74131.28431.4648
B2D*-aDE/rand/1/bin19.17570.123619.84740.039919.3810.074818.01170.3944
Mesh6.5960.72586.51540.266.39560.31992.75580.229
Time-optimal6.59140.74136.53210.27136.41160.31772.75930.2315
Robust6.67090.79046.51090.26116.19470.25932.74110.2276
B2D*-bDE/rand/1/bin16.99540.400518.13360.114417.71870.134513.19930.8169
Mesh5.57290.85285.23040.32394.87880.25442.21920.3461
Time-optimal5.58250.8365.2330.32274.8680.24632.22380.3472
Robust5.64790.9815.2150.32664.75780.23592.18180.334
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guzmán-Gaspar, J.-Y.; Mezura-Montes, E.; Domínguez-Isidro, S. Differential Evolution in Robust Optimization Over Time Using a Survival Time Approach. Math. Comput. Appl. 2020, 25, 72. https://0-doi-org.brum.beds.ac.uk/10.3390/mca25040072

AMA Style

Guzmán-Gaspar J-Y, Mezura-Montes E, Domínguez-Isidro S. Differential Evolution in Robust Optimization Over Time Using a Survival Time Approach. Mathematical and Computational Applications. 2020; 25(4):72. https://0-doi-org.brum.beds.ac.uk/10.3390/mca25040072

Chicago/Turabian Style

Guzmán-Gaspar, José-Yaír, Efrén Mezura-Montes, and Saúl Domínguez-Isidro. 2020. "Differential Evolution in Robust Optimization Over Time Using a Survival Time Approach" Mathematical and Computational Applications 25, no. 4: 72. https://0-doi-org.brum.beds.ac.uk/10.3390/mca25040072

Article Metrics

Back to TopTop