%定義遺傳算法參數
Dim=20;               %變量維數
NIND=20;              %個體數目(Number of individuals)
Preci=20;             %變量的二進制位數(Precision of variables)
MAXGEN=100;            %最大遺傳代數(Maximum number of generations)
GGAP=0.8;             %代溝(Generation gap)
SEL_F='sus';          %選擇函數名
XOV_F='xovsp';        %重組函數名
MUT_F='mut';          %變異函數名
OBJ_F='objdopi';      %目標函數名
FieldDR=feval(OBJ_F,[],1);                    %計算目標函數值
%建立區域描述器(Build field descriptor)
FieldDD=[rep([Preci],[1,Dim]);FieldDR;rep([1;0;1;1],[1,Dim])];
Chrom=crtbp(NIND, Dim*Preci);                 %創建初始種群
gen=0;
Best=NaN*ones(MAXGEN,1);                      %最優解初值
while gen<MAXGEN                              %最大循環次數
    ObjV=feval(OBJ_F,bs2rv(Chrom,FieldDD));   %計算目標函數值
    Best(gen+1)=min(ObjV);                    %最優解
    plot(log10(Best),'bo');
    FitnV=ranking(ObjV);                      %分配適應度值(Assign fitness values)
    SelCh=select(SEL_F,Chrom,FitnV,GGAP);     %選擇
    SelCh=recombin(XOV_F,SelCh);              %重組
    SelCh=mutate(MUT_F,SelCh);                %變異
    Chrom=reins(Chrom,SelCh);                 %重插入
    gen=gen+1;
end
grid;
xlabel('迭代次數');ylabel('目標函數值(取對數)');