【正文】
for i=1:D1 eval=eval+(sol(i)^210*cos(2*pi*sol(i))+10)。)%適應(yīng)度函數(shù)%function eval=fitness(x,D)sol=x。)TmaxyPbmaxPbminPbavrdisp(39。 endendtocdisp(39。for z=1:U if Pbmaxy(z) Pbmax=y(z)。%U次中的最差值和最好值Pbmax=y(1)。%計(jì)算平均適應(yīng)最優(yōu)值sum=Pb+sum。 end end t=t+1。 end if cost(i)=Pb Pb=cost(i)。 if(score=cost(i)) x(i,1:D)=trial(1:D)。 end if trial(k)b1(k) trial(k)=b1(k)。 else trial(k)=x(i,k)。 end %對(duì)每個(gè)個(gè)體進(jìn)行交叉操作,得試驗(yàn)個(gè)體 jrand=floor(rand*D+1)。c~=b break。 if c~=iamp。b~=a break