function s1=select(v0,f,population)
%輪盤賭選擇

matrixsum=sum(f);
roulette=cumsum(f);%求向量累加
choiseP=roulette/matrixsum;%得選擇概率

for i=1:population
    
    tempP=rand(1);
    
        if (tempP<choiseP(1))
            selected(i,:)=v0(1,:);
           else
              for j=2:length(choiseP)
                 if ((tempP>=choiseP(j-1))&(tempP<=choiseP(j)))
                    selected(i,:)=v0(j,:);
                    break;
                 end
           end
        end
  end
s1=selected;