遺傳算法matlab代碼

1、functionyouhuafunD=codeN=50%Tunablemaxgen=50%Tunablecrossrate=0.5%Tunablemuterate=0.08%Tunablegeneration=1num=length(D)fatherr=rint(numN3)sce=zeros(maxgenN)whilegeneration=maxgenind=rperm(N2)2%隨機(jī)配對(duì)交叉A=fatherr(:ind(1:(N2)

2、2))B=fatherr(:ind((N2)21:end))%多點(diǎn)交叉rnd=r(num(N2)2)ind=rndtmp=A(ind)A(ind)=B(ind)B(ind)=tmp%%兩點(diǎn)交叉%fkk=1:(N2)2%rndtmp=rint(11num)1%tmp=A(1:rndtmpkk)%A(1:rndtmpkk)=B(1:rndtmpkk)%B(1:rndtmpkk)=tmp%endfatherr=[fatherr(:1:2)AB

3、]%變異rnd=r(numN)ind=rnd[mn]=size(ind)tmp=rint(mn2)1tmp(:1:2)=0fatherr=tmpfatherrfatherr=mod(fatherr3)%fatherr(ind)=tmp%評(píng)價(jià)、選擇sceN=scefun(fatherrD)%求得N個(gè)個(gè)體的評(píng)價(jià)函數(shù)sce(generation:)=sceN[sceStsceind]=st(sceN)sumsce=cumsum(sceSt)s

4、umsce=sumsce.sumsce(end)childind(1:2)=sceind(end1:end)fk=3:Ntmprnd=rF0index=find(FF4k==ite)if~isempty(F0index)tmpMat=F3(F0index)tmpSco=sum(tmpMat)SceBin(ite)=mod(tmpSco300)endendScek(k)=sum(SceBin)endSceN=SceNScek遺傳算法實(shí)例:

5、%下面舉例說(shuō)明遺傳算法%%求下列函數(shù)的最大值%%f(x)=10sin(5x)7cos(4x)x∈[010]%%將x的值用一個(gè)10位的二值形式表示為二值問(wèn)題，一個(gè)10位的二值數(shù)提供的分辨率是每為(100)(2^101)≈0.01。%%將變量域[010]離散化為二值域[01023]x=010b1023其中b是[01023]中的一個(gè)二值數(shù)。%%%%%%%%編程%%2.1初始化(編碼)%initpop.m函數(shù)的功能是實(shí)現(xiàn)群體的初始化，popsi

6、ze表示群體的大小，chromlength表示染色體的長(zhǎng)度(二值數(shù)的長(zhǎng)度)，%長(zhǎng)度大小取決于變量的二進(jìn)制編碼的長(zhǎng)度(在本例中取10位)。%遺傳算法子程序%Name:initpop.m%初始化functionpop=initpop(popsizechromlength)pop=round(r(popsizechromlength))%r隨機(jī)產(chǎn)生每個(gè)單元為01行數(shù)為popsize，列數(shù)為chromlength的矩陣，%roud對(duì)矩陣的每個(gè)單