function NewChrIx = sus(FitnV,Nsel);
% FitnV-種群中每個染色體的適值
% Nsel-選擇的染色體數目

[Nind,ans] = size(FitnV);
% 進行隨機采樣
cumfit = cumsum(FitnV);
trials = cumfit(Nind) / Nsel * (rand + (0:Nsel-1)');
Mf = cumfit(:, ones(1, Nsel));
Mt = trials(:, ones(1, Nind))';
[NewChrIx, ans] = find(Mt < Mf & [ zeros(1, Nsel); Mf(1:Nind-1, :) ] <= Mt);

%得到新的種群
[ans, shuf] = sort(rand(Nsel, 1));
NewChrIx = NewChrIx(shuf);