三分钟学会遗传算法


三分钟学会遗传算法

文章插图
遗传算法
此节介绍最著名的遗传算法(GA) 。遗传算法属于进化算法,基本思想是取自“物竞天泽、适者生存”的进化法则 。简单来说,遗传算法就是将问题编码成为染色体,然后经过不断选择、交叉、变异等操作来更新染色体的编码并进行迭代,每次迭代保留上一代好的染色体,丢弃差的染色体,最终达到满足目标的最终染色体 。整个流程由下图构成(手写,见谅 -_-!!)
三分钟学会遗传算法

文章插图
流程图
步骤由以下几步构成:
编码(coding)——首先初始化及编码 。在此步,根据问题或者目标函数(objective function)构成解数据(solutions),在遗传算法中,该解数据就被称为染色体(chromosome) 。值得一提的是,遗传算法为多解(population based)算法,所以会有多条染色体 。初始化中会随机生成N条染色体,, 这里表示染色体包含了n条 。其中,这里表示第i条染色体由d维数值构成 。GA会以这个N个数据作为初始点开始进行进化 。
【三分钟学会遗传算法】评估适应度(evaluate fitness)——这一步用染色体来进行目标函数运算,染色体的好坏将被指明 。
选择(selection)——从当前染色体中挑选出优良的个体,以一定概率使他们成为父代进行交叉或者变异操作,他们的优秀基因后代得到保留 。物竞天择这里得以体现 。
交叉(crossover)——父代的两个两个染色体,通过互换染色体构成新的染色体 。例如下图,父亲母亲各提供两个基因给我 。这样我既保留了父母的基于,同时又有自己的特性 。
三分钟学会遗传算法

文章插图
交叉
变异(mutation)——以一定概率使基因发生突变 。该算子一般以较低概率发生 。如下图所示:
三分钟学会遗传算法

文章插图
变异
下面我们将一步一步为各位呈现如何用matlab编写一个简单的GA算法 。
本问题为实数最小化minimization问题 。我们需要在解空间内找到最小值或近似最小值,此处我们使用sphere函数作为目标函数(读者可以自行修改为其他的目标函数) 。
三分钟学会遗传算法

文章插图
sphere function
  • 初始化:在这一步中,我们将在给定问题空间内生成随机解,代码如下:
% %% 初始化% % 输入:chromes_size,dim维数,lb下界,ub上界% % 输出:chromes新种群function chromes=init_chromes(chromes_size,dim,lb,ub) % 上下界中随机生成染色体 chromes = rand(chromes_size,dim)*(ub-lb)+lb;end
  • 选择:选择是从当前代中挑选优秀的染色体保留以繁殖下一代 。我们这里采取的方法是俄罗斯轮盘选择方式 。谁的解优,谁获得选中的概率越高 。首先,我们需要先求出各染色体的fitness倒数 。
  • ,然后求出各染色体的比重,比重越大,被保留机会越大 。。代码如下:
%% 选择%俄罗斯轮盘赌function [newchromes,newfitness] = selection(chromes,fitness) weights = 1./fitness; % fitness倒数 totalfit=sum(weights); % 累加所有weights totalf = weights./totalfit; %求出各染色体比重 index = []; for i = 1:size(chromes,1) % 循环选出较优染色体 pick = rand; while pick == 0 pick = rand; end for j = 1:size(chromes,1) pick = pick - totalf(j); if pick<0 index = [index j]; break end end end newchromes =chromes(index,:); newfitness = fitness(index);end
  • 交叉:此步会随机选取两个选择过后的染色体作为父代,从两个染色体中各截取一部分基因生成新的染色体,代码如下:
%% 交叉function newchromes = crossover(chromes,pc) % 生成一个新的解 newchromes = ones(size(chromes)); for i = 1:size(chromes,1) % 随机选取两个染色体 parents=randperm(10,2); %随机选取一个位置 pos = round(rand*size(chromes,2)); if(rand<pc) % 交叉生成新的解 newchromes(i,:)=[chromes(parents(1),1:pos) chromes(parents(2),pos+1:size(chromes,2))]; else newchromes(i,:)=chromes(i,:); end endend
  • 变异:以一各小概率生成随机变异一个gene,代码如下:
% 变异function newchromes= muatation(chromes,pm,lb,ub) for i = 1:size(chromes,1) newchromes(i,:) = chromes(i,:); if (rand<pm) %随机选取一位 pos = ceil(rand*size(chromes,2)); %变异 newchromes(i,pos)= rand*(ub-lb)+lb; end endend
  • 主函数,首先初始化各参数,然后进行迭代,当满足终止条件停止:
% 清除workspace,清屏clearclc % 染色体数量chromes_size = 20;% 问题维数dim = 10;% 交叉概率pc =0.9;% 变异概率pm = 0.2;% 问题上下边界lb = -1;ub = 1;% 循环次数maxiter = 1000;% 目标方程namefunc= 'objfun';fd = str2func(namefunc);?% 初始化chromes = init_chromes(chromes_size,dim,lb,ub);% 求个染色体fitness for i = 1:chromes_size fitness(i)=feval(fd,chromes(i,:)); end% 求出最优解 [bestfitness bestindex]=min(fitness); bestchrome = chromes(bestindex,:); % 主循环for iter=1:maxiter % 选择 [chromes,newfitness] = selection(chromes,fitness); % 交叉 chromes= crossover(chromes,pc); % 变异 chromes= muatation(chromes,pm,lb,ub); % 更新最优 for i = 1:chromes_size fitness(i)=feval(fd,chromes(i,:)); if fitness(i)<bestfitness bestfitness= fitness(i) bestchrome = chromes(i,:); end end trace(iter)=bestfitness; end?% 绘制fitness曲线plot(trace)title('Fitness curve')xlabel('Iterations')ylabel('Fitness value')


推荐阅读