遗传算法用matlab求函数极大值.docVIP

遗传算法用matlab求函数极大值 题目： 寻找f(x)=x2,,当x在0~31区间的最大值。 源程序： %遗传算法求解函数最大值 %本程序用到了英国谢菲尔德大学(Sheffield)开发的工具箱GATBX，该工具箱比matlab自带的GATOOL使用更加灵活，但在编写程序方面稍微复杂一些 Close all; Clear all; figure(1); fplot(variable*variable,[0,31]); %画出函数曲线 %以下定义遗传算法参数 GTSM=40; %定义个体数目 ZDYCDS=20; %定义最大遗传代数 EJZWS=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(x)=x2 while genZDYCDS %进行循环控制，当当前代数小于定义的最大遗传代数时，继续循环，直至代数等于最大遗传代数 FitnV=ranking(-ObjV); %分配适应度值 SelCh=select(sus, Chrom, FitnV, DG); %选择，即对个体按照他们的适配值进行复制 SelCh=recombin(xovsp, SelCh, 0.7); %交叉，即首先将复制产生的匹配池中的成员随机两两匹配，再进行交叉繁殖 SelCh=mut(SelCh); %变异，以一个很小的概率随机地改变一个个体串位的值 variable=bs2rv(SelCh, FieldD); %子代个体的十进制转换 ObjVSel=variable*variable; %计算子代的目标函数值 [Chrom ObjV]=reins(Chrom, SelCh, 1, 1, ObjV, ObjVSel); %再插入子代的新种群，其中Chrom为包含当前种群个体的矩阵，SelCh为包好当前种群后代的矩阵 variable=bs2rv(Chrom, FieldD); %十进制转换 gen=gen+1; %代数计数器增加 %输出最优解及其序号，并在目标函数图像中标出，Y为最优解,I为种群的 %序号 [Y, I]=max(ObjV);hold on; %求出其最大目标函数值 plot(variable(I), Y, bo); trace(1, gen)=max(ObjV); %遗传算法性能跟踪 trace(2, gen)=sum(ObjV)/length(ObjV); end %以下为结果显示部分，通过上面计算出的数值进行绘图 variable=bs2rv(Chrom, FieldD); %最优个体进行十进制转换 hold on, grid; plot(variable,ObjV,b*); %将结果画出 三、运行结果： 由图可见该函数为单调递增函数，即当X=31时，该取得最大值f(x)max=961。