1. Johngo学长首页
  2. 技术点详细复盘专题
  3. 人工智能

python遗传算法之geatpy学习

• 人工智能 • 阅读 100

😻今天我们来学习python中的遗传算法的使用，我们这里使用的是geatpy的包进行学习，本博客主要从geatpy中的各种数据结构一步一步进行学习，请大家耐心看完。

🐤其实以前也学习过遗传算法，但是主要使用matlab进行编程的，后面觉得matlab太麻烦了，还是使用python方便些，于是开始继续学习。

目录

首先是安装geatpy，使用pip3命令进行安装即可：

pip3 install geatpy

出现如下提示即安装成功：

python遗传算法之geatpy学习
  1. geatpy的基础数据结构

geatpy中的大部分数据都是都是使用numpy的数组进行存储和计算的，下面我将介绍遗传算法中的概念是如何用numpy数据表示，以及行和列的含义。

2.1 种群染色体

遗传算法中最重要的就是个体的染色体表示，在geatpy中种群染色体用Chrom表示，这是一个二维数组，其中每一行对应一个个体的染色体编码，Chrom的结构如下：其中lind表示编码的长度，Nind表示的是种群的规模（个体数量）。

python遗传算法之geatpy学习

; 2.2 种群表现型

种群表现型是是指种群染色体矩阵Chrom经过解码后得到的基因表现型矩阵Phen，每一行对应一个个体，每一列对应一个决策变量，Phen的结构如下：其中Nvar表示变量的个数

python遗传算法之geatpy学习
Phen的值与采用的解码方式有关。Geatpy提供二进制/格雷码编码转十进制整数或实数的解码方式。另外，在Geatpy也可以使用不需要解码的”实值编码”种群，这种种群的染色体的每一位就对应着决策变量的实际值，即这种编码方式下Phen等价Chrom。

2.3 目标函数值

Geatpy采用numpy的array类型矩阵来存储种群的目标函数值。一般命名为ObjV，每一行对应每一个个体，因此它拥有与Chrom相同的行数；每一列对应一个目标函数，因此对于单目标函数，ObjV会只有1列；而对于多目标函数，ObjV会有多列, ObjV的表示形式如下：

python遗传算法之geatpy学习

; 2.4 个体适应度

Geatpy采用列向量来存储种群个体适应度（适应度函数计算而来）。一般命名为FitnV，它同样是numpy的array类型，每一行对应种群矩阵的每一个个体。因此它拥有与Chrom相同的行数，FitnV的格式如下：

python遗传算法之geatpy学习
注意：Geatpy中的适应度遵循”最小适应度为0″的约定。

2.5 违反约束程度矩阵

Geatpy采用numpy的array类型的矩阵CV(Constraint Violation Value)来存储种群个体违反各个约束条件的程度。命名为CV，它的每一行对应种群的每一个个体，因此它拥有与Chrom相同的行数；每一列对应一个约束条件，因此若有一个约束条件，那么CV矩阵就会只有一列，如有多个约束条件，CV矩阵就会有多列。如果设有num个约束，则CV矩阵的结构如下图所示：

python遗传算法之geatpy学习
CV矩阵的某个元素若小于或等于0，则表示该元素对应的个体满足对应的约束条件。若大于0，则表示违反约束条件，在大于0的条件下值越大，该个体违反该约束的程度就越高。Geatpy提供两种处理约束条件的方法，一种是罚函数法，另一种是可行性法则。在使用可行性法则处理约束条件时，需要用到CV矩阵。

; 2.6 译码矩阵

所谓的译码矩阵，只是用来描述种群染色体特征的矩阵，如染色体中的每一位元素所表达的决策变量的范围、是否包含范围的边界、采用二进制还是格雷码、是否使用对数刻度、染色体解码后所代表的决策变量的是连续型变量还是离散型变量等等。

在只使用工具箱的库函数而不使用Geatpy提供的面向对象的进化算法框架时，译码矩阵可以单独使用。若采用Geatpy提供的面向对象的进化算法框架时，译码矩阵可以与一个存储着种群染色体编码方式的字符串Encoding来配合使用。

目前Geatpy中有三种Encoding，分别为：

  • BG：(二进制/格雷码)
  • RI：（(实整数编码，即实数和整数的混合编码）
  • P：（排列编码，即染色体每一位的元素都是互异)

注：‘RI’和’P’编码的染色体都不需要解码，染色体上的每一位本身就代表着决策变量的真实值，因此”实整数编码”和”排列编码”可统称为”实值编码”

以BG编码为例，我们展示一下编译矩阵FieldD。FieldD的结构如下：

python遗传算法之geatpy学习
其中，lens,lb,ub,codes,scales,lbin,ubin,varTypes都是行向量，其长度等于决策变量的个数。
  • lens：代表以条染色体中，每个子染色体的长度。
  • lb：代表每个变量的上界
  • ub：代表每个变量的下界
  • codes：代表染色体字串用的编码方式，[1,0,1]代表第一个变量用的格雷编码，第二个变量用的二进制编码，第3个变量用的格雷编码。
  • scales：指明每个子串用的是算术刻度还是对数刻度。scales[i] = 0为算术刻度，scales[i] = 1为对数刻度（对数刻度很少用，可以忽略。）
  • lbin：代表变量上界是否包含其范围边界。0代表不包含，1代表包含。'[ ‘和 ‘（’ 的区别
  • ubin：代表变量下界是否包含其范围边界。0代表不包含，1代表包含。
  • varTypes：代表决策变量的类型，元素为0表示对应位置的决策变量是连续型变量；1表示对应的是离散型变量。

例如：有以下一个译码矩阵

python遗传算法之geatpy学习
它表示待解码的种群染色体矩阵Chrom解码后可以表示成3个决策变量，每个决策变量的取值范围分别是[1,10], [2,9], [3,15]。其中第一第二个变量采用的是二进制编码，第三个变量采用的是格雷编码，且第一、第三个决策变量为连续型变量；第二个为离散型变量。 

import geatpy as ea
Phen = ea.bs2ri(Chrom, FieldD)

2.7 进化追踪器

在使用Geatpy进行进化算法编程时，常常建立一个进化追踪器(如pop_trace)来记录种群在进化的过程中各代的最优个体，尤其是采用无精英保留机制时，进化追踪器帮助我们记录种群在进化过程中的最优个体。待进化完成后，再从进化追踪器中挑选出”历史最优”的个体。这种进化记录器有多种，其中一种是numpy的array类型的，结构如下：其中MAXGEN是种群进化的代数（迭代次数）。

python遗传算法之geatpy学习
trace的每一列代表不同的指标，比如第一列记录各代种群的最佳目标函数值，第二列记录各代种群的平均目标函数值…trace的每一行对应每一代，如第一行代表第一代，第二行代表第二代…另外一种进化记录器是一个列表，列表中的每一个元素都是一个拥有相同数据类型的数据。比如在Geatpy的面向对象进化算法框架中的pop_trace，它是一个列表，列表中的每一个元素都是历代的种群对象。

; 3. geatpy的种群结构

3.1 Population类

在Geatpy提供的面向对象进化算法框架中，种群类(Population)是一个存储着与种群个体相关信息的类。它有以下基本属性：

  • sizes : int -种群规模，即种群的个体数目。
  • ChromNum : int -染色体的数目，即每个个体有多少条染色体。
  • Encoding : str -染色体编码方式。
  • Field : array -译码矩阵，可以是FieldD或FieldDR。
  • Chrom : array -种群染色体矩阵，每一行对应一个个体的一条染色体。
  • Lind : int -种群染色体长度。
  • ObjV : array -种群目标函数值矩阵。
  • FitnV : array -种群个体适应度列向量。
  • CV : array -种群个体违反约束条件程度的矩阵。
  • Phen : array -种群表现型矩阵。

可以直接对种群对象进行提取个体、个体合并等操作，比如pop1和pop2是两个种群对象，则通过语句”pop3 = pop1 + pop2″，即可把两个种群的个体合并，得到一个新的种群。在合并的过程中，实际上是把种群的各个属性进行合并，然后用合并的数据来生成一个新的种群(详见Population.py)。又比如执行语句”pop3 = pop1[[0]]”，可以把种群的第0号个体抽取出来，得到一个新的只有一个个体的种群对象pop3。值得注意的是，种群的这种个体抽取操作要求下标必须为列表或是Numpy array类型的行向量，不能是标量(详见Population.py)

3.2 PsyPopulation类

PsyPopulation类是Population的子类，它提供Population类所不支持的多染色体混合编码。它有以下基本属性：

  • sizes : int -种群规模，即种群的个体数目。
  • ChromNum : int -染色体的数目，即每个个体有多少条染色体。
  • Encodings : list -存储各染色体编码方式的列表。
  • Fields : list -存储各染色体对应的译码矩阵的列表。
  • Chroms : list -存储种群各染色体矩阵的列表。
  • Linds : list -存储种群各染色体长度的列表。
  • ObjV : array -种群目标函数值矩阵。
  • FitnV : array -种群个体适应度列向量。
  • CV : array -种群个体违反约束条件程度的矩阵。
  • Phen : array -种群表现型矩阵。

可见PsyPopulation类基本与Population类一样，不同之处是采用Linds、Encodings、Fields和Chroms分别存储多个Lind、Encoding、Field和Chrom。

PsyPopulation类的对象往往与带”psy”字样的进化算法模板配合使用，以实现多染色体混合编码的进化优化。

  1. 求解标准测试函数——McCormick函数

遗传算法求解以下函数的最小值: m i n f ( x , y ) = s i n ( x + y ) + ( x − y ) 2 − 1.5 x + 2.5 y + 1 minf(x, y) =sin(x+y) + (x−y)2−1.5x+ 2.5y+ 1 m i n f (x ,y )=s i n (x +y )+(x −y )2 −1 .5 x +2 .5 y +1

代码实现：


import numpy as np
import geatpy as ea
import time
"""============================目标函数============================"""
def aim(Phen):
    x1 = Phen[:, [0]]
    x2 = Phen[:, [1]]
    return np.sin(x1 + x2) + (x1 - x2) ** 2 - 1.5 * x1 + 2.5 * x2+1
"""============================变量设置============================"""
x1 = [-1.5, 4]
x2 = [-3, 4]
b1 = [1, 1]
b2 = [1, 1]

ranges=np.vstack([x1, x2]).T

borders=np.vstack([b1, b2]).T
varTypes = np.array([0, 0])
"""==========================染色体编码设置========================="""
Encoding ='BG'
codes = [1, 1]
precisions =[6, 6]
scales = [0, 0]
FieldD =ea.crtfld(Encoding,varTypes,ranges,borders,precisions,codes,scales)

"""=========================遗传算法参数设置========================"""
NIND     = 20
MAXGEN   = 100
maxormins = np.array([1])
selectStyle ='sus'
recStyle ='xovdp'
mutStyle ='mutbin'
Lind =int(np.sum(FieldD[0, :]))
pc= 0.9
pm= 1/Lind
obj_trace = np.zeros((MAXGEN, 2))
var_trace = np.zeros((MAXGEN, Lind))

"""=========================开始遗传算法进化========================"""
start_time = time.time()
Chrom = ea.crtpc(Encoding,NIND, FieldD)
variable = ea.bs2ri(Chrom, FieldD)
ObjV = aim(variable)
best_ind = np.argmin(ObjV)

for gen in range(MAXGEN):
    FitnV = ea.ranking(maxormins * ObjV)
    SelCh = Chrom[ea.selecting(selectStyle,FitnV,NIND-1),:]
    SelCh = ea.recombin(recStyle, SelCh, pc)
    SelCh = ea.mutate(mutStyle, Encoding, SelCh, pm)

    Chrom = np.vstack([Chrom[best_ind, :], SelCh])
    Phen = ea.bs2ri(Chrom, FieldD)
    ObjV = aim(Phen)

    best_ind = np.argmin(ObjV)
    obj_trace[gen,0]=np.sum(ObjV)/ObjV.shape[0]
    obj_trace[gen,1]=ObjV[best_ind]
    var_trace[gen,:]=Chrom[best_ind,:]

    end_time = time.time()
    ea.trcplot(obj_trace, [['种群个体平均目标函数值','种群最优个体目标函数值']])

"""============================输出结果============================"""
best_gen = np.argmin(obj_trace[:, [1]])
print('最优解的目标函数值：', obj_trace[best_gen, 1])
variable = ea.bs2ri(var_trace[[best_gen], :], FieldD)
print('最优解的决策变量值为：')
for i in range(variable.shape[1]):
    print('x'+str(i)+'=',variable[0, i])
    print('用时：', end_time - start_time,'秒')

效果图：

python遗传算法之geatpy学习
结果如下：
python遗传算法之geatpy学习

5.参考文章

链接: geatpy说明文档.

Original: https://blog.csdn.net/weixin_48077303/article/details/124079269
Author: 柳小葱
Title: python遗传算法之geatpy学习

原创文章受到原创版权保护。转载请注明出处：https://www.johngo689.com/673496/

转载文章受原作者版权保护。转载请注明原作者出处！

赞 (0)
0

【自取】最近整理的，有需要可以领取学习：



Linux核心资料大放送~

全栈面试题汇总（持续更新&可下载）

一个提高学习100%效率的工具！

【超详细】深度学习面试题目！

LeetCode Python刷题答案下载！

LeetCode Java版刷题答案下载！

LeetCode C++ 版本，抓紧保存！

LeetCode GO语言 刷题答案下载！

大家都在看

亲爱的 Coder【最近整理，可免费获取】👉 最新必读书单  | 👏 面试题下载  | 🌎 免费的AI知识星球