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