The Differential Evolution (DE) algorithm is one of the most successful evolutionary computation techniques. However, its structure is not trivially translatable in terms of mathematical transformations that describe its population dynamics. In this work, analytical expressions are developed for the probability of enhancement of individuals after each application of a mutation operator followed by a crossover operation, assuming a population distributed radially around the optimum for the sphere objective function, considering the DE/rand/1/bin and the DE/rand/1/exp algorithm versions. These expressions are validated by numerical experiments. Considering quadratic functions given by f ( x ) = x T D T D x and populations distributed according to the linear transformation D - 1 of a radially distributed population, it is also shown that the expressions still hold in the cases when f ( x ) is separable ( D is diagonal) and when D is any nonsingular matrix and the crossover rate is C r = 1 . 0 . The expressions are employed for the analysis of DE population dynamics. The analysis is extended to more complex situations, reaching rather precise predictions of the effect of problem dimension and of the choice of algorithm parameters.

This paper proposes a local search optimizer that, employed as an additional operator in multiobjective evolutionary techniques, can help to find more precise estimates of the Pareto-optimal surface with a smaller cost of function evaluation. The new operator employs quadratic approximations of the objective functions and constraints, which are built using only the function samples already produced by the usual evolutionary algorithm function evaluations. The local search phase consists of solving the auxiliary multiobjective quadratic optimization problem defined from the quadratic approximations, scalarized via a goal attainment formulation using an LMI solver. As the determination of the new approximated solutions is performed without the need of any additional function evaluation, the proposed methodology is suitable for costly black-box optimization problems.