27.06.2013 Views

Evolution and Optimum Seeking

Evolution and Optimum Seeking

Evolution and Optimum Seeking

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

240 Summary <strong>and</strong> Outlook<br />

this, which in certain cases yields a considerable improvement in the rate of progress. It<br />

can be achieved either by separate variation of the st<strong>and</strong>ard deviations i for i = 1(1)n,<br />

by recombination alone, or, even better, by both measures together. Whereas in the two<br />

membered scheme, in which (unless the (0)<br />

i are initially given di erentvalues) the contour<br />

lines of equiprobable steps are circles, or hyperspherical surfaces, they are now ellipses or<br />

hyperellipsoids that can extend or contract along the coordinate directions following the<br />

n-dimensional normal distribution of the set of n r<strong>and</strong>om components zi for i = 1(1)n :<br />

1<br />

w(z) =<br />

(2 ) n 2<br />

nQ<br />

i=1<br />

i<br />

This is not yet, however, the most general form of a normal distribution, which is rather:<br />

p<br />

Det A<br />

w(z) =<br />

(2 ) n exp<br />

2<br />

; 1<br />

2 (z ; )T A (z ; )<br />

The expectation value vector can be regarded as a deterministic part of the r<strong>and</strong>om step<br />

z. However, the comparison made by Schrack <strong>and</strong>Borowski (1972) between the r<strong>and</strong>om<br />

strategies of Schumer-Steiglitz <strong>and</strong> Matyas shows that even an ingenious learning scheme<br />

for adapting to the local conditions only improves the convergence in special cases. A<br />

much more important feature seems to be the step length adaptation. It is now possible<br />

for the elements of the matrix A to be chosen so as to give the ellipsoid of variation any<br />

desired orientation in the space. Its axes, the regression directions of the r<strong>and</strong>om vector,<br />

only coincide with the coordinate axes if A is a diagonal matrix. In that case the old<br />

scheme is recovered whereby thevariances ii or the 2<br />

i reappear as diagonal elements of<br />

the inverse matrix A ;1 . If, however, the other elements, the covariances ij = ji are nonzero,<br />

the ellipsoids are rotated in the space. The r<strong>and</strong>om components zi become mutually<br />

dependent, or correlated. The simplest kind of correlation is linear, which is the only case<br />

to yield hyperellipsoids as surfaces of constant step probability. Instead of just n strategy<br />

parameters i one would now have tovary n<br />

2 (n + 1) di erent quantities ij. Although in<br />

principle the multimembered evolution strategy allows an arbitrary number of strategy<br />

variables to be included in the mutation-selection process, in practice the adaptation of<br />

so many parameters could take too long <strong>and</strong> cancel out the advantage of more degrees of<br />

freedom. Furthermore, the ij must satisfy certain compatibility conditions (Sylvester's<br />

criteria, see Faddejew <strong>and</strong> Faddejewa, 1973) to ensure an orthogonal coordinate system<br />

exp<br />

; 1<br />

2<br />

or a positive de nite matrix A. In the simplest case, n =2,with<br />

there is only one condition:<br />

<strong>and</strong> the quantity de ned by<br />

12 =<br />

A ;1 =<br />

"<br />

11 12<br />

21 22<br />

2<br />

12 = 2<br />

21 < 11 22 = 2<br />

1<br />

#<br />

nX<br />

i=1<br />

2<br />

2<br />

zi<br />

12 q<br />

( 1 2) ;1 < 12 < 1<br />

i<br />

2 !

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!