数学建模十类经典算法（2）

百年孤独

2、最优化理论的三大非经典算法
/ _; J# e( s) L8 e& s+ g这十几年来最优化理论有了飞速发展，模拟退火法、神经网络、遗传算法这三类算法发展很快。近几年的赛题越来越复杂，很多问题没有什么很好的模型可以借鉴，于是这三类算法很多时候可以派上用场，比如：97 年A 题的模拟退火算法，00 年B 题的神经网络分类算法，象01 年B 题这种难题也可以使用神经网络，还有美国竞赛89 年A 题也和BP 算法有关系，当时是86 年刚提出BP 算法，89 年就考了，说明赛题可能是当今前沿科技的抽象体现。目前算法最佳的是遗传算法。
  Q9 i  w. n! I* ?" ~2 s$ l# e
4 Z$ Q& [: }! ]" w/ G! I; x遗传算法简介：
  R: y2 D8 \5 ]* O9 P- N遗传算法是一类借鉴生物界自然选择和自然遗传机制的随机化搜索算法，由美国J.Holland教授提出，其主要特点是群体搜索策略和群体中个体之间的信息交换，搜索不依赖于梯度信息。它尤其适用于传统搜索方法难于解决的复杂和非线性问题，可广泛用于组合优化、机器学习、自适应控制、规划设计和人工生命等领域，是21世纪有关智能计算中的关键技术之一。
- L! K4 A# M4 j$ P) O- a$ D' p3 ~在人工智能领域中，有不少问题需要在复杂和庞大的搜索空间中寻找最优解或准最优解。象货郎担问题和规划问题等组合优化问题就是典型的例子。在求解此类问题时，若不能利用问题固有知识来缩小搜索空间则会产生搜索的组合爆炸。
1 t) U9 {: B, H9 A; }  }! p
因此，研究能在搜索过程中自动获取和积累有关搜索空间的知识，并自适应地控制搜索过程，从而得到最优解地通用搜索方法一直是令人瞩目地课题。遗传算法就是这种特别有效地算法。生物的进化是一个奇妙的优化过程，它通过选择淘汰，突然变异，基因遗传等规律产生适应环境变化的优良物种。遗传算法是根据生物进化思想而启发得出的一种全局优化算法。尽管遗传算法本身在理论和应用方法上仍有许多待进一步研究地问题，但它已在很多领域地应用中展现了其特色和魅力。
+ u4 Z, \7 u+ t! m7 T& b8 k遗传算法的基本概念
遗传算法的基本思想是基于Darwin进化论和Mendel的遗传学说的。
; Z% X8 J& I5 @* s  A1 K5 x0 GDarwin进化论最重要的是适者生存原理。它认为每一物种在发展中越来越适应环境。物种每个个体的基本特征由后代所继承，但后代又会产生一些异于父代的新变化。在环境变化时，只有那些能适应环境的个体特征方能保留下来。
/ H" u0 X3 l+ g2 t, n. H8 DMendel遗传学说最重要的是基因遗传原理。它认为遗传以密码方式存在细胞中，并以基因形式包含在染色体内。每个基因有特殊的位置并控制某种特殊性质；所以，每个基因产生的个体对环境具有某种适应性。基因突变和基因杂交可产生更适应于环境的后代。经过存优去劣的自然淘汰，适应性高的基因结构得以保存下来。
由于遗传算法是由进化论和遗传学机理而产生的直接搜索优化方法；故而在这个算法中要用到各种进化和遗传学的概念。这些概念如下：
' Y. \/ ]3 k+ ^* y- ?0 _! o一、串(String)
它是个体(Individual)的形式，在算法中为二进制串，并且对应于遗传学中的染色体(Chromosome)。
# x% k: j3 f; l* u" _二、群体(Population)
个体的**称为群体，串是群体的元素
三、群体大小(Population Size)
6 L" v# ?: d; L" B8 P! U在群体中个体的数量称为群体的大小。
! j+ W3 n5 C/ V3 q1 O4 }四、基因(Gene)
基因是串中的元素，基因用于表示个体的特征。例如有一个串S＝1011，则其中的1，0，1，1这4个元素分别称为基因。它们的值称为等位基因(Alletes)。
+ G- V) a2 o/ ~! f五 、基因位置(Gene Position)
+ x* g7 P" Y0 J" F一个基因在串中的位置称为基因位置，有时也简称基因位。基因位置由串的左向右计算，例如在串S＝1101中，0的基因位置是3。基因位置对应于遗传学中的地点(Locus)。
; p( K: R& H/ W, J4 T2 h
+ U6 `. J5 q4 z  V六、基因特征值(Gene Feature)
在用串表示整数时，基因的特征值与二进制数的权一致；例如在串S=1011中，基因位置3中的1，它的基因特征值为2；基因位置1中的1，它的基因特征值为8。
1 Z5 E0 X6 _" H7 R: _七、串结构空间SS
在串中，基因任意组合所构成的串的**。基因操作是在结构空间中进行的。串结构空间对应于遗传学中的基因型(Genotype)的**。
八、参数空间SP
( x' d# h7 X; h+ C* [+ Y9 q这是串空间在物理系统中的映射，它对应于遗传学中的表现型(Phenotype)的**。
% \, C# _7 y( X* l1 Q九、非线性
8 X1 R( `: w& d0 u5 @  w它对应遗传学中的异位显性(Epistasis)
* r% x& x9 U8 r# b3 q- L十、适应度(Fitness)
表示某一个体对于环境的适应程度。遗传算法的原理
) E, k) d# ]; V- u4 V1 O遗传算法GA把问题的解表示成“染色体”，在算法中也即是以二进制编码的串。并且，在执行遗传算法之前，给出一群“染色体”，也即是假设解。然后，把这些假设解置于问题的“环境”中，并按适者生存的原则，从中选择出较适应环境的“染色体”进行复制，再通过交叉，变异过程产生更适应环境的新一代“染色体”群。这样，一代一代地进化，最后就会收敛到最适应环境的一个“染色体”上，它就是问题的最优解。
一、遗传算法的目的
典型的遗传算法CGA(Canonical Genetic Algorithm)通常用于解决下面这一类的静态最优化问题：
0 y1 P" ?# Q9 k( ^* A2 M2 |) }8 I考虑对于一群长度为L的二进制编码bi，i＝1，2，…，n；有
bi∈{0,1}
& `0 U, J- d7 [给定目标函数f，有f(bi)，并且
0<f(bi)<∞
" S0 G( Q; o! l) C  {4 T3 t同时
f(bi)≠f(bi+1)
求满足下式
( F  a! O: \2 _4 ~' xmax{f(bi)|bi∈{0,1}的bi。
很明显，遗传算法是一种最优化方法，它通过进化和遗传机理，从给出的原始解群中，不断进化产生新的解，最后收敛到一个特定的串bi处，即求出最优解。
3 U* j4 h  L2 l+ U$ r  |: p
二、遗传算法的基本原理
1 w6 ^+ z# E4 d0 M长度为L的n个二进制串bi(i＝1，2，…，n)组成了遗传算法的初解群，也称为初始群体。在每个串中，每个二进制位就是个体染色体的基因。根据进化术语，对群体执行的操作有三种：
1．选择(Selection)
1 {8 K& Q+ w' U" ?这是从群体中选择出较适应环境的个体。这些选中的个体用于繁殖下一代。故有时也称这一操作为再生(Reproduction)。由于在选择用于繁殖下一代的个体时，是根据个体对环境的适应度而决定其繁殖量的，故而有时也称为非均匀再生(differential reproduction)。
2．交叉(Crossover)
& W6 c9 U0 b; h这是在选中用于繁殖下一代的个体中，对两个不同的个体的相同位置的基因进行交换，从而产生新的个体。
3．变异(Mutation)
这是在选中的个体中，对个体中的某些基因执行异向转化。在串bi中，如果某位基因为1，产生变异时就是把它变成0；反亦反之。三、遗传算法的步骤
1．初始化
选择一个群体，即选择一个串或个体的**bi，i=1，2，...n。这个初始的群体也就是问题假设解的**。一般取n＝30-160。
通常以随机方法产生串或个体的**bi,i＝1，2，...n。问题的最优解将通过这些初始假设解进化而求出。
2．选择
根据适者生存原则选择下一代的个体。在选择时，以适应度为选择原则。适应度准则体现了适者生存，不适应者淘汰的自然法则。
6 R& u2 m  i7 K; \8 p给出目标函数f，则f(bi)称为个体bi的适应度。以
% Y9 \# H: h# e- j
: P& _7 M6 }9 G; H6 D% q/ l
, r* }" H5 g/ k为选中bi为下一代个体的次数。

显然：
(1)适应度较高的个体，繁殖下一代的数目较多。
(2)适应度较小的个体，繁殖下一代的数目较少；甚至被淘汰。
" y7 u* _6 P2 }/ [1 v! }5 O! z  \这样，就产生了对环境适应能力较强的后代。对于问题求解角度来讲，就是选择出和最优解较接近的中间解。
0 C' M- }/ i  F3 ^选择的方法有：
适应度比例法
$ w6 p( A$ {( C期望值法
排位次法
精华保存法

1 ^! L& F8 X) H: d; a; @6 j1 h/ R" f3．交叉
: W9 ?: D( O& S7 a对于选中用于繁殖下一代的个体，随机地选择两个个体的相同位置，按交叉概率P。在选中的位置实行交换。这个过程反映了随机信息交换；目的在于产生新的基因组合，也即产生新的个体。交叉时，可实行单点交叉或多点交叉。
1 m7 s' t7 w& I+ _& s9 d0 [; d
+ K, [+ l# W/ p9 f例如：有个体
. k- B( R0 n: h, r( qS1=100101
S2=010111
, a: S+ T4 o1 W2 ^6 Z# u2 y选择它们的左边3位进行交叉操作，则有
S1=010101
. f9 T% x; Z1 K- V) IS2=100111
! H& n( R/ ^  E; o一般而言，交叉概率P，取值为0.25—0.75。4．变异
根据生物遗传中基因变异的原理，以变异概率Pm对某些个体的某些位执行变异。在变异时，对执行变异的串的对应位求反，即把1变为0，把0变为1。变异概率Pm与生物变异极小的情况一致，所以，Pm的取值较小，一般取0.01-0.2。
- I  m. }) I. h0 K: ?
例如：
有个体S＝101011。
' G7 g, _0 _' |4 F1 h! m' G对其的第1，4位置的基因进行变异，则有
5 B0 Y( ?: B* h+ [2 }S'=001111
6 N1 s8 @) I1 z6 @& W& k" O& @单靠变异不能在求解中得到好处。但是，它能保证算法过程不会产生无法进化的单一群体。因为在所有的个体一样时，交叉是无法产生新的个体的，这时只能靠变异产生新的个体。也就是说，变异增加了全局优化的特质。
' Y6 G" K& {. [+ L5．全局最优收敛(Convergence to the global optimum)
当最优个体的适应度达到给定的阀值，或者最优个体的适应度和群体适应度不再上升时，则算法的迭代过程收敛、算法结束。否则，用经过选择、交叉、变异所得到的新一代群体取代上一代群体，并返回到第2步即选择操作处继续循环执行。
遗传算法基本处理流程图如下：

二、遗传算法的应用关键
遗传算法在应用中最关键的问题有如下3个
8 f; D* [5 c: ^2 A2 ]1．串的编码方式
$ f% Z9 H5 C, c) y4 t6 O6 q这本质是问题编码。一般把问题的各种参数用二进制编码，构成子串；然后把子串拼接构成“染色体”串。串长度及编码形式对算法收敛影响极大。
* c  O4 C, E; t" R; n3 D0 f2．适应函数的确定
! x' @, @( a" X" A) c& G. ~适应函数(fitness function)也称对象函数(object function)，这是问题求解品质的测量函数；往往也称为问题的“环境”。一般可以把问题的模型函数作为对象函数；但有时需要另行构造。
% Y' {* \2 P0 P5 ]! b3．遗传算法自身参数设定
5 _' o, ]! t7 O6 j1 @遗传算法自身参数有3个，即群体大小n、交叉概率Pc和变异概率Pm。
群体大小n太小时难以求出最优解，太大则增长收敛时间。一般n＝30-160。交叉概率Pc太小时难以向前搜索，太大则容易破坏高适应值的结构。一般取Pc=0.25-0.75。变异概率Pm太小时难以产生新的基因结构，太大使遗传算法成了单纯的随机搜索。一般取Pm＝0．01—0．2。
" V! w' w6 X3 W4 c5 s0 qmatlab遗传算法工具箱函数及实例讲解
$ L% u  @" x! \3 B; f核心函数：
(1)function [pop]=initializega(num,bounds,eevalFN,eevalOps,options)--初始种群的生成函数
【输出参数】
pop--生成的初始种群
" p* Q1 {3 x1 t【输入参数】
num--种群中的个体数目
* G' a4 d# e. S) Ebounds--代表变量的上下界的矩阵
: i# K- H& D8 |7 n3 V0 beevalFN--适应度函数
9 v& X) {) H  qeevalOps--传递给适应度函数的参数
options--选择编码形式(浮点编码或是二进制编码)[precision F_or_B],如
9 O+ J" M4 K& A! q8 f# Qprecision--变量进行二进制编码时指定的精度
0 `( h! ^- J' m1 hF_or_B--为1时选择浮点编码，否则为二进制编码,由precision指定精度)
2)function [x,endPop,bPop,traceInfo] = ga(bounds,evalFN,evalOps,startPop,opts,...
% u7 g! y: L4 Q. G0 ?termFN,termOps,selectFN,selectOps,xOverFNs,xOverOps,mutFNs,mutOps)--遗传算法函数
2 S' E+ c, d- Q* R4 r+ B* S【输出参数】
6 E" t; J8 O1 bx--求得的最优解
endPop--最终得到的种群
bPop--最优种群的一个搜索轨迹
4 D4 x) ~* L; p1 s" `【输入参数】
/ v4 Y* K. S$ n1 y, j) k! ebounds--代表变量上下界的矩阵
2 k) y, M5 O& S, W( s  WevalFN--适应度函数
& O4 ?7 x( `& jevalOps--传递给适应度函数的参数
% \' _5 ^: ^9 E1 o4 l2 MstartPop-初始种群
opts[epsilon prob_ops display]--opts(1:2)等同于initializega的options参数，第三个参数控制是否输出，一般为0。如[1e-6 1 0]
termFN--终止函数的名称,如[‘maxGenTerm’]
termOps--传递给终止函数的参数,如[100]
9 Q$ I, F' _/ L* w, CselectFN--选择函数的名称,如[‘normGeomSelect’]
selectOps--传递给选择函数的参数,如[0.08]
xOverFNs--交叉函数名称表，以空格分开，如['arithXover heuristicXover simpleXover']
xOverOps--传递给交叉函数的参数表，如[2 0;2 3;2 0]
9 R; a4 j5 `  b% L; _$ Z( r7 q  MmutFNs--变异函数表，如['boundaryMutation multiNonUnifMutation nonUnifMutation unifMutation']
' Y) b* f7 h9 h( \mutOps--传递给交叉函数的参数表,如[4 0 0;6 100 3;4 100 3;4 0 0]
( l4 W$ `* q( B( r% B【问题】求f(x)=x+10*sin(5x)+7*cos(4x)的最大值，其中0<=x<=9
+ K& ]7 n: T' `1 G4 z" h. j& E9 b7 t3 b【分析】选择二进制编码，种群中的个体数目为10，二进制编码长度为20，交叉概率为0.95,变异概率为0.08
【程序清单】
%编写目标函数
function[sol,eval]=fitness(sol,options)
' R6 n0 a& t7 Ox=sol(1);
  ~+ h8 I5 {, E% c+ d) t( k2 n/ G, Zeval=x+10*sin(5*x)+7*cos(4*x);
%把上述函数存储为fitness.m文件并放在工作目录下
' x# m/ D  s' \4 d7 s& [
. H9 M" T, F7 [' M3 ZinitPop=initializega(10,[0 9],'fitness');%生成初始种群，大小为10
[x endPop,bPop,trace]=ga([0 9],'fitness',[],initPop,[1e-6 1 1],'maxGenTerm',25,'normGeomSelect',...
[0.08],['arithXover'],[2],'nonUnifMutation',[2 25 3]) %25次遗传迭代
运算结果为：x = 7.8562 24.8553(当x为7.8562时，f（x）取最大值24.8553)

5 V& N$ G8 C5 E+ q: U注：1、遗传算法一般用来取得近似最优解，而不是最优解。
; v6 u( T5 ~' Z7 d  C2、matlab工具箱函数必须放在工作目录下