%定義遺傳算法參數
NIND=40;                    %個體數目(Number of individuals)
MAXGEN=400;                 %最大遺傳代數(Maximum number of generations)
GGAP=0.9;                   %代溝(Generation gap)
trace=zeros(MAXGEN,2);      %遺傳算法性能跟蹤初始值
BaseV=crtbase(15,8);
Chrom=crtbp(NIND, BaseV)+ones(NIND,15);    %初始種群
gen=0;
ObjV=targetalloc(Chrom);                   %計算初始種群函數值
while gen<MAXGEN
    FitnV=ranking(-ObjV);                  %分配適應度值(Assign fitness values)
    SelCh=select('sus',Chrom,FitnV,GGAP);               %選擇
    SelCh=recombin('xovsp',SelCh,0.7);                  %重組
    f=rep([1;8],[1,15]);
    SelCh=mutbga(SelCh, f);SelCh=fix(SelCh);            %變異
    ObjVSel=targetalloc(SelCh);                         %計算子代目標函數值
    [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);   %重插入
    gen=gen+1;
    trace(gen,1)=max(ObjV);                             %遺傳算法性能跟蹤
    trace(gen,2)=sum(ObjV)/length(ObjV);
end
[Y, I]=max(ObjV);Chrom(I,:),Y                           %最優解及其目標函數值
plot(trace(:,1),'-.');hold on;
plot(trace(:,2));grid
legend('解的變化','種群均值的變化')