ga_main.m

基本的遗传算法matlab语言实现

%**********************%主函数*****************************************
function main()
global chrom lchrom oldpop newpop varible fitness  popsize  sumfitness %定义全局变量
global pcross pmutation temp bestfit maxfit gen bestgen
global maxgen po pp mp np
lchrom=12;         %染色体长度
popsize=200;        %种群大小
pcross=0.6;        %交叉概率
pmutation=0.05;     %变异概率
maxgen=100;        %最大代数
po=0.1;            %淘汰概率
pp=0.1;            %保护概率
mp=floor(pp*popsize);  %保护的个数
np=floor(po*popsize);  %淘汰的个数
%
initpop;              % 初始种群
%
for gen=1:maxgen
    generation;
end
%
best;
bestfit     % 最佳个体适应度值输出
bestgen    % 最佳个体所在代数输出
gen=1:maxgen;
plot(gen,maxfit(1,gen));  % 进化曲线
%
%**********************  产生初始种群  ************************************ 
%
%初始群体长生的很好，可以参考
function initpop()
global lchrom oldpop popsize chrom
for i=1:popsize
    chrom=rand(1,lchrom);    %  lchrom=12  染色体长度
    for j=1:lchrom
        if chrom(1,j)<0.5
            chrom(1,j)=0;
       else 
           chrom(1,j)=1;
        end
    end
    oldpop(i,1:lchrom)=chrom;    
end
%
%*************************%产生新一代个体**********************************
%
function generation()
objfun;                  %计算适应度值
pp_po;                   %保护淘汰操作
select;                  %选择操作
crossover;
mutation;
%
%************************%计算适应度值************************************
%
function objfun()
global lchrom oldpop fitness popsize chrom maxfit gen varible
a=0;
%b=100;
b=30;
for i=1:popsize
    chrom=oldpop(i,:);
    c=decimal(chrom);
    varible(1,i)=a+c*(b-a)/(2.^lchrom-1);  %对应变量值
    fitness(1,i)=varible(1,i)*varible(1,i);
%    fitness(1,i)=20*cos(0.25*varible(1,i))-12*sin(0.33*varible(1,i))+40;  %个体适应度函数值
end
%
lsort;                          % 个体排序
%
maxfit(1,gen)=max(fitness);     %求最大适应度值
%
%************************%二进制转十进制**********************************
%
function c=decimal(chrom)
global lchrom popsize
c=0;
for j=1:lchrom
    c=c+chrom(1,j)*2.^(lchrom-j);
end
%
%*************************%个体排序****************************************
%
function lsort()
global popsize fitness oldpop
for i=1:popsize
    j=i+1;
    while j<=popsize
        if fitness(1,i)>fitness(1,j)
            tf=fitness(1,i);
            tc=oldpop(i,:);
            fitness(1,i)=fitness(1,j);
            oldpop(i,:)=oldpop(j,:);
            fitness(1,j)=tf;
            oldpop(j,:)=tc;
        end
        j=j+1;
    end
end
%**************************************************************************
%
%*************************%保护/淘汰操作************************************
function pp_po()
global popsize oldpop np
i=np+1;                      % np=floor(po*popsize);  %淘汰的个数
j=1;
while i<=popsize
    toldpop(j,:)=oldpop(i,:);
    j=j+1;
    i=i+1;
end
for i=1:(popsize-np)
    oldpop(i,:)=toldpop(i,:);
end
%**************************************************************************
%*************************%转轮法选择操作**********************************
function select()
global fitness popsize sumfitness oldpop temp mp np
sumfitness=0;
for i=1:(popsize-np-mp)
    sumfitness=sumfitness+fitness(1,i);
end
for i=1:(popsize-mp-np)
    p(1,i)=fitness(1,i)/sumfitness;   % 个体染色体的选择概率
end
q=cumsum(p);                          % 个体染色体的累积概率
b=sort(rand(1,(popsize-mp)));         % mp保护个体数
j=1;
k=1;
while j<=(popsize-mp)
    if b(1,j)<q(1,k)
       temp(j,:)=oldpop(k,:);
       j=j+1;
   else
       k=k+1;
   end
end
for i=(popsize-mp+1):popsize
    j=popsize-np-mp+1;
    temp(i,:)=oldpop(j,:);
    j=j+1;
end
%**************************************************************************
%**************************%交叉操作***************************************
function crossover()
global temp popsize pcross lchrom mp
n=floor(pcross*(popsize-mp));           %交叉发生的次数
if rem(n,2)~=0
    n=n+1;
end
j=1;
m=0;
for i=1:(popsize-mp)
    p=rand;
    if p<pcross
        parent(j,:)=temp(i,:);      %产生两个父代
        k(1,j)=i;
        j=j+1;
        m=m+1;
        if (j==3)&(m<=n)
            pos=round(rand*(lchrom-1))+1;
            for i=1:pos
                child1(1,i)=parent(1,i);
                child2(1,i)=parent(2,i);
            end
            for i=(pos+1):lchrom
                child1(1,i)=parent(2,i);
                child2(1,i)=parent(1,i);
            end
            i=k(1,1);
            j=k(1,2);
            temp(i,:)=child1(1,:);
            temp(j,:)=child2(1,:);
            j=1;
        end
       %j=1;
    end
end
%**************************************************************************
%****************************%变异操作*************************************
function mutation()
global popsize lchrom pmutation temp newpop oldpop mp 
m=lchrom*(popsize-mp);               %总的基因数
n=round(pmutation*m);           %变异发生的次数
for i=1:n
    k=round(rand*(m-1))+1;
    j=ceil(k/lchrom);
    l=rem(k,lchrom);
    if l==0
       temp(j,lchrom)=~temp(j,lchrom);
   else
       temp(j,l)=~temp(j,l);
   end
end
for i=1:popsize
    newpop(i,:)=temp(i,:);     %产生新的个体
    oldpop(i,:)=newpop(i,:);
end
%**************************************************************************
%*********************%最佳个体********************************************
function best()
global maxfit bestfit gen maxgen bestgen
bestfit=maxfit(1,1);
gen=2;
while gen<=maxgen
      if maxfit(1,1)<maxfit(1,gen)
         bestfit=maxfit(1,gen);
         bestgen=gen;
     end
     gen=gen+1;
 end
%**************************************************************************