遗传算法实验六演示文稿

遗传算法实验六演示文稿第一页，共七十三页。遗传算法实验六第二页，共七十三页。§5.1简单一元函数优化实例第三页，共七十三页。程序简化代码如下：NIND=40;%定义个体数目MAXGEN=25;%定义最大遗传代数PRECI=20;%编码长度GGAP=0.9;%代沟FieldD=[20;-1;2;1;0;1;1];%区域描述器Chrom=crtbp(NIND,PRECI);%产生初始种群gen=0;%代计数器variable=bs2rv(Chrom,FieldD);%初始种群的十进制转换ObjV=variable.*sin(10*pi*variable)+2.0;%计算初始种群目标函数值whilegen<MAXGEN+1FitnV=ranking(-ObjV);%分配适应度值SelCh=select('sus',Chrom,FitnV,GGAP);%选择SelCh=recombin('xovsp',SelCh,0.7);%重组SelCh=mut(SelCh); %变异variable=bs2rv(SelCh,FieldD);%子代个体的十进制转换ObjVSel=variable.*sin(10*pi*variable)+2.0;%计算子代的目标函数值

[ChromObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);%重插入子代到种群Y=max(ObjV)%当前种群最优解gen=gen+1;%代计数器增加end第四页，共七十三页。figure(1);%画出函数曲线fplot('variable.*sin(10*pi*variable)+2.0',[-1,2]);%定义遗传算法参数NIND=40;%群体中个体数目MAXGEN=25;%最大遗传代数PRECI=20;%变量的二进制位数具有性能跟踪和图像输出功能的程序代码第五页，共七十三页。%代沟GGAP=0.9;%寻优结果的初始值trace=zeros(2,MAXGEN);%建立区域描述器FieldD=[20;-1;2;1;0;1;1];%生成初始种群Chrom=crtbp(NIND,PRECI);具有性能跟踪和图像输出功能的程序代码第六页，共七十三页。%代计数器gen=0;%计算初始种群的十进制转换variable=bs2rv(Chrom,FieldD);%计算目标函数值 ObjV=variable.*sin(10*pi*variable)+2.0;具有性能跟踪和图像输出功能的程序代码第七页，共七十三页。whilegen<MAXGEN%分配适应度值FitnV=ranking(-ObjV);%选择SelCh=select('sus',Chrom,FitnV,GGAP);%重组 SelCh=recombin('xovsp',SelCh,0.7);%变异SelCh=mut(SelCh);具有性能跟踪和图像输出功能的程序代码第八页，共七十三页。%子代个体的十进制转换variable=bs2rv(SelCh,FieldD);%计算子代的目标函数值ObjVSel=variable.*sin(10*pi*variable)+2.0;%重插入子代的新种群[ChromObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); variable=bs2rv(Chrom,FieldD);具有性能跟踪和图像输出功能的程序代码第九页，共七十三页。%代计数器增加gen=gen+1;%输出最优解及其序号，并在目标函数图像%中标出，Y为最优解,I为种群的序号[Y,I]=max(ObjV);holdon;plot(variable(I),Y,'bo'); %遗传算法性能跟踪trace(1,gen)=max(ObjV);%每一代最优解trace(2,gen)=sum(ObjV)/length(ObjV);end具有性能跟踪和图像输出功能的程序代码第十页，共七十三页。%最优个体的十进制转换variable=bs2rv(Chrom,FieldD);holdon,gridon;plot(variable,ObjV,'b*');figure(2);plot(trace(1,:));holdon;plot(trace(2,:),'-.');gridofflegend('解的变化','种群均值的变化')具有性能跟踪和图像输出功能的程序代码第十一页，共七十三页。程序运行追踪结果：代序寻优结果自变量函数值12.02923.842322.21043.781932.35773.806342.48813.81552.64013.8072第十二页，共七十三页。程序运行追踪结果62.61333.755572.86583.807382.94003.804392.91203.8321102.86773.8321112.69733.8320122.79633.8326132.73173.8318第十三页，共七十三页。程序运行追踪结果142.77393.8318152.65213.8274162.78443.8274172.51583.8274182.77683.8285192.79073.8277202.80803.8379第十四页，共七十三页。程序运行追踪结果213.09473.8489223.15253.8496232.91473.8496242.91443.8493253.02593.8493经过25次迭代后最优解及种群均值的变化如下图：

第十五页，共七十三页。第十六页，共七十三页。第十七页，共七十三页。例题中用到的一些绘图函数说明1、绘制函数图像函数—fplot调用格式：FPLOT(FUN,LIMS)FPLOT(FUN,LIMS,TOL)FPLOT(FUN,LIMS,N)FPLOT(FUN,LIMS,'LineSpec')参数说明：LIMS=[XMINXMAX]或LIMS=[XMINXMAXYMINYMAX]——给出变量范围第十八页，共七十三页。TOL—可接受的相对误差，默认值为2e-3,即0.2%N—限定绘制的函数至少有N+1个点，默认值为1.最大步长限定为(1/N)*(XMAX-XMIN).LineSpec—指定线型FUN—要绘制图像的函数例如：fplot('variable.*sin(10*pi*variable)+2.0',[-1,2]);第十九页，共七十三页。2、定义零矩阵—ZEROSZEROS(N)—产生N阶零矩阵ZEROS(M,N)或ZEROS([M,N]))—产生M×N阶零矩阵；ZEROS(M,N,P,...)orZEROS([MNP...]))—产生M×N×P×...阶零矩阵；ZEROS(SIZE(A)))—产生与矩阵A一样大小的零矩阵。3、求最大值函数—MAXMAX(X)—X为向量时，返回向量X的最大元素；X为矩阵时，返回一个行向量，包含矩阵X的每列的最大元素.第二十页，共七十三页。4、绘图函数——PLOTPLOT(X,Y)—绘制以向量X为横坐标，向量Y为纵坐标的线图.如果X或Y是一个矩阵,则绘制多线图；如果X是一个标量，而Y是一个向量,则绘声绘色制的是length(Y)个不连续点。PLOT(Y)—以矩阵Y的行标为横坐标，每一列为纵坐标绘制图形。如果Y是一个复矩阵，PLOT(Y)相当于PLOT(real(Y),imag(Y)).PLOT(X,Y,S)—指定线型绘图，S为一字符串，指定绘图方式第二十一页，共七十三页。代表字符颜色代表符号线型C青色-实线M洋红--虚线Y黄色：点连线R红色-.点划线G绿色none无线B兰色W白色k黑色线型和颜色第二十二页，共七十三页。标记符绘图方式标记符绘图方式+十字号^反勾号o小圆圈v勾号*星号>大于号.小黑点<小于号X叉号pentagram五角星Square小正方形Hexagram六角星diamond菱形符号none无标记数据点标记字符第二十三页，共七十三页。5、建立图形窗口函数—FIGUREFIGURE(H)—使句柄为H的图形窗口为当前图形；如果图形窗口不存在，则建立一个句柄为H的图形窗口。GCF—返回当前图形窗口的句柄。6、设置网格线—GRIDGRIDON—给坐标系添加网格线GRIDOFF—去除坐标系中的网格线第二十四页，共七十三页。§5.2多元单峰函数优化实例第二十五页，共七十三页。简化程序代码如下：NIND=40; %个体数目MAXGEN=500; %最大遗传代数NVAR=20; %变量的维数PRECI=20; %编码长度GGAP=0.9; %代沟FieldD=[rep([PRECI],[1,NVAR]);rep([-512;512],…[1,NVAR]);rep([1;0;1;1],[1,NVAR])];%建立区域描述器Chrom=crtbp(NIND,NVAR*PRECI);%创建初始种群gen=0; %代计数器ObjV=objfun1(bs2rv(Chrom,FieldD));%计算初始种群个体的目标函数值whilegen<MAXGEN %迭代FitnV=ranking(ObjV); %分配适应度值SelCh=select('sus',Chrom,FitnV,GGAP);%选择SelCh=recombin('xovsp',SelCh,0.7); %重组SelCh=mut(SelCh); %变异ObjVSel=objfun1(bs2rv(SelCh,FieldD)); %计算子代目标函数值[ChromObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); %重插入gen=gen+1; %代计数器增加trace(gen,1)=min(ObjV); %遗传算法性能跟踪trace(gen,2)=sum(ObjV)/length(ObjV);end[Y,I]=min(ObjV)X=bs2rv(Chrom,FieldD);X(I,:)第二十六页，共七十三页。前5个语句进行的初始化工作，包括确定群体中个体数目，最大进化代数，编码长度，代沟，建立区域描述器等。第二十七页，共七十三页。1、矩阵复制函数REP功能用来复制一个矩阵调用格式：MatOut=rep(MatIn,REPN);参数说明MatIn—输入矩阵MatOut—输出矩阵REPN—一个二维向量，其中REPN(1)指定垂直方向复制次数；REPN(2)指定水平方向复制次数程序中的函数说明第二十八页，共七十三页。例如：>>MatIn=[123]MatIn=123>>REPN=[12]REPN=12>>MatOut=rep(MatIn,REPN)MatOut=123123程序中的函数说明第二十九页，共七十三页。例如：>>MatIn=[123]MatIn=123>>REPN=[21]REPN=21>>MatOut=rep(MatIn,REPN)MatOut=123123程序中的函数说明第三十页，共七十三页。例如：>>MatIn=[123]MatIn=123>>REPN=[32]REPN=32>>MatOut=rep(MatIn,REPN)MatOut=123123123123123123程序中的函数说明第三十一页，共七十三页。具有性能跟的程序代码如下：%定义遗传算法参数NIND=40; %个体数目MAXGEN=500;%最大遗传代数NVAR=20;%变量的维数PRECI=20; %编码长度GGAP=0.9;%代沟trace=zeros(MAXGEN,2);%建立区域描述器FieldD=[rep([PRECI],[1,NVAR]);rep([-512;512],…[1,NVAR]);rep([1;0;1;1],[1,NVAR])];%创建初始种群Chrom=crtbp(NIND,NVAR*PRECI);%代计数器gen=0; %计算初始种群个体的目标函数值ObjV=objfun1(bs2rv(Chrom,FieldD));第三十二页，共七十三页。whilegen<MAXGEN %迭代FitnV=ranking(ObjV); %分配适应度值SelCh=select('sus',Chrom,FitnV,GGAP); %选择SelCh=recombin('xovsp',SelCh,0.7); %重组SelCh=mut(SelCh); %变异

%计算子代目标函数值ObjVSel=objfun1(bs2rv(SelCh,FieldD));%重插入[ChromObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);%代计数器增加 gen=gen+1;%遗传算法性能跟踪 trace(gen,1)=min(ObjV);trace(gen,2)=sum(ObjV)/length(ObjV);end第三十三页，共七十三页。plot(trace(:,1));holdon;plot(trace(:,2),'-.');grid;legend('种群均值的变化','解的变化')%输出最优解及其对应的20个自变量的十进制值,Y为最优解,%

I为种群的序号[Y,I]=min(ObjV)X=bs2rv(Chrom,FieldD);X(I,:)第三十四页，共七十三页。functionObjVal=objfun1(Chrom,switch1);%目标函数OBJFUN1.M%调用格式:ObjVal=objfun1(Chrom,switch1)%输入参数:%Chrom-当前种群%switch1-如果Chrom==[]则switch1==1%并返回边界；如果switch1==2则返回标题；%如果switch1==3则返回全局最小值%输出变量：%ObjVal-各个体的目标函数值目标函数第三十五页，共七十三页。ifNind==0ifswitch1==2ObjVal=['DEJONGfunction1-'int2str(Dim)];elseifswitch1==3ObjVal=0;elseObjVal=100*[-5.12;5.12];ObjVal=ObjVal(1:2,ones(Dim,1));endelseifNvar==DimObjVal=sum((Chrom.*Chrom)')';elseerror('sizeofmatrixChromisnotcorrectforfunctionevaluation');end第三十六页，共七十三页。第三十七页，共七十三页。程序运行所求得的最优解为：y=1.0320第三十八页，共七十三页。§5.3多元多峰函数优化实例第三十九页，共七十三页。Shubert函数的图像为第四十页，共七十三页。目标函数functionz=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)); 第四十六页，共七十三页。whilegen<MAXGENFitnV=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));[ChromObjV]=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);holdon;plot(ObjV,'b*');grid;endend第四十八页，共七十三页。figure(3);clf;plot(trace(:,1));holdon;plot(trace(:,2),'-.');gridlegend('解的变化','种群均值的变化')第四十九页，共七十三页。第五十页，共七十三页。第五十一页，共七十三页。§5.4在多目标优化中的应用第五十二页，共七十三页。5.4.1多目标优化的概念第五十三页，共七十三页。5.4.1多目标优化的概念第五十四页，共七十三页。5.4.1多目标优化的概念第五十五页，共七十三页。5.4.1多目标优化的概念第五十六页，共七十三页。5.4.2多目标优化问题的遗传算法第五十七页，共七十三页。5.4.2多目标优化问题的遗传算法第五十八页，共七十三页。5.4.2多目标优化问题的遗传算法第五十九页，共七十三页。5.4.2多目标优化问题的遗传算法第六十页，共七十三页。5.4.2多目标优化问题的遗传算法第六十一页，共七十三页。5.4.2多目标优化问题的遗传算法第六十二页，共七十三页。5.4.2多目标优化问题的遗传算法第六十三页，共七十三页。5.4.3应用举例第六十四页，共七十三页。5.4.2多目标优化问题的遗传算法NIND=100;%个体数目MAXGEN=50;%最大遗传代数NVAR=2;%变量个数PRECI=20;%变量的二进制位数GGAP=0.9;%代沟trace1=[];trace2=[];trace3=[];%性能跟踪%建立区域描述器FieldD=[rep([PRECI],[1,NVAR]);[1,1;4,2];…rep([1;0;1;1],[1,NVAR])];第六十五页，共七十三页。5.4.2多目标优化问题的遗传算法Chrom=crtbp(NIND,NVAR*PRECI);%初始种群v=bs2rv(Chrom,FieldD);%初始种群十进制转换gen=1;whilegen<MAXGEN[NIND,N]=size(Chrom);M=fix(NIND/2);ObjV1=f1(v(1:M,:));%分组后第一目标函数值FitnV1=ranking(ObjV1);%分配适应度值第六十六页，共七十三页。5.4.2多目标优化问题的遗传算法SelCh1=select('sus',Chrom(1:M,:),FitnV1,GGAP);ObjV2=f2(v(M+1:NIND,:));%分组后第二目标函数值FitnV2=ranking(ObjV2);SelCh2=select('sus',Chr