matlab

clear;
S1=[ ones(1,5)*(-500);
    ones(1,5)*500  ];
popsize_sch=30;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
S2=[ ones(1,2)*(-100);
    ones(1,2)*100  ];


jedinci_sch=genrpop(popsize_sch,S1);
jedinci_kvad=genrpop(popsize_sch,S2);


posledna_generacia=200;
amp=ones(1,10)*1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

fit_matica_kvad=1:30;

for i=1:posledna_generacia
       fit_matica=schwefel(jedinci_sch);
       nova_generacia=selbest(jedinci_sch,fit_matica,[3 2 1]);
        ostatny=selsus(jedinci_sch,fit_matica,24);
        krizenie=crossov(ostatny,2,0);
        mutac=mutx(krizenie,0.1,S1);
        mutacia=muta(mutac,0.1,amp,S1);
        jedinci_sch=[nova_generacia;mutacia];

        graf(i)=min(fit_matica);
end
plot(graf);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%




for i=1:posledna_generacia
       fit_matica_kvad=schwefel(jedinci_sch);
       nova_generacia=selbest(jedinci_sch,fit_matica,[3 2 1]);
        ostatny=selsus(jedinci_sch,fit_matica,24);
        krizenie=crossov(ostatny,2,0);
        mutac=mutx(krizenie,0.1,S1);
        mutacia=muta(mutac,0.1,amp,S1);
        jedinci_sch=[nova_generacia;mutacia];

        graf(i)=min(fit_matica);
end


plot(graf);
hold on;
fit_matica