本帖最后由 梦里花111 于 2013-9-25 22:28 编辑 1 L& I* W# k" p" V9 q+ Q' a' }1 R! Y
在工程实践中,经常会接触到一些比较“新颖”的算法或理论,比如:等。这些算法或理论都有一些共同的特性(比如模拟自然过) s7 A1 o; n* j! O
程),通称为“智能算法”。它们在解决一些复杂的工程问题时大有用武之地。 7 q- n# [7 Q* r7 Q, l1 N( ~! K & A4 O! x# W$ E& [( Y3 _
, C% W }% A L0 {, B5 {! h 这些算法都有什么含义?首先给出个局部搜索,模拟退火,遗传算法,禁忌搜索的形象比喻: 6 g. _2 ~9 |- k5 J - T* ?( R; z/ c9 z N
/ \5 }' \( Z; T+ T 为了找出地球上最高的山,一群有志气的兔子们开始想办法。 ' Y7 t6 v3 `4 Z: ^' h 1 E. x* Z7 O6 p0 [+ X
$ O) }* O6 s& y6 ~4 O# c 1.兔子朝着比现在高的地方跳去。他们找到了不远处的最高山峰。但是这座山不一定是珠穆朗玛峰。这就是局部搜索,它不能保证局部最优值就是全局最优值。9 s+ W* x/ o! v
. e+ d: j$ l' I9 m A 2 E2 f8 D2 F% a8 h; {
2.兔子喝醉了。他随机地跳了很长时间。这期间,它可能走向高处,也可能踏入平地。但是,他渐渐清醒了并朝最高方向跳去。这就是模拟退火。( Y( H5 n" K, G9 `' f& k% X
7 e6 F0 }- _( P% b1 _/ [% ?2 K, o ' t0 Y% s- m5 I
3.兔子们吃了失忆药片,并被发射到太空,然后随机落到了地球上的某些地方。他们不知道自己的使命是什么。但是,如果你过几年就杀死一部分海拔低的兔子,多产的' Q( L# p: D9 j8 t; c. L3 i
兔子们自己就会找到珠穆朗玛峰。这就是遗传算法。. @" K- O# F' Q$ l
% h" [6 @6 A' w9 q$ d: u
4 Y6 d: w. E% M, f# X
4.兔子们知道一个兔的力量是渺小的。他们互相转告着,哪里的山已经找过,并且找过的每一座山他们都留下一只兔子做记号。他们制定了下一步去哪里寻找的策略。这 4 D2 w' l& M* F; c& l3 V就是禁忌搜索。 1 e' d; b& x Z+ ^ P y- b " U4 Q- W- G w 2 ]$ h2 ~6 Y/ P6 a
智能优化算法的概述 2 Z- s; x. \) E# n : u6 x( u! t3 ~9 m* {# O % J/ N3 P5 @# l6 Q) B
智能优化算法要解决的一般是最优化问题。最优化问题可以分为(1)求解一个函数中,使得函数值最小的自变量取值的函数优化问题和(2)在一个解空间里面,寻找最 : T% r4 w0 y, l4 J! D: F6 ^优解,使目标函数值最小的组合优化问题。典型的组合优化问题有:旅行商问题(Traveling - P; t/ @0 C3 h; p F; z' Q8 x
Salesman Problem,TSP),加工调度问题(Scheduling 1 `" |* g5 D3 ~- P4 W% k: V) n$ { Problem),0-1背包问题(Knapsack " G. `$ ~2 r, k! }3 c Problem),以及装箱问题(Bin Packing Problem)等。 / M+ S8 d* q8 A0 O F: s) j 优化算法有很多,经典算法包括:有线性规划,动态规划等;改进型局部搜索算法包括爬山法,最速下降法等,本文介绍的模拟退火、遗传算法以及禁忌搜索称作指导性搜4 c0 H* t; X. E
索法。而神经网络,混沌搜索则属于系统动态演化方法。# v' G* L& p4 |
' N8 `) O4 y/ W- n: J4 f 3 N3 t/ O% t, F5 }: ?$ d; s( I 优化思想里面经常提到邻域函数,它的作用是指出如何由当前解得到一个(组)新解。其具体实现方式要根据具体问题分析来定。8 B! h1 E: E. @& F
! w& ?% s4 c0 R$ C* x: I + {7 C& m+ X" y( l
一般而言,局部搜索就是基于贪婪思想利用邻域函数进行搜索,若找到一个比现有值更优的解就弃前者而取后者。但是,它一般只可以得到“局部极小解”,就是说,可能; q, e. C5 u6 ~* W
这只兔子登“登泰山而小天下”,但是却没有找到珠穆朗玛峰。而模拟退火,遗传算法,禁忌搜索,神经网络等从不同的角度和策略实现了改进,取得较好的“全局最小解 ”。 # A( m$ }* _- u: i! n( a; D8 Z / b5 H1 q8 l0 Y- V
" n5 W% l/ a. M. o, i( R) a6 E
模拟退火算法(Simulated Annealing,SA)# a6 [% n. W7 }2 S5 o P, E
$ k( J# Y8 n& t2 j$ H + V% R! k' V2 [; O0 Y; o$ w
模拟退火算法的依据是固体物质退火过程和组合优化问题之间的相似性。物质在加热的时候,粒子间的布朗运动增强,到达一定强度后,固体物质转化为液态,这个时候再5 p8 R, A+ j6 s B
进行退火,粒子热运动减弱,并逐渐趋于有序,最后达到稳定。 ( X5 H. K" ^9 C ' A& t6 w, M# M9 M ^ & i: y% j# ^0 l% M 模拟退火的解不再像局部搜索那样最后的结果依赖初始点。它引入了一个接受概率p。如果新的点(设为pn)的目标函数f(pn)更好,则p=1,表示选取新点;否* o" h! k, @% z
则,接受概率p是当前点(设为pc)的目标函数f(pc),新点的目标函数f(pn)以及另一个控制参数“温度”T的函数。也就是说,模拟退火没有像局部搜索那* k1 a% P7 H# b7 o
样每次都贪婪地寻找比现在好的点,目标函数差一点的点也有可能接受进来。随着算法的执行,系统温度T逐渐降低,最后终止于某个低温,在该温度下,系统不再接受变( X9 {- E& k# Y- B+ I
化。 0 v1 E8 x- `& B; [ - x$ d+ H5 \* g9 u y7 t7 k
1 ^4 {6 L) m, N/ p* M: G9 F# u& M
模拟退火的典型特征是除了接受目标函数的改进外,还接受一个衰减极限,当T较大时,接受较大的衰减,当T逐渐变小时,接受较小的衰减,当T为0时,就不再接受衰 " G/ ~' W/ P4 @" i; Z! ^, x减。这一特征意味着模拟退火与局部搜索相反,它能避开局部极小,并且还保持了局部搜索的通用性和简单性。 / p4 o$ ]$ ^" J% \/ \: k 在物理上,先加热,让分子间互相碰撞,变成无序状态,内能加大,然后降温,最后的分子次序反而会更有序,内能比没有加热前更小。就像那只兔子,它喝醉后,对比较* v9 A1 P. h; d0 q- W% C
近的山峰视而不见,迷迷糊糊地跳一大圈子,反而更有可能找到珠峰。 " x$ \" i4 f, p) a' s& r $ N. |& ]+ e& I, r' V$ W0 B+ S
$ [, u: U3 ]6 | P$ [
值得注意的是,当T为0时,模拟退火就成为局部搜索的一个特例。 3 u! ]9 J! u. m5 }: M: f 6 b" z; }$ W1 W; i: E4 `+ r1 u! [
' s, l2 x! T1 e) g- r& N
模拟退火的伪码表达: : Y( s" C; v3 E" N8 N procedure simulated annealing 8 O. v4 y1 G1 N4 N begin F$ m& b9 j( [7 h. b t:=0; " L' l. M" c; T* W
initialize temperature T # n) ^+ A1 y+ ^2 p9 r. L$ ^ select a current string vc at random; # C. o8 x& F, k evaluate vc; . k9 J" p* m- D! X5 D, l1 K
repeat 3 S0 \# l Z' g8 [: M0 k repeat + G3 w, W: P, S' o select a new string vn in the neighborhood of vc; (1) : V' r4 g3 G3 u0 u' A4 W* ^
if f(vc) then vc:=vn; 1 U# n1 t+ w/ i0 [: l: P: g- x4 F
else if random [0,1] then vc:=vn; ) W; {" T1 Z6 V/ G
until (termination-condition) (3) . {* E. N2 r9 L5 E) U- P
T:=g(T,t); (4) 0 y3 P5 i2 i0 Y+ B T:=t+1; 4 F) ?) |/ ~+ v% w: Z; Q until (stop-criterion) (5) * q3 F2 t: O3 q2 `) q9 A
end; : G( S( h+ k$ F4 }: I! t0 q / b% \! ?8 x1 L6 D; }
1 u0 M* V6 W7 A; H. I
上面的程序中,关键的是(1)新状态产生函数,(2)新状态接受函数,(3)抽样稳定准则,(4)退温函数,(5)退火结束准则 3 T* J/ S8 }* y0 Q% m' b' a (简称三函数两准则)是直接影响优化结果的主要环节。虽然实验结果证明初始值对于最后的结果没有影响,但是初温越高,得到高质量解的概率越大。所以,应该尽量选 7 `( F$ R5 y. X1 N, F N8 {取比较高的初温。7 T1 j3 m( l! Y/ B) S, b% I. q) m
6 f5 {7 l7 v6 h4 ^, d 3 m8 D' w- c' S& G$ l; w6 u4 k
上面关键环节的选取策略:/ s% u* V: ?- Y. T# ~
(1)状态产生函数:候选解由当前解的邻域函数决定,可以取互换,插入,逆序等操作产生,然后根据概率分布方式选取新的解,概率可以取均匀分布、正态分布、高斯' h1 H* M1 L# I, b
分布、柯西分布等。 & ?- O" K( c" x7 o; A+ w& a + t, }* e: j5 g. P8 k6 v/ p1 K ' s5 f. L! S& T7 ?6 C2 g
(2)状态接受函数:这个环节最关键,但是,实验表明,何种接受函数对于最后结果影响不大。所以,一般选取min ) i; k O& Z# b- \0 B
[1, exp ((f (vn)-f (vc))/T)]。 " r9 `; T) }5 } (3)抽样稳定准则:一般常用的有:检验目标函数的均值是否稳定;连续若干步的目标值变化较小;规定一定的步数; 7 A8 J1 ~ p5 ` . m% n0 i3 ~* ]" Y+ L& [; a
8 e M1 J$ T9 Q (4)退温函数:如果要求温度必须按照一定的比率下降,SA算法可以采用 + B" k. _0 V, w6 ?4 x ,但是温度下降很慢;快速SA中,一般采用 ! s# w8 S! q' V/ ]0 |2 H, O 。目前,经常用的是 , 是一个不断变化的值。* R; [5 S+ Z: e$ V& w
(5)退火结束准则:一般有:设置终止温度;设置迭代次数;搜索到的最优值连续多次保持不变;检验系统熵是否稳定。; N u* N* J3 S( Y7 i$ K$ ?
% n: S/ |2 u1 H# G ) |% ^8 O2 G3 {. _% h k
遗传算法(Genetic Algorithm, GA) ( _" P" U+ D( k! |/ ?/ d. y4 j & C/ _8 S% k1 x0 p2 [/ w3 }
7 h3 m% r0 ~7 G) [* C2 J6 l3 W
“物竞天择,适者生存”,是进化论的基本思想。遗传算法就是模拟自然界想做的事。遗传算法可以很好地用于优化问题,若把它看作对自然过程高度理想化的模拟,更能) p i9 s4 r/ g* V. r
显出它本身的优雅——虽然生存竞争是残酷的。( w+ E+ Y5 z8 |3 G
% p! I) Y& C1 o
: ~( [% y: S$ ^ ?* U% I 遗传算法以一种群体中的所有个体为对象,并利用随机化技术指导对一个被编码的参数空间进行高效搜索。其中,选择、交叉和变异构成了遗传算法的遗传操作;参数编码 _, c% B' D" P& ?& ^1 f2 X、初始群体的设定、适应度函数的设计、遗传操作设计、控制参数设定五个要素组成了遗传算法的核心内容。, |+ G8 h6 l/ U% a+ M/ u3 i
作为一种新的全局优化搜索算法,遗传算法以其简单通用、健壮性强、适于并行处理以及高效、实用等显著特点,在各个领域得到了广泛应用,取得了良好效果,并逐渐成 % f& Z4 n \% {2 U为重要的智能算法之一。 ' s. ^3 @+ i4 y: i7 W% `2 ` ( |; ~4 t, g+ N% ^8 U6 y5 N! M 5 y& `) F2 ~9 U* \
遗传算法的伪码:! T$ A, d/ m. U. y( s! s X2 X
5 `7 Z$ F3 s9 \" b# B 6 b( ^; q! V+ M- G$ Q: G
procedure genetic algorithm 8 e$ h k+ v1 Y6 x1 C begin 8 {. E5 e' R+ Q+ I( w- C$ g
initialize a group and evaluate the fitness value ; (1) - Q: ]! I5 I7 b2 o while not convergent (2) + P+ K c( `5 y7 A$ S4 Y) i begin : c; G. x _; W- k* {! D
select; (3)0 x* I# ~* @; I0 |6 L
if random[0,1] crossover; (4) 4 \0 P) x' w3 H1 S- t, ` if random (0,1) mutation; (5) + A5 }6 W/ N/ f9 L3 d+ R" R end; 3 D2 |/ a8 q; c- K1 }
end 9 n1 W; @5 K9 h
上述程序中有五个重要的环节:: t1 `6 U# ]0 t, Q, Z, A! w
(1)编码和初始群体的生成:GA在进行搜索之前先将解空间的解数据表示成遗传空间的基因型串结构数据,这些串结构数据的不同组合便构成了不同的点。然后随机产 0 D5 d3 _5 e3 m生N个初始串结构数据,每个串结构数据称为一个个体,$ m% c. @9 {- o* r! h6 q7 c3 ~- s
N个体构成了一个群体。GA以这N个串结构数据作为初始点开始迭代。) h2 X/ J; C F* t$ x
$ }1 q: T' w4 V1 o5 S; ]
. \# C! h- ~5 L8 P V. E 比如,旅行商问题中,可以把商人走过的路径进行编码,也可以对整个图矩阵进行编码。编码方式依赖于问题怎样描述比较好解决。初始群体也应该选取适当,如果选取的& d) H, Y7 z& i8 J' `; {
过小则杂交优势不明显,算法性能很差(数量上占了优势的老鼠进化能力比老虎强),群体选取太大则计算量太大。7 ?+ H( P8 N/ x/ K& J: ~+ i
, o2 r' D- E- `
7 X1 `. [; c- [, g7 O0 ?+ p (2)检查算法收敛准则是否满足,控制算法是否结束。可以采用判断与最优解的适配度或者定一个迭代次数来达到。0 O9 ^, B* @$ E. e/ X" K6 F: T9 l
' B0 o5 z8 Y7 p 5 u; W( |! K4 b2 q* u$ h/ E! c9 @, m (3)适应性值评估检测和选择:适应性函数表明个体或解的优劣性,在程序的开始也应该评价适应性,以便和以后的做比较。不同的问题,适应性函数的定义方式也不同: w) k6 q3 `2 h+ C+ Z
。根据适应性的好坏,进行选择。选择的目的是为了从当前群体中选出优良的个体,使它们有机会作为父代为下一代繁殖子孙。遗传算法通过选择过程体现这一思想,进行 9 w& |$ }/ ~$ c% a" r选择的原则是适应性强的个体为下一代贡献一个或多个后代的概率大。选择实现了达尔文的适者生存原则。 ! t D* a, `3 t0 x 5 N. |1 o6 l1 g) p ; j. h; {1 H7 ^7 C7 R) G0 h0 B; g
(4)杂交:按照杂交概率(pc)进行杂交。杂交操作是遗传算法中最主要的遗传操作。通过杂交操作可以得到新一代个体,新个体组合了其父辈个体的特性。杂交体现; ?' N* }& f+ L; t' D6 r
了信息交换的思想。' n4 r/ C% R! A
/ A3 H% {# S4 X0 Z
0 [8 @: n7 O$ O! W* p% n" u5 N 可以选定一个点对染色体串进行互换,插入,逆序等杂交,也可以随机选取几个点杂交。杂交概率如果太大,种群更新快,但是高适应性的个体很容易被淹没,概率小了搜 , w, L' [: ^5 d/ R5 d索会停滞。 % h" O3 e; p' O; ~ D H + e+ j5 x( I9 r* U4 v o % u# b4 p0 }* B: F5 v9 D& L E (5)变异:按照变异概率(pm)进行变异。变异首先在群体中随机选择一个个体,对于选中的个体以一定的概率随机地改变串结构数据中某个串的值。同生物界一样, GA中变异发生的概率很低。变异为新个体的产生提供了机会。) k& U' ?7 F3 h, j
0 x/ }4 m! e' S5 L0 d
7 o& \1 g7 B7 ^& X 变异可以防止有效基因的缺损造成的进化停滞。比较低的变异概率就已经可以让基因不断变更,太大了会陷入随机搜索。想一下,生物界每一代都和上一代差距很大,会是; T/ V( z& d8 c( x$ h' |
怎样的可怕情形。 0 T/ X- }5 N R0 W 4 `& n" K) a* r- \% ~ Y! f 4 y7 v4 s9 E! w. O% \. [
就像自然界的变异适和任何物种一样,对变量进行了编码的遗传算法没有考虑函数本身是否可导,是否连续等性质,所以适用性很强;并且,它开始就对一个种群进行操作 1 v" |2 A2 j2 w; x) x,隐含了并行性,也容易找到“全局最优解”。 % X. X5 V7 m. U- M6 i p% ^ " R; \+ Q4 i3 ~& H0 ^" {3 r 3 J7 W7 }+ i+ S/ O/ n 禁忌搜索算法(Tabu Search,TS) : K& O4 t. X) O5 z8 M( T1 d - r' Z# h; S7 F8 @
* N% [7 ?# t: Y9 D
为了找到“全局最优解”,就不应该执着于某一个特定的区域。局部搜索的缺点就是太贪婪地对某一个局部区域以及其邻域搜索,导致一叶障目,不见泰山。禁忌搜索就是 9 `+ v3 N: e2 c9 @4 @* Y' k对于找到的一部分局部最优解,有意识地避开它(但不是完全隔绝),从而获得更多的搜索区间。兔子们找到了泰山,它们之中的一只就会留守在这里,其他的再去别的地 6 B6 z; k1 f' U' j8 L. j方寻找。就这样,一大圈后,把找到的几个山峰一比较,珠穆朗玛峰脱颖而出。 x: h; X v/ j
% ~2 f) ]1 M) {/ B
. R# q; K( U; t: ] 当兔子们再寻找的时候,一般地会有意识地避开泰山,因为他们知道,这里已经找过,并且有一只兔子在那里看着了。这就是禁忌搜索中“禁忌表(tabu 3 o6 M3 ~9 b# D d/ _
list)”的含义。那只留在泰山的兔子一般不会就安家在那里了,它会在一定时间后重新回到找最高峰的大军,因为这个时候已经有了许多新的消息,泰山毕竟也有一 6 H5 Q8 ?& A' w$ S" ^6 z个不错的高度,需要重新考虑,这个归队时间,在禁忌搜索里面叫做“禁忌长度(tabu 6 M8 M" y- e2 [$ k" a
length)”;如果在搜索的过程中,留守泰山的兔子还没有归队,但是找到的地方全是华北平原等比较低的地方,兔子们就不得不再次考虑选中泰山,也就是说,当 o: a4 U+ s* v X8 ` J一个有兔子留守的地方优越性太突出,超过了“best # _& F/ v, F) k7 }( L7 Q to % w2 M! P ~4 h! ?
far”的状态,就可以不顾及有没有兔子留守,都把这个地方考虑进来,这就叫“特赦准则(aspiration * Y' Y8 n4 U$ g, `# m
criterion)”。这三个概念是禁忌搜索和一般搜索准则最不同的地方,算法的优化也关键在这里。 / D' t% {0 _% W6 I4 i & k+ l5 G- U0 Q2 f4 g9 k+ M
" O$ q7 g$ Q/ C4 u 伪码表达: % y0 m' k+ L }0 I. E6 x! | procedure tabu search; 7 s) g8 u& V# A2 Y# ?$ |
begin 0 Y* P; N7 l' i+ i' T initialize a string vc at random,clear up the tabu list; 8 R/ n- S/ A& ?% g1 ]
cur:=vc; 1 x/ E$ _' m" V% Q
repeat , X4 z8 Q# @, _2 k0 A! n select a new string vn in the neighborhood of vc; $ y) c- p4 p- _
if va>best_to_far then {va is a string in the tabu list} / _: a6 m. |' E/ N: W: D" X% u) q9 V0 a begin : s" _3 M# j. A' o0 v9 @ cur:=va; 9 e: @# e; s% Q* @& J
let va take place of the oldest string in the tabu list; 1 X) ~, o6 J z+ z. B) y best_to_far:=va; : ?9 }; l; T$ b# s% ]1 K$ R end else 9 o4 N* O2 Y0 c% E0 a( C Q( T
begin ( ?$ U, D& E0 J9 U4 [0 t* R
cur:=vn; ( s) k0 M1 I, c* W& |; A! L let vn take place of the oldest string in the tabu list; $ D. b9 |% `' k4 }% B2 h- w end; 6 e7 J/ W+ e+ s$ H until (termination-condition); ) F" }- e/ G0 A0 N
end; 2 o' W/ y- \+ Y9 d. V- m% Y9 S : T, {! U0 D: S* J& u
/ i, q8 V4 O b9 W. W F
以上程序中有关键的几点: E- B; l; z9 O: ~% y4 |0 [
(1)禁忌对象:可以选取当前的值(cur)作为禁忌对象放进tabu 3 W- D' j0 P1 n9 G* \ list,也可以把和当然值在同一“等高线”上的都放进tabu : B) y8 ?4 X* M/ d list。" n% I8 U' r1 e) w* `
(2)为了降低计算量,禁忌长度和禁忌表的集合不宜太大,但是禁忌长度太小容易循环搜索,禁忌表太小容易陷入“局部极优解”。 / C$ R5 Q2 g$ E/ h9 n' [ ) f$ B/ p: n& g3 a: q* Y% P4 i
4 m3 G$ h, M( e i (3)上述程序段中对best_to_far的操作是直接赋值为最优的“解禁候选解”,但是有时候会出现没有大于best_to_far的,候选解也全部被禁的 “死锁”状态,这个时候,就应该对候选解中最佳的进行解禁,以能够继续下去。2 v& |; V0 f8 n
2 u& \3 V. D/ ?3 s* v# I : H: m9 Y- {0 ?* v
(4)终止准则:和模拟退火,遗传算法差不多,常用的有:给定一个迭代步数;设定与估计的最优解的距离小于某个范围时,就终止搜索;当与最优解的距离连续若干步 " z; b! K- f6 c8 }保持不变时,终止搜索;7 b9 n; Q- Q1 C4 K
, |* Q. }, h# |, g $ p7 `/ x. W3 F) ` 禁忌搜索是对人类思维过程本身的一种模拟,它通过对一些局部最优解的禁忌(也可以说是记忆)达到接纳一部分较差解,从而跳出局部搜索的目的。& c- z1 @ c( Q# H$ O
8 c: x* N6 i* ?) i
R; S }) ^/ b' z9 D, w 人工神经网络(Artificial Neural Network,ANN) ' ]+ U$ W; W/ K 6 g+ m" j) ]0 y& g. O $ |! l( B8 I9 S* c0 _3 S
神经网络从名字就知道是对人脑的模拟。它的神经元结构,它的构成与作用方式都是在模仿人脑,但是也仅仅是粗糙的模仿,远没有达到完美的地步。和冯·诺依曼机不同! A5 ~. ^8 p7 b0 c5 P/ V- n h
,神经网络计算非数字,非精确,高度并行,并且有自学习功能。 # h& d- n( g8 @; x " V" y9 h, W; R% i