欢迎来到天天文库
浏览记录
ID:13007529
大小:65.50 KB
页数:4页
时间:2018-07-20
《遗传算法求函数极大值(matlab实现)》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、遗传算法求函数最大值(matlab实现)一、题目:寻找f(x)=x2,,当x在0~31区间的最大值。二、源程序:%遗传算法求解函数最大值%本程序用到了英国谢菲尔德大学(Sheffield)开发的工具箱GATBX,该工具箱比matlab自带的GATOOL使用更加灵活,但在编写程序方面稍微复杂一些Closeall;Clearall;figure(1);fplot('variable*variable',[0,31]);%画出函数曲线%以下定义遗传算法参数GTSM=40;%定义个体数目ZDYCDS=20;%定义最大遗传代数EJZWS=
2、5;%定义变量的二进制位数DG=0.9;%定义代沟trace=zeros(2,ZDYCDS);%最优结果的初始值FieldD=[5;-1;2;1;0;1;1];%定义区域描述器的各个参数%以下为遗传算法基本操作部分,包括创建初始种群、复制、交叉和变异Chrom=crtbp(GTSM,EJZWS);%创建初始种群,即生成给定规模的二进制种群和结构gen=0;%定义代数计数器初始值variable=bs2rv(Chrom,FieldD);%对生成的初始种群进行十进制转换ObjV=variable*variable;%计算目标函数值f
3、(x)=x2whilegen4、s2rv(SelCh,FieldD);%子代个体的十进制转换ObjVSel=variable*variable;%计算子代的目标函数值[ChromObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);%再插入子代的新种群,其中Chrom为包含当前种群个体的矩阵,SelCh为包好当前种群后代的矩阵variable=bs2rv(Chrom,FieldD);%十进制转换gen=gen+1;%代数计数器增加%输出最优解及其序号,并在目标函数图像中标出,Y为最优解,I为种群的%序号[Y,I]=max(ObjV5、);holdon;%求出其最大目标函数值plot(variable(I),Y,'bo');trace(1,gen)=max(ObjV);%遗传算法性能跟踪trace(2,gen)=sum(ObjV)/length(ObjV);end%以下为结果显示部分,通过上面计算出的数值进行绘图variable=bs2rv(Chrom,FieldD);%最优个体进行十进制转换holdon,grid;plot(variable,ObjV,'b*');%将结果画出三、运行结果:由图可见该函数为单调递增函数,即当X=31时,该取得最大值f(x)ma6、x=961。
4、s2rv(SelCh,FieldD);%子代个体的十进制转换ObjVSel=variable*variable;%计算子代的目标函数值[ChromObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);%再插入子代的新种群,其中Chrom为包含当前种群个体的矩阵,SelCh为包好当前种群后代的矩阵variable=bs2rv(Chrom,FieldD);%十进制转换gen=gen+1;%代数计数器增加%输出最优解及其序号,并在目标函数图像中标出,Y为最优解,I为种群的%序号[Y,I]=max(ObjV
5、);holdon;%求出其最大目标函数值plot(variable(I),Y,'bo');trace(1,gen)=max(ObjV);%遗传算法性能跟踪trace(2,gen)=sum(ObjV)/length(ObjV);end%以下为结果显示部分,通过上面计算出的数值进行绘图variable=bs2rv(Chrom,FieldD);%最优个体进行十进制转换holdon,grid;plot(variable,ObjV,'b*');%将结果画出三、运行结果:由图可见该函数为单调递增函数,即当X=31时,该取得最大值f(x)ma
6、x=961。
此文档下载收益归作者所有