function [gt]=mutation(gt,pm) %均勻變異
%每條染色體獲得的變異概率是均等的  隨機產生兩點進行交換
[m,n]=size(gt);
ran=rand(1,m);  % 給每條染色體隨機分配一個0-1均勻分布的數
[x,mu]=find(ran<pm);   % 確定哪些染色體要變異
for k=1:length(mu)
    position=rand(1,2)*n;
    position=ceil(position);
    var=gt(mu(k),position(1));
    gt(mu(k),position(1))=gt(mu(k),position(2));
    gt(mu(k),position(2))=var;
end