clear
clc
tic;
gen=500;N=65;L=22;n=10;pm=0.1;fat=0.2;
P=init(N,L);
xmin=0;xmax=9;
f='x+10*sin(x.*5)+7*cos(x.*4)';
vfx=[];vmfit=[];
it=0;
while it<gen
    it=it+1;
    x=decode(P(:,1:22),xmin,xmax);
    fit=eval(f);
    if it==1
        figure(1);fplot(f,[xmin,xmax]);grid on;hold on;
        plot(x,fit,'k*');title('(a)抗體的初始位置分布圖');xlabel('x');ylabel('y');
    end
    if it==gen
        figure(2);fplot(f,[xmin,xmax]);grid on;hold on;
        plot(x,fit,'k*');title('(b)抗體的最終位置分布圖');xlabel('x');ylabel('y');
    end
    T=[];
    [a,ind]=sort(fit);
    valx=x(ind(end-n+1:end));
    fx=a(end-n+1:end);
    vfx=[vfx fx(end)];
    [T,pcs]=reprod(n,fat,N,ind,P,T);
    T=Hypermutation(T,L,pm,xmax,xmin);
    T(pcs,:)=P(fliplr(ind(end-n+1:end)),:);
     x=decode(T(:,1:22),xmin,xmax);
     fit=eval(f);pcs=[0 pcs];vmfit=[vmfit mean(fit)];
     for i=1:n
         [out(i),bcs(i)]=max(fit(pcs(i)+1:pcs(i+1)));
         bcs(i)=bcs(i)+pcs(i);
     end
     P(fliplr(ind(end-n+1:end)),:)=T(bcs,:);
end
 [aa,bb]=max(fit)
 x(bb)
 toc   
 
 
 function ab=init(n,s)
ab=2.*rand(n,s)-1;
ab=hardlim(ab);  

function x=decode(v,xmin,xmax)
v=fliplr(v);
s=size(v);
aux=0:1:21;aux=ones(s(1),1)*aux;
x1=sum((v.*2.^aux)');
x=xmin+(xmax-xmin)*x1./4194303;

function T=Hypermutation(T,L,pm,xmax,xmin);
M=rand(size(T,1),L)<=pm;
M=T-2.*(T.*M)+M;
k=round(log(10*(xmax-xmin)));
k=1;
T(:,k:L)=M(:,k:L);

function [T,pcs]=reprod(n,fat,N,ind,P,T)
if n==1
    cs=N;
    T=ones(N,1)*P(ind(1),:);
else
    for i=1:n
        %cs(i)=ceil(fat*N/i);
        cs(i)=round(fat*N);
        pcs(i)=sum(cs);
        T=[T;ones(cs(i),1)*P(ind(end-i+1),:)];
    end
end