%用遺傳算法進(jìn)行簡(jiǎn)單函數(shù)的優(yōu)化
clear

bn=22; %個(gè)體串長(zhǎng)度
inn=50; %初始種群大小
gnmax=200;  %最大代數(shù)
pc=0.75; %交叉概率
pm=0.05; %變異概率

%產(chǎn)生初始種群
s=round(rand(inn,bn));

%計(jì)算適應(yīng)度,返回適應(yīng)度f(wàn)和累積概率p
[f,p]=objf(s);  

gn=1;
while gn<gnmax+1
   for j=1:2:inn
      
      %選擇操作
      seln=sel(s,p);
      
      %交叉操作
      scro=cro(s,seln,pc);
      scnew(j,:)=scro(1,:);
      scnew(j+1,:)=scro(2,:);
      
      %變異操作
      smnew(j,:)=mut(scnew(j,:),pm);
      smnew(j+1,:)=mut(scnew(j+1,:),pm);
   end
   s=smnew;  %產(chǎn)生了新的種群
   
   %計(jì)算新種群的適應(yīng)度   
   [f,p]=objf(s);
   
   %記錄當(dāng)前代最好和平均的適應(yīng)度
   [fmax,nmax]=max(f);
   fmean=mean(f);
   ymax(gn)=fmax;
   ymean(gn)=fmean;
   %記錄當(dāng)前代的最佳個(gè)體
   x=n2to10(s(nmax,:));
   xx=-1.0+x*3/(power(2,bn)-1);
   xmax(gn)=xx;
   
   gn=gn+1
end
gn=gn-1;

%繪制曲線
subplot(2,1,1);
plot(1:gn,[ymax;ymean]);
title('歷代適應(yīng)度變化','fonts',10);
legend('最大適應(yīng)度','平均適應(yīng)度');
string1=['最終適應(yīng)度',num2str(ymax(gn))];
gtext(string1);
subplot(2,1,2);
plot(1:gn,xmax,'r-');
legend('自變量');
string2=['最終自變量',num2str(xmax(gn))];
gtext(string2);