Populations Can Be Essential in Tracking Dynamic Optima.

Real-world optimisation problems are often dynamic. Previously good solutions must be updated or replaced due to changes in objectives and constraints. It is often claimed that evolutionary algorithms are particularly suitable for dynamic optimisation because a large population can contain different solutions that may be useful in the future. However, rigorous theoretical demonstrations for how populations in dynamic optimisation can be essential are sparse and restricted to special cases. This paper provides theoretical explanations of how populations can be essential in evolutionary dynamic optimisation in a general and natural setting. We describe a natural class of dynamic optimisation problems where a sufficiently large population is necessary to keep track of moving optima reliably. We establish a relationship between the population-size and the probability that the algorithm loses track of the optimum.

Introduction

In a classical optimisation setting, so-called static optimisation, the focus is usually directed to finding an optimal or a high quality solution as fast as possible. In real-world optimisation, problem specific data may change over time, thus previously good solutions can lose their quality and must be updated or replaced. Automatic optimal control is a typical illustration of these situations, e. g. parameters of a machine can be set optimally under ideal conditions of a factory but they need to be adapted to changes in the real environment upon deployment. Dynamic optimisation is an area of research that is concerned with such optimisation problems that change over time. A specific characteristic is that it does not only focus on locating an optimal solution but also on tracking a moving optimum (see [7] for a definition).

It is often suggested that Evolutionary Algorithms (EAs), especially the ones with populations, are suitable for dynamic optimisation because a large population can contain different solutions which could be useful in the future [22]. However, theoretical demonstrations for how populations in dynamic optimisation can be essential are sparse and restricted to special cases. The ability of a very simple EA without a population, the EA, to track a target bitstring in a OneMax-like function is analysed in [4, 20]. The analysis has recently been extended from bitstrings to larger alphabets [11]. The influence of magnitude and frequency of changes on the efficiency of the EA in optimising a specifically designed function was investigated in [19], showing that some dynamic optimisation problems become easier with higher frequency of change. The analysis of the EA that uses a larger offspring population but still not a real population on a simple lattice problem is presented in [9]. The efficiencies of specific diversity mechanisms when using an actual population were compared in [18]. This was done for a specific example function (introduced by [19]) and considering low frequency of changes. It was shown in [12] that a Min-Max Ant System (MMAS) can beat the EA in a deterministic dynamic environment. The comparison was later extended to general alphabets and to the EA that preserves genotype diversity [15]. With that particular setting of the EA, the size of the alphabets defines a threshold on the parent population size so that the algorithm is able to track and reach the optimal solution in polynomial time. The result was also extended to the single-destination Shortest Path Problem [14]. Comparisons were also made between EAs and Artificial Immune System (AIS) on a OneMax-like problem with the dynamic being periodic [10].

Considering the existing analyses we can in summary note two shortcomings that leave the impression that important fundamental questions about dynamic optimisation are still not answered satisfactorily. One shortcoming is the concentration on simple evolutionary algorithms and other search heuristics that do not make use of an actual population. Clearly, the advantages of a population-based approach cannot be explored and explained this way. The other is that many studies consider very complex dynamic environments that make it hard to see the principal and fundamental issues. Therefore, the fundamental question why even a simple population without complicated diversity mechanisms can be helpful in dynamic environments requires more attention.

Motivated by the above facts, we will use a simple argument considering a very general class of dynamic functions to show that a population is essential to keep track of the optimal region. We define our function class on the most often used search space, bit strings of a fixed length. However, it is not difficult to extend the function class to be defined for any finite search space and any unary mutation operator . The class is called -stable on with respect to , where n is the required number of bits to specify a search point of and c and are positive constants independent of n. The function class is only restricted by the probability of recovering the optimal region via the mutation operator (see Definition 4). The definition of the function class does not refer explicitly to other function characteristics, such as the topology or the fitness values of the set of optimal solutions, or the distribution of fitness values of the set of non-optimal solutions.

We will use the Moving Hamming Ball function from [2] as an illustrative example over the search space and with respect to the bitwise mutation operator. We also use this specific function to argue that an approach based on a single individual, such as the EA, is inefficient in tracking the optimal region in spite of being equipped with the same mutation operator. On the other hand, we show that a population-based algorithm with a sufficiently large population can efficiently track the moving optimal region of any dynamic function of the class defined for any given finite search space.

The remainder of the paper is organised as follows. The next section first gives a formal description of dynamic optimisation and efficient tracking, then the class of dynamic functions that we consider is described with an example function. Next, we consider the EA and RLS on the function class and provide an analysis to serve as an example how search heuristics based on single solutions are not able to track the optimal solutions over time. The efficiency of population-based algorithms is then explained by proving a positive result about their performance. Here, we use the setting of non-elitist populations and show that, with a sufficient selective pressure, the ability of the population to track the moving optimal region is overwhelmingly high with respect to the population size. On the other hand, as a consequence of a fair comparison to a single-individual approach, the population must not be too big in order to capture the frequency of changes. Finally, we summarise, conclude and point out directions for future research.

The paper uses the following notation and terminology. For any positive integer n, define . The natural logarithm is denoted by , and the logarithm to the base 2 is denoted by . The Hamming distance is denoted by and the Iverson bracket is denoted by . We use to denote the indicator function of a set A, i.e. if , and 0 otherwise. For a given bitstring , the Hamming ball around x with radius r is denoted by . The bitstring containing n one-bits and no zero-bits is denoted . An event is said to occur with overwhelmingly high probability (w. o. p.) with respect to a parameter n, if the probability of the event is bounded from below by .

A General Class of Dynamic Functions

Before defining the class of -stable dynamic functions which will be studied in this paper, we first formalise our notion of dynamic optimisation, and we define what we mean when saying that a dynamic search heuristic tracks a moving optimal region efficiently.

A Formal Description of Dynamic Optimisation

We focus on optimisation of pseudo-Boolean functions with discrete-time dynamics, as formalised below. Note that our formalisation can be generalised to any finite search space , e.g. replacing with , and our results for population-based algorithms also hold for this generalisation.

Definition 1

A dynamic function F is a random sequence of functions where for all . The optimal regions associated with F is the sequence , where .

The perhaps simplest, non-trivial example of a dynamic function is a periodic function which deterministically alternates between two functions, say and , such that and for all . We will consider more complex dynamic functions, where the sequence of functions is random and non-periodic. Although the sequence of functions in a dynamic function is random, each individual function is deterministic, i.e., we do not consider dynamic optimisation with noisy functions.

In this paper, we do not make any assumption about the changes of the function and the speed of the algorithm. It has been pointed out that it is important to consider the relationship between the speed of the execution platform where the algorithm runs and the speed of change of the function because this has significant influence on the performance [10]. Almost all studies assume that the function cannot change within one generation of the algorithm. The only exception we are aware of is a paper by Branke and Wang [1] who analyse a (1, 2) EA. We follow this idea but consider a much broader class of algorithms.

When applying a search heuristic to a dynamic function, we therefore have to consider two time lines: the first is associated with the evolution of the dynamic function, and the second is associated with the search points generated by the heuristic. Following the convention from black-box complexity [6], we assume that the function evaluations are the most expensive operations, for sake of the analysis becoming the basic time steps of an algorithm. The time consumed by all other operations, such as sampling an individual or applying a mutation operator, is assumed to be negligible. We connect the two time-lines by assuming that every time the heuristic queries a search point, the time-line of the dynamic function increases by one. We allow dynamic search heuristics some flexibility in that search points can be queried not only with respect to the most recent function , but also with respect to past functions. For example, the individuals in a population can be evaluated with respect to one particular time. We also assume that the decisions made by the search heuristic does not influence the dynamic of the function. The dynamic optimisation-scenario we have described is summarised in the following definition.

Definition 2

A dynamic search heuristic is an algorithm which given a search history of elements in , selects a search point and an evaluation time , and evaluates .

An element in a search history describes the search point queried by the algorithm in step t, the time point with which the search point is evaluated, and the corresponding function value . We can now formalise the notion of efficient tracking of optima.

Definition 3

A search heuristic is said to efficiently track the optima of a dynamic function F if there exist and constants such thatwhere is the sequence of search points queried by the heuristic, and is the sequence of optimal search points of function F.

Informally, Definition 3 means that the algorithm queries optimal search points frequently. More precisely, within every sub-interval of length within the exponentially long time interval from to , a constant fraction of the queried search points are optimal. Note that the optimality of a search point is defined with respect to the query time, and regardless of the function with which the algorithm evaluates the search point. The constraint on the length of sub-intervals guarantees that the time between generation of two optimal search points is bounded from above by a polynomial. It is clear from the definition that an algorithm is inefficient if with a sufficiently high probability, e.g. at least constant, it loses track of the optimal region and does not recover it within a polynomial number of steps.

A Class of Stable Dynamic Functions

The class of -stable dynamic functions with respect to a variation operator is defined as follows.

Definition 4

Let be any unary variation operator, and , . If there exist constants such that with probability at least , the optimal regions of a function F satisfy for all time points t and with , and for all search points ,then F is called -stable with respect to .

Definition 4 covers a large class of dynamic optimisation functions for any given pair of parameters . The optimal regions over time can take many shapes, including disconnected pieces over as long as the distances between them and the cardinality of the intersections allow the probabilistic condition to hold. Figure 1 illustrates the required condition.

Fig. 1

Illustration of a -stable dynamic function, in which any search point in the optimal region of time t can be mutated into the optimal region of time with probability at least

Given an operator , we focus on the class of -stable functions where c and are positive constants. We will show that a population-based algorithm with a sufficiently large population and a sufficiently strong selection pressure can track the optimal region of any function in the class efficiently. The next section gives an example function of the class for being bitwise mutation and explains how it fits within the framework of -stable function. We will then use the example function to argue that algorithms that base their search on a single individual, such as the EA, can be inefficient.

Example of a Stable Pseudo-Boolean Function

We consider the Moving Hamming Ball function as described in [2]. The static version of the function has the following form.

Definition 5

The Hamming Ball function around a target bitstring and a radius r is defined as,

It suffices to change in sequential steps to create a dynamic version from the static one. We use the following dynamic setting for the function: the points in time when the target is changed are determined by a sequence of random variables drawn from a Poisson distribution.

Definition 6

Let be a sequence of random variables independently sampled from a Poisson distribution with parameter , , be some integer in [n], and be a sequence of bitstrings generated byThe Moving Hamming Ball () function with parameters r, , and is defined as

The function fits within the stability framework of Definition 4 with respect to the bitwise mutation operator . This variation operator, which has a parameter , flips each position in the bitstring independently with probability . Hence, the probability of mutating a bitstring into a bitstring is

Lemma 7

For all positive constants d, and , the function is -stable with respect to the bitwise mutation operator with parameter .

Proof

For any given time t, let X be the random variable associated with the number of time steps in the future that the target bitstring will be changed. If we pick , then it is clear that within the next time steps, there will be more than one change of the target bitstring if and only if . It follows from Lemmas 18 and 19 thatIt suffices to pick the constant so that with a probability of at least , there is at most one change to the target function within the next time steps. Under that condition, it holds for all and for all , that for some .

In the case that , e.g. the target does not move or it moves closer to x, it suffices to not flip any of the bits. For any constant , it holds for all thatIn the case that , it suffices to recover the bits among the mismatched ones, soNote that , henceand is -stable with respect to .

It is not difficult to see that the stability condition of the function class still holds with the following relaxations:We will not consider these relaxations as they are not required to distinguish between the effectiveness of single-individual and population-based evolutionary algorithms.

The fitness of the solutions inside the Hamming ball changes when the target string moves,

The fitness of the solutions outside the current Hamming Ball can be distributed differently, as long as they are less than the current optimal fitness,

The moving step is relaxed to be sampled from any discrete distribution over in each change of the target bitstring.

Algorithms

We will compare the performance of population-based and single-individual based evolutionary algorithms. In this section we first define these classes of algorithms.

We are considering dynamic optimisation problems and, as mentioned in the introduction and discussed by Jansen and Zarges [10], it is important to clarify how the algorithms deal with change of the fitness functions, in particular if this happens during one generation. In this paper, we consider algorithms that make use of consistent comparisons when applying on a dynamic function: when an algorithm has to make fitness comparisons on a set of solutions, it will first make a static copy of the dynamic function and evaluate the solutions on this copy. This approach corresponds to an implementation where the necessary data to evaluate the optimisation function is collected before evaluating a set of solutions. Meanwhile the real function may have changed more or less depending on the number of solutions in the set.

We first consider the single-individual approach described in Algorithm 1. The algorithm keeps a current search point x. In each iteration, it produces a new candidate solution , and compares it with the current search point using the same function. Hence, static copies of the dynamic function are made in every two time steps. This corresponds to a frequent update of the dynamic function. We let be the bitwise mutation operator described in Sect. 2.3, and obtain the well-known EA [5]. However, the result can be easily generalised to other mutation operators, such as the single-bit flip operator used in the RLS algorithm.

We are mostly interested in the influence of the population size, designated by the parameter , on the ability of a population-based algorithm to track the moving optimal region. We focus on the non-elitist setting as described in Algorithm 2. The algorithm uses a unary variation operator denoted by , no crossover operator, and an unspecified selection mechanism . The selection mechanism is any random operator that given a population P and access to a fitness function returns one of the individuals in P. By specifying different and , Algorithm 2 can instantiate a large number of population-based search heuristics, such as the () EA. The number of search points produced in each round is the only parameter that appears in the description of Algorithm 2. The () EA fits within this framework by making sure that the selection in line 4 only takes into account the best of the search points created in the last round.

The algorithm maintains a population of individuals which during one generation (steps 2–6) is replaced by a newly created population of the same size. As for the EA, we assume that the initial population is contained in the first optimal region . Each individual in the next population is created by first making a copy x of one parent individual which is selected from the current population (step 4 selection), then modifying the copy using operator (step 5, mutation).

When selecting individuals, the algorithm must take into account that multiple changes to the fitness function can occur during one generation. Here, we assume that the algorithm makes a static copy of the fitness function at the beginning of each generation, i.e. at time . The selection mechanism compares all individuals in a generation using the static copy. Note that if the population size is too large with respect to the problem parameter (which controls the frequency of change of the dynamic function), then the optimal region may change several times between two consecutive generations. Hence, the population size should not be too large. However, we will show in the next section that a sufficiently large population is also essential to keep the population within the optimal region.

The result for populations will be first shown for any finite search space and any mutation operator because the class of dynamic function is defined with respect to the operator . Then we will use over as a specific example.

Although Algorithm 2 can use any selection mechanism , we are looking for choices of that allows the algorithm to track optima efficiently. Formally, applied on finite populations of size is represented by the transition matrix , where represents the probability of selecting individual P(i), i.e. the i-th individual, of P. We also write , e.g. in the algorithm, to express that x is sampled from the distribution over P given by . We use to denote the ith best individual of P, or the so-called -ranked individual. Similar to [3, 13], we characterise by the cumulative selection probability.

Definition 8

([13]). Given a fitness function , the cumulative selection probability associated with selection mechanism is defined on f for all and a by

Informally, is the probability of selecting an individual with fitness at least as good as that of the -ranked individual, assuming that P is sorted according to fitness values. We are interested in a lower bound function of . Most often this lower bound is independent of P, in which case we simply write it as a function of only, i. e. as .

Performance Analysis

In this section, we first show that single-individual approaches are inefficient in tracking moving optima on at least one example function of the class, precisely on . Then we prove a general result that an appropriately parameterised population-based algorithms can efficiently track the moving optima of any function in the class.

Inefficiency of a Single Individual

In this section, we will show that the EA spends an exponential fraction of its time outside the optimal region of a function, for a sufficiently small constant b, any constant and any . That is to say the algorithm is inefficient even in tracking such a stable function.

To prove such a result, we have to analyse the behaviour of the algorithm both inside and outside the moving Hamming ball: We assume that the algorithm starts at the center of the first optimal region and show that after some initial time, whenever the center of the ball moves, there is a constant probability that the EA will memorise a search point outside of the new ball; Whenever the algorithm is outside of the optimal region there is also a constant probability that the memorised search point will drift away from the optimal region (eventually get lost), before an optimal solution inside the new ball is discovered. Since the changes to the function happens within an expected polynomial number steps, we can conclude that with a high probability, the EA only spends a polynomial number of time steps inside the moving Hamming ball.

We start with the first argument, the behaviour of the algorithm inside the Hamming ball. We notice that the changes induced by the dynamics of the fitness function strongly drag the target away from the current memorised search point, however this does not happen in every iteration. In every iteration, the changes by mutation drive the memorised solution away from the center of the current Hamming ball, but the elitist selection also keeps the memorised solution from falling outside. We have the following analysis of the drift.

We consider the process , where is the Hamming distance to the border of the optimal region of at time t, i.e. . The process starts with , e.g. exactly at the center of the Hamming ball. Given , define , then is the drift towards the border at time t and where .

First of all, the dynamic now only kicks in every cn time steps in expectation. Also, the contributing drift is positive. For example, if the dynamic kicks in, let Z be the number of bits being corrected by the dynamic, then we have , and the contributing drift is for any .

We now compute the drift by mutation at time t and where . Let X and Y be the number of bits being corrected and being messed up respectively by the mutation, so , and the two variables are independent. Note that for all integers , and , it holdsThus for , stochastically dominates and we also havefor any constant . Therefore, with any constant , we have that . Hence, for we have at least a constant drift away from the centerThe only position where we have a drift toward the center is the one at the border, e.g. . However, this is not a strong drift. When the target bitstring does not move in the next iteration, we have , then the negative drift is no more thanIn summary, we get the drift by mutation: It is then suggested that the equilibrium state of the memorised search point is around the border. Furthermore, we can quantify the expected fraction of time that the search point is found at the border, using the following tool.

Lemma 9

Given a stochastic process over a state space , and two constants such thatthen for all

, and

,

Define . For all , it holdsIt follows thatFinally, since

The following lemma considers non-negative, integer-valued stochastic processes with positive drift at most in state 0, and negative drift at least elsewhere. It provides a lower bound on the probability of such a process being in state 0 after some time.

Lemma 10

Let be any stochastic process with support in for some fixed , which satisfies the properties of Lemma 9 for some . Then for any random variable which is independent of , it holds

Choose , and define where and , i.e., we consider the process from a random starting point . Due to independence between , and , we haveLemma 9 applied to now implies

We now show that once the EA has lost track of the optimal region it will take a long time to recover. We assume that the objective function is with radius , i.e. for some constant (note that b can depend on n). The first step in this proof is to show that with not too small probability the EA ends up far away (more specifically, in a linear distance) from the Hamming ball before recovering it.

Lemma 11

Given , let be a sequence of random bit strings such that and for some , i.e. for some . For any , define . It holds that where for a not too large constant .

We begin with considering another random sequence where for each the point is created by flipping one randomly selected bit in . Let be defined as but with respect to instead of .

Let , i. e., the probability to enter the Hamming ball before reaching distance d given the process is started with Hamming distance x. Note that, for symmetry reasons, is well defined, i. e., the probability does only depend on the Hamming distance x and not the specific choice of .

By definition of we have for and for . For all other values of x, i. e., for we haveby definition of the sequence because with probability the Hamming distance to the centre of the Hamming ball o is increased by 1 and with the remaining probability x / n it is decreased by 1. If we pessimistically assume that the probability to move towards the Hamming ball is always equal to we obtain an upper bound on and are in the situation of the gambler’s ruin problem with initial funds and , , and and the probability to be ruinedgives an upper bound on the probability to enter the Hamming ball before reaching distance d when starting with Hamming distance x to the centre of the Hamming ball. We consider the probability q(r, d, x) for different values of r, d and x. We are interested in the results for and consider for this where we chose the constant such that for some positive constant . It is clear that due to the upper bound on r such a constant exists.

It is not difficult to see that . For this is and the best bound we can obtain. For we need to be more precise.

We begin with the case and . For this setting we consider and know that holds. Now we consider and see that holds.

Finally, for the case , we also consider and know that holds. Now we consider and see that holds.

Together, we have for all values of r and . Since the sequence corresponds to ‘local mutations’ this proves the claim for random local search. We generalise the statement to the EA in the following way. We can express a standard bit mutation as a process where first a random number is chosen and then k bits are selected uniformly at random to be flipped. The case does not flip any bit and can be ignored. The case is covered by the analysis for RLS. For larger we observe that such a step is very similar to a sequence of k steps where exactly 1 bit is flipped. The difference does not change the limits we considered above. Since in one standard bit mutation k bits flip with probability we can ignore steps where bits flips since they contribute too little to change the asymptotic result.

Theorem 12

On with any constants , and the EA with mutation rate 1 / n will spend only an exponentially small proportion of its time in optimal regions.

We assume the process starts inside the first Hamming ball, and consider the process as as described before in Lemma 9. For the standard mutation, we have the drifts according to Eqs. 1 and 2 are and . Applying Lemma 10 gives that after time steps, whenever the center of the Hamming ball is moved, it holds that . Conditioned on this event, the probability that the dynamic move so that is where , e.g. the number of bit positions being corrected by the dynamic (see the definition of Z before Lemma 9). Therefore, and by Markov’s inequality and it holds thatWhen the is outside of the current Hamming ball, it follows from Lemma 11 that there is a probability of at least that the EA reaches linear Hamming distance to the Hamming ball before finding its way back to it. Application of the negative drift theorem [17] yields that the probability to find the way back into the optimal region within steps is for a sufficiently small constant .

We have just show that in every time steps, whenever a change occurs to the target bitstring there is a probability of at least that the EA will lose track of the optimal region where and are constants. Applying this argument n times (a change occurs approximately every cn time steps), we conclude that with an overwhelmingly high probability, the EA will spend no more than time steps within the optimal region of the function.

Efficiency of Non-elitist, Population-Based Algorithms

Theorem 13, which is the main result in this section, gives conditions under which the non-elitist, population-based Algorithm 2 tracks the optimal regions of dynamic functions efficiently. We show that these conditions can be satisfied for the moving Hamming-balls function for any constant .

Theorem 13

If there are constants and such thatthen Algorithm 2 initialised with tracks the optima of F efficiently.

F is a -stable dynamic function wrt. with , and

satisfies for all ,

Condition 1 of the theorem requires that the optimal region of the function does not move too much relatively to the variation operator during one generation. The population size is a parameter of the algorithm which can be chosen freely. So if the function is -stable, then the first condition can be satisfied by setting population size . Condition 2 requires that the selection mechanism induces a sufficiently high selective pressure. Note that increasingly high selective pressure is required for decreasing values of , where is the probability of recovering the optimal search region via mutation (see Definition 4).

The central argument in the analysis is illustrated in Fig. 2. It follows from the stability-assumption that any search point in can be mutated into for any with probability at least . Hence, if the algorithm selects a search point in with probability , then the offspring belongs to with probability at least . This argument is invoked in both of the two steps of the analysis.

Fig. 2

Illustration of Lemmas 14 and 15

Lemma 14

Assume that conditions 1 and 2 of Theorem 13 hold. Then for any , , if , then any offspring in generation belongs to with probability at least .

The algorithm produces an individual in if the algorithm selects an individual in and mutates this individual into . The probability of this event is .

Lemma 15, which is the first step of the analysis, implies that in every generation , a large fraction of the population belongs to . This can be shown inductively by arguing using Lemma 14 that if many individuals in belong to , then whp. many individuals in belong to . Knowing that many individuals in belong to for every generation gives us some control on the dynamics of the population. However it does not imply that the dynamic performance measure in Definition 3 is satisfied because the individuals in may not necessarily have been optimal when they were generated. A second step in the analysis is therefore required, showing that if sufficiently many individuals in population belong to , then many offspring in generation were optimal at the time they were generated. This second step is contained in the proof of Theorem 13.

Lemma 15

Assume that conditions 1 and 2 of Theorem 13 hold. Then for any generation , if , then

By Lemma 14, any offspring in generation belongs to independently with probability . Hence, by a Chernoff bound, the probability that less than offspring belongs to is .

We are now in position to prove the main result of this section.

Proof of Theorem 13

We say that generation fails if and . By Lemma 15 and a union bound, the probability that any of the first generations fails is , assuming that is a sufficiently small constant. By Lemma 2 and assuming no failure, any individual with belongs to the optimal region with probability at least . By the definition of the algorithm, individuals within the same generation are produced independently. During any time interval where , at least individuals are produced in the same generation, and hence independently. It therefore holds by a Chernoff bound that for any time interval with ,The theorem now follows by taking into account the failure probability with a union bound, and choosing the parameters and in Definition 3.

Theorem 13 implies that with a sufficiently slow dynamic, e.g. for any constant , the population-based algorithm can efficiently track the moving optima of the function, given that induces a sufficiently strong selective pressure. We now show that given any constant , it is possible to parameterise many selection mechanisms so that they satisfy this requirement on . The selection mechanisms are:The following theorem shows how these selection mechanisms can be parameterised to satisfy the second requirement of Theorem 13, and hence ensure that Algorithm 2 tracks the moving optima of any -stable function with respect to the mutation operator .

In k-tournament selection, k individuals are sampled uniformly at random with replacement from the population, and the fittest of these individuals is returned.

In -selection, parents are sampled uniformly at random among the fittest individuals in the population.

A function is a ranking function [8] if for all , and . In ranking selection with ranking function , the probability of selecting individuals ranked or better is . Linear ranking selection uses for some . Exponential ranking selection uses for some .

Theorem 16

For any constant , let F be any -stable function wrt. for . If there is a constant such that Algorithm 2 initialised with , and selection mechanism eitherthen the algorithm tracks the optima of F efficiently.

k-tournament selection with ,

-selection with ,

linear ranking selection with , or

exponential ranking selection with ,

The result follows from Theorem 13 if we can show that there exist constants and such that for all . The results for k-tournament, -selection and linear ranking follow from Lemmas 5, 6 and 7 from [13] with in place of . For exponential ranking, we notice thatthe result then follows similarly to k-tournament as in the proof of Lemma 5 in [13] with in place of k (Equations (3) and (4) in [13] show that , then the constants and are shown to exist given the condition on k).

Finally, we apply Theorem 16 to show that population-based EAs can track the optima of the example Moving Hamming Ball function efficiently. Note that the parameter in linear ranking selection can only take values in the interval (1, 2]. The conditions of Theorem 16 can therefore only be satisfied if , i.e., the optimal regions can only change slightly. For the last part of the paper, we therefore exclude linear ranking selection.

Corollary 17

For any constants , , and , Algorithm 2 with the bitwise mutation operator for , with population size , and selection mechanism eithercan efficiently track the moving optima of .

k-tournament selection with ,

-selection with ,

exponential ranking selection with .

It follows from Lemma 7 that for any constant , is -stable with respect to the mutation operator . Since for a sufficiently small , the function is also -stable. The result then follows by applying Theorem 16.

Conclusion

This paper has considered the frequently stated intuition that evolutionary algorithms maintaining a population of diverse solutions can be more resilient to dynamic changes in the objective function than algorithms maintaining single solutions. We have described a general class of fitness functions where population-based evolutionary algorithms outperform single-individual evolutionary algorithms. We have proved that for this function class, single-individual approaches, such as the EA and RLS, have a constant risk of losing the optimal solution region at any given time. Moreover, these single-individual algorithms not only lose the optimal region with constant probability, but are also likely to drift away from the optimal region subsequently.

On the other hand, assuming a not too high frequency of change, we describe sufficient conditions such that a non-elitist population-based evolutionary algorithm will remain within the optimal region with overwhelmingly high probability. Our analysis covers a range of the most commonly used selection mechanisms, and we provide appropriate parameter settings for each of them. Furthermore, the success of the population-based evolutionary algorithm does not rely on an explicit diversity mechanism. Our analysis gives further explanations of how and why populations can be essential and widely used in dynamic optimisation.

As future work, we would like to investigate further the influence of population settings within this class of dynamic functions, such as elitist populations, the necessary condition for the population size with respect to the frequency and magnitude of changes and how a population could rebuild itself after losing a few optimal solutions.