function [pops]=changes(cpop,bounds,len,p)
%基因突變函數
%function [pops]=changes(pop,bounds,len,p)
%pop        種群數目
%bounds     邊界約束
%len        每個變量的編碼長度
%           如len為[4 3 3];表示有三個變量，第一個變量的二進制編碼長度為4,依次類推
%p          突變概率
%pops       返回突變后的基因
%p1         基因突變數目

if isempty(p)
    p=0.01;
end
[n,m]=size(cpop);
pop=cpop;
p1=round(sum(len)*n*p);
k=0;q=[];v=[];
while(k<p1)
    k=k+1;
    q(k)=round(rand*(sum(len)*n-1))+1;
    for i=1:k-1
        if q(k)==q(i)
            q(k)=[];
            k=k-1;
        end
    end
end

for i=1:n
[B(i,:),len]=B2F(pop(i,:),bounds);
end

v=reshape(B,1,n*sum(len));

for i=1:p1
    if v(q(i))==0
        v(q(i))=1;
    else
        v(q(i))=0;
    end
end
v=reshape(v,n,sum(len));
for i=1:n
   pop(i,:)=F2B(v(i,:),bounds,len);
end
pops=pop
cpop