figure(1);
fplot('variable.*sin(10*pi*variable)+2.0',[-1,2]);   %畫(huà)出函數(shù)曲線(xiàn)
%定義遺傳算法參數(shù)
NIND=40;        %個(gè)體數(shù)目(Number of individuals)
MAXGEN=25;      %最大遺傳代數(shù)(Maximum number of generations)
PRECI=20;       %變量的二進(jìn)制位數(shù)(Precision of variables)
GGAP=0.9;       %代溝(Generation gap)
trace=zeros(2, MAXGEN);                        %尋優(yōu)結(jié)果的初始值
FieldD=[20;-1;2;1;0;1;1];                      %區(qū)域描述器(Build field descriptor)
Chrom=crtbp(NIND, PRECI);                      %初始種群
gen=0;                                         %代計(jì)數(shù)器
variable=bs2rv(Chrom, FieldD);                 %計(jì)算初始種群的十進(jìn)制轉(zhuǎn)換
ObjV=variable.*sin(10*pi*variable)+2.0;        %計(jì)算目標(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);                                      %變異
   variable=bs2rv(SelCh, FieldD);                         %子代個(gè)體的十進(jìn)制轉(zhuǎn)換
   ObjVSel=variable.*sin(10*pi*variable)+2.0;             %計(jì)算子代的目標(biāo)函數(shù)值
   [Chrom ObjV]=reins(Chrom, SelCh, 1, 1, ObjV, ObjVSel); %重插入子代的新種群
   variable=bs2rv(Chrom, FieldD);
   gen=gen+1;                                             %代計(jì)數(shù)器增加
   %輸出最優(yōu)解及其序號(hào)，并在目標(biāo)函數(shù)圖像中標(biāo)出，Y為最優(yōu)解,I為種群的序號(hào)
   [Y, I]=max(ObjV);hold on;
   plot(variable(I), Y, 'bo');
   trace(1, gen)=max(ObjV);                               %遺傳算法性能跟蹤
   trace(2, gen)=sum(ObjV)/length(ObjV);
end
variable=bs2rv(Chrom, FieldD);                            %最優(yōu)個(gè)體的十進(jìn)制轉(zhuǎn)換
hold on, grid;
plot(variable,ObjV,'b*');
figure(2);
plot(trace(1,:));
hold on;
plot(trace(2,:),'-.');grid
legend('解的變化','種群均值的變化')