NIND=100;                              %個體數(shù)目(Number of individuals)
MAXGEN=50;                             %最大遺傳代數(shù)(Maximum number of generations)
NVAR=2;                                %變量個數(shù)
PRECI=20;                              %變量的二進(jìn)制位數(shù)(Precision of variables)
GGAP=0.9;                              %代溝(Generation gap)
trace1=[];trace2=[];trace3=[];         %性能跟蹤
%建立區(qū)域描述器(Build field descriptor)
FieldD=[rep([PRECI],[1,NVAR]);[1,1;4,2];rep([1;0;1;1],[1,NVAR])];
Chrom=crtbp(NIND,NVAR*PRECI);          %初始種群
v=bs2rv(Chrom,FieldD);                 %初始種群十進(jìn)制轉(zhuǎn)換
gen=1;
while gen<MAXGEN
    ObjV1=f1(v(1:NIND/2,:));           %分組后第一目標(biāo)函數(shù)值
    FitnV1=ranking(ObjV1);             %分配適應(yīng)度值(Assign fitness values)
    SelCh1=select('sus',Chrom(1:NIND/2,:),FitnV1,GGAP);                 %選擇
    ObjV2=f2(v(NIND/2+1:NIND,:));      %分組后第二目標(biāo)函數(shù)值
    FitnV2=ranking(ObjV2);
    SelCh2=select('sus',Chrom(NIND/2+1:NIND,:),FitnV2,GGAP);            %選擇
    ObjV=[ObjV1;ObjV2];
    SelCh=[SelCh1;SelCh2];             %合并
    SelCh=recombin('xovsp',SelCh,0.7); %重組
    SelCh=mut(SelCh);
    v=bs2rv(SelCh,FieldD);
    ObjVSel1=f1(v(1:NIND/2*GGAP,:));
    ObjVSel2=f2(v(NIND/2*GGAP+1:NIND*GGAP,:));
    ObjVSel=[ObjVSel1;ObjVSel2];       %計(jì)算子代的兩個目標(biāo)函數(shù)值后合并在一起
    [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);                    %重插入子代的新種群
    v=bs2rv(Chrom,FieldD);
    
    trace1(gen,1)=min(f1(v));
    trace1(gen,2)=sum(f1(v))/length(f1(v));
    trace2(gen,1)=min(f2(v));
    trace2(gen,2)=sum(f2(v))/length(f2(v));
    trace3(gen,1)=min(f1(v)+f2(v));
    trace3(gen,2)=sum(f1(v))/length(f1(v))+sum(f2(v))/length(f2(v));
    gen=gen+1;
end
figure(1);clf;
plot(trace1(:,1));hold on;plot(trace1(:,2),'-.');
plot(trace1(:,1),'.');plot(trace1(:,2),'.');grid;
legend('解的變化','種群均值的變化')
xlabel('迭代次數(shù)');ylabel('第一目標(biāo)函數(shù)值');
figure(2);clf;
plot(trace2(:,1));hold on;
plot(trace2(:,2),'-.');
plot(trace2(:,1),'.');
plot(trace2(:,2),'.');grid;
legend('解的變化','種群均值的變化')
xlabel('迭代次數(shù)');ylabel('第二目標(biāo)函數(shù)值');
figure(3);clf;
plot(trace3(:,1));hold on;
plot(trace3(:,2),'-.');
plot(trace3(:,1),'.');
plot(trace3(:,2),'.');grid;
legend('解的變化','種群均值的變化')
xlabel('迭代次數(shù)');ylabel('目標(biāo)函數(shù)值之和');
figure(4);clf;plot(f1(v));hold on;
plot(f2(v),'r-.');grid;