%變異
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