ga_demo.m

patten regnization source从1-14章能运行

% 蜡傈磊 舅绊府硫阑 捞侩茄 鉴雀 魄概盔 巩力(TSP)狼 秦搬
% N: 档矫狼 荐. 档矫甸篮 1俊辑 N栏肺 扼骇傅. 
clear all
N=input('鉴雀且 档矫狼 荐: ');
%======================================================================
% N俺狼 档矫俊 措茄 TSP 巩力 积己 
pos=rand(N,2);  % 窜困 沥荤阿屈郴俊 档矫狼 困摹甫 烙狼肺 汲沥 
distance=dist(pos,pos); %档矫埃狼 芭府 青纺  
figure(1),clf
plot(pos(:,1),pos(:,2),'o') %档矫 困摹 钎矫
axis([0 1 0 1]), axis square, hold on
%弊覆惑俊 档矫 锅龋 扼骇傅
for i=1:N
   text(pos(i,1)+.03,pos(i,2),int2str(i));
end
hold off
%========================================================================
% 蜡傈磊 沤祸 
disp('********** 蜡傈磊 舅绊府硫*****************');
ngen=input('柳拳矫懦 技措荐: ');
ngpool=input('# 蜡傈磊 钱(gene pool)郴狼 俺眉(chromosoms)荐: '); % 蜡傈磊 钱狼 农扁
gpool=zeros(ngpool,N); % 蜡傈磊 钱(gene pool)
% 蜡傈磊 钱 檬扁拳 
for i=1:ngpool, 
   gpool(i,:)=[1 randomize([2:N]')'];
end
oldstep=0;
oldcostmin=99999;
costmin=N, tourmin=zeros(1,N); cost=zeros(1,ngpool); 
oldtourmin=gpool(2,:);
%========================================================================
for step=1:ngen, 
   % 1窜拌. 泅犁 俺眉狼 利钦档 乞啊, 鉴雀 芭府啊 利钦档啊 登哥 陋阑荐废 亮促.  
   for i=1:ngpool,
      cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
   end
   % 泅犁 弥惑狼 菩胶甫 扁废 
   [costmin,idx]=min(cost);
   tourmin=gpool(idx,:);
   % 2窜拌. 背瞒(crossover)客 倒楷函捞(mutation) 楷魂 荐青 
   [csort,ridx]=sort(cost); % 厚侩捞 利篮巴俊辑 奴巴栏肺 沥纺茄促. 
   
   for i=2:ngpool/2 
       sameidx=[gpool(2*i-1,:)==gpool(2*i,:)]; 
       diffidx=find(sameidx==0); 
       if length(diffidx)<=2 
           % 倒楷函捞积己
           gpool(2*i,:)=[1 randomize([2:N]')']; 
       else
           if i%3 == 0 % 何盒背眉背硅
             [gpool(2*i-1,:),gpool(2*i,:)]=PMX(gpool(ridx(2*i-1),:),gpool(ridx(2*i),:));         
           end
           if i%3 == 1  % 鉴辑背硅
             [gpool(2*i-1,:),gpool(2*i,:)]=OX(gpool(ridx(2*i-1),:),gpool(ridx(2*i),:));
           end
           if i%3 == 2 % 林扁背硅
             [gpool(2*i-1,:),gpool(2*i,:)]=CX(gpool(ridx(2*i-1),:),gpool(ridx(2*i),:));
           end
       end 
   end 
   
   %pause;
   %弥利鉴雀 风飘啊 函窍绰 何盒父 免仿秦林扁 困茄 贸府 
   sameidx=[oldtourmin==tourmin]; % 鞍篮 镑俊绰 #1肺 付胶农甫 竞快绊                            
   diffidx=find(sameidx==0); % 促弗 牢郸胶甫 茫酒辑 
   if length(diffidx) ~= 0 & costmin < oldcostmin     
       figure;
       plot(pos(:,1),pos(:,2),'o'),axis([0 1 0 1]), axis square, hold on 
       trip=pos(tourmin',:); trip=[trip; trip(1,:)]; 
       plot(trip(:,1),trip(:,2),'-');
       title(['技措:' num2str(step) ',  蜡傈磊 八祸, 厚侩:' num2str(costmin)])
       hold off;
       drawnow
       disp(['技措' num2str(step) '::弥利鉴雀:'  num2str(tourmin)  '::厚侩:'  num2str(costmin) ]);
      oldcostmin = costmin;
      oldtourmin = tourmin;
      oldstep=step;
   end
end

disp(['厚侩 窃荐 乞啊 : ' int2str(ngen*2+ngpool) ' times!'])
disp(['弥家 鉴雀 芭府 : ' num2str(oldcostmin)])
disp('弥利 鉴雀 : ') 
disp(num2str(oldtourmin))

figure
plot(pos(:,1),pos(:,2),'o'),axis([0 1 0 1]), axis square, hold on
trip=pos(oldtourmin',:); trip=[trip; trip(1,:)];
plot(trip(:,1),trip(:,2),'-')
title(['技措:' num2str(oldstep) ',  蜡傈磊 八祸, 厚侩:' num2str(oldcostmin)])
hold off