mutation.m

来自「遗传算法得一种编译程序,不是很正统」· M 代码 · 共 19 行

function [newpop]=mutation(pop,pm)
[px,py]=size(pop);
newpop=ones(size(pop));
for i=1:px
    if (rand<pm)
        mpoint=round(rand*py);
    if  mpoint<=0
        mpoint=1;
    end
    newpop(i)=pop(i);
    if any(newpop(i,mpoint))==0
        newpop(i,mpoint)=1;
    else
        newpop(i,mpoint)=0;
    end
    else
        newpop(i)=pop(i);
    end
end