crossgens.m

来自「求解TSP问题的遗传算法matlab程序」· M 代码 · 共 37 行

function newroad=crossgens(road1,road2,road3,Map)
n=length(road1);
road=[road1;road2;road3];
point=road1(1);
newroad=[point];
for k=2:n
  current=[];
  next=[];
  edge=[];
  for i=1:3
      current(i)=find(road(i,:)==point);
      if current(i)==n
         next(i)=1;
      else
         next(i)=current(i)+1;
      end
      if sum(newroad==road(i,next(i)))==0 %next point already appear
     
        edge(i)=distance(Map(point,:),Map(road(i,next(i)),:));
      else
        edge(i)=inf;
      end
 end
 [minedge,i]=min(edge);
 if minedge<inf
     point=road(i,next(i));
 else
     for j=1:n
         if sum(newroad==j)==0 %j not appear
             point=j;
             break;
         end
     end
 end
 newroad=[newroad,point];
end