crossover.m

遗传算法的代码,可以仿照这个学习一下,希望有所帮助

%--------------------------------------------------
%                    交叉算子             
%--------------------------------------------------

function new_gen = crossover(old_gen,pc)

[m,n] = size(old_gen);
[nouse,mating] = sort(rand(m,1),1);
mat_gen = old_gen(mating,:);
pairs = size(mat_gen,1) / 2;
bits = size(mat_gen,2);
cpairs = rand(pairs,1) < pc;
cpoints = randint(pairs,1,[1,bits]);
cpoints = cpairs.*cpoints;
for i = 1:pairs
    new_gen([2*i-1 2*i],:) = [mat_gen([2*i-1 2*i],1:cpoints(i)) mat_gen([2*i 2*i-1],cpoints(i)+1:bits)];
end