%定義遺傳算法參數
NIND=40;               %個體數目(Numbe of individuals)
MAXGEN=500;            %最大遺傳代數(Maximum number of generations)
NVAR=20;               %變量的維數
PRECI=20;              %變量的二進制位數(Precision of variables)
GGAP=0.9;              %代溝(Generation gap)
trace=zeros(MAXGEN, 2);
%建立區域描述器(Build field descriptor)
FieldD=[rep([PRECI],[1,NVAR]);rep([-512;512],[1, NVAR]);rep([1;0;1;1],[1,NVAR])];
Chrom=crtbp(NIND, NVAR*PRECI);                       %創建初始種群
gen=0;                                               %代計數器
ObjV=objfun1(bs2rv(Chrom, FieldD));                  %計算初始種群個體的目標函數值
while gen<MAXGEN                                     %迭代
    FitnV=ranking(ObjV);                             %分配適應度值(Assign fitness values)
    SelCh=select('sus', Chrom, FitnV, GGAP);         %選擇
    SelCh=recombin('xovsp', SelCh, 0.7);             %重組
    SelCh=mut(SelCh);                                %變異
    ObjVSel=objfun1(bs2rv(SelCh, FieldD));           %計算子代目標函數值 
    [Chrom ObjV]=reins(Chrom, SelCh, 1, 1, ObjV, ObjVSel);     %重插入
    gen=gen+1;                                                 %代計數器增加
    trace(gen, 1)=min(ObjV);                                   %遺傳算法性能跟蹤
    trace(gen, 2)=sum(ObjV)/length(ObjV);
end
plot(trace(:,1));hold on;
plot(trace(:,2),'-.');grid;
legend(' 種群均值的變化','解的變化')
%輸出最優解及其對應的20個自變量的十進制值,Y為最優解,I為種群的序號
[Y, I]=min(ObjV)
X=bs2rv(Chrom, FieldD);
X(I,:)