遗传算法优化BP神经网络_遗传算法权值权值-CSDN博客

本文链接：https://blog.csdn.net/u014271114/article/details/53445701

遗传算法原理

遗传算法（Genetic Algorithms）是1962年由美国Michigan大学Holland教授提出的模拟自然界遗传机制和生物进化论而形成的一种并行随机搜索最优化方法。它把自然界中的“优胜劣汰，适者生存”的生物进化原理引入优化参数形成的编码串联群体中，按照所选择的适应度函数并通过遗传中的选择、交叉和变异对个体进行筛选，是适应值好的个体被保留，适应度差的个体倍淘汰，新的种群集成上一代的信息，有由于上一代。这样反复循环，直到满足条件。遗传算法基本的操作分为：

选择操作
选择操作是指从旧的群体中以一定的概率选择个体到新的群体中，个体被选中的概率是跟适应度值有关，个体适应度值越好，被选中的概率越大。
交叉操作
交叉操作是指从个体中选择两个个体，通过两个染色体的狡猾组合，来查省心的优秀个体。交叉过程为从种群中任选两个染色体，随机选择一点或多点染色体未知进行交换。交叉操作如图所示：
变异操作
变异操作是指从种群中任选一个个体，选择染色体的一点进行变异以产生更优秀的个体。变异操作如图所示:

遗传算法具有高效启发式搜索、并行计算等特点，目前已经应用在函数优化、组合优化以及生产调度等方面。

遗传算法的基本要素

遗传算法的基本要素包括染色体编码方法、适应度函数、遗传操作和运行参数。

其中染色体编码方法是指个体的编码方法，目前主要包括二进制法、实数法等。二进制法是指把个体编码成一个二进制串，实数法是指把个体编码成一个实数串。

适应度函数是指根据进化目标编写的计算个体适应度值得函数，通过适应度函数计算每一个个体的适应度值，提供给选择算子进行选择。

遗传操作是指选择操作、交叉操作、和便已操作。

运行参数是遗传算法在初始化时确定的参数，主要包括群体大小 $M$ ，遗传代数 $G$ ，交叉概率 $P_c$ 和变异概率 $P_m$

算法流程

遗传算法优化BP神经网络分为BP神经网络结构确定、遗传算法优化和BP神经网络预测3个部分。其中，BP神经网络的结构确定部分根据拟合函数输入输出参数个数确定BP神经网络结构，进而确定遗传算法个体的长度。遗传算法优化使用遗传算法优化BP神经网络的权值和阈值，种群中的每个个体都包含了一个网络所有的权值和阈值，个体通过适应度函数计算个体的适应度值，遗传算法通过选择、交叉和变异操作找到最优适应度值对应个体。BP神经网络预测用遗传算法得到最优个体对网络初始权值和阈值赋值，网络经训练后预测函数输出。

算法流程

遗传算法的实现

种群初始化
个体编码方法为实数编码，每隔个体均为一个实数串，由输入层和隐含层连接权值、隐含层阈值、隐含层与输出层连接权值以及输出层阈值4个部分组成。个体包含了神经网络全部的权值和阈值，在网路结构一直的情况下，就可以构成一个结构、权值、阈值确定的神经网络。
适应度函数
根据个体得到BP神经网络的初始权值和阈值，用训练数据训练BP神经网络后预测系统输出，把预测输出个期望的输出之间的误差绝对值和 $E$ 作为个体适应度值F计算公式为
$F = k (\sum i = 1 n a s b (y i - o i))$ $F=k(\sum_{i=1}^n asb(y_i-o_i))$
式中，n为网络输出节点数； $y_i$ 为BP神经网络第i个节点的期望输出 $o_i$ 为第i个节点的实际输出；k为系数。
选择操作
遗传算法选择操作有轮盘赌法、锦标赛法等多种方法。选择轮盘赌法时，即基于适应度比例的选择策略，每个个体i的选择概率 $p_i$ 为：
$f i = k F i$ $f_i=\frac{k}{F_i}$ $p i = f i \sum N j = 1 f j$ $p_i=\frac{f_i}{\sum_{j=1}^N f_j}$
式中， $F_i$ 为个体i的适应度值，由于适应度值越小越好，所以在个体选择前对适应度求倒数；k为系数；N为种群个体数目。
交叉操作
由于个体采用实数编码，所以交叉操作方法采用实数交叉法，第k个染色体 $a_k$ 和第l个染色体a_l在j为的交叉操作方法如下：
$a k j = a k j (1 - b) + a l j b$ $a_{kj}=a_{kj}(1-b)+a_{lj}b$ $a l j = a l j (1 - b) + a k j b$ $a_{lj}=a_{lj}(1-b)+a_{kj}b$ 式中，b是[0,1]之间的随机数。
变异操作
选取第i个个体的第j个基因 $a_{ij}$ 进行变异变异操作方法如下：

$a i j = {a i j + (a i j - a m a x) * f (g) r > 0.5 a i j + (a m i n - a i j) * f (g) r < = 0.5$ $a_{ij}= \begin{cases} a_{ij}+(a_{ij}-a_{max})*f(g) \qquad r>0.5 \\ a_{ij}+(a_{min}-a_{ij})*f(g) \qquad r<=0.5 \\ \end{cases}$ 式中， $a_{max}$ 为基因 $a_{ij}$ 的上界； $a_{min}$ 为基因的下界； $f(g)=r_2(1-\frac{g}{G_{max}})^2$ $r_2$ 是一个随机数；g为当前迭代次数； $G_max$ 是最大的进化次数；r为[0,1]间的随机数。