遗传算法Matlab实例讲解PPT课件.pptx

遗传算法（GA）实例讲解求函数f(x)=10*sin(5x)+7*cos(4x)，x∈[0,10]的最大值：1.

函数最大值f(x)=10*sin(5x)+7*cos(4x)的最大值f1(x)=-f(x)-10*sin(5x)-7*cos(4x)的最小值[x,fval]=fminbnd(‘-10*sin(5*x)-7*cos(4*x)’,0,10)x=1.5708fval=-17.0000x=7.8540fval=-17.0000f(x)=10*sin(5x)+7*cos(4x)，x∈[0,10]的最大值为：x1=1.5708x2=7.8540F(x1)=f(x2)=17.00002.

主函数Matlab仿真主程序将x的值用一个10位的二值形式表示为二值问题（编码），一个10位的二值数提供的分辨率是(10-0)/(2^10-1)≈0.01。3.

群体初始化

将x的值用一个10位的二值形式表示为二值问题（编码），一个10位的二值数提供的分辨率是每为(10-0)/(2^10-1)≈0.01。4.

计算目标函数值function[objvalue]=calobjvalue(pop)temp1=decodechrom(pop,1,10);x=temp1*10/1023;objvalue=10*sin(5*x)+7*cos(4*x);end5.

计算个体适应值functionfitvalue=calfitvalue(objvalue)globalCmin;Cmin=0;[px,py]=size(objvalue);fori=1:pxifobjvalue(i)+Cmin0temp=Cmin+objvalue(i);elsetemp=0.0;endfitvalue(i)=temp;endfitvalue=fitvalue;end6.

选择function[newpop]=selection(pop,fitvalue)totalfit=sum(fitvalue);fitvalue=fitvalue/totalfit;fitvalue=cumsum(fitvalue);[px,py]=size(pop);ms=sort(rand(px,1));fitin=1;newin=1;whilenewin=pxif(ms(newin))fitvalue(fitin)newpop(newin)=pop(fitin);newin=newin+1;elsefitin=fitin+1;endendend7.

交叉function[newpop]=crossover(pop,pc)[px,py]=size(pop);newpop=ones(size(pop));fori=1:2:px-1if(randpc)cpoint=round(rand*py);newpop(i,:)=[pop(i,1:cpoint),pop(i+1,cpoint+1:py)];newpop(i+1,:)=[pop(i+1,1:cpoint),pop(i,cpoint+1:py)];elsenewpop(i,:)=pop(i);newpop(i+1,:)=pop(i+1);endendend8.

变异function[newpop]=mutation(pop,pm)[px,py]=size(pop);fori=1:pxif(randpm)mpoint=round(rand*py);ifmpoint=0mpoint=1;endnewpop(i)=pop(i);ifany(newpop(i,mpoint)==0)newpop(i,mpoint)=1;elsenewpop(i,mpoint)=0;endelsenewpop(i)=pop(i);endendend9.

计算最大适应值function[bestindividual,bestfit]=best(pop,fitvalue)[px,py]=size(pop);bestindividual=pop(1,:);bestfit=fitvalue