CCC / лаби / лаб5 / Текстовый документ

clc; clear;
n=2;
X=zeros(n,5);
Xmin=-3;
Xmax=1;
for i=1:n
X(i,:)=Xmin+(Xmax-Xmin)*rand(1,5);
end;
for i=1:n
F1(i)=0;
for j=1:5
F1(i)=F1(i)+j*X(i,j)^2;
end;
end;
Z1=[X,F1'];
tic;
p=10;
Pmut=0.59;
best=[];
mut_level=0.01;
while toc<1
for i=1:n
for j=1:p
Y(p*(i-1)+j,:)=X(i,:);
end;
end;

for i=1:n*p
for j=1:5
if rand(1,1)<Pmut
Y(i,j)=Y(i,j)+mut_level*randn(1,1);
end;
end;
end;
for i=1:n*p
F(i)=0;
for j=1:5
F(i)=F(i)+j*Y(i,j)^2;
end;
end;
Z=[Y,F'];
Rezult=[Z1;Z];
[g,q]=size(Rezult);
Rezult=sortrows(Rezult,q);
for i=1:n
for j=1:5
X(i,j)=Rezult(i,j);
end;
end;
best=[best;Rezult(1,q)];

end;
disp(' koren raven ');
X(1,:)
disp(' celevaja funkcija ravna ');
min(best)
semilogy(best);