function NewChrom = mut(OldChrom,Pm,BaseV)
% OldChrom-當代種群的染色體
% Pm-變異概率
% BaseV-染色體基

[Nind, Lind] = size(OldChrom) ;
if nargin < 2, Pm = 0.7/Lind ; end
if isnan(Pm), Pm = 0.7/Lind; end
if (nargin < 3), BaseV = crtbase(Lind);  end
if (isnan(BaseV)), BaseV = crtbase(Lind);  end
if (isempty(BaseV)), BaseV = crtbase(Lind);  end
if (nargin == 3) & (Lind ~= length(BaseV))
   error('OldChrom and BaseV are incompatible'), end
% 生成變異屏蔽矩陣
BaseM = BaseV(ones(Nind,1),:) ;

% 進行變異操作
NewChrom = rem(OldChrom+(rand(Nind,Lind)<Pm).*ceil(rand(Nind,Lind).*(BaseM-1)),BaseM);