Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat procedury krzyżowania jednopunktowego

Innym bardzo podobnym operatorem jest krzyżowanie dwupunktowe, w którym to chromosomy dzielimy na trzy części a nie na dwie części jak w przypadku krzyżowania jednopunktowego. Podobnie jak w krzyżowaniu jednopunktowym losujemy dwóch osobników X1 i X2z populacji do krzyżowania z prawdopodobieństwem pk i losujemy im punkty przecięcia g1i g2 chromosomu - wektora binarnego. Podzielone na trzy części chromosomy wymieniają się tylko środkowym fragmentem kodu genetycznego.

Z=X₁¹, X_g¹₁, X_g²₁₊₁, ..., X_g²₂, X_g¹₂₊₁, ..., X_n¹

Schemat procedury krzyżowania dwupunktowego

Rys. 3. Krzyżowanie dwupunktowe

Legenda:

X₁, X₂-osobniki rodzice X⁰₁, X₂⁰ -osobniki dzieci X_a, X_b -kod który ulega krzyżowaniu

g i g’-losowo wybrane punkty podziału

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat procedury krzyżowania dwupunktowego

Istnieje jeszcze kilka różnych metod krzyżowania wielopunktowego, w którym chromosom dzielimy na 3, 4 lub inne części, a wymiana

fragmentu genu jest dowolna w zależności od naszych upodobań. Oprócz tego jest jeszcze szereg innych krzyżowań i wciąż powstają nowe pomysły na zoptymalizowanie operatora krzyżowania.

W przypadku reprezentacji zmiennopozycyjnej w krzyżowaniu

arytmetycznym zamiast wymieniać geny, jak w przypadku krzyżowania binarnego uśrednia się wartości z dwóch elementów w wybranych punktach. Operator zmiennopozycyjny krzyżowania arytmetycznego ma taką własność, że każdy element nowego chromosomu pozostaje nadal w dziedzinie. Krzyżowanie arytmetyczne można stosować albo do wybranych elementów dwóch wektorów, albo do całych wektorów. Przy zbiorach nie wypukłych potomek powstaje na ”krawędzi” zbioru albo należy szukać gdzieś indziej pary do skrzyżowania lub zastosować algorytm naprawy.

Schemat procedury krzyżowania dwupunktowego

Krzyżowanie arytmetyczne działające na dwóch argumentach definiuje się jako liniową kombinację dwóch wektorów. Jeżeli w takim krzyżowaniu mają być brane pod uwagę para rodziców x₁ix₂to potomkowie ich mają postać:

x₁⁰ = ax₁+ (1 − a)x₂ , x₂⁰ = ax₂+ (1 − a)x₂ W operatorze tym używa się losowej wartości a∈[0,1], co zawsze gwarantuje, że są to punkty dopuszczalne, czyli należące do dziedziny przeszukiwanej przestrzeni zadania. Jeśli a=¹₂, to takie krzyżowanie nosi nazwę gwarantowanego średniego krzyżowania, o którym szerzej w pracy Davisa L. ”Adapting Operator Probabilities in Genetic Algoritms”.

Trzecią metodą krzyżowania dla liczb zmiennopozycyjnych, którą warto przedstawić jest krzyżowanie heurystyczne, o której szerzej z kolei w pracy Wrighta A. H. ” Genetic Algorithms for Real Parameter

Optimization”.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat procedury krzyżowania dwupunktowego

Do określenia kierunku poszukiwań operator ten używa wartości funkcji celu, tworzy tylko jednego potomka i może on w ogóle nie utworzyć potomka. Operator ten tworzy jednego potomka z3 z dwojga rodziców z1

i z2według wzoru:

z3= r (z2 - z1) + z2

gdzie r jest liczbą losową leżącą między 0 a 1, zaś rodzic z₂nie jest gorszy od z₁, to znaczy f (z₂) f (z₁) dla maksymalizacji i f (z₂) ¬ f (z1) dla zadania minimalizacji. Jeśli operator utworzy potomka, który jest niedopuszczalny, to w takim przypadku generuje inną wartość losową r i tworzy innego potomka. Jednak, gdy po q próbach operator nie jest w stanie utworzyć osobnika dopuszczalnego, to nie dochodzi do utworzenia nowego potomka. Operator ten jest bardzo silnym narzędziem, gdyż znakomicie obsadza osobników w optimum lokalnym i przeszukuje przestrzeń rozwiązań w obiecujących kierunkach. Oczywiście istnieje jeszcze szereg innych rodzajów i podejść do operatorów krzyżowania dla kodowania rzeczywistego - zmiennopozycyjnego.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Operatory mutacji

W przypadku mutacji w reprezentacji binarnej sprawa jest bardzo prosta wystarczy np. zanegować jeden bit w rozwiązaniu, aby otrzymać zupełnie nowego osobnika. Oczywiście mutacja może być bardziej urozmaicona (negacja losowej liczby bitów, odwracanie kolejności losowej liczby bitów, przesunięcie losowej liczby bitów i inne.). Należy jednak pamiętać, że operować ona może tylko na jednym rozwiązaniu. Mutowanie polega na negacji jednego z genów danego osobnika. Proces mutowania bieżącej populacji polega na:

mutowania.

Wybór osobnika odbywa się w następujący sposób: Pobieramy 1 osobnika bieżącej populacji, po czym losujemy progowo bez zwracania sprawdzając, czy osobnik może być zmutowany. Jeżeli tak, to losujemy z rozkładem

równomiernym q, wartość losową z przedziału: < left(gen), right(gen) >; -długość genu, określający gen poddania negacji. Po zmutowaniu pobieramy następnego osobnika i sprawdzamy, czy może on być zmutowany i powtarzamy ten cykl aż do wyczerpania osobników do mutowania.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Operatory mutacji

wyborze osobników, których należy zmutować,

przeprowadzaniu po kolei na każdym wybranym osobniku operacji mutowania.

Operatory mutacji

wyborze osobników, których należy zmutować,

przeprowadzaniu po kolei na każdym wybranym osobniku operacji mutowania.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat procedury krzyżowania dwupunktowego

Rys. 4. Mutacja jednopunktowa

W mutacji wielopunktowej stworzenie nowego chromosomu odbywa się za pomocą wygenerowania maski, która określa, które pozycje bitu zamienić. Działanie maski mutacji powoduje inwersję bitów kodu x na pozycjach, na których występują jedynki w kodzie maski i. Maska mutacji i = (a₀, a₁, . . . , a_{l −1}) losowana jest niezależnie dla każdej operacji, bez względu na wartość x, przy czym prawdopodobieństwo wylosowania jedynki jest stałe, równe prawdopodobieństwu mutacji p_m.

Schemat procedury krzyżowania dwupunktowego

Rys. 4. Mutacja jednopunktowa

W mutacji wielopunktowej stworzenie nowego chromosomu odbywa się za pomocą wygenerowania maski, która określa, które pozycje bitu zamienić.

Działanie maski mutacji powoduje inwersję bitów kodu x na pozycjach, na których występują jedynki w kodzie maski i. Maska mutacji i = (a0, a₁, . . . , a_{l −1}) losowana jest niezależnie dla każdej operacji, bez względu na wartość x, przy czym prawdopodobieństwo wylosowania jedynki jest stałe, równe prawdopodobieństwu mutacji p_m.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat mutacji wielopunktowej

Rys. 5. Mutacja wielopunktowa

Istnieje jeszcze wiele innych metod mutacji binarnych takich jak np.: mutacja brzegowa podobna do równomiernej lub mutacja

nierównomierna. W przypadku zmiennopozycyjnej (rzeczywistej) reprezentacji mutacja może polegać np. na wpisaniu do losowego pola struktury losowej wartości, która jest zgodna z cechami osobnika.

Schemat mutacji wielopunktowej

Rys. 5. Mutacja wielopunktowa

Istnieje jeszcze wiele innych metod mutacji binarnych takich jak np.:

mutacja brzegowa podobna do równomiernej lub mutacja nierównomierna. W przypadku zmiennopozycyjnej (rzeczywistej) reprezentacji mutacja może polegać np. na wpisaniu do losowego pola struktury losowej wartości, która jest zgodna z cechami osobnika.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat mutacji wielopunktowej

Dla każdego chromosomu generujemy liczbę losową q z zakresu [0,1].

Jeżeli q < prawdopodobieństwo mutacji to wybieramy chromosom do mutowania. Mutację dokonuje się poprzez dodanie wektora x’= x + m, gdzie:

m ∈ R_n m = (m₁,...,m_n).

Liczba mi jest losowana niezależnie dla i=1,...,n z rozkładem N(0,?), o N-wymiarowym rozkładzie normalnym i wartości oczekiwanej 0.

Rys. Wektor mutacji

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat mutacji wielopunktowej

Musimy pamiętać o tym, aby nie wypaść poza zbiór dopuszczalny.

Jeśli taka sytuacja się zdarzy to możemy zastosować mutację brzegową tzn. przyjąć punkt przecięcia wektora z końcowym obszarem zbioru dopuszczalnego jako rozwiązanie. Mutacja brzegowa prowadzi do zbieżności populacji, co w niektórych rozpatrywanych problemach jest korzystne.

Rolą mutacji jest nie tyle zmiana jakościowa chromosomów, co wprowadzenie nowego materiału genetycznego do danej populacji.

zależności od rodzaju problemu różne są funkcje sprawdzające przystosowanie osobników. W przypadku problemu

f(xi)=eval(code(xi))

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Schemat mutacji wielopunktowej

Musimy pamiętać o tym, aby nie wypaść poza zbiór dopuszczalny.

Rolą mutacji jest nie tyle zmiana jakościowa chromosomów, co wprowadzenie nowego materiału genetycznego do danej populacji.

Ocena osobnika i funkcja przystosowania

Ocena osobnika to badanie dobroci poszczególnych osobników. W zależności od rodzaju problemu różne są funkcje sprawdzające przystosowanie osobników. W przypadku problemu

f(xi)=eval(code(xi))

Ocena osobnika i funkcja przystosowania

znajdowania ekstremum funkcji dobrocią jest po prostu wartość funkcji w zadanym x. W przypadku rozwiązywania problemu komiwojażera jest to długość trasy, jaka jest reprezentowana przez danego osobnika. W czasie obliczania funkcji przystosowania ze zbioru osobników wydzielane są genotypy, które później są dekodowane a wartości, czyli fenotypy są przypisywane dla każdego osobnika. Dla każdego fenotypu wyliczana jest wartość funkcji przystosowania (ang. fitness) i przypisywana dla każdego osobnika. Otrzymany wynik powinien być proporcjonalny do użyteczności lub zdolności osobnika reprezentowanego przez ten chromosom. W metodach obliczeń ewolucyjnych funkcja oceny jest jedynym połączeniem zadania z algorytmem. Funkcja oceny szereguje osobników w populacji.

Lepsze osobniki mają większe szanse przeżycia i reprodukcji. Dlatego podstawową rzeczą jest określenie funkcji oceny, która charakteryzuje zadanie w doskonały sposób.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Ocena osobnika i funkcja przystosowania

Wraz ze wzrostem wymiarowości zadania zwiększa się stopień komplikacji rozwiązywanego problemu. Podwajają się automatycznie liczby różnych wartości w algorytmie genetycznym a nawet w niektórych przypadkach dochodzi do zaniku minimów lokalnych wraz z rosnącą liczbą wymiarów.

Intuicyjnie funkcja przystosowania, nazywana funkcją celu lub przez biologów - przystosowania do środowiska, stanowi pewien miernik zysku użyteczności lub innej wielkości, którą chcielibyśmy z maksymalizować. W przyrodzie przystosowanie, czyli potomstwo, które dorasta do poziomu reprodukcyjnego, jest ostatecznym i jedynym kryterium życia.

Metody selekcji

Krok ten jest esencją całej genetyki. W tym miejscu tworzona jest nowa populacja na podstawie już istniejącej. W zależności od wartości funkcji oceny (obliczanej w poprzednim kroku) dany osobnik ma większe, (gdy jest ’dobry’) lub mniejsze, (gdy jest ’słaby’) szanse na znalezienie się w kolejnym pokoleniu. Zadaniem reprodukcji jest wybranie rodziców dla operacji genetycznych prowadzących do utworzenia populacji potomnej O^t.

Operację selekcji można zrealizować na wiele sposobów. Jednym z najprostszych jest symulacja odpowiednio wykalibrowanej tarczy obrotowej rys. 7, każdemu ciągowi kodowemu (chromosomowi)

odpowiada sektor o rozmiarze proporcjonalnym do przystosowania. tzw.

ruletki.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji

Koło ruletki polega na n krotnym losowaniu (n - liczba osobników w populacji) ze starej populacji osobników, które zostaną przepisane do nowej populacji. Oczywiście wszystkie osobniki mają różne

prawdopodobieństwa wylosowania. Prawdopodobieństwo to jest liczone z następującego wzoru:

Wartosc przystosowania danego osobnika suma wartosci przystosowania wszystkich osobnikow

Powyższy wzór jest poprawny tylko wtedy, gdy maksymalizujemy funkcję oceny. Gdybyśmy ją minimalizowali można zastosować następujący wzór:

Wartosc najgorszego osobnika−Wartosc danego osobnika+1 suma wartosci wszystkich osobnikow +1

Wzór ten odwraca minimalizację na maksymalizację. Mamy trzy osobniki o następujących wartościach przystosowania dla maksymalizacji: 5 , 1 , 2

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Odpowiadające im wartości prawdopodobieństwa są zatem równe:

Pierwszy

1.Dla każdego ciągu kodowego obliczyć wartość przystosowania eval(vi) , vi (i = 1, ..., rozmiar populacji),

2. Obliczyć całkowite przystosowanie populacji F (suma po wszystkich wartościach funkcji przystosowania) ,

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Odpowiadające im wartości prawdopodobieństwa są zatem równe:

Pierwszy

1.Dla każdego ciągu kodowego obliczyć wartość przystosowania eval(vi) , vi (i = 1, ..., rozmiar populacji),

2. Obliczyć całkowite przystosowanie populacji F (suma po wszystkich wartościach funkcji przystosowania) ,

Metody selekcji

3. Obliczyć prawdopodobieństwo selekcji pi każdego chromosomu vi

4. Obliczyć dystrybuantę qi dla każdego chromosomu vi . qi =

j =1

Ciągi kodowe losujemy w następujący sposób (Obrót ruletką):

1. Wygeneruj r z przedziału [0..1]

2. Jeśli r < q₁, to wybieramy chromosom v₁, gdy q_{i −1}< r ¬ q_i wybieramy v₁ (2 ¬ i ¬ rozmiar populacji).

Weźmy dowolne 4 ciągi 5-bitowe (rys. poniżej). Niech każdy z nich reprezentuje liczbę całkowitą zapisaną w systemie dwójkowym, natomiast funkcja celu niech będzie kwadratem liczby reprezentowanej przez chromosom.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji

Rys. 6. Zasada ruletki

Rys. 7. Podział ze względu na przystosowanie osobnika

Metody selekcji

Reprodukcji dokonujemy czterokrotnie uruchamiając ruletkę. W naszym przykładzie chromosom nr 1 ma wskaźnik przystosowania 169, co odpowiada 14,4% sumarycznej wartości. Na ciąg ten przypada zatem 14,4% obwodu koła, tak więc w każdej próbie prawdopodobieństwo otrzymania ciągu nr 1 wynosi 0,144. Za każdym razem, gdy jest potrzebny nowy potomek, uruchamiamy ruletkę wybierając jednego z kandydatów do reprodukcji. Dokonujemy tego w następujący sposób:

losujemy liczbę z przedziału [0 ; 1], zaś chromosomom przypisujemy określone rozłączne podprzedziały tego przedziału (np. ciąg nr 1 - (0 ; 0,144), ciąg nr 2 - (0,144 ; 0,636), gdzie 0,636 = 0,144 + 49,2 itd.).

Dzięki temu lepiej przystosowane ciągi kodowe wprowadzają większą liczbę potomków do następnego pokolenia. Proces ten jest kontynuowany dopóki, dopóty wszystkie miejsca w nowej populacji nie zostaną

zapełnione. Oczywiście pewne chromosomy będą wybrane więcej niż raz.

Jest to zgodne z twierdzeniem o schematach, które mówi, że z

najlepszych chromosomów powstaje więcej kopii, ze średnich tyle samo , a najgorsze wymierają.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji

Istnieje kilka sposobów obliczania ’szansy’ poszczególnych osobników.

Selekcja metodą rankingu linowego jest bardzo podobna do selekcji metodą koła ruletki.

Modyfikacja polega jedynie na zmianie funkcji określającej

prawdopodobieństwo wyboru danego osobnika. Przed przystąpieniem do tej selekcji należy nadać każdemu z osobników pewną wartość

(przystosowanie) zależną od jego położenia na liście posortowanej względem wartości funkcji oceny. Jeśli chcemy maksymalizować to wartości powinny być posortowane rosnąco, w przypadku minimalizacji wartości powinny być posortowane malejąco.

Metody selekcji

Istnieje kilka sposobów obliczania ’szansy’ poszczególnych osobników.

Selekcja metodą rankingu linowego jest bardzo podobna do selekcji metodą koła ruletki.

rozm pop

i =1

eval (vi)

Modyfikacja polega jedynie na zmianie funkcji określającej

prawdopodobieństwo wyboru danego osobnika. Przed przystąpieniem do tej selekcji należy nadać każdemu z osobników pewną wartość

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji

Aby obliczyć prawdopodobieństwo wybrania każdego osobnika można skorzystać ze wzoru:

Prawdopodobienstwo = przystosowanie

suma przystosowania wszystkich osobnikow

pi = ^{eval (v}_F ⁱ⁾

Mając takie osobniki: 5, 1, 2 wartości prawdopodobieństw wyglądałyby następująco:

Metody selekcji

Aby obliczyć prawdopodobieństwo wybrania każdego osobnika można skorzystać ze wzoru:

Prawdopodobienstwo = przystosowanie

suma przystosowania wszystkich osobnikow

pi = ^{eval (v}_F ⁱ⁾

Mając takie osobniki: 5, 1, 2 wartości prawdopodobieństw wyglądałyby następująco:

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji

Ten sposób wyliczania prawdopodobieństw zmniejsza przewagę, jaką mają najlepsze rozwiązania, gdy ich przewaga jest bardzo duża i zwiększa przewagę, gdy jest ona bardzo mała.

Metoda turniejowa jest zupełnie różna od powyższych i polega na losowym wyborze z całej populacji kilku osobników (jest to tzw. grupa turniejowa), a później z tej grupy wybierany jest osobnik najlepiej przystosowany i on przepisywany jest do nowo tworzonej populacji.

Losowanie grup turniejowych oraz wybieranie z nich najlepszego osobnika należy powtórzyć aż do utworzenia całej nowej populacji. Selekcja metodą turniejową jest pozbawiona niedogodności metody koła ruletki, gdzie wymagana jest maksymalizacja funkcji oceny, w turnieju ważna jest jedynie informacja o łepszym”jednego rozwiązania nad innym.

Φ(X^q) Φ(Y ), dla każdego Y ∈ Q

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji gdzie:

Q - populacja wybrana z populacji potomnej do przeprowadzenia turnieju.

X^q - zwycięzca turnieju.

Proces wyboru osobników do populacji Q może być realizowany dwoma sposobami:

Losowanie osobnika do turnieju ze zwracaniem jest mniej restrykcyjne, ponieważ pozwala na reprodukcję dowolnego osobnika z niezerowym prawdopodobieństwem. Natomiast losowanie bez zwracania najgorszych osobników nie będzie w ogóle reprodukować.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji gdzie:

Q - populacja wybrana z populacji potomnej do przeprowadzenia turnieju.

X^q - zwycięzca turnieju.

Proces wyboru osobników do populacji Q może być realizowany dwoma sposobami:

w losowaniu bez zwracania

w losowaniu ze zwracaniem.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji gdzie:

Q - populacja wybrana z populacji potomnej do przeprowadzenia turnieju.

X^q - zwycięzca turnieju.

Proces wyboru osobników do populacji Q może być realizowany dwoma sposobami:

w losowaniu bez zwracania w losowaniu ze zwracaniem.

prawdopodobieństwem. Natomiast losowanie bez zwracania najgorszych osobników nie będzie w ogóle reprodukować.

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Metody selekcji gdzie:

Q - populacja wybrana z populacji potomnej do przeprowadzenia turnieju.

X^q - zwycięzca turnieju.

Proces wyboru osobników do populacji Q może być realizowany dwoma sposobami:

w losowaniu bez zwracania w losowaniu ze zwracaniem.

Projekt

Celem projektu jest przebadanie stopnia przydatności algorytmów ewolucyjnych do znajdowania optimów lokalnych i globalnych funkcji dwóch zmiennych.

Rys.9. Badana funkcja nr 1

Techniki ewolucyjne - algorytm genetyczny i nie tylko

Projekt

Rys.9. Badana funkcja nr 2

W dokumencie Techniki ewolucyjne - algorytm genetyczny i nie tylko (Stron 22-56)