crossover.m

来自「ga三中操作算子的改进」· M 代码 · 共 23 行

function [newpop,newquanzhi]=crossover(pop,pc,quanzhi)
%交叉
[px,py]=size(pop);
newpop=zeros(size(pop));
newquanzhi=zeros(size(quanzhi));
for i=1:2:px-1
    if (rand<pc)
       cpoint=round(rand*py);
       if cpoint<=0
          cpoint=1;
       end
       newpop(i,:)=[pop(i,1:cpoint) pop(i+1,cpoint+1:py)];
       newpop(i+1,:)=[pop(i+1,1:cpoint) pop(i,cpoint+1:py)];
       newquanzhi(i,:)=[quanzhi(i,1:cpoint) quanzhi(i+1,cpoint+1:py)];
       newquanzhi(i+1,:)=[quanzhi(i+1,1:cpoint) quanzhi(i,cpoint+1:py)];
    else
       newpop(i,:)=pop(i,:);
       newpop(i+1,:)=pop(i+1,:);
       newquanzhi(i,:)=quanzhi(i,:);
       newquanzhi(i+1,:)=quanzhi(i+1,:);       
    end
end