利用遗传算法求函数的极大值.doc

利用遗传算法求函数的极大值 该函数有两个局部极大值点，分别是f(2.048,-2.048)=3897.7342和f(2.048,-2.048)=3905.9262，其中，后者为全局最大点。可以分别用二进制编码和十进制编码遗传算法求函数极大值 遗传算法二进制编码求函数极大值程序 %Generic Algorithm for function f(x1,x2) optimum clear all; close all; %Parameters参数 Size=80; %群体大小 G=100; %终止进化代数 CodeL=10; %代码长度 umax=2.048; umin=-2.048; E=round(rand(Size,2*CodeL)); %Initial Code 最初代码 %Main Program主程序 for k=1:1:G time(k)=k; for s=1:1:Size m=E(s,:); y1=0;y2=0; %X对应的十进制代码 %Uncoding m1=m(1:1:CodeL); for i=1:1:CodeL y1=y1+m1(i)*2^(i-1); %将y1转换为十进制数 end x1=(umax-umin)*y1/1023+umin; m2=m(CodeL+1:1:2*CodeL); for i=1:1:CodeL y2=y2+m2(i)*2^(i-1); %将y2转换为十进制数 end x2=(umax-umin)*y2/1023+umin; %求x对应的十进制数 F(s)=100*(x1^2-x2)^2+(1-x1)^2; %个体适应度函数 end Ji=1./F; %个体适应度函数的倒数 %****** Step 1 : Evaluate BestJ ****** BestJ(k)=min(Ji); fi=F; %Fitness Function 适应函数 [Oderfi,Indexfi]=sort(fi); %Arranging fi small to bigger Bestfi=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 kk=1; for i=1:1:Size for j=1:1:fi_S(i) %Select and Reproduce TempE(kk,:)=E(Indexfi(i),:); kk=kk+1; %kk is used to reproduce end end %************ Step 3 : Crossover Operation ************ pc=0.60; %交叉概率 n=ceil(20*rand); %种群大小 fo