VPF::[Математика][Pascal]Генетический алгоритм

v4v4n

Дата 2.5.2009, 19:01 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 45
Регистрация: 26.9.2008

Репутация: 2
Всего: 4

Есть задача:
Селекция - это оператор, посредством которого индивиды (т.е. хромосомы) выбираются для спаривания и порождения потомков. Индивиды с более высокой пригодностью должны выбираться с большей вероятностью. Существует большое число различных моделей селекции, некоторые из которых не имеют биологических аналогов. Большинство схем селекций, используемых в ГА, создают промежуточную популяцию и затем выбирают из нее случайным образом пары индивидов для скрещивания...
...
Ранговая селекция
Индивиды сортируются (ранжируются) на основе их пригодности таким образом, чтобы f(i)>f(j) для i > j.
Затем каждому индивиду назначается вероятность pi быть отобранным, взятая из заданного распределения с ограничением
∑pi = 1
Типичные распределения:
Линейное: pi = a * i + b(a < 0).
Отрицательное экспоненциальное: pi = a * exp(b * i + c).

реализую процедуру ранговой селекции на паскале.

Как вычислить коэффициенты a, b, c ? И как потом ими пользоваться для отбора индивидов?

v4v4n

Дата 2.5.2009, 20:08 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 45
Регистрация: 26.9.2008

Репутация: 2
Всего: 4

Тупо, но работает smile

Код


procedure select(popsize:integer; var pop:TPopulation);
var i : integer;
    newpop : TPopulation;
    a,b : float;

    procedure Sort;
    var ind : TIndividual;
        sorted : boolean;
        i : integer;
    Begin
        sorted:=false;
        While not sorted do begin
            sorted:=true;
            for i:=1 to popsize-1 do
                if pop[i].fitness<pop[i+1].fitness then begin
                   ind:=pop[i];
                   pop[i]:=pop[i+1];
                   pop[i+1]:=ind;
                   sorted:=false;
                end;
        end;
    End;
    function selection:integer;
    var i:integer; z, summ : float;
    Begin
        z:=(1-random); summ:=0; i:=0;
        repeat
           inc(i);
           summ:=summ+a*i+b;
        until summ>=z;
        selection:=i;
    End;
Begin
    a:=2/(sqr(popsize)+popsize);    b:=0;
    sort;
    for i:=1 to popsize do newpop[i]:=pop[selection];
    pop:=newpop;
End;

Подскажите рационально получилось? Или это вообще на ранговую селекцию не тянет?

Это сообщение отредактировал(а) v4v4n - 2.5.2009, 20:17

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »