% clc;
% clear all
data=load('burma14.txt');%eil51 Att48  burma14 st70 Oliver30 bayg29 ch150 ch130 pr226
cityCoor=[data(:,2) data(:,3)];%³ÇÊÐ×ø±ê¾ØÕó*
nMax=500;                     %½ø»¯´ÎÊý
indiNumber=20;                 %¸öÌåÊýÄ¿
Nc=10;                        %´ÎÊý
n=size(cityCoor,1);            %³ÇÊÐÊýÄ¿*
cityDist=zeros(n,n);           %³ÇÊоàÀë¾ØÕó*
Gbest=zeros(1,nMax);      %¸÷´úµÄ×îÓÅÖµ*
for i=1:n
    for j=1:n
        if i~=j
            cityDist(i,j)=((cityCoor(i,1)-cityCoor(j,1))^2+...
                (cityCoor(i,2)-cityCoor(j,2))^2)^0.5;
        end
        cityDist(j,i)=cityDist(i,j);
    end
end
individual=zeros(indiNumber,n);
for i=1:indiNumber
    individual(i,:)=randperm(n); %³õʼ»¯ÖÖȺ  
end
individual;d
indiFit=fitness(individual,cityCoor,cityDist);%¼ÆËãÖÖȺÊÊÓ¦¶È
[value,index]=min(indiFit);
Gbest=value;
deta=zeros(indiNumber,n);
tourGbest=individual(index,:);         %µ±Ç°ÈºÌå×îÓÅ
%%
for K=1:nMax
          for i=1:indiNumber
           indiFit(i)=fitness(individual(i,:),cityCoor,cityDist);%¼ÆËãÖÖȺÊÊÓ¦¶È
           [value,index]=min(indiFit) ;
           Gbest(K)=value;
           tourGbest=individual(index,:);
          end  
%% ²½³¤µü´ú
   for N=1:Nc
%         deta=randint(indiNumber,n);
        for i=1:indiNumber      
                r=sqrt(sum((individual(i,1)-individual(i,2)).^2));
                if (fitness(individual(i,1))>fitness(individual(i,2)))%¼ÆËãÖÖȺÊÊÓ¦¶È
                          gamma=1;betamin=0.2;beta0=1; beta=(beta0-betamin)*exp(-gamma*r.^2)+betamin;
                         tmpf=0.5*(rand(1,d)-0.5).*n;
                      individual1(i,1)=individual(i,1).*(1-beta)+individual(i,2).*beta+tmpf
                end         
             individual1=DupCK(individual1,i,n);  %µ÷ÓòéÖغ¯Êý     
        end
   end
   %%

end
   %% ½á¹û×÷ͼ
   [Bm,Bma]=min(Gbest)
plot(1:nMax,Gbest,'--c')
title('Ñ°ÓŹý³Ì')
xlabel('µü´ú´ÎÊý')
ylabel('ÊÊÓ¦¶ÈÖµ')
hold on



function indiFit=fitness(individual,cityCoor,cityDist)
%% ¸Ãº¯ÊýÓÃÓÚ¼ÆËã¸öÌåÊÊÓ¦¶ÈÖµ
%x           input     ¸öÌå
%cityCoor    input     ³ÇÊÐ×ø±ê
%cityDist    input     ³ÇÊоàÀë
%indiFit     output    ¸öÌåÊÊÓ¦¶ÈÖµ

m=size(individual,1);
n=size(cityCoor,1);
indiFit=zeros(m,1);
for i=1:m
    for j=1:n-1
    indiFit(i)=indiFit(i)+cityDist(individual(i,j),individual(i,j+1));
    end
    indiFit(i)=indiFit(i)+cityDist(individual(i,1),individual(i,n));
end

%%  ²é¿´Öظ´µÄºÍÒÅ©µÄ
function individual=DupCK(individual,i,n)

ww=tabulate(individual(i,:));
%ͳ¼ÆÊýÖµÐÍÊý×éÖи÷ÔªËسöÏÖµÄƵÂÊ£¬ÆµÊý
[w,t,m]=unique(individual(i,:));
%tÌåÏÖwÖÐÔªËØÔÚÔ­ÏòÁ¿ÖеÄλÖã¬mÌåÏÖÔ­ÏòÁ¿ÔÚwÖеÄλÖÃ
            ll=sort(t);
            x=tabulate(ll);

            l1=length(ww(:,2));
              if l1<n
                  for k=l1+1:n
                      ww(k,1)=k;
                      ww(k,2)=0;
                   end
              end

            l=length(x(:,2));
              if l<n
                  for k=l+1:n
                      x(k,1)=k;
                      x(k,2)=0;
                   end
              end
    r=0;r1=0;

        for s=1:n
                if x(s,2)==0
                    r=r+1;
                   y(r)=x(s,1) ;      %ÔÚÕâЩλÖÃÉÏÖظ´µÄ
                end
                if ww(s,2)==0
                    r1=r1+1;
                    w1(r1)=ww(s,1);   %ûÓгöÏֵijÇÊУ¬¼´Â©µôµÄ³ÇÊбàºÅ
                end
        end
       for r2=1:r1
           individual(i,y(r2))=w1(r2);  
       end

2Â¥: Originally posted by somomo91 at 2017-07-03 06:18:57
ûÓÐÊäÈëÊý¾Ý£¬Ã»·¨ÔËÐгÌÐòµÄ
¿ÉÒÔÌṩһ¸ö¼òµ¥µÄÑù±¾

2Â¥: Originally posted by somomo91 at 2017-07-03 06:18:57
ûÓÐÊäÈëÊý¾Ý£¬Ã»·¨ÔËÐгÌÐòµÄ
¿ÉÒÔÌṩһ¸ö¼òµ¥µÄÑù±¾

for N=1:Nc
         deta=randint(indiNumber,n,[-a,a]);
        for i=1:indiNumber  
            for j=1:indiNumber
                r=sqrt(sum((individual(i,:)-individual(j,:)).^2));
                if (indiFit(individual(i,:))<indiFit(individual(j,:)))%¼ÆËãÖÖȺÊÊÓ¦¶È
                          gamma=1;betamin=0.2;beta0=1; beta=(beta0-betamin)*exp(-gamma*r.^2)+betamin;
                         tmpf=0.5*(rand(1,d)-0.5).*n;
                      individual(i,:)=individual(i,:).*(1-beta)+individual(j,:).*beta+deta;
                      while individual(i,:)>n
                          individual(i,:)=individual(i,:)-n;
                   end
                    while individual(i,:)<1
                          individual(i,:)=individual1(i,:)+n;
                    end  
                end         
               individual=DupCK(individual,i,n);  %µ÷ÓòéÖغ¯Êý
            end
        end
%    end