GA.zip_GA优化_优化算法

clear; clc; NIND=40;%一代中的个体数目 MAXGEN=100;%最大遗传代数 NVAR=2; PRECI=25;%变量的二进制位数 GGAP=0.7;%相邻代之间的交配概率 trace=zeros(MAXGEN,2);%寻优结果的初始值全为0 FieldD=[25,25;0.2,0.2;1,1;0,0;0,0;1,1;1,1];%区域描述器 Chrom=crtbp(NIND,NVAR*PRECI);%初始种群的建立 gen=0; %代计数器 x=bs2rv(Chrom,FieldD); ObjV=peaks1(x(:,1),x(:,2)); while gen<MAXGEN, FitnV=ranking(-ObjV);%分配适应度的值 SelCh=select('sus',Chrom,FitnV,GGAP);%选择 SelCh=recombin('xovsp',SelCh,0.7);%重组 SelCh=mut(SelCh);%变异 x=bs2rv(SelCh,FieldD);%子代个体的十进制转换 ObjVSel=peaks1(x(:,1),x(:,2));%计算子代的目标函数 [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);%重新插入子代的新种群 gen=gen+1; %输出最优解及其对应的10个变量的十进制值 [Y,I]=max(ObjV); Y,bs2rv(Chrom(I,:),FieldD) %Y;bs2rv(Chrom(I,:),FieldD); X=bs2rv(Chrom(I,:),FieldD); Z0=peaks1(X(1,1),X(1,2)) %拟合坐标在原函数值 %zx2=N1dx2(X(1,1),X(1,2)) %拟合坐标二阶导 %zy2=N1dy2(X(1,1),X(1,2)) % zxy=N1dxdy(X(1,1),X(1,2)) % K=gauss(X(1,1),X(1,2)) %拟合坐标值在拟合函数中的高斯曲率 % d3=juli3(X(1,1),X(1,2)) % d5=juli5(X(1,1),X(1,2)) %d10=juli10(X(1,1),X(1,2)) % plot(I,Y,'bo'); % for i=1:length(I) % text(I,Y,['(',num2str(I),',',num2str(Y),')']) %在(x,y)的每一点处添加文字，文字内容为(x,y),其中每一点的x,y用num2str把数字转化成字符串 % end %trace (gen,1)=min(ObjV); %trace (gen,2)=sum(ObjV)/length(ObjV); % if(gen==MAXGEN) %figure(2); % plot(ObjV);hold on; %plot(ObjV,'b*');grid; %end end variable=bs2rv(Chrom,FieldD);%最优个体的十进制转换 hold on grid; plot(variable',ObjV','b*'); figure(3);clf; plot (trace (:,1)); hold on; plot(trace (:,2),'r-.');grid; legend('解的变化','种群均值的变化');