数学建模社区-数学中国

标题: 原来数模的十大经典算法可以这样理解哦 [打印本页]

作者: 梦里花111    时间: 2013-9-25 22:26
标题: 原来数模的十大经典算法可以这样理解哦
本帖最后由 梦里花111 于 2013-9-25 22:28 编辑
: W+ r! s8 Y8 R
; X8 f3 j8 K) {( f! P' |4 A在工程实践中,经常会接触到一些比较“新颖”的算法或理论,比如:等。这些算法或理论都有一些共同的特性(比如模拟自然过
& M& L* _. b# U5 L$ N: p5 ~程),通称为“智能算法”。它们在解决一些复杂的工程问题时大有用武之地。

! k6 i% G  T' |2 g" S5 O  

: U5 F: r( A1 x  
) w0 J8 S! s4 M  这些算法都有什么含义?首先给出个局部搜索,模拟退火,遗传算法,禁忌搜索的形象比喻:$ m; K/ j# ?$ g8 Z0 A
  

7 R" _9 x3 e: K% b8 L9 U, b    W9 h, d8 M' q- H  C" \
  为了找出地球上最高的山,一群有志气的兔子们开始想办法。2 {1 V$ a! y0 e
  
, l  F- a9 i7 r1 G
  
: q& b) u  R+ U* {- V  1.兔子朝着比现在高的地方跳去。他们找到了不远处的最高山峰。但是这座山不一定是珠穆朗玛峰。这就是局部搜索,它不能保证局部最优值就是全局最优值。
% b. B9 C* q- u) y. a# [  

* S9 S& E; H( Y% [7 t, C% r  W1 e9 P  2 `. C5 p! [: e
  2.兔子喝醉了。他随机地跳了很长时间。这期间,它可能走向高处,也可能踏入平地。但是,他渐渐清醒了并朝最高方向跳去。这就是模拟退火。
7 ?( X, i6 @' E+ z  

0 |% o/ _5 }1 }8 Z& \) d  
# r7 c$ J( G( a1 @, b- `. Q  3.兔子们吃了失忆药片,并被发射到太空,然后随机落到了地球上的某些地方。他们不知道自己的使命是什么。但是,如果你过几年就杀死一部分海拔低的兔子,多产的8 W/ Y& I8 n7 s$ S
兔子们自己就会找到珠穆朗玛峰。这就是遗传算法。
+ e7 H: g1 T/ V' |- |+ N
  

- ?4 z( M' j3 |2 F7 r) F  + d( m$ j0 P) G/ \; K1 u- A
  4.兔子们知道一个兔的力量是渺小的。他们互相转告着,哪里的山已经找过,并且找过的每一座山他们都留下一只兔子做记号。他们制定了下一步去哪里寻找的策略。这
* t7 Q! C2 F' @) Q  t, Z# @就是禁忌搜索。
; r% s8 I% p$ J9 z* a6 J' G% Q6 ]
  
' ?, k1 S% _* Z6 j% U8 g
  
. j, ]9 T: E9 ^, K  智能优化算法的概述
5 U* t6 S8 R2 [" U5 {  

$ q9 p; w1 j- q' @  5 J! H0 o7 ~, r2 i+ y# X4 M8 W
  智能优化算法要解决的一般是最优化问题。最优化问题可以分为(1)求解一个函数中,使得函数值最小的自变量取值的函数优化问题和(2)在一个解空间里面,寻找最$ V; R7 k% e( X2 s1 N8 {3 w
优解,使目标函数值最小的组合优化问题。典型的组合优化问题有:旅行商问题(
Traveling ( }" Y4 T4 W6 l! b
  Salesman Problem,TSP),加工调度问题(
Scheduling 7 V/ m: B  [# s7 z" N& r
  Problem),0-1背包问题(
Knapsack * ]  h0 J+ I9 c/ {- h
  Problem),以及装箱问题(Bin Packing Problem)等。
, Z0 v1 W( G, k' O+ ^) d! H/ k
  优化算法有很多,经典算法包括:有线性规划,动态规划等;改进型局部搜索算法包括爬山法,最速下降法等,本文介绍的模拟退火、遗传算法以及禁忌搜索称作指导性搜
3 f! L: L5 {& D& W7 T; l% @索法。而神经网络,混沌搜索则属于系统动态演化方法。

. z2 s9 i# [% \4 O  

1 Z. i1 [" E/ C0 c/ w0 h  - P& l9 `3 ~4 O* a8 U* K( G
  优化思想里面经常提到邻域函数,它的作用是指出如何由当前解得到一个(组)新解。其具体实现方式要根据具体问题分析来定。
/ n9 ^2 V- M, j7 |* v' ]/ k+ o  

1 ^8 Y1 |; t% A8 K6 D  
$ j1 S7 R6 x: S. L  一般而言,局部搜索就是基于贪婪思想利用邻域函数进行搜索,若找到一个比现有值更优的解就弃前者而取后者。但是,它一般只可以得到“局部极小解”,就是说,可能
7 }; }5 Y) y: ?% t这只兔子登“登泰山而小天下”,但是却没有找到珠穆朗玛峰。而模拟退火,遗传算法,禁忌搜索,神经网络等从不同的角度和策略实现了改进,取得较好的“全局最小解 ”。
. X& c  x" G; K2 D  H1 y# H
  
* f1 f0 J; `- S4 o! X9 m) f$ Q
  
7 g) K6 x8 I% s% i$ H  模拟退火算法(Simulated Annealing,SA)* L: v/ x) k2 U$ c
  
$ W+ d# k. G, g# @5 G0 N
  
9 a" b; n: }# b  模拟退火算法的依据是固体物质退火过程和组合优化问题之间的相似性。物质在加热的时候,粒子间的布朗运动增强,到达一定强度后,固体物质转化为液态,这个时候再8 I3 a) C( U1 H% @9 W$ C
进行退火,粒子热运动减弱,并逐渐趋于有序,最后达到稳定。
# @4 W2 A" @+ M; T4 A3 a
  

3 A& {/ P8 O! H- q& ^- l    m$ J) _  \3 ~: l/ M" q# c' z
  模拟退火的解不再像局部搜索那样最后的结果依赖初始点。它引入了一个接受概率p。如果新的点(设为pn)的目标函数f(pn)更好,则p=1,表示选取新点;否
6 W; L: d  Q1 l  J则,接受概率p是当前点(设为pc)的目标函数f(pc),新点的目标函数f(pn)以及另一个控制参数“温度”T的函数。也就是说,模拟退火没有像局部搜索那
. Q  C6 B* m4 k样每次都贪婪地寻找比现在好的点,目标函数差一点的点也有可能接受进来。随着算法的执行,系统温度T逐渐降低,最后终止于某个低温,在该温度下,系统不再接受变
! S: R/ T' \9 {! E3 q$ M2 B. ~7 Q5 v化。

- k: K/ |3 `$ Z  
# y3 u9 E2 J& T
  
7 p: X+ Z+ Z3 L1 f6 f  b6 y% D) T7 |  模拟退火的典型特征是除了接受目标函数的改进外,还接受一个衰减极限,当T较大时,接受较大的衰减,当T逐渐变小时,接受较小的衰减,当T为0时,就不再接受衰
6 e1 [+ ^, L" G3 c减。这一特征意味着模拟退火与局部搜索相反,它能避开局部极小,并且还保持了局部搜索的通用性和简单性。
. }2 d" N4 `. j2 b' D/ H9 y
  在物理上,先加热,让分子间互相碰撞,变成无序状态,内能加大,然后降温,最后的分子次序反而会更有序,内能比没有加热前更小。就像那只兔子,它喝醉后,对比较6 y1 V. w4 `7 t6 G& H# t
近的山峰视而不见,迷迷糊糊地跳一大圈子,反而更有可能找到珠峰。
6 Z4 E) j0 h! G" B
  
, z; Y5 d- \1 s% U
  
+ g( p! d2 y& w1 @' U* x  值得注意的是,当T为0时,模拟退火就成为局部搜索的一个特例。* ^7 z% E6 I  R$ h
  
/ V& s3 l- w9 G6 H  u
  
$ _. ^: C5 c; [, a0 Q  模拟退火的伪码表达:! W  C4 _6 S, d8 D4 u6 ~# ?0 P
  
procedure simulated annealing
/ ~" j. g- `, q. p! n2 u  
begin % `) Q! s% }6 ^9 r8 `
  
t:=0;
6 j$ p6 K( ^# w' R' v  
initialize temperature T / Y2 I5 H2 |9 q' j
  
select a current string vc at random; / j4 a! y" n. W) p/ ?
  
evaluate vc;
8 h" [/ ?9 y5 \. C6 B$ ]  
repeat
2 W; V0 r& n  ?: h. ?, m  
repeat ! F- _/ U0 d7 V. o. l* q3 I
  
select a new string vn in the neighborhood of vc; (1)
" R- D( C! @; m. s( F* l3 F& u+ L; @   if f(vc)  
then vc:=vn; 9 Q: p: s) ~+ }& D" X# e
   else if random [0,1]   
then vc:=vn; % k6 C+ F, R( l7 W4 ?
  
until (termination-condition) (3) , O. s  Z9 M  h* k) G
  
T:=g(T,t); (4) 5 }( v9 R9 Y' J. w4 I
  
T:=t+1;
+ o  f' C7 S/ D  
until (stop-criterion) (5)
) [4 v7 Z' m- E" K  
end; 0 I5 s- C* j7 O. }' Y
  
8 Y9 u) D- a1 v- V
  
0 \& s1 v7 j! V5 P  上面的程序中,关键的是(1)新状态产生函数,(2)新状态接受函数,(3)抽样稳定准则,(4)退温函数,(5)退火结束准则" t6 G+ |0 _: c$ o7 S9 l9 d
  (简称三函数两准则)是直接影响优化结果的主要环节。虽然实验结果证明初始值对于最后的结果没有影响,但是初温越高,得到高质量解的概率越大。所以,应该尽量选6 C1 M; q1 X- y$ p7 v
取比较高的初温。
. a( _- T8 [$ E% [6 o
  

" M0 F, G: K9 r" z" R  - j& ^3 L6 Q: W, E" N$ U
  上面关键环节的选取策略:
/ {9 V( ~" P; l% h  v3 }  (1)状态产生函数:候选解由当前解的邻域函数决定,可以取互换,插入,逆序等操作产生,然后根据概率分布方式选取新的解,概率可以取均匀分布、正态分布、高斯
1 U0 g' J9 B0 S# r7 N7 A分布、柯西分布等。

* p8 m, c. w2 z  

5 ~5 i8 d( M0 J' t$ E. Y3 n  4 [* m9 _  ^# T, q  q# c( [2 y+ }
  (2)状态接受函数:这个环节最关键,但是,实验表明,何种接受函数对于最后结果影响不大。所以,一般选取min
2 I7 L. D; \$ s1 x& }2 K  K( |3 @  [1, exp ((f (vn)-f (vc))/T)]。
8 W! K' G7 s& J' N3 S1 c
  (3)抽样稳定准则:一般常用的有:检验目标函数的均值是否稳定;连续若干步的目标值变化较小;规定一定的步数;
. T" U) V6 |( y2 x" ]! J3 T
  
; j2 [: c: v. p% W
  
7 ~3 w9 R) i7 Y5 M9 J  (4)退温函数:如果要求温度必须按照一定的比率下降,SA算法可以采用2 T4 B) S6 r( S3 \
  ,但是温度下降很慢;快速SA中,一般采用

! K- \( c5 z! a& b9 \  。目前,经常用的是 , 是一个不断变化的值。

7 u# n2 ~4 N6 ?+ L  (5)退火结束准则:一般有:设置终止温度;设置迭代次数;搜索到的最优值连续多次保持不变;检验系统熵是否稳定。
6 U( M/ a! X1 I2 D- _7 u9 p
  
$ M* o  w# R, ~1 B4 @2 E
  * h4 N( R% I2 b
  为了保证有比较优的解,算法往往采取慢降温、多抽样、以及把“终止温度”设的比较低等方式,导致算法运行时间比较长,这也是模拟退火的最大缺点。人喝醉了酒办起! @3 ?& P* R5 F# x
事来都不利索,何况兔子?

& i# K6 ~* o% @  z' J; p. z- }  
' P4 B1 G; E# H/ F
  / m* A$ `2 M* Z8 H8 K& b
  遗传算法(Genetic Algorithm, GA)% t" j$ r5 f& }& ^8 [0 O3 h7 k4 n
  

/ ?7 d- U$ b* d( C& ~- f9 m8 u2 v/ M  - ]2 \3 u) J# @5 o& b& |& [
  “物竞天择,适者生存”,是进化论的基本思想。遗传算法就是模拟自然界想做的事。遗传算法可以很好地用于优化问题,若把它看作对自然过程高度理想化的模拟,更能
8 X# J% i) Z' I# _6 `1 g显出它本身的优雅——虽然生存竞争是残酷的。

) ]) \- Q# q- u- q  

8 f. [, A1 u2 |# o. L2 ]  
4 F+ E4 b& @" m) U8 U  遗传算法以一种群体中的所有个体为对象,并利用随机化技术指导对一个被编码的参数空间进行高效搜索。其中,选择、交叉和变异构成了遗传算法的遗传操作;参数编码
9 \. o; ~( y) W! q! S7 {+ U、初始群体的设定、适应度函数的设计、遗传操作设计、控制参数设定五个要素组成了遗传算法的核心内容。

0 {" }, s- T6 v: E$ Q# e6 Z* W  作为一种新的全局优化搜索算法,遗传算法以其简单通用、健壮性强、适于并行处理以及高效、实用等显著特点,在各个领域得到了广泛应用,取得了良好效果,并逐渐成
7 @  o% T/ \2 m2 v% P8 o3 ~为重要的智能算法之一。
2 s( |, |# \+ ]- F/ {1 K% E8 X
  
* j. Z7 X9 m' [- [
  9 G) d( o4 o- n0 \7 s2 X) q6 P
  遗传算法的伪码:  I/ g  {  l: b3 ]+ R1 P
  
( W8 c/ A% k# M1 h" l
  ) z( \. t7 ~; H2 z: S
  procedure genetic algorithm
# }4 N- S* K* R# W$ x! S3 @( M  
begin
( ~/ l- N2 T9 p6 z   initialize a group and evaluate the fitness value ; (1)

% y( i+ g/ ]% Q   while not convergent (2)
" H6 I% U" B5 t( B' H
  
begin
9 F) [* c0 s% f# h$ U- Z: U   select; (3)

$ Y& I9 }7 @5 O, t. `& I# S   if random[0,1]   crossover; (4)

3 F# u8 B- k; w! f2 C; O3 P   if random (0,1)   mutation; (5)
+ e7 c* K9 R/ w7 J8 C# M* v
  
end;
1 y* ?( `$ E1 r* ?$ y# Q: w  
end
% z5 {' ^! X; }: V  上述程序中有五个重要的环节:

# E# Z* V: Q* B' ?3 v  (1)编码和初始群体的生成:GA在进行搜索之前先将解空间的解数据表示成遗传空间的基因型串结构数据,这些串结构数据的不同组合便构成了不同的点。然后随机产4 \( \/ m5 d% T7 x& b
生N个初始串结构数据,每个串结构数据称为一个个体,

2 l) M; @  N3 a  N个体构成了一个群体。GA以这N个串结构数据作为初始点开始迭代。
# N5 ~/ o7 k$ o8 d
  

) F7 U) {, O$ z" A  $ d% P( i3 Z* n7 ~$ D; p% W
  比如,旅行商问题中,可以把商人走过的路径进行编码,也可以对整个图矩阵进行编码。编码方式依赖于问题怎样描述比较好解决。初始群体也应该选取适当,如果选取的
- R# h9 U/ k4 b% x" g3 \% U. L过小则杂交优势不明显,算法性能很差(数量上占了优势的老鼠进化能力比老虎强),群体选取太大则计算量太大。

9 |9 u& D: a6 h; w5 B+ X1 Y  

, d5 T: u! u" ~% ~5 A  
& `( n" T& n+ z: V  (2)检查算法收敛准则是否满足,控制算法是否结束。可以采用判断与最优解的适配度或者定一个迭代次数来达到。$ h0 r" p4 u4 ~1 Y* p! A" V5 t! M  |
  
" ?1 O6 W2 Y, N( ?* D
  4 E* K9 G: W' [5 j9 u0 x
  (3)适应性值评估检测和选择:适应性函数表明个体或解的优劣性,在程序的开始也应该评价适应性,以便和以后的做比较。不同的问题,适应性函数的定义方式也不同7 ?, ^& {% @! i; Y
。根据适应性的好坏,进行选择。选择的目的是为了从当前群体中选出优良的个体,使它们有机会作为父代为下一代繁殖子孙。遗传算法通过选择过程体现这一思想,进行+ Y% R, U* q$ j
选择的原则是适应性强的个体为下一代贡献一个或多个后代的概率大。选择实现了达尔文的适者生存原则。

/ y) s; [- u0 Q, q. \2 s6 H. z  
: v! ]& z% r1 Y$ s. t3 C% ?
  
9 @: z8 _/ t+ x  (4)杂交:按照杂交概率(pc)进行杂交。杂交操作是遗传算法中最主要的遗传操作。通过杂交操作可以得到新一代个体,新个体组合了其父辈个体的特性。杂交体现
+ ]- x* Q/ D' r/ h2 `了信息交换的思想。
2 r9 D& f2 m/ j$ K( R/ `
  
& {# G* @6 N7 D. o8 O- @/ v
  
: n( W1 k. Z+ r( I# H) y  可以选定一个点对染色体串进行互换,插入,逆序等杂交,也可以随机选取几个点杂交。杂交概率如果太大,种群更新快,但是高适应性的个体很容易被淹没,概率小了搜: T4 y4 x. x/ g( u  K  G
索会停滞。

8 S, q( X6 y" L7 I  _  

  Q$ Q9 @. ~6 r; K0 Z  6 |, s2 Z' R% i
  (5)变异:按照变异概率(pm)进行变异。变异首先在群体中随机选择一个个体,对于选中的个体以一定的概率随机地改变串结构数据中某个串的值。同生物界一样, GA中变异发生的概率很低。变异为新个体的产生提供了机会。, Y+ [! y! a! P, L
  

' j3 U! F2 W( k& K5 U  
1 C8 q  x& @' O1 f3 n  变异可以防止有效基因的缺损造成的进化停滞。比较低的变异概率就已经可以让基因不断变更,太大了会陷入随机搜索。想一下,生物界每一代都和上一代差距很大,会是
, n5 Y% R- T3 z怎样的可怕情形。
0 a! O$ s+ e) [. m: e
  
4 @9 ~) k, K3 O
  
6 L/ U7 T- K$ l+ C( P  就像自然界的变异适和任何物种一样,对变量进行了编码的遗传算法没有考虑函数本身是否可导,是否连续等性质,所以适用性很强;并且,它开始就对一个种群进行操作
% L* S+ v( h( Q6 U,隐含了并行性,也容易找到“全局最优解”。
; f4 A1 a9 G1 e) ~3 I% R
  

3 z( }# A  q2 q, C2 Z  
2 z) a8 L8 s* b( B  禁忌搜索算法(Tabu Search,TS)
5 X3 Y$ O9 \3 H4 T  
9 _2 T' r9 E$ I  x8 Q" ?* y
  
5 K+ j9 c9 e( M' L# [  为了找到“全局最优解”,就不应该执着于某一个特定的区域。局部搜索的缺点就是太贪婪地对某一个局部区域以及其邻域搜索,导致一叶障目,不见泰山。禁忌搜索就是
. h: X5 S/ R8 N0 L对于找到的一部分局部最优解,有意识地避开它(但不是完全隔绝),从而获得更多的搜索区间。兔子们找到了泰山,它们之中的一只就会留守在这里,其他的再去别的地: B/ K9 k, N1 t1 p/ R6 q' f2 M
方寻找。就这样,一大圈后,把找到的几个山峰一比较,珠穆朗玛峰脱颖而出。
; G! C- z6 R6 Z$ L- y: L) U
  
' {+ [" V* g2 u5 d! a+ E4 y2 B
  
; u' Q, Q0 u) C$ f  当兔子们再寻找的时候,一般地会有意识地避开泰山,因为他们知道,这里已经找过,并且有一只兔子在那里看着了。这就是禁忌搜索中“禁忌表(tabu
( O6 r. c. R2 k# l% Z  list)”的含义。那只留在泰山的兔子一般不会就安家在那里了,它会在一定时间后重新回到找最高峰的大军,因为这个时候已经有了许多新的消息,泰山毕竟也有一
" H+ m* v, G; z4 t8 S3 `个不错的高度,需要重新考虑,这个归队时间,在禁忌搜索里面叫做“禁忌长度(
tabu
' P. ?, Q2 r* Q+ ]' b7 X2 h) Q4 o  length)”;如果在搜索的过程中,留守泰山的兔子还没有归队,但是找到的地方全是华北平原等比较低的地方,兔子们就不得不再次考虑选中泰山,也就是说,当! T( |& E6 Z# [+ ?; x
一个有兔子留守的地方优越性太突出,超过了
“best
9 m7 ~; F. Y$ g' `) S2 e8 z# y  
to
/ K, m* m4 \' k/ ]8 \; S  far”的状态,就可以不顾及有没有兔子留守,都把这个地方考虑进来,这就叫“特赦准则(
aspiration ' v5 W3 ]0 i# @' I$ a, _; O( G. E
  criterion)”。这三个概念是禁忌搜索和一般搜索准则最不同的地方,算法的优化也关键在这里。

/ y- k; b2 U2 e. a  l  

. p2 N; {* v/ G& P, T# k4 ~  
9 y2 H+ }* I' C& R1 H# ?9 Z6 w) W  伪码表达:( [9 {3 z( _3 t) k
  
procedure tabu search; ' i4 Q5 ?: l) x' n+ T; D* x
  
begin
- a. F; M8 X  ?2 E. R4 ^  
initialize a string vc at random,clear up the tabu list;
0 i& F5 {( |* ~: }/ X0 n  
cur:=vc;
- n$ G+ M* o1 v" P+ W3 K' V  
repeat
; U0 x+ C$ U2 `6 k  
select a new string vn in the neighborhood of vc;
: o' k# E! J8 n  
if va>best_to_far then {va is a string in the tabu list} + J+ y4 P( k2 o% R3 L- G( ~
  
begin
3 ^* i8 O5 ]$ Z0 l$ B+ y& t& ~  
cur:=va; 0 B6 c2 Q$ o8 V8 p" f% D$ J6 @% E) o
  
let va take place of the oldest string in the tabu list; 4 c) E* t) P. W! e5 b
  
best_to_far:=va; $ q% h9 c* D( y
  
end else
  N6 }& F. W% n6 P0 N0 P% f  
begin
. p5 [+ e$ W# p- R* z0 w  
cur:=vn; ! Y+ p2 ^: {2 `3 Z
  
let vn take place of the oldest string in the tabu list; 8 F* }- A% k. t6 L% W; w
  
end;
' F" K  Y& J, o$ x2 q  
until (termination-condition);
4 Q8 ?5 Q: }1 j8 K  
end; 4 H: C4 G* O% @
  

" F. |  a  Q4 u. r- ]7 h, p& X  * S' \1 v9 G" a+ e$ f
  以上程序中有关键的几点:
  d* K/ a3 N0 }4 n* r; S: c% ^  (1)禁忌对象:可以选取当前的值(cur)作为禁忌对象放进
tabu / U5 `5 I" ^% U+ w
  list,也可以把和当然值在同一“等高线”上的都放进
tabu
  j- V9 X: ~! G. b) ]0 e  list。
8 ]. M- z; @& m, n3 V" ]: d
  (2)为了降低计算量,禁忌长度和禁忌表的集合不宜太大,但是禁忌长度太小容易循环搜索,禁忌表太小容易陷入“局部极优解”。

8 m6 P+ P/ [3 }# m1 ~1 k6 k- ~% G  

0 y  z9 u8 J* M8 F4 k) O' L" O  
  N% f2 d$ E# e% ?  N3 N6 p  (3)上述程序段中对best_to_far的操作是直接赋值为最优的“解禁候选解”,但是有时候会出现没有大于best_to_far的,候选解也全部被禁的 “死锁”状态,这个时候,就应该对候选解中最佳的进行解禁,以能够继续下去。
3 B1 n. C& {7 O( G! k$ s+ S. N* f  
+ C2 ^4 e' a! j4 o
  
& X9 G& A0 E9 ?; T  (4)终止准则:和模拟退火,遗传算法差不多,常用的有:给定一个迭代步数;设定与估计的最优解的距离小于某个范围时,就终止搜索;当与最优解的距离连续若干步" p+ Q# m( E: p+ v+ s
保持不变时,终止搜索;
; H2 l# K/ v' j. i
  

6 t, A) X( [- e4 N9 f* A) ]  
1 m  ]9 _2 A/ c9 X% u2 Z3 y  禁忌搜索是对人类思维过程本身的一种模拟,它通过对一些局部最优解的禁忌(也可以说是记忆)达到接纳一部分较差解,从而跳出局部搜索的目的。$ U9 ^! T5 A! B- V; V5 \% C' H  q
  

+ y  d$ g8 m: Z+ k$ Z9 r; q9 y  7 k( m8 W9 A! ^' P/ ]3 d% q4 b4 f
  人工神经网络(Artificial Neural Network,ANN)
' |) K+ u3 r4 i  
8 F& h2 N/ S% u( k  v/ C8 U
    O( B7 [. b$ U1 D
  神经网络从名字就知道是对人脑的模拟。它的神经元结构,它的构成与作用方式都是在模仿人脑,但是也仅仅是粗糙的模仿,远没有达到完美的地步。和冯·诺依曼机不同8 h( s5 J: v7 X' \4 X
,神经网络计算非数字,非精确,高度并行,并且有自学习功能。
- ~$ B! ?9 l# Q2 K9 d
  
* n) h% d8 S' h! U9 Q
  ; |6 ~1 {6 O6 O  v
  生命科学中,神经细胞一般称作神经元,它是整个神经结构的最基本单位。每个神经细胞就像一条胳膊,其中像手掌的地方含有细胞核,称作细胞体,像手指的称作树突,. F1 d- x% R& G* c3 N+ @. ^: G
是信息的输入通路,像手臂的称作轴突,是信息的输出通路;神经元之间错综复杂地连在一起,互相之间传递信号,而传递的信号可以导致神经元电位的变化,一旦电位高
7 t1 w0 r$ G1 G! F8 o2 c出一定值,就会引起神经元的激发,此神经元就会通过轴突传出电信号。

: ~+ g: W1 i; K  

1 f5 T; W) g$ Y+ O5 V  + h4 l2 u. v" c7 D* d
  而如果要用计算机模仿生物神经,就需要人工的神经网络有三个要素:(1)形式定义人工神经元;(2)给出人工神经元的连接方式,或者说给出网络结构;(3)给出% W  S) \# ^. m
人工神经元之间信号强度的定义。

9 \2 v, R; i( h+ P$ h  
0 V, k6 p- ?2 k) Q/ F) s
  
2 H5 r( j' B9 l7 X4 P/ F  历史上第一个人工神经网络模型称作M-P模型,非常简单:0 t. D9 g2 q3 o- m: a, J
  

# P8 @- a+ V2 X4 `  
/ t. W* V% T8 J! V; J  其中,
5 t+ r) O3 }8 b* u6 s  表示神经元i在t时刻的状态,为1表示激发态,为0表示抑制态;

6 H& U* _" K+ W. h% l4 l+ l  是神经元i和j之间的连接强度;

, D5 U0 K) X0 D! s' X4 V  表示神经元i的阈值,超过这个值神经元才能激发。

" S% k( W( ^: N4 h' L& ]0 ?  这个模型是最简单的神经元模型。但是功能已经非常强大:此模型的发明人McCulloch和Pitts已经证明,不考虑速度和实现的复杂性,它可以完成当前数字
2 g* V& q! O% T! \5 F. q计算机的任何工作。

4 O$ ]) W- Y8 R  
. ?# O0 u) V1 W9 k$ F" ]
    V- l; M/ B) n# V* U5 {! z
  以上这个M-P模型仅仅是一层的网络,如果从对一个平面进行分割的方面来考虑的话,M-P网络只能把一个平面分成个半平面,却不能够选取特定的一部分。而解决的
8 [* _. u4 M# y& L$ W6 w办法就是“多层前向网路”。
9 m5 r8 x5 s: }
  
$ I; ]8 S+ W  D* e3 g
  
1 Q) T! ~; V0 f5 s0 J   图2
2 y! k6 f3 F" [( Z  图2是多层前向网络的示意图。最下面的

8 d. y6 f3 y  `% V  称作输入层,最上面一层称作输出层,任何一个中间层都接受来自前一层的所有输入,加工后传入后一层。每一层的神经元之间没有联系,输入输出层之间也没有直接联系
& q7 D; O; J( j" m,并且仅仅是单向联系,没有反馈。这样的网络被称作“多层前向网络”。数据在输入后,经过每一层的加权,最后输出结果。

6 v0 y* a/ ^; S5 f& n0 E  

' o3 T2 `1 V1 n: s  ! e# W. }2 ~% b; y# E' o9 w
   图3
7 p3 s* q& M: J: X% G% s  如图3,用可覆盖面来说明多层网络的功能:单层网络只能把平面分成两部分,双层网络就可以分割任意凸域,多层网络则可以分割任意区域。
  R- \8 v# Y+ A, B+ z
  为了让这种网络有合适的权值,必须给网络一定的激励,让它自己学习,调整。一种方法称作“向后传播算法(
Back : N4 b' q# \: T3 z# \8 W
  Propagation,BP)”,其基本思想是考察最后输出解和理想解的差异,调整权值,并把这种调整从输出层开始向后推演,经过中间层,达到输入层。
$ ?) a' ^4 S/ F; x( `( r
  

9 H/ ~2 g6 u8 o% j5 h  " R7 H4 S5 X5 u0 b" O4 A$ {& X
  可见,神经网络是通过学习来达到解决问题的目的,学习没有改变单个神经元的结构和工作方式,单个神经元的特性和要解决的问题之间也没有直接联系,这里学习的作用; C/ w7 O! k9 u) ~0 H
是根据神经元之间激励与抑制的关系,改变它们的作用强度。学习样本中的任何样品的信息都包含在网络的每个权值之中。
0 K+ P" ]; }' t1 ~
  
4 K9 Q2 ~. f) p( n7 B0 A7 @  `
  , |6 ?+ d8 o1 ~" }* r: I) W7 |# [
  BP算法中有考察输出解和理想解差异的过程,假设差距为w,则调整权值的目的就是为了使得w最小化。这就又包含了前文所说的“最小值”问题。一般的BP算法采用
! V, p/ K5 D/ p- ]- g的是局部搜索,比如最速下降法,牛顿法等,当然如果想要得到全局最优解,可以采用模拟退火,遗传算法等。当前向网络采用模拟退火算法作为学习方法的时候,一般成4 H$ V9 }* v9 ~" e" X' A4 B
为“波尔兹曼网络”,属于随机性神经网络。

2 t3 p6 F% J! \8 b. Y! x  

8 [# L) d3 `' h0 @% T  5 N! p# [9 z& X9 W: V3 O9 v/ o
  在学习BP算法学习的过程中,需要已经有一部分确定的值作为理想输出,这就好像中学生在学习的时候,有老师的监督。如果没有了监督,人工神经网络该怎么学习?  G+ T, y  D7 v. T/ Y( P
  
; t3 d9 K9 }( l/ z2 E
  # s" C0 D: c. L! ?$ P4 O
  就像没有了宏观调控,自由的市场引入了竞争一样,有一种学习方法称作“无监督有竞争的学习”。在输入神经元i的若干个神经元之间开展竞争,竞争之后,只有一个神
( ?8 F" s- D; l" c经元为1,其他均为0,而对于失败的神经元,调整使得向对竞争有利的方向移动,则最终也可能在一次竞争中胜利;
2 ^% C' W9 F  A% E
  

) ~! m1 `3 J5 N+ o9 \- _  $ E7 z% u- X% q( {/ \9 N5 I
  人工神经网络还有反馈网络如Hopfield网络,它的神经元的信号传递方向是双向的,并且引入一个能量函数,通过神经元之间不断地相互影响,能量函数值不断下
* r/ L/ K% M( h* J6 B7 `降,最后能给出一个能量比较低的解。这个思想和模拟退火差不多。
/ v4 C. I, S1 C5 d7 |# B7 `# H
  

4 s* H- z) B; F' b- @2 d  % R9 n; @# T) y/ C  x: f* r( A
  人工神经网络应用到算法上时,其正确率和速度与软件的实现联系不大,关键的是它自身的不断学习。这种思想已经和冯·诺依曼模型很不一样。# w) ~7 N8 b# V  ]/ j& \
  

7 n; T% |' X& S# I4 q, a  
5 w! e/ l5 K; y+ v  总结
# o% D1 V. w( v+ u  模拟退火,遗传算法,禁忌搜索,神经网络在解决全局最优解的问题上有着独到的优点,并且,它们有一个共同的特点:都是模拟了自然过程。模拟退火思路源于物理学中+ q; a. `- C' c* J4 n
固体物质的退火过程,遗传算法借鉴了自然界优胜劣汰的进化思想,禁忌搜索模拟了人类有记忆过程的智力过程,神经网络更是直接模拟了人脑。
# b4 D1 x# v, i" D$ G
  
, {  k  E: _- R% E, h! j# D
  + g, ]6 K8 _: X, [
  
, |- g: F2 ^; S  它们之间的联系也非常紧密,比如模拟退火和遗传算法为神经网络提供更优良的学习算法提供了思路。把它们有机地综合在一起,取长补短,性能将更加优良。
( |/ \9 q0 \( L  
  D  I9 G; e. M: F. w
  这几种智能算法有别于一般的按照图灵机进行精确计算的程序,尤其是人工神经网络,是对计算机模型的一种新的诠释,跳出了冯·诺依曼机的圈子,按照这种思想来设计$ b! E# u1 V/ A; @- L
的计算机有着广阔的发展前景

% F* R" K, F5 p) J' l禁忌搜索算法(Tabu Search,TS)
# @0 y* F6 I) {
为了找到“全局最优解”,就不应该执着于某一个特定的区域。局部搜索的缺点就是太贪婪地对某一个局部区域以及其邻域搜索,导致一叶障目,不见泰山。禁忌搜索就是对于找到的一部分局部最优解,有意识地避开它(但不是完全隔绝),从而获得更多的搜索区间。兔子们找到了泰山,它们之中的一只就会留守在这里,其他的再去别的地方寻找。就这样,一大圈后,把找到的几个山峰一比较,珠穆朗玛峰脱颖而出。
2 G3 T8 H, i, _* J1 Q0 Y0 {当兔子们再寻找的时候,一般地会有意识地避开泰山,因为他们知道,这里已经找过,并且有一只兔子在那里看着了。这就是禁忌搜索中“禁忌表(tabu list)”的含义。那只留在泰山的兔子一般不会就安家在那里了,它会在一定时间后重新回到找最高峰的大军,因为这个时候已经有了许多新的消息,泰山毕竟也有一个不错的高度,需要重新考虑,这个归队时间,在禁忌搜索里面叫做“禁忌长度(tabu length)”;如果在搜索的过程中,留守泰山的兔子还没有归队,但是找到的地方全是华北平原等比较低的地方,兔子们就不得不再次考虑选中泰山,也就是说,当一个有兔子留守的地方优越性太突出,超过了“best to far”的状态,就可以不顾及有没有兔子留守,都把这个地方考虑进来,这就叫“特赦准则(aspiration criterion)”。这三个概念是禁忌搜索和一般搜索准则最不同的地方,算法的优化也关键在这里。
5 e/ N8 E! A4 Y" \伪码表达:
. i/ ?: K( F: C& kprocedure tabu search;4 f) @4 Y6 A' h' L7 _0 V! N
begin% o/ d4 ]4 n7 n6 h1 L: y
initialize a string vc at random,clear up the tabu list;
  \9 h; S% e' `- _cur:=vc;+ ]) @8 l+ c' E+ q' }
repeat$ ~- Y- Z+ T" ?! W
select a new string vn in the neighborhood of vc;
2 {! z+ n% s8 V5 W5 C2 Z! Rif va>best_to_far then {va is a string in the tabu list}
9 z9 O0 A3 s1 S# V# ~: m# i# }begin
8 G3 ~# f& \- |cur:=va;$ J$ C' l* l+ |# [
let va take place of the oldest string in the tabu list;( o& V+ h& h$ I  Y0 v- W. f0 W) k
best_to_far:=va;
8 {  X2 G, p1 `end else- @2 s* v" C' Y$ D, {
begin/ \7 d+ M1 y" H1 Z" `1 G9 @/ T
cur:=vn;
# r% ]) h# q4 H9 Y' E$ Llet vn take place of the oldest string in the tabu list;) h1 O5 k# O7 {% O3 U
end;
" X; h, x6 p: H( P5 wuntil (termination-condition);! a+ }9 f4 q% p1 Z/ r# G8 y2 e
end;9 }$ y2 y: ?' ]4 _
) x& I- T, ?! M) R
以上程序中有关键的几点:
4 M7 i+ ^! r; Z, v(1)禁忌对象:可以选取当前的值(cur)作为禁忌对象放进tabu list,也可以把和当然值在同一“等高线”上的都放进tabu list。, r7 ?) {7 y5 |& O' g
(2)为了降低计算量,禁忌长度和禁忌表的集合不宜太大,但是禁忌长度太小容易循环搜索,禁忌表太小容易陷入“局部极优解”。( ~9 f% m* ?- t7 r3 f7 O
(3)上述程序段中对best_to_far的操作是直接赋值为最优的“解禁候选解”,但是有时候会出现没有大于best_to_far的,候选解也全部被禁的“死锁”状态,这个时候,就应该对候选解中最佳的进行解禁,以能够继续下去。" P! Q4 B4 ~' z1 Z" X, N" K7 r* ]
(4)终止准则:和模拟退火,遗传算法差不多,常用的有:给定一个迭代步数;设定与估计的最优解的距离小于某个范围时,就终止搜索;当与最优解的距离连续若干步保持不变时,终止搜索;
, i" ?' T3 P. d$ B# i禁忌搜索是对人类思维过程本身的一种模拟,它通过对一些局部最优解的禁忌(也可以说是记忆)达到接纳一部分较差解,从而跳出局部搜索的目的。+ |% C3 Z. W  S, D* [+ Q* L  K& a/ A
人工神经网络(Artificial Neural Network,ANN). e% w$ I* c& C0 D3 X; s) l7 j
神经网络从名字就知道是对人脑的模拟。它的神经元结构,它的构成与作用方式都是在模仿人脑,但是也仅仅是粗糙的模仿,远没有达到完美的地步。和冯·诺依曼机不同,神经网络计算非数字,非精确,高度并行,并且有自学习功能。
3 q0 @7 ^2 _& p3 B) A6 S/ q生命科学中,神经细胞一般称作神经元,它是整个神经结构的最基本单位。每个神经细胞就像一条胳膊,其中像手掌的地方含有细胞核,称作细胞体,像手指的称作树突,是信息的输入通路,像手臂的称作轴突,是信息的输出通路;神经元之间错综复杂地连在一起,互相之间传递信号,而传递的信号可以导致神经元电位的变化,一旦电位高出一定值,就会引起神经元的激发,此神经元就会通过轴突传出电信号。
( i2 ~8 N- J0 _: L# R! s而如果要用计算机模仿生物神经,就需要人工的神经网络有三个要素:(1)形式定义人工神经元;(2)给出人工神经元的连接方式,或者说给出网络结构;(3)给出人工神经元之间信号强度的定义。, I' I; O; B8 _7 V. T! @9 ~4 @& h
历史上第一个人工神经网络模型称作M-P模型,非常简单:5 M  \" z, V9 }0 S" \
其中, 表示神经元i在t时刻的状态,为1表示激发态,为0表示抑制态; 是神经元i和j之间的连接强度; 表示神经元i的阈值,超过这个值神经元才能激发。$ ~" v* R" z. u/ I
这个模型是最简单的神经元模型。但是功能已经非常强大:此模型的发明人McCulloch和Pitts已经证明,不考虑速度和实现的复杂性,它可以完成当前数字计算机的任何工作。
. u) i. e% r: ~  r" ?) Z7 Y+ ?+ u以上这个M-P模型仅仅是一层的网络,如果从对一个平面进行分割的方面来考虑的话,M-P网络只能把一个平面分成个半平面,却不能够选取特定的一部分。而解决的办法就是“多层前向网路”。& T3 b$ v4 R. u, ~" q: G) ?, N
图2
" C; |" f/ s( s' Q5 {- D图2 是多层前向网络的示意图。最下面的称作输入层,最上面一层称作输出层,任何一个中间层都接受来自前一层的所有输入,加工后传入后一层。每一层的神经元之间没有联系,输入输出层之间也没有直接联系,并且仅仅是单向联系,没有反馈。这样的网络被称作“多层前向网络”。数据在输入后,经过每一层的加权,最后输出结果。
$ q# Z) n1 N& M2 u
图39 y; \; a+ p0 @) V2 A, e7 \
如图3,用可覆盖面来说明多层网络的功能:单层网络只能把平面分成两部分,双层网络就可以分割任意凸域,多层网络则可以分割任意区域。

6 C' U& c4 `" A6 ]% f" t为了让这种网络有合适的权值,必须给网络一定的激励,让它自己学习,调整。一种方法称作“向后传播算法(Back Propagation,BP)”,其基本思想是考察最后输出解和理想解的差异,调整权值,并把这种调整从输出层开始向后推演,经过中间层,达到输入层。
3 f" m! |' I0 @+ B+ V/ a$ h& R可见,神经网络是通过学习来达到解决问题的目的,学习没有改变单个神经元的结构和工作方式,单个神经元的特性和要解决的问题之间也没有直接联系,这里学习的作用是根据神经元之间激励与抑制的关系,改变它们的作用强度。学习样本中的任何样品的信息都包含在网络的每个权值之中。
9 E; m! I+ a% ]  i0 ^' ]BP算法中有考察输出解和理想解差异的过程,假设差距为w,则调整权值的目的就是为了使得w最小化。这就又包含了前文所说的“最小值”问题。一般的BP算法采用的是局部搜索,比如最速下降法,牛顿法等,当然如果想要得到全局最优解,可以采用模拟退火,遗传算法等。当前向网络采用模拟退火算法作为学习方法的时候,一般成为“波尔兹曼网络”,属于随机性神经网络。
  C$ k7 x2 d2 x1 _在学习BP算法学习的过程中,需要已经有一部分确定的值作为理想输出,这就好像中学生在学习的时候,有老师的监督。如果没有了监督,人工神经网络该怎么学习?$ u& B8 C* j# [  A, l( @' ~
就像没有了宏观调控,自由的市场引入了竞争一样,有一种学习方法称作“无监督有竞争的学习”。在输入神经元i的若干个神经元之间开展竞争,竞争之后,只有一个神经元为1,其他均为0,而对于失败的神经元,调整使得向对竞争有利的方向移动,则最终也可能在一次竞争中胜利;" ]7 n/ r! S8 H% Z
人工神经网络还有反馈网络如Hopfield网络,它的神经元的信号传递方向是双向的,并且引入一个能量函数,通过神经元之间不断地相互影响,能量函数值不断下降,最后能给出一个能量比较低的解。这个思想和模拟退火差不多。
  [' |3 U3 W# r( y4 a人工神经网络应用到算法上时,其正确率和速度与软件的实现联系不大,关键的是它自身的不断学习。这种思想已经和冯·诺依曼模型很不一样。. f  m: T% c3 ~$ o, h* k( t
总结' p8 x4 l+ @0 a* {9 u9 x
模拟退火,遗传算法,禁忌搜索,神经网络在解决全局最优解的问题上有着独到的优点,并且,它们有一个共同的特点:都是模拟了自然过程。模拟退火思路源于物理学中固体物质的退火过程,遗传算法借鉴了自然界优胜劣汰的进化思想,禁忌搜索模拟了人类有记忆过程的智力过程,神经网络更是直接模拟了人脑。
4 u8 c% h# o9 p+ E4 y# i) P, ^它们之间的联系也非常紧密,比如模拟退火和遗传算法为神经网络提供更优良的学习算法提供了思路。把它们有机地综合在一起,取长补短,性能将更加优良。- r+ M: L' r2 @( u, \( J, m0 c& {
这几种智能算法有别于一般的按照图灵机进行精确计算的程序,尤其是人工神经网络,是对计算机模型的一种新的诠释,跳出了冯·诺依曼机的圈子,按照这种思想来设计的计算机有着广阔的发展前景
0 ]" Y. y8 S& k6 l7 W8 J2 l
作者: 花映碧水    时间: 2013-9-27 01:09

作者: shuang0813    时间: 2013-12-18 21:04
嘿嘿~~有点意思!
作者: 梦里花111    时间: 2013-12-19 08:44
有意思吧。




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5