《遗传算法求解函数最大值(matlab》由会员分享,可在线阅读,更多相关《遗传算法求解函数最大值(matlab(5页珍藏版)》请在人人文库网上搜索。
1、遗传算法求解函数F(x1,x2)=100*(x12-x2)2+(1-x1)2;的最大值(MATLAB)%Generic Algorithm for function f(x1,x2) optimum (最大值)clear all;close all;%ParametersSize=80; G=100; CodeL=10;umax=2.048;umin=-2.048;E=round(rand(Size,2*CodeL); %Initial Code%Main Programfor k=1:1:Gtime(k)=k;for s=1:1:Sizem=E(s,:);y1=0;y2=0;%Uncodin。
2、gm1=m(1:1:CodeL);for i=1:1:CodeLy1=y1+m1(i)*2(i-1);endx1=(umax-umin)*y1/1023+umin;m2=m(CodeL+1:1:2*CodeL);for i=1:1:CodeLy2=y2+m2(i)*2(i-1);endx2=(umax-umin)*y2/1023+umin;F(s)=100*(x12-x2)2+(1-x1)2;endJi=1./(F+1);%* Step 1 : Evaluate BestJ *BestJ(k)=min(Ji);fi=F; %Fitness FunctionOderfi,Indexfi=sort。
3、(fi); %Arranging fi small to biggerBestfi=Oderfi(Size); %Let Bestfi=max(fi)BestS=E(Indexfi(Size),:); %Let BestS=E(m), m is the Indexfi belong to max(fi)bfi(k)=Bestfi;%* Step 2 : Select and Reproduct Operation*fi_sum=sum(fi);fi_Size=(Oderfi/fi_sum)*Size;fi_S=floor(fi_Size); %Selecting Bigger fi value。
4、 (取整)kk=1;for i=1:1:Sizefor j=1:1:fi_S(i) %Select and Reproduce TempE(kk,:)=E(Indexfi(i),:); kk=kk+1; %kk is used to reproduceendend%* Step 3 : Crossover Operation *pc=0.60;n=ceil(20*rand);for i=1:2:(Size-1)temp=rand;if pctemp %Crossover Conditionfor j=n:1:20TempE(i,j)=E(i+1,j);TempE(i+1,j)=E(i,j);e。
5、ndendendTempE(Size,:)=BestS;E=TempE;%* Step 4: Mutation Operation *%pm=0.001;%pm=0.001-1:1:Size*(0.001)/Size; %Bigger fi, smaller Pm%pm=0.0; %No mutationpm=0.1; %Big mutationfor i=1:1:Sizefor j=1:1:2*CodeLtemp=rand;if pmtemp %Mutation Conditionif TempE(i,j)=0TempE(i,j)=1;elseTempE(i,j)=0;endendenden。
6、d%Guarantee TempPop(30,:) is the code belong to the best individual(max(fi)TempE(Size,:)=BestS;E=TempE;endMax_Value=BestfiBestSx1x2figure(1);plot(time,BestJ); xlabel(Times);ylabel(Best J);figure(2);plot(time,bfi);xlabel(times);ylabel(Best F);轐行绒果Map_Falue =3.9059e+ 03BestS =0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 01 = -2&0480x2 =-2.0080。