function cro=cros(nind,old,n,pc)
a=randint(1,1,9);
b=rand(1);
for i=1:2:nind-1
if a==0|a==1
cro(i,:)=old(i,:);
cro(i+1,:)=old(i+1,:);
elseif a>=2&b>=pc
A=[]; B=[];
for k=1:n
for j=a+1:n
if old(i+1,k)==old(i,j)
A=[A,old(i+1,k)];
end
end
end
cro(i,:)=[old(i,1:a),A];
for m=1:n
for l=a+1:n
if old(i,m)==old(i+1,l)
B=[B,old(i,m)];
end
end
end
cro(i+1,:)=[old(i+1,1:a),B];
elseif a>=2&b cro(i,:)=old(i,:);
cro(i+1,:)=old(i+1,:);
end
end
if rem(nind,2)==0 %求余
cro;
else
cro=[cro;old(nind,:)];
end
cro;
%这个函数是交叉函数,nind指种群规模,即个体数目,
%old指可行解矩阵,即个体解矩阵,
%n指个体长度,即指派问题的任务数,
%pc指交叉概率。
%cros(6,[1 2 3 4 ;3 2 1 4;1 3 2 4;3 1 4 2;4 3 2 1;4 2 3 1],4,1)
%cros(6,[1 2 3 4 ;3 2 1 4;1 3 2 4;3 1 4 2;4 3 2 1;4 2 3 1],4,0.6)