遗传算法实验六完整版.ppt

精选文档 程序运行所求得的最优解为：y=1.0320 精选文档 §5.3多元多峰函数优化实例 精选文档 * Shubert函数的图像为 精选文档 * 目标函数 function z=shubert(x,y) z=((1*cos((1+1)*x+1))+(2*cos((2+1)*x+2))+… (3*cos((3+1)*x+3))+(4*cos((4+1)*x+4))+… (5*cos((5+1)*x+5))).*((1*cos((1+1)*y+1))+… (2*cos((2+1)*y+2))+(3*cos((3+1)*y+3))+… (4*cos((4+1)*y+4))+(5*cos((5+1)*y+5))); 精选文档 * 1、将绘图区域划分为矩形网格—MESHGRID 功能： 将向量x,y指定的区域转化为矩形X，Y 调用格式： [X,Y] = MESHGRID(x,y) 2、 3维图形输出函数 SURF(X,Y,Z)—着色表面图 MESH(X,Y,Z)—网线图 例如： [X,Y] = meshgrid(-2:.2:2, -2:.2:2); Z = X .* exp(-X.^2 - Y.^2); surf(X,Y,Z); MESH(X,Y,Z); 程序中的函数说明 精选文档 精选文档 精选文档 程序代码如下： [x1,x2]=meshgrid(-10:.1:10); %画出Shubert函数图像figure(1);mesh(x1,x2,shubert(x1,x2)); %定义遗传算法参数 NIND=40; %个体数目 MAXGEN=50; %最大遗传代数 NVAR=2; %变量数目 PRECI=25; %变量的二进制位数 GGAP=0.9; %代沟 精选文档 %建立区域描述器 FieldD=[rep([PRECI],[1,NVAR]);… rep([-10;10],[1,NVAR]);rep([1;0;1;1],[1,NVAR])]; %创建初始种群 Chrom=crtbp(NIND, NVAR*PRECI); gen=0; %遗传算法性能跟踪初始值 trace=zeros(MAXGEN, 2); %初始种群十进制转换 x=bs2rv(Chrom, FieldD); %初始种群的目标函数值 ObjV=Shubert(x(:,1),x(:,2)); 精选文档 while genMAXGEN FitnV=ranking(ObjV); %分配适应度值 SelCh=select(sus,Chrom,FitnV,GGAP); SelCh=recombin(xovsp,SelCh,0.7); SelCh=mut(SelCh); x=bs2rv(SelCh,FieldD); ObjVSel=Shubert(x(:,1),x(:,2)); [Chrom ObjV]=reins(Chrom,SelCh,1,1,… ObjV,ObjVSel); gen=gen+1; 精选文档 [Y, I]=min(ObjV); %输出每一次的最优解及其对应的自变量值 Y,bs2rv(Chrom(I,:),FieldD) %遗传算法性能跟踪 trace(gen,1)=min(ObjV); trace(gen,2)=sum(ObjV)/length(ObjV); %迭代数为50时画出目标函数值分布图 if(gen==50) figure(2); plot(ObjV);hold on; plot(ObjV,b*);grid; end end 精选文档 figure(3);clf; plot(trace(:,1));hold on; plot(trace(:,2),-.);grid legend(解的变化,种群均值的变化) 精选文档 精选文档 精选文档 §5.4 在多目标优化中的应用 精选文档 5.4.1多目标优化的概念 精选文档 figure(1); %画出函数曲线 fplot(variable.*sin(10