my_sga_mutation.asv

多种群的链式智能体遗传算法

% 2.6 变异
% 变异(mutation)，基因的突变普遍存在于生物的进化过程中。变异是指父代中的每个个体的每一位都以概率 pm 翻转，即由“1”变为“0”，
% 或由“0”变为“1”。遗传算法的变异特性可以使求解过程随机地搜索到解可能存在的整个空间，因此可以在一定程度上求得全局最优解。
%遗传算法子程序
%Name: mutation.m
%变异
function [newpop]=sga_mutation(popsize,chromlength,pop,pm,up_range)
%函数说明
%入口参数：pop         当前处理的种群,
%          popsize     种群大小,
%          chromlength 基因位长度,
%          up_range    自变量的相对上限,
%          down_range  自变量的相对下限
%出口参数：best_fit    各代最有适应度值的集合
%功能说明：初始种群pop通过遗传交叉变异在给定定义域中寻优，在遗传最大代数dai之前得到最优值就自动退出遗传，
%          判定最优值的标准：前后两代间最优值之间差值连续30次小于 阈值0.00001
%          如果达到遗传最大代数也没有达到退出标准，则强制结束遗传。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% popsize=10;
% pop=[1023,980,2133,726,1946,657,1702,1224,1522,807;
%      2025,2565,1773,2340,2354,1144,2061,240,1573,825]
% chromlength=12;
% pm=0.7;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


 pop1=pop(1,:);                %将种群矩阵拆分
 pop1;
 pop2=pop(2,:);
 pop2;
 
 pop1=dec2bin(pop1,chromlength);  %%各自转为二进制
 pop1;
 pop2=dec2bin(pop2,chromlength);
 pop2;
 
newpop1=pop1;
newpop2=pop2;
for i=1:popsize;

   if(rand<pm) %产生一随机数与变异概率比较
     mpoint=round(rand*chromlength);
%      mpoint
     if mpoint<=0
        mpoint=1;
     end

    if newpop1(i,mpoint)=='0'
       newpop1(i,mpoint)='1';
    else
       newpop1(i,mpoint)='0';
    end
   end  
   
    if(rand<pm) %产生一随机数与变异概率比较
     mpoint=round(rand*chromlength);
%      mpoint
     if mpoint<=0
        mpoint=1;
     end  
     
     if newpop2(i,mpoint)=='0'
        newpop2(i,mpoint)='1';
     else
        newpop2(i,mpoint)='0';
     end
    end

end

newpop1=bin2dec(newpop1);
newpop2=bin2dec(newpop2);


%如果产生的个体超过了相对定义域  则取上限
for i=1:popsize
    if newpop1(i)>=up_range
%         pop_init=randperm(20000); 
%         newpop1(i)=newpop1(i)-3095;  
        newpop1(i)=rand;
    end
    
    if newpop2(i)>=up_range
%         newpop2(i)=newpop2(i)-3095;       
%         pop_init=randperm(20000); 
        newpop2(i)=rand;       
    end
    
end

newpop(1,:)=newpop1';
newpop(2,:)=newpop2';
newpop;