Matlab 贪婪交叉算子遗传算法

贪婪交叉算子
function r=crossover(p,n,k,x,y)
r1=randperm(k,1);
r2=randperm(k,1);
while r1r2
r2=randperm(k,1);
end
l1=p(r1,:);
l2=p(r2,:);
r=randperm(n);
for p1=2:n
r(p1)=0;
end
for p1=2:n
m1=find(l1r(p1-1));
m2=find(l2r(p1-1));
if m1length(l1)
m1=1;
else
m1=m1+1;
end
if m2length(l2)
m2=1;
else
m2=m2+1;
end
if distance(l1(m1),r(p1-1),x,y)<=distance(l2(m2),r(p1-1),x,y)
r(p1)=l1(m1);
else
r(p1)=l2(m2);
end
m1=find(l1r(p1-1));
m2=find(l2==r(p1-1));
l1(m1)=[];
l2(m2)=[];
end
end