matlab遺傳算法工具箱函數(shù)及實(shí)例講解_百度文庫(kù)

1、matlab遺傳算法工具箱函數(shù)及實(shí)例講解核心函數(shù)： (1function pop=initializega(num,bounds,eevalFN,eevalOps,options-初始種群的生成函數(shù) 【輸出參數(shù)】 pop-生成的初始種群 【輸入?yún)?shù)】 num-種群中的個(gè)體數(shù)目 bounds-代表變量的上下界的矩陣 eevalFN-適應(yīng)度函數(shù) eevalOps-傳遞給適應(yīng)度函數(shù)的參數(shù) options-選擇編碼形式(浮點(diǎn)編碼或是二進(jìn)制編碼precision F_or_B,如 precision-變量進(jìn)行二進(jìn)制編碼時(shí)指定的精度 F_or_B-為1時(shí)選擇浮點(diǎn)編碼，否則為二進(jìn)制編碼,由precision

2、指定精度 (2function x,endPop,bPop,traceInfo = ga(bounds,evalFN,evalOps,startPop,opts,. termFN,termOps,selectFN,selectOps,xOverFNs,xOverOps,mutFNs,mutOps-遺傳算法函數(shù) 【輸出參數(shù)】 x-求得的最優(yōu)解 endPop-最終得到的種群 bPop-最優(yōu)種群的一個(gè)搜索軌跡 【輸入?yún)?shù)】 bounds-代表變量上下界的矩陣 evalFN-適應(yīng)度函數(shù) evalOps-傳遞給適應(yīng)度函數(shù)的參數(shù) startPop-初始種群 optsepsilon prob_ops dis

3、play-opts(1:2等同于initializega的options參數(shù)，第三個(gè)參數(shù)控制是否輸出，一般為0。如1e-6 1 0 termFN-終止函數(shù)的名稱(chēng),如maxGenTerm termOps-傳遞個(gè)終止函數(shù)的參數(shù),如100 selectFN-選擇函數(shù)的名稱(chēng),如normGeomSelect selectOps-傳遞個(gè)選擇函數(shù)的參數(shù),如0.08 xOverFNs-交叉函數(shù)名稱(chēng)表，以空格分開(kāi)，如arithXover heuristicXover simpleXover xOverOps-傳遞給交叉函數(shù)的參數(shù)表，如2 0;2 3;2 0 mutFNs-變異函數(shù)表，如boundaryMutat

4、ion multiNonUnifMutation nonUnifMutation unifMutation mutOps-傳遞給交叉函數(shù)的參數(shù)表,如4 0 0;6 100 3;4 100 3;4 0 0 matlab遺傳算法工具箱附件【注意】matlab工具箱函數(shù)必須放在工作目錄下 【問(wèn)題】求f(x=x+10*sin(5x+7*cos(4x的最大值，其中0=x=9 【分析】選擇二進(jìn)制編碼，種群中的個(gè)體數(shù)目為10，二進(jìn)制編碼長(zhǎng)度為20，交叉概率為0.95,變異概率為0.08 【程序清單】 %編寫(xiě)目標(biāo)函數(shù) functionsol,eval=fitness(sol,options x=sol(1;

5、eval=x+10*sin(5*x+7*cos(4*x; %把上述函數(shù)存儲(chǔ)為fitness.m文件并放在工作目錄下 initPop=initializega(10,0 9,fitness;%生成初始種群，大小為10 x endPop,bPop,trace=ga(0 9,fitness,initPop,1e-6 1 1,maxGenTerm,25,normGeomSelect,. 0.08,arithXover,2,nonUnifMutation,2 25 3 %25次遺傳迭代 運(yùn)算結(jié)果為：x = 7.8562 24.8553(當(dāng)x為7.8562時(shí)，f（x）取最大值24.8553 注：遺傳算法一

6、般用來(lái)取得近似最優(yōu)解，而不是最優(yōu)解。另外遺傳算法的收斂性跟其初始值有關(guān)，大家運(yùn)行上面的命令所得到的借過(guò)可能跟我的借過(guò)不同或是差別很大。但多執(zhí)行幾次上面的命令（隨即取不同的初始群體）一定可以得到近似最優(yōu)解。遺傳算法實(shí)例2 【問(wèn)題】在5=Xi=5,i=1,2區(qū)間內(nèi)，求解 f(x1,x2=-20*exp(-0.2*sqrt(0.5*(x1.2+x2.2-exp(0.5*(cos(2*pi*x1+cos(2*pi*x2+22.71282的最小值。 【分析】種群大小10，最大代數(shù)1000，變異率0.1,交叉率0.3 【程序清單】 源函數(shù)的matlab代碼 function eval=f(sol numv

7、=size(sol,2; x=sol(1:numv; eval=-20*exp(-0.2*sqrt(sum(x.2/numv-exp(sum(cos(2*pi*x/numv+22.71282; %適應(yīng)度函數(shù)的matlab代碼 function sol,eval=fitness(sol,options numv=size(sol,2-1; x=sol(1:numv; eval=f(x; eval=-eval; %遺傳算法的matlab代碼 bounds=ones(2,1*-5 5; p,endPop,bestSols,trace=ga(bounds,fitness 注：前兩個(gè)文件存儲(chǔ)為m文件并放

8、在工作目錄下，運(yùn)行結(jié)果為 p = 0.0000 -0.0000 0.00551.【分析】選擇二進(jìn)制編碼，種群中的個(gè)體數(shù)目為10，二進(jìn)制編碼長(zhǎng)度為20，交叉概率為0.95,變異概率為0.08 initPop=initializega(10,0 9,fitness;%生成初始種群，大小為10x endPop,bPop,trace=ga(0 9,fitness,initPop,1e-6 1 1,maxGenTerm,25,normGeomSelect,.0.08,arithXover,2,nonUnifMutation,2 25 3 %25次遺傳迭代addoil2.【分析】種群大小10，最大代數(shù)10

9、00，變異率0.1,交叉率0.3p,endPop,bestSols,trace=ga(bounds,fitness看了一下，發(fā)現(xiàn)和在你的代碼里沒(méi)有體現(xiàn)出來(lái)？evalops是傳遞給適應(yīng)度函數(shù)的參數(shù)，opts是二進(jìn)制編碼的精度，termops是選擇maxGenTerm結(jié)束函數(shù)時(shí)傳遞個(gè)maxGenTerm 的參數(shù)，即遺傳代數(shù)。xoverops是傳遞給交叉函數(shù)的參數(shù)。mutops是傳遞給變異函數(shù)的參數(shù)，具體含義我也沒(méi)弄懂，我覺(jué)得有點(diǎn)怪。多目標(biāo)優(yōu)化的處理 目前我見(jiàn)到過(guò)提到的有這樣一些方法： 1.小生境算法； 2.轉(zhuǎn)化為單目標(biāo)算法（實(shí)質(zhì)上就是對(duì)各個(gè)目標(biāo)進(jìn)行加權(quán)求和轉(zhuǎn)換為單目標(biāo)）； 3.共享函數(shù)法等. .

10、上面他講的matlab中的工具箱,我沒(méi)有下載 估計(jì)也試gaot v5 我大致看過(guò)這個(gè)工具箱 如果你想直接應(yīng)用這個(gè)工具箱 好像是不可以的 但是,如果你只是先自己設(shè)計(jì)好一些處理方法 譬如,用小生境,或者是加權(quán)轉(zhuǎn)換 你還是可以利用其中的很多函數(shù)的, 譬如: 編碼的, 交叉的,變異的. 你可以在自己設(shè)計(jì)的程序中適當(dāng)?shù)卣{(diào)用它們問(wèn)：我們可以控制matlab的啟動(dòng)嗎？ 答： 在toolboxlocal目錄下有兩個(gè)文件，matlabrc.m，startup.m文件，它們關(guān)系matlab 的啟動(dòng)。Matlab啟動(dòng)是自動(dòng)執(zhí)行這兩個(gè)文件，只要修改這兩個(gè)文件，我們就可以控制ma tlab的啟動(dòng)了。 startup.m

11、不是安裝時(shí)就有的，原型是startupsav.m 在執(zhí)行matlab.exe還可以添加參數(shù) /automatin 啟動(dòng)作為一個(gè)自動(dòng)化服務(wù)器 /logfile logfilename 自動(dòng)把matlab輸出寫(xiě)到指定的日志文件 /minimize 最小化啟動(dòng) /nosplash 不顯示扉屏 /r M_file 自動(dòng)執(zhí)行指定的M文件 /regserver 注冊(cè)ActiveX條目 /unregserver 刪除ActiveX條目 這些參數(shù)對(duì)在vb中用shell命令啟動(dòng)matlab時(shí)是很有用的。 在執(zhí)行matlab.exe還可以添加參數(shù) /automatin 啟動(dòng)作為一個(gè)自動(dòng)化服務(wù)器 /logfile

12、logfilename 自動(dòng)把matlab輸出寫(xiě)到指定的日志文件 /minimize 最小化啟動(dòng) /nosplash 不顯示扉屏 /r M_file 自動(dòng)執(zhí)行指定的M文件 /regserver 注冊(cè)ActiveX條目 /unregserver 刪除ActiveX條目 這些參數(shù)對(duì)在vb中用shell命令啟動(dòng)matlab時(shí)是很有用的。 另： 在matlab退出時(shí)，會(huì)自動(dòng)執(zhí)行toolboxlocal目錄下finish.m文件。addoil wrote: 【注意】matlab工具箱函數(shù)必須放在工作目錄下 【問(wèn)題】求f(x=x+10*sin(5x+7*cos(4x的最大值，其中0=x In C:MATL

13、AB6p5toolboxmatlabgeneralpath.m at line 116 In C:MATLAB6p5toolboxgastartup.m at line 1 In C:MATLAB6p5toolboxlocalmatlabrc.m at line 19 不知是什么原因，請(qǐng)問(wèn)如何解決，盼回復(fù)！ 哥們： 這個(gè)問(wèn)題我也遇到過(guò)，解決辦法很簡(jiǎn)單。把遺傳算法工具箱中的startup.m的語(yǔ)句屏蔽掉就可以了。如果問(wèn)題解決了，版主給加分。引用 使用道具 報(bào)告 回復(fù) TOP wzchu1979 該用戶(hù)已被刪除 42# 大 中 小 發(fā)表于 2006-8-12 10:45 只看該作者 如何如何控制交

14、叉概率和變異概率?各位大俠: 小弟請(qǐng)教三個(gè)問(wèn)題. matlab遺傳算法工具箱函數(shù)是如何如何控制交叉概率和變異概率? xOverFNs交叉函數(shù)為什么會(huì)有多個(gè),例如缺省的是arithXover heuristicXover simpleXover ? xOverOps參數(shù)是如何起作用的?引用 使用道具 報(bào)告 回復(fù) TOP yxz32 新手會(huì)員帖子2精華0積分0貢獻(xiàn)積分0 閱讀權(quán)限10 個(gè)人空間 發(fā)短消息 加為好友 當(dāng)前離線 43# 大 中 小 發(fā)表于 2006-8-28 19:54 只看該作者 我是一個(gè)新手，看了帖子上的程序，可是無(wú)法通過(guò)?。扛呤謧兘o解答一下吧！%源函數(shù)的matlab代碼 func

15、tion eval=f(sol numv=size(sol,2; x=sol(1:numv; eval=-20*exp(-0.2*sqrt(sum(x.2/numv-exp(sum(cos(2*pi*x/numv+22.71282; %適應(yīng)度函數(shù)的matlab代碼 function sol,eval=fitness(sol,options numv=size(sol,2-1; x=sol(1:numv; eval=f(x; eval=-eval; %遺傳算法的matlab代碼 bounds=ones(2,1*-5 5; p,endPop,bestSols,trace=ga(bounds,fit

16、ness ? Error: File: E:matlabworkzz4.m Line: 8 Column: 49Missing operator, comma, or semicolon.UID185946精華0在線時(shí)間3 小時(shí)注冊(cè)時(shí)間2006-8-27最后登錄2007-1-22查看詳細(xì)資料引用 使用道具 報(bào)告 回復(fù) TOP kying 該用戶(hù)已被刪除 44# 大 中 小 發(fā)表于 2006-10-14 21:04 只看該作者 我也正需要這方面的資料，實(shí)在太好了，謝謝了。引用 使用道具 報(bào)告 回復(fù) TOP beveron88 新手會(huì)員帖子32精華0積分0貢獻(xiàn)積分12 閱讀權(quán)限10 個(gè)人空間 發(fā)短

17、消息 加為好友 當(dāng)前離線 45# 大 中 小 發(fā)表于 2007-7-25 15:54 只看該作者 謝謝樓主熱心分享本帖最近評(píng)分記錄 bainhome 貢獻(xiàn)積分 -5 灌水 2007-7-25 17:39 UID183531精華0在線時(shí)間47 小時(shí)注冊(cè)時(shí)間2006-8-16最后登錄2008-5-15查看詳細(xì)資料引用 使用道具 報(bào)告 回復(fù) TOP chrismon 新手會(huì)員帖子12精華0積分0貢獻(xiàn)積分-82 閱讀權(quán)限10 個(gè)人空間 發(fā)短消息 加為好友 當(dāng)前離線 46# 大 中 小 發(fā)表于 2007-7-26 21:23 只看該作者 好 謝謝 學(xué)習(xí)本帖最近評(píng)分記錄 dingzixiu 貢獻(xiàn)積分 -2

18、 惡意灌水 2007-7-26 22:07 UID70706精華0在線時(shí)間4 小時(shí)注冊(cè)時(shí)間2004-9-21最后登錄2007-12-7查看詳細(xì)資料引用 使用道具 報(bào)告 回復(fù) TOP pyyyzl007 新手會(huì)員帖子8精華0積分0貢獻(xiàn)積分3 閱讀權(quán)限10 個(gè)人空間 發(fā)短消息 加為好友 當(dāng)前離線 47# 大 中 小 發(fā)表于 2007-7-27 21:34 只看該作者 function pop = initializega(num, bounds, evalFN,evalOps,optionsif nargin5options=1e-6 1;endif nargin4evalOps=;endif any(evalFN48 %Not a .m fileif options(2=1 %Float GA estr=x=pop(i,1; pop(i,xZomeLength=, evalFN ;else %Binary GA estr=x=b2f(pop(i,:,bounds,bits; pop(i,xZomeLength=, evalFN ; endelse %A .m fileif options(2=1 %Float GA estr= pop(i,: pop(i,xZomeLength= eval