mytry4.m

自己编的遗传算法

clear;
figure(1);
fplot('5776/variabl+variabl',[50,100]);   %画出函数曲线
%定义遗传算法参数
NIND=50;        %个体数目(Number of individuals)
MAXGEN=25;      %最大遗传代数(Maximum number of generations)
PRECI=20;       %变量的二进制位数(Precision of variables)
MIGGEN=20;        %每20代迁移个体
MIGR=0.2;              %迁移率
GGAP=0.8;       %代沟(Generation gap)
trace=zeros(2, MAXGEN);                        %寻优结果的初始值
FieldD=[20;50;100;0;0;1;1];                      %区域描述器(Build field descriptor)
Chrom=crtbp(NIND, PRECI);                      %初始种群
gen=0;                                         %代计数器
variable=bs2rv(Chrom, FieldD);                 %计算初始种群的十进制转换
ObjV=5776./variable+variable;        %计算目标函数值

while gen<MAXGEN
   FitnV=ranking(-ObjV);                                  %分配适应度值(Assign fitness values)         
   SelCh=select('sus', Chrom, FitnV, GGAP);               %选择
   SelCh=recombin('xovsp', SelCh, 0.7);                   %重组
   SelCh=mut(SelCh);                                      %变异
 
   variable=bs2rv(SelCh, FieldD);                         %子代个体的十进制转换
   ObjVSel=5776./variable+variable;            %计算子代的目标函数值
   [Chrom ObjV]=reins(Chrom, SelCh, 1, 1, ObjV, ObjVSel); %重插入子代的新种群
   variable=bs2rv(Chrom, FieldD);
   gen=gen+1;                                             %代计数器增加
   %输出最优解及其序号，并在目标函数图像中标出，Y为最优解,I为种群的序号
   [Y, I]=min(ObjV);hold on;
   plot(variable(I), Y, 'bo');
       Y,bs2rv(Chrom(I,:),FieldD);
     
   trace(1, gen)=min(ObjV);                               %遗传算法性能跟踪
   trace(2, gen)=sum(ObjV)/length(ObjV);
   
   
      %在子种群之间迁移个体
       if(rem(gen,MIGGEN)==0)
        [Chrom, ObjV]=migrate(Chrom,1,[MIGR, 1, 1],ObjV);
       end
    
       
end
min(ObjV)
variable=bs2rv(Chrom, FieldD);                            %最优个体的十进制转换
hold on, grid;
plot(variable,ObjV,'b*');
figure(2);
plot(trace(1,:));
hold on;
plot(trace(2,:),'-.');grid
legend('解的变化','种群均值的变化')