




已阅读5页,还剩5页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 产品持续维护承诺函9篇范文
- 2025年制造业高端装备研发能力考核试卷-机器人与自动化医疗手术机器人避障算法
- 难点解析人教版八年级物理上册第5章透镜及其应用综合练习练习题(详解)
- 2.2025年互联网与信息技术继续教育云原生架构设计与优化之云原生环境下深度学习模型部署优化考核试卷
- 2025年房地产项目土地闲置处置法律实务考核试卷
- 考点解析人教版八年级上册物理光现象《平面镜成像》综合测试试题(含解析)
- 在动手体验的过程中领悟概念的本质
- 开放题在小学数学“学讲”课堂的教学实践
- 考点解析-人教版八年级物理上册第5章透镜及其应用-透镜专题练习试题(含详细解析)
- 解析卷人教版八年级上册物理光现象《光的反射》专题攻克试卷(含答案解析)
- 老年人安全知识宣传方案
- 沉香合作协议合同范例
- 象棋入门课件教学
- 街道办消防安全知识培训课件
- 学生外出安全协议书模板学生外出安全协议书
- 初中班会 10月31日世界勤俭日宣传介绍主题班会 课件
- 锅炉使用单位安全员题库
- 五年级上册体育与健康教案
- 医院培训课件:《中医拔罐类技术相关性感染预防与控制指南》
- FZT 73012-2017 文胸行业标准
- 新概念第二册课文和单词
评论
0/150
提交评论