clc;
clear?all;
figure（1）;
fplot（‘variable.*sin（10*pi*variable）+2.0‘[-12]）;???%畫出函數曲線

tic?????????????%計算程序運行時間，和toc連用，tic表示開始，toc表示結束
NIND=40;
MAXGEN=25;
PRECI=10;
GGAP=0.9;
trace=zeros（2MAXGEN）;
FieldD=[PRECI;-1;2;1;0;1;1];???%該處包含編碼長度、編碼范圍等信息
Chrom=crtbp（NINDPRECI）;
gen=0;
variable=bs2rv（ChromFieldD）;
ObjV=variable.*sin（10*pi*variable）+2.0;
while?gen????FitnV=ranking（-ObjV）;
????SelCh=select（‘sus‘ChromFitnVGGAP）;
????SelCh=recombin（‘xovsp‘SelCh0.7）;
????SelCh=mut（SelCh）;
????variable=bs2rv（SelChFieldD）;
????ObjVSel=variable.*sin（10*pi*variable）+2.0;
????[Chrom?ObjV]=reins（ChromSelCh11ObjVObjVSel）;
????gen=gen+1;
????[YI]=max（ObjVSel）;hold?on;
????plot（variable（I）Y‘ro‘）;
????trace（1gen）=max（ObjVSel）;
????trace（2gen）=sum（ObjVSel）/length（ObjVSel）;
end
t=toc

variable=bs2rv（ChromFieldD）;
hold?ongrid;
plot（variable‘ObjV‘‘b*‘）;
figure（2）;
plot（trace（1:）‘）;
hold?on;
plot（trace（2:）‘‘-.‘）;grid;
%legend（‘解的變化‘‘種群均值的變化‘）;