简单的遗传算法MATLAB实现

遗传算法是对达尔文生物进化理论的简单模拟，其遵循“适者生存”、“优胜略汰”的原理。遗传算法模拟一个人工种群的进化过程，并且通过选择、杂交以及变异等机制，种群经过若干代以后，总是达到最优（或近最优）的状态。

自从遗传算法被提出以来，其得到了广泛的应用，特别是在函数优化、生产调度、模式识别、神经网络、自适应控制等领域，遗传算法更是发挥了重大的作用，大大提高了问题求解的效率。遗传算法也是当前“软计算”领域的重要研究课题。

本文首先结合MATLAB对遗传算法实现过程进行详细的分析，然后通过1个实际的函数优化案例对其应用进行探讨。

1. 遗传算法实现过程

现实生活中很多问题都可以转换为函数优化问题，所以本文将以函数优化问题作为背景，对GA的实现过程进行探讨。大部分函数优化问题都可以写成求最大值或者最小值的形式，为了不是一般性，我们可以将所有求最优值的情况都转换成求最大值的形式，例如，求函数f(x)的最大值，

若是求函数f(x)的最小值，可以将其转换成

g(x)=-f(x)，然后求g(x)的最大值，

这里x可以是一个变量，也可是是一个由k个变量组成的向量, x=(x1, x2, …, xk)。每个xi, i=1,2,…,k, 其定义域为Di，Di=[ai, bi]。

一般规定f(x)在其定义域内只取正值，若不满足，可以将其转换成以下形式，

其中C是一个正常数。 1.1 编码与解码

要实现遗传算法首先需要弄清楚如何对求解问题进行编码和解码。对于函数优化问题，一般来说，有两种编码方式，一是实数编码，一是二进制编码，两者各有优缺点，二进制编码具有稳定性高、种群多样性大等优点，但是需要的存储空间大，需要解码过程并且难以理解；而实数编码直接用实数表示基因，容易理解并且不要解码过程，但是容易过早收敛，从而陷入局部最优。本文以最常用的二进制编码为例，说明遗传编码的过程。

从遗传算法求解的过程来看，需要处理好两个空间的问题，一个是编码空间，另一个是解空间，如下图所示

从解空间到编码空间的映射过程成为编码过程；从编码空间到解空间的映射过程成为解码过程。下面就以求解一个简单的一维函数f(x) = -(x-1)^2+4, x的取值范围为[-1,3]最大值为例，来说明编码及解码过程。 编码：

在编码之前需要确定求解的精度，在这里，我们设定求解的精度为小数点后四位，即1e-4。这样可以将每个自变量xi的解空间划分为

个等分。以上面

这个函数为例，即可以将x的解空间划分为(3-(-1))*1e+4=40000个等分。使ni满足

，这里ni表示使上式成立的最小整数，

即表示自变量xi的基因串的长度。因为2<40000<2，这里ni取16。例如0000110110000101就表示一个解空间中的基因串。表示所有自变量x=(x1, x2, …, xk)的二进制串的总长度称为一个染色体(Chromosome)的

15

16

长度或者一个个体(Individual)的长度，

。编码过程一般在实现遗传算法之前需

要指定。 解码：

解码即将编码空间中的基因串翻译成解空间中的自变量的实际值的过程。对于二进制编码而言，每个二进制基因串都可以这样翻译成一个十进制实数值，

。例如基因串

0000110110000101，可以翻译为

，这里二进制基因串转变成十

进制是从左至右进行的。 1.2 初始化种群

在开始遗传算法迭代过程之前，需要对种群进行初始化。设种群大小为pop_size，每个染色体或个体的长度为chromo_size，种群的大小决定了种群的多样性，而染色体的长度则是由前述的编码过程决定的。一般随机生成初始种群，但是如果知道种群的实际分布，也可以按照此分布来生成初始种群。假设生成的初始种群为(v1, v2, …, vpop_size)。 1.3 选择操作

选择操作即从前代种群中选择个体到下一代种群的过程。一般根据个体适应度的分布来选择个体。以初始种群(v1, v2, …, vpop_size)为例，假设每个个体的适应度为(fitness(v1), fitness(v2),…, fitness(vpop_size))，一般适应度可以按照解码的过程进

行计算。以轮盘赌的方式选择个体，如下图

随机转动一下轮盘，当轮盘停止转动时，若指针指向某个个体，则该个体被选中。很明显，具有较高适应度的个体比具有较低适应度的个体更有机会被选中。但是这种选择具有随机性，在选择的过程中可能会丢失掉比较好的个体，所以可以使用精英机制，将前代最优个体直接选到下一代中。

轮盘赌选择具体算法如下（这里假定种群中个体是按照适应度从小到大进行排列的，如果不是，可以按照某种排序算法对种群个体进行重排）： Selection Algorithm

var pop, pop_new;/*pop为前代种群，pop_new为下一代种群*/

varfitness_value,

fitness_table;/*fitness_value为种群的适应度，fitness_table为种群累积适应度*/ fori=1:pop_size

r = rand*fitness_table(pop_size);/*随机生成一个随机数，在0和总适应度之间，因为

fitness_table(pop_size)为最后一个个体的累积适应度，即为总适应度*/ first = 1; last = pop_size;

mid = round((last+first)/2); idx = -1;

/*下面按照排中法选择个体*/ while (first <= last) && (idx == -1) if r >fitness_table(mid) first = mid;

elseif r

mid = round((last+first)/2); if (last - first) == 1 idx = last; break; end if end while

for j=1:chromo_size pop_new(i,j)=pop(idx,j); end for end for

/*是否精英选择*/ if elitism

p = pop_size-1; else

p = pop_size; end if fori=1:p

for j=1:chromo_size

pop(i,j) = pop_new(i,j);/*若是精英选择，则只将pop_new前pop_size-1个个体赋给pop，最后一个为前代最优个体保留*/

end for end for 1.3 交叉操作

交叉操作是对任意两个个体进行的（在这里我们实现的算法是直接对相邻的两个个体进行的）。随机选择两个个体，如下图所示

然后随机生成一个实数0<=r<=1, 如果r

if(rand

cross_pos = round(rand * chromo_size);/*交叉位置*/

if or (cross_pos == 0, cross_pos == 1) continue;/*若交叉位置为0或1，则不进行交叉*/ end if

for j=cross_pos:chromo_size

pop(i,j)<->pop(i+1,j);/*交换*/ end for end if end for 1.4 变异操作

变异操作是对单个个体进行的。首先生成一个随机实数0<=r<=1, 如果r

如个体需要进行变异操作，首先需要确定变异位置(rand*chromo_size)，若为0则不进行变异，否则则

对该位置的二进制数字进行变异：1变成0, 0变成1.（当然也可以选择多点进行变异） 单点变异的具体算法描述如下： Mutation algorithm fori=1:pop_size

if rand

mutate_pos = round(rand*chromo_size);/*变异位置*/

ifmutate_pos == 0

continue;/*若变异位置为0，则不进行变异*/ end if

pop(i,mutate_pos) = 1 - pop(i, mutate_pos);/*将变异位置上的数字至反*/ end if end for

1.5 遗传算法流程

遗传算法计算流程图如下图所示

end end cleari; clear j; clear k; clear min; clear temp; clear temp1;

选择操作： %轮盘赌选择操作 %pop_size: 种群大小 %chromo_size: 染色体长度 %cross_rate: 是否精英选择

function selection(pop_size, chromo_size, elitism) global pop;

globalfitness_table;

fori=1:pop_size

r = rand * fitness_table(pop_size); first = 1; last = pop_size;

mid = round((last+first)/2); idx = -1;

while (first <= last) && (idx == -1) if r >fitness_table(mid) first = mid;

elseif r

mid = round((last+first)/2); if (last - first) == 1 idx = last; break; end end

for j=1:chromo_size pop_new(i,j)=pop(idx,j); end end if elitism

p = pop_size-1; else

p = pop_size; end fori=1:p

for j=1:chromo_size pop(i,j) = pop_new(i,j); end end cleari; clear j; clearpop_new; clear first; clear last; clearidx; clear mid;

交叉操作： %单点交叉操作 %pop_size: 种群大小 %chromo_size: 染色体长度 %cross_rate: 交叉概率

function crossover(pop_size, chromo_size, cross_rate) global pop; fori=1:2:pop_size if(rand

cross_pos = round(rand * chromo_size); if or (cross_pos == 0, cross_pos == 1) continue; end

for j=cross_pos:chromo_size temp = pop(i,j); pop(i,j) = pop(i+1,j); pop(i+1,j) = temp; end end

end cleari; clear j; clear temp; clearcross_pos;

变异操作： %单点变异操作 %pop_size: 种群大小 %chromo_size: 染色体长度 %cross_rate: 变异概率

function mutation(pop_size, chromo_size, mutate_rate) global pop;

fori=1:pop_size if rand

mutate_pos = round(rand*chromo_size); ifmutate_pos == 0 continue;

end

pop(i,mutate_pos) = 1 - pop(i, mutate_pos); end end cleari;

clearmutate_pos; 打印算法迭代过程： %打印算法迭代过程

%generation_size: 迭代次数

functionplotGA(generation_size) globalfitness_avg; x = 1:1:generation_size; y = fitness_avg; plot(x,y) 算法主函数： %遗传算法主函数 %pop_size: 输入种群大小 %chromo_size: 输入染色体长度 %generation_size: 输入迭代次数

%cross_rate: 输入交叉概率 %cross_rate: 输入变异概率 %elitism: 输入是否精英选择 %m: 输出最佳个体 %n: 输出最佳适应度 %p: 输出最佳个体出现代 %q: 输出最佳个体自变量值

function [m,n,p,q] = GeneticAlgorithm(pop_size, chromo_size, generation_size, cross_rate, mutate_rate, elitism)

global G ; %当前代

global fitness_value;%当前代适应度矩阵 global best_fitness;%历代最佳适应值 global fitness_avg;%历代平均适应值矩阵 global best_individual;%历代最佳个体 global best_generation;%最佳个体出现代

fitness_avg = zeros(generation_size,1);

disp \

fitness_value(pop_size) = 0.; best_fitness = 0.; best_generation = 0;

initilize(pop_size, chromo_size);%初始化 for G=1:generation_size

fitness(pop_size, chromo_size); %计算适应度

rank(pop_size, chromo_size); %对个体按适应度大小进行排序

selection(pop_size, chromo_size, elitism);%选择操作

crossover(pop_size, chromo_size, cross_rate);%交叉操作

mutation(pop_size, chromo_size, mutate_rate);%变异操作 end

plotGA(generation_size);%打印算法迭代过程 m = best_individual;%获得最佳个体 n = best_fitness;%获得最佳适应度

p = best_generation;%获得最佳个体出现代

%获得最佳个体变量值，对不同的优化目标，此处需要改写 q = 0.;

for j=1:chromo_size ifbest_individual(j) == 1 q = q+2^(j-1); end end

q = -1+q*(3.-(-1.))/(2^chromo_size-1); cleari; clear j;

2. 案例研究

对上一节中的函数进行优化，设置遗传算法相关参数，程序如下

functionrun_ga()

elitism = true;%选择精英操作 pop_size = 20;%种群大小

chromo_size = 16;%染色体大小 generation_size = 200;%迭代次数 cross_rate = 0.6;%交叉概率 mutate_rate = 0.01;%变异概率

[m,n,p,q] = GeneticAlgorithm(pop_size, chromo_size, generation_size, cross_rate, mutate_rate,elitism); disp \最优个体\m

disp \最优适应度\n

disp \最优个体对应自变量值\q

disp \得到最优结果的代数\p clear;

结果如下： \最优个体\

本文来源：https://www.bwwdw.com/article/54s6.html