mutation1.m

请大家看一看, 我编的这个用遗传算法求 f(x)=xsin(10pi*x)+2.0 x为-1到2区间的值

%mutation 遗传算法中变异操作的实现
function newpop=mutatiaon1(oldpop,pm)
[px,py]=size(oldpop);
for i=1:px
    if rand<pm
     msite=round(rand*py);
     if msite<=0
         msite=1;
     end
     if any(oldpop(i,msite))==0
         oldpop(i,msite)=1;
     else
         oldpop(i,msite)=0;
     end
   end
end
 newpop=oldpop;