newroute.asv

TSP问题,即旅行商问题,matlab环境下可以运行

function [x,changex,l]=newroute(x0,l0)
global customer;
global depot;
global demand;


x0=[1 3 0
    1 10 0
   3 1 8 
   3  6 4
   3 2 7 
   3 5 9 ];
l0=[1 1 2 2 2 2];
demand=[142 85 385
       163 175 430
        87 96 540
        63 57 324
        204 217 238
        130 165 296
        180 210 340
        93 114 286
        170 163 392
        78 103 310
        50 148 0
        34 82 0
        185 198 0]; 
    
NumberRoute=length(l0);  
a=ceil(NumberRoute*rand);  
  %从a上选择顾客  
c1_index=ceil(l0(a)*rand);
c1=x0(a,c1_index+1);
%将c1路线a删除
if l0(a)==1
    x1(a,1)=x0(a,1);
    x1(a,2)=0;
else
   x1(a,1:c1_index)=x0(a,1:c1_index);
   x1(a,c1_index+1:l0(a))=x0(a,c1_index+2:l0(a)+1);
end
l1(a)=l0(a)-1;
for j=1:depot
    alldis(j)=distance(demand,c1,depot+customer);
end
[dis,n]=min(alldis);

%增一新路

for v=1:NumberRoute
    if v~=a 
        x1(v,1:l0(v)+1)=x0(v,1:l0(v)+1);
        l1(v)=l0(v);
    end
end


x1(NumberRoute+1,1:2)=[n c1];
l1(NumberRoute+1)=1;
NumberRoute=NumberRoute+1;
if l1(a)==0
    x(1:a-1,:)=x1(1:a-1,:);
    x(a:NumberRoute-1,:)=x1(a+1:NumberRoute,:);
    l(1:a-1)=l1(1:a-1);
    l(a:NumberRoute-1)=l1(a+1:NumberRoute);
else
    x=x1;
    l=l1;
end


changex=[5 c1 0];