- 在线时间
- 18 小时
- 最后登录
- 2012-9-22
- 注册时间
- 2012-7-20
- 听众数
- 7
- 收听数
- 0
- 能力
- 0 分
- 体力
- 1995 点
- 威望
- 0 点
- 阅读权限
- 50
- 积分
- 653
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 105
- 主题
- 7
- 精华
- 0
- 分享
- 0
- 好友
- 18
升级   13.25% TA的每日心情 | 开心 2012-9-22 23:13 |
|---|
签到天数: 44 天 [LV.5]常住居民I
- 自我介绍
- 数学爱好者
 群组: 学术交流A |
数学建模十大算法漫谈
2 T5 m$ m$ M! r9 d4 e6 S6 e: L
0 m7 c- w8 _1 G/ d作者:July 二零一一年一月二十九日
; Y: B) _1 d: r* _. X/ Z, L本文参考:* j+ k5 T" S4 q6 f3 Z+ P v4 f: Y
I、 细数二十世纪最伟大的十大算法 [译者:本人July] s3 k0 q, H! a# V3 L
II、 本BLOG内 经典算法研究系列
Q2 J: c4 {; u5 o" ?( ?0 bIII、维基百科 Z, N# u. s7 `+ ^' N3 Q
------------------------------------------
2 I+ F( Y& W2 ^( J: T. b9 k博主说明:. [. _" S+ ~: [2 M' Q2 T
1、此数学建模十大算法依据网上的一份榜单而写,本文对此十大算法作一一简单介绍。
8 k' S) A! x" i- e' E$ F6 W这只是一份榜单而已,数学建模中还有很多的算法,未一一囊括。欢迎读者提供更多的好的算法。$ x p l a$ }& h
2、在具体阐述每一算法的应用时,除了列出常见的应用之外,
; ^: {- U+ p& d, X) W) B3 n同时,还会具体结合数学建模竞赛一一阐述。
" Z' q2 B- Y T/ e/ N7 e毕竟,此十大算法,在数学建模竞赛中有着无比广泛而重要的应用。8 T5 a+ i. ]" K! }* n, y2 M5 b
且,凡是标着“某某年某国某题”,即是那一年某个国家的数学建模竞赛原题。
( N" C; U0 A/ x3、此十大算法,在一些经典的算法设计书籍上,无过多阐述。$ m9 |6 x% c. I$ v) x0 m& a
若要具体细致的深入研究,还得请参考国内或国际上关于此十大算法的优秀论文。( s+ e( |& C* h) O! A
谢谢。/ b1 A4 j/ z" H' z& k
; K; K. v, M6 Q- l# o% ]
& \! e T% a$ a* ]/ ^一、蒙特卡罗算法- |; S0 q3 _" m0 P" w
1946年,美国拉斯阿莫斯国家实验室的三位科学家John von Neumann,Stan Ulam 和 Nick Metropolis
2 w# Z6 M8 a* x) e1 ]0 @) \7 f: E共同发明了,蒙特卡罗方法。5 O4 _# U0 ^6 `9 \# b' X
2 Y2 V9 T3 X, @; |* Z
此算法被评为20世纪最伟大的十大算法之一,详情,请参见我的博文:4 R) U# t# N- g' S# V
http://blog.csdn.net/v_JULY_v/archive/2011/01/10/6127953.aspx; m/ J( h" F8 U9 Q
. x* q+ i% P( g: F: u' j4 d蒙特卡罗方法(Monte Carlo method),又称随机抽样或统计模拟方法,是一种以概率统计理论为指导
8 P/ w5 g" D o+ |$ A的一类非常重要的数值计算方法。此方法使用随机数(或更常见的伪随机数)来解决很多计算问题的方
; ^: b1 Q" [% T+ H5 ]- I' y法。: z" p( u: Q( S! b+ x, X
6 F3 L+ n. w8 V. G: _由于传统的经验方法由于不能逼近真实的物理过程,很难得到满意的结果,而蒙特卡罗方法由于能够真1 w+ h# [* k+ S( {% i- H {
实地模拟实际物理过程,故解决问题与实际非常符合,可以得到很圆满的结果。
8 L" m/ j, V9 @5 u蒙特卡罗方法的基本原理及思想如下:" d& U% v; ^# ]1 z3 d+ D
当所求解问题是某种随机事件出现的概率,或者是某个随机变量的期望值时,通过某种“实验”的方法
9 y2 [/ _/ @2 o$ t3 t* v$ ~,以这种事件出现的频率估计这一随机事件的概率,或者得到这个随机变量的某些数字特征,并将其作
4 u4 r/ H2 a$ ^1 Z$ o1 q为问题的解。: B, b! h2 z, Z5 d$ |$ H
0 C. | `; \4 K m/ R有一个例子可以使你比较直观地了解蒙特卡洛方法:
% m- U, d4 q i8 a) ^% }3 I; |9 {假设我们要计算一个不规则图形的面积,那么图形的不规则程度和分析性计算(比如,积分)的复杂程
- Y/ J- U7 v7 F% h9 G度是成正比的。蒙特卡洛方法是怎么计算的呢?假想你有一袋豆子,把豆子均匀地朝这个图形上撒,然, Y* c) M& L5 S
后数这个图形之中有多少颗豆子,这个豆子的数目就是图形的面积。当你的豆子越小,撒的越多的时候
) D" p/ L- v8 [,结果就越精确。4 ]+ Q. ^3 p! W6 c0 L4 c# F
在这里我们要假定豆子都在一个平面上,相互之间没有重叠。
& Y* x# U; l- z/ ?7 k: J3 B
7 y2 g1 I( L' {( _: M& i. n蒙特卡罗方法通过抓住事物运动的几何数量和几何特征,利用数学方法来加以模拟,即进行一种数字模8 Q; x) r7 F; `' K3 z2 u* Q7 t
拟实验。它是以一个概率模型为基础,按照这个模型所描绘的过程,通过模拟实验的结果,作为问题的
/ [; Y) m, I4 f9 @近似解。" n. ^ q0 E0 k4 \6 @ m
6 m- { e5 |9 ?* j; I
蒙特卡罗方法与一般计算方法有很大区别,一般计算方法对于解决多维或因素复杂的问题非常困难,而0 [" Y. K* z; _% A/ I. I
蒙特卡罗方法对于解决这方面的问题却比较简单。其特点如下: 6 ]' b3 G4 N* g
I、 直接追踪粒子,物理思路清晰,易于理解。
" f4 ~' s! L1 @+ QII、 采用随机抽样的方法,较真切的模拟粒子输运的过程,反映了统计涨落的规律。1 \# G7 g1 p8 E" t
III、不受系统多维、多因素等复杂性的限制,是解决复杂系统粒子输运问题的好方法。
4 c% {8 U* ~8 B* L3 }! P等等。% p+ M3 P o7 M% `* p! x8 |
此算法,日后还会在本BLOG 内详细阐述。
$ f. ^4 c( Y7 o. D 2 ~0 W* q( {4 b8 B
- C- C% `8 c% g' r$ \二、数据拟合、参数估计、插值等数据处理算法
7 a( g0 G l9 ^$ k; P我们通常会遇到大量的数据需要处理, 而处理数据的关键就在于这些算法,通常使用Matlab作为工具。7 ]% _$ F+ X% \4 y6 N% x! N
数据拟合在数学建模比赛中中有应用,与图形处理有关的问题很多与拟合有关系,一个例子就是98年数3 v( w9 R- m7 V4 K4 X( J1 {
学建模美国赛A题,生物组织切片的三维插值处理,94年A题逢山开路,山体海拔高度的插值计算,还有
& u6 c- `/ s' G吵的沸沸扬扬可能会考的“非典”问题也要用到数据拟合算法,观察数据的走向进行处理。2 L/ x* t7 T8 G4 a- h8 |& G
2 T# _# h N+ G& z) [) E) q1 ` K5 }
此类问题在 MATLAB 中有很多现成的函数可以调用,熟悉MATLAB,这些方法都能游刃有余的用好。
. O) o7 L1 W( e1 f
2 \" s% j1 p/ a. B3 ?% v
0 _' a T. h8 z- ]三、线性规划、整数规划、多元规划、二次规划等规划类问题' \3 H* i6 ?+ C# j
数学建模竞赛中很多问题都和数学规划有关,可以说不少的模型都可以归结为一组不等式作为约束条件+ n! c# T5 F5 e" H! U0 o" e9 x
、几个函数表达式作为目标函数的问题,遇到这类问题,求解就是关键了,比如98年B题,用很多不等式
$ `( b6 ?# F p完全可以把问题刻画清楚,因此列举出规划后用 Lindo 、 Lingo 等软件来进行解决比较方便,所以还3 X" \7 \- }: Y
需要熟悉这两个软件。
/ p) x. v4 K: o8 Z5 i" b , x3 x& W0 d# k1 X0 {9 [0 _" h6 V
. @% U1 \6 g6 a1 ~- i四、图论算法
: _: m2 J; [# L+ b' F这类问题算法有很多,2 H$ D, K& y' W) k
包括: Dijkstra 、 Floyd 、 Prim 、 Bellman-Ford ,最大流,二分匹配等问题。
$ H7 S( h6 d) n. S % E* R0 F. b' S
关于此类图论算法,可参考Introduction to Algorithms--算法导论,关于图算法的第22章-第26章。
- G+ Z6 G$ A+ W! O同时,本BLOG内经典算法研究系列,对Dijkstra算法有所简单描述,1 M: H- S( z; e, Y# Z
-----------" |. C( u$ _" g3 j8 j& r/ I/ X% i
经典算法研究系列:二、Dijkstra 算法初探; D- r3 Y/ m' E0 J( C5 \6 @2 Z
http://blog.csdn.net/v_JULY_v/archive/2010/12/24/6096981.aspx
" F: } K+ M3 T2 G5 J更多,请关注本BLOG 日后更新的博文。
9 ]# ?# {+ e$ |0 r4 Z9 W
( K$ @3 [# s$ v" N; i! p6 r0 c6 s2 G* ]6 d5 j3 L
五、动态规划、回溯搜索、分治算法、分支定界等计算机算法
3 w' G9 ?8 |- I( m+ O0 {在数学建模竞赛中,如:92 年B题用分枝定界法, 97年B题是典型的动态规划问题,. k( a" {7 ~) r. X: ^, f& P
此外 98 年 B 题体现了分治算法。0 Q, Y1 o$ P/ `% U
! m% {7 h$ ~0 t/ a9 T1 M
这方面问题和 ACM 程序设计竞赛中的问题类似,
) h: {5 q. v1 u推荐看一下算法导论,与《计算机算法设计与分析》(电子工业出版社)等与计算机算法有关的书。+ X( u3 J" f" H @/ ]
# k) K7 D* ?8 [" C( F& T5 H3 c1 F5 W( V& g# o
六、最优化理论的三大经典算法:模拟退火法、神经网络、遗传算法
" e% K1 |6 K% \& y这十几年来最优化理论有了飞速发展,模拟退火法、神经网络、遗传算法这三类算法发展很快。
0 W0 d) V" G9 A* C在数学建模竞赛中:比如97年A题的模拟退火算法,00年B题的神经网络分类算法,01年B题这种难题也可
4 R- G, Q7 e8 _5 e7 C0 n以使用神经网络,还有美国竞赛89年A题也和 BP 算法有关系,当时是86年刚提出BP算法,89年就考了,
l. R7 l$ C: w9 O, w. C说明赛题可能是当今前沿科技的抽象体现。 , K2 N# y2 k6 O7 h; v
03 年 B 题伽马刀问题也是目前研究的课题,目前算法最佳的是遗传算法。, t1 r, F% T! B+ K
3 u- ^4 L" J; Q" z( I1 |. x6 ~/ w另,本人对人工智能非常感兴趣,遗传算法已在本BLOG内有所阐述,敬请参见。+ G) Q1 ?$ P3 z2 M9 Y
----------
4 V7 x. I3 j7 K7 B经典算法研究系列:七、深入浅出遗传算法,透析GA本质( t$ R) S3 l( G( Q/ u0 {0 L
http://blog.csdn.net/v_JULY_v/archive/2011/01/12/6132775.aspx
) y% m' Z+ P+ G1 K$ y ) D# m) a$ Q8 y! [* i& c
其它俩大算法,模拟退火法,与神经网络,也定会在本BLOG内日后的博文更新中,详细阐述。
- K, T7 \4 r% J7 K4 L2 O ! O! {$ d! c3 G2 s5 U0 f( a* Z. Q
, @& c" E0 K }& X/ w
七、网格算法和穷举法# ?% _3 ^5 [: F& E
网格算法和穷举法一样,只是网格法是连续问题的穷举。0 F) N3 L5 v& S" w$ l, \
比如要求在 N 个变量情况下的最优化问题,那么对这些变量可取的空间进行采点,' n, C4 q# M& U" f
比如在 [ a; b ] 区间内取 M +1 个点,就是 a; a +( b ? a ) =M; a +2 ¢ ( b ? a ) =M ; …;b
5 Y5 b+ ^% g8 L# _. ~那么这样循环就需要进行 ( M + 1) N 次运算,所以计算量很大。6 q/ w3 K% b& [; P; G( h. G
( r* s @3 ^' }' h0 d: u) j3 Q在数学建模竞赛中:比如 97 年 A 题、 99 年 B 题都可以用网格法搜索,这种方法最好在运算速度较
; s* P( `5 @" Y( s/ m8 Q) N) A快的计算机中进行,还有要用高级语言来做,最好不要用 MATLAB 做网格,否则会算很久。6 |+ C: [' g" X# L
穷举法大家都熟悉,自不用多说了。 3 y3 ~! S/ D9 Z3 v8 C
a5 @9 K2 A# W( ~; c4 m! x" t5 m. S+ \. d4 C3 d3 L
八、一些连续离散化方法& D$ K- p: L X. t4 T
大部分物理问题的编程解决,都和这种方法有一定的联系。物理问题是反映我们生活在一个连续的世界9 T# N: _5 ^, d* B) {
中,计算机只能处理离散的量,所以需要对连续量进行离散处理。
/ Q8 ^: A. C- `" V
, \$ W2 |( E. I) s% c& Y6 p这种方法应用很广,而且和上面的很多算法有关。
+ y" S! D! k% ^事实上,网格算法、蒙特卡罗算法、模拟退火都用了这个思想。
9 A9 V) c6 k* [+ r) O; F7 c0 H+ l0 D% D, f
1 Z5 S, ^1 t6 l1 Q; c# r9 C1 l- C/ }- s$ o7 G0 t
九、数值分析算法
( U9 H4 z1 `8 S$ z i) X2 [数值分析(numerical analysis),是数学的一个分支,主要研究连续数学(区别于离散数学)问题的
0 k( `* u8 H- w6 |* J* \算法。* j# g0 i- D# r. z6 ^
如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比 如方程组求解、矩阵运算、$ \6 ?/ D6 L/ x% P& Y
函数积分等算法就需要额外编写库函数进行调用。. V( F) p4 U- E9 U# J$ R
这类算法是针对高级语言而专门设的,如果你用的是 MATLAB 、 Mathematica ,大可不必准备,+ a5 T; I2 Z- x
因为像数值分析中有很多函数一般的数学软件是具备的。
1 S" i* R; r, J7 Y- V) i 1 `: P' O- S3 [1 v* M8 F' r
2 v9 r$ R6 {; _. @十、图象处理算法4 a3 J" q# ?) i8 X& b
在数学建模竞赛中:比如01 年 A 题中需要你会读 BMP 图象、美国赛 98 年 A 题需要你知道三维插值8 V% ]; l0 ^( O, _" o5 F: T1 e
计算, 03 年 B 题要求更高,不但需要编程计算还要进行处理,而数模论文中也有很多图片需要展示,
I6 t; A( A- i7 b% Z因此图象处理就是关键。做好这类问题,重要的是把MATLAB 学好,特别是图象处理的部分。
7 ?' b& A' G0 B2 Q0 U # A- J+ ]2 v; N
此数学建模十大算法的程序源码打包,请于此处下载:
, f+ k' F$ M3 S# F5 K" bhttp://download.csdn.net/source/30073363 N' V# F: H" V. ~7 H
/ A6 }, k6 E& s2 @! z5 X8 b! L
本人对算法,尤其感兴趣,且日渐愈浓,* Y9 l9 Q6 U3 P1 K
日后,更多的、好的、经典实用算法将会在本BLOG内有所详细而细致入微的阐述与深入研究。4 R9 z. m- `* y% V. l$ _
完。
; L& D( s H4 y7 d + @) `0 b$ \) Q: p7 ^$ V
7 j. r4 F( r5 A% S% L, i
作者声明:
3 B! c4 H* Z# C; `5 W+ Q6 n0 k; v本人July对本博客所有任何文章、内容和资料享有版权,+ R* x4 b5 d" G% B
转载请注明作者本人July及出处。谢谢。二零一一年一月二十九日。 |
zan
|