matlab遗传算法工具箱函数及实例讲解.doc

matlab遗传算法工具箱函数及实例讲解核心函数： (1)function pop=initializega(num,bounds,eevalFN,eevalOps,options)-初始种群的生成函数 【输出参数】 pop-生成的初始种群 【输入参数】 num-种群中的个体数目 bounds-代表变量的上下界的矩阵 eevalFN-适应度函数 eevalOps-传递给适应度函数的参数 options-选择编码形式(浮点编码或是二进制编码)precision F_or_B,如 precision-变量进行二进制编码时指定的精度 F_or_B-为1时选择浮点编码，否则为二进制编码,由precision指定精度) (2)function x,endPop,bPop,traceInfo = ga(bounds,evalFN,evalOps,startPop,opts,. termFN,termOps,selectFN,selectOps,xOverFNs,xOverOps,mutFNs,mutOps)-遗传算法函数 【输出参数】 x-求得的最优解 endPop-最终得到的种群 bPop-最优种群的一个搜索轨迹 【输入参数】 bounds-代表变量上下界的矩阵 evalFN-适应度函数 evalOps-传递给适应度函数的参数 startPop-初始种群 optsepsilon prob_ops display-opts(1:2)等同于initializega的options参数，第三个参数控制是否输出，一般为0。如1e-6 1 0 termFN-终止函数的名称,如maxGenTerm termOps-传递个终止函数的参数,如100 selectFN-选择函数的名称,如normGeomSelect selectOps-传递个选择函数的参数,如0.08 xOverFNs-交叉函数名称表，以空格分开，如arithXover heuristicXover simpleXover xOverOps-传递给交叉函数的参数表，如2 0;2 3;2 0 mutFNs-变异函数表，如boundaryMutation multiNonUnifMutation nonUnifMutation unifMutation mutOps-传递给交叉函数的参数表,如4 0 0;6 100 3;4 100 3;4 0 0 matlab遗传算法工具箱附件【注意】matlab工具箱函数必须放在工作目录下 【问题】求f(x)=x+10*sin(5x)+7*cos(4x)的最大值，其中0<=x<=9 【分析】选择二进制编码，种群中的个体数目为10，二进制编码长度为20，交叉概率为0.95,变异概率为0.08 【程序清单】 %编写目标函数 functionsol,eval=fitness(sol,options) x=sol(1); eval=x+10*sin(5*x)+7*cos(4*x); %把上述函数存储为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次遗传迭代 运算结果为：x = 7.8562 24.8553(当x为7.8562时，f（x）取最大值24.8553) 注：遗传算法一般用来取得近似最优解，而不是最优解。另外遗传算法的收敛性跟其初始值有关，大家运行上面的命令所得到的借过可能跟我的借过不同或是差别很大。但多执行几次上面的命令（随即取不同的初始群体）一定可以得到近似最优解。遗传算法实例2 【问题】在5<=Xi<=5,i=1,2区间内，求解 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，最大代数1000，变异率0.1,交叉率0.3 【程序清单】 源函数的matlab代码 function 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; %适应度函数的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) 注：前两个文件存储为m文件并放在工作目录下，运行结果为 p = 0.0000 -0.0000 0.00551.【分析】选择二进制编码，种群中的个体数目为10，二进制编码长度为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，最大代数1000，变异率0.1,交叉率0.3p,endPop,bestSols,trace=ga(bounds,fitness)看了一下，发现和在你的代码里没有体现出来？问得好啊，snow_man_0。这些东西都是我从树上抄下来的，书上没有讲清楚，而且我也觉得有些参数有问题。 至于第一个例子，选择二进制编码和种群中的个体数为10是没有问题的，二进制的编码长度为20是通过精度（1e-6） 来确定的。我也觉得交叉概率0.95和变异概率0.08没有在程序中体现出来，程序只所以能够正确运行，我想可能是采用了默认的交叉概率0.6以及变异概 率0.05的原因吧。程序中的参数arithXover,2,nonUnifMutation,2 25 3) 我也觉得奇怪，大家一起看看源代码，好好研究研究啊。evalops是传递给适应度函数的参数，opts是二进制编码的精度，termops是选择maxGenTerm结束函数时传递个maxGenTerm 的参数，即遗传代数。xoverops是传递给交叉函数的参数。mutops是传递给变异函数的参数，具体含义我也没弄懂，我觉得有点怪。多目标优化的处理 目前我见到过提到的有这样一些方法： 1.小生境算法； 2.转化为单目标算法（实质上就是对各个目标进行加权求和转换为单目标）； 3.共享函数法等. . 上面他讲的matlab中的工具箱,我没有下载 估计也试gaot v5 我大致看过这个工具箱 如果你想直接应用这个工具箱 好像是不可以的 但是,如果你只是先自己设计好一些处理方法 譬如,用小生境,或者是加权转换 你还是可以利用其中的很多函数的, 譬如: 编码的, 交叉的,变异的. 你可以在自己设计的程序中适当地调用它们问：我们可以控制matlab的启动吗？ 答： 在toolboxlocal目录下有两个文件，matlabrc.m，startup.m文件，它们关系matlab 的启动。Matlab启动是自动执行这两个文件，只要修改这两个文件，我们就可以控制ma tlab的启动了。 startup.m不是安装时就有的，原型是startupsav.m 在执行matlab.exe还可以添加参数 /automatin 启动作为一个自动化服务器 /logfile logfilename 自动把matlab输出写到指定的日志文件 /minimize 最小化启动 /nosplash 不显示扉屏 /r M_file 自动执行指定的M文件 /regserver 注册ActiveX条目 /unregserver 删除ActiveX条目 这些参数对在vb中用shell命令启动matlab时是很有用的。 在执行matlab.exe还可以添加参数 /automatin 启动作为一个自动化服务器 /logfile logfilename 自动把matlab输出写到指定的日志文件 /minimize 最小化启动 /nosplash 不显示扉屏 /r M_file 自动执行指定的M文件 /regserver 注册ActiveX条目 /unregserver 删除ActiveX条目 这些参数对在vb中用shell命令启动matlab时是很有用的。 另： 在matlab退出时，会自动执行toolboxlocal目录下finish.m文件。addoil wrote: 【注意】matlab工具箱函数必须放在工作目录下 【问题】求f(x)=x+10*sin(5x)+7*cos(4x)的最大值，其中0<=x<=9 【分析】选择二进制编码，种群中的个体数目为10，二进制编码长度为20，交叉概率为0.95,变异概率为0.08 【程序清单】 %编写目标函数 functionsol,eval=fitness(sol,options) x=sol(1); eval=x+10*sin(5*x)+7*cos(4*x); %把上述函数存储为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次遗传迭代 运算借过为：x = 7.8562 24.8553(当x为7.8562时，f（x）取最大值24.8553) 注：遗传算法一般用来取得近似最优解，而不是最优解。另外遗传算法的收敛性跟其初始值有关，大家运行上面的命令所得到的借过可能跟我的借过不同或是差别很大。但多执行几次上面的命令（随即取不同的初始群体）一定可以得到近似最优解。请教addoil： （1）这里面编写目标函数那3句为什么要另外存成一个文件？这3句与下面的程序存成一个文件不可以吗？ （2）目标函数中参数sol的含义是什么？为什么在函数定义行“”的两边都有这个sol参数，如何理解？ 谢谢！sol是当前个体，若你的求解空间是n维，则sol是n+1个元素的行向量。一般通过sol来编目标函数,sol(i)就是第i个自变量。第n1个元素好像就是优化的目标吧。addoil 大侠，我下载了你所提供的遗传算法工具箱，能够使用了，但是在command 框中总是出现下面的出错提示： Using Toolbox Path Cache. Type help toolbox_path_cache for more info. To get started, select MATLAB Help from the Help menu. Warning: Name is nonexistent or not a directory: c:MATLAB6p5toolboxgagaotv5. > In C:MATLAB6p5toolboxmatlabgeneralpath.m at line 116 In C:MATLAB6p5toolboxgastartup.m at line 1 In C:MATLAB6p5toolboxlocalmatlabrc.m at line 19 不知是什么原因，请问如何解决，盼回复！ 哥们： 这个问题我也遇到过，解决办法很简单。把遗传算法工具箱中的startup.m的语句屏蔽掉就可以了。如果问题解决了，版主给加分。引用 使用道具 报告 回复 TOP wzchu1979 该用户已被删除 42# 大 中 小 发表于 2006-8-12 10:45 只看该作者 如何如何控制交叉概率和变异概率?各位大侠: 小弟请教三个问题. matlab遗传算法工具箱函数是如何如何控制交叉概率和变异概率? xOverFNs交叉函数为什么会有多个,例如缺省的是arithXover heuristicXover simpleXover ? xOverOps参数是如何起作用的?引用 使用道具 报告 回复 TOP yxz32 新手会员帖子2精华0积分0贡献积分0 阅读权限10 个人空间 发短消息 加为好友 当前离线 43# 大 中 小 发表于 2006-8-28 19:54 只看该作者 我是一个新手，看了帖子上的程序，可是无法通过啊？高手们给解答一下吧！%源函数的matlab代码 function 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; %适应度函数的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) ? Error: File: E:matlabworkzz4.m Line: 8 Column: 49Missing operator, comma, or semicolon.UID185946精华0在线时间3 小时注册时间2006-8-27最后登录2007-1-22查看详细资料引用 使用道具 报告 回复 TOP kying 该用户已被删除 44# 大 中 小 发表于 2006-10-14 21:04 只看该作者 我也正需要这方面的资料，实在太好了，谢谢了。引用 使用道具 报告 回复 TOP beveron88 新手会员帖子32精华0积分0贡献积分12 阅读权限10 个人空间 发短消息 加为好友 当前离线 45# 大 中 小 发表于 2007-7-25 15:54 只看该作者 谢谢楼主热心分享本帖最近评分记录 bainhome 贡献积分 -5 灌水 2007-7-25 17:39 UID183531精华0在线时间47 小时注册时间2006-8-16最后登录2008-5-15查看详细资料引用 使用道具 报告 回复 TOP chrismon 新手会员帖子12精华0积分0贡献积分-82 阅读权限10 个人空间 发短消息 加为好友 当前离线 46# 大 中 小 发表于 2007-7-26 21:23 只看该作者 好 谢谢 学习本帖最近评分记录 dingzixiu 贡献积分 -2 恶意灌水 2007-7-26 22:07 UID70706精华0在线时间4 小时注册时间2004-9-21最后登录2007-12-7查看详细资料引用 使用道具 报告 回复 TOP pyyyzl007 新手会员帖子8精华0积分0贡献积分3 阅读权限10 个人空间 发短消息 加为好友 当前离线 47# 大 中 小 发表于 2007-7-27 21:34 只看该作者 function pop = initializega(num, bounds, evalFN,evalOps,options)if 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); p