- 在线时间
- 18 小时
- 最后登录
- 2012-9-22
- 注册时间
- 2012-7-20
- 听众数
- 7
- 收听数
- 0
- 能力
- 0 分
- 体力
- 1963 点
- 威望
- 0 点
- 阅读权限
- 50
- 积分
- 643
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 105
- 主题
- 7
- 精华
- 0
- 分享
- 0
- 好友
- 18
升级   10.75% TA的每日心情 | 开心 2012-9-22 23:13 |
---|
签到天数: 44 天 [LV.5]常住居民I
- 自我介绍
- 数学爱好者
 群组: 学术交流A |
数学建模十大算法漫谈
) a1 \* ^! Y7 [' h4 T
5 ^5 E2 K4 E3 V' C* E作者:July 二零一一年一月二十九日# f- _: J) X" ]: a! k
本文参考:. }( a- |, Q$ y1 S1 p
I、 细数二十世纪最伟大的十大算法 [译者:本人July]
* X6 b, K! A+ n1 I% S7 III、 本BLOG内 经典算法研究系列6 z }% |: `6 p( d4 y+ b
III、维基百科$ l7 {3 J0 |% ~ `& Y+ _: o/ r
------------------------------------------& K4 U, g2 \# z7 |; \& ]
博主说明:/ k. D. I% k: ~) {; m6 f5 X
1、此数学建模十大算法依据网上的一份榜单而写,本文对此十大算法作一一简单介绍。- V! D x. l. v9 j) u& o
这只是一份榜单而已,数学建模中还有很多的算法,未一一囊括。欢迎读者提供更多的好的算法。
% M2 `. H; A/ i. D3 A2、在具体阐述每一算法的应用时,除了列出常见的应用之外,7 q8 B' p# [* N& Z1 x
同时,还会具体结合数学建模竞赛一一阐述。
) p5 q6 I9 w' g! ~, T& S- B+ L Y毕竟,此十大算法,在数学建模竞赛中有着无比广泛而重要的应用。
' ^8 v6 E8 [, d/ G+ k/ X且,凡是标着“某某年某国某题”,即是那一年某个国家的数学建模竞赛原题。
1 b$ E( r0 ?; X. Y" A/ a. v3、此十大算法,在一些经典的算法设计书籍上,无过多阐述。
/ u" u& p0 k# D% l9 f) x5 r若要具体细致的深入研究,还得请参考国内或国际上关于此十大算法的优秀论文。
- A1 T2 F: A6 k, E! {6 y4 v7 L谢谢。
9 X/ i" ~$ P' o
4 F; ?, O1 P' h7 g% C9 \3 f) M" j 9 l; r: \7 }% g
一、蒙特卡罗算法2 H/ j! M/ Q8 z2 x# V) U3 A# `
1946年,美国拉斯阿莫斯国家实验室的三位科学家John von Neumann,Stan Ulam 和 Nick Metropolis
3 B2 ?* @ B- ^# ?2 R( W共同发明了,蒙特卡罗方法。2 n. L4 k+ h* q. S$ E
+ E8 `6 W" r9 N+ f此算法被评为20世纪最伟大的十大算法之一,详情,请参见我的博文:4 {4 o% e8 `% E5 b! [/ i( Q
http://blog.csdn.net/v_JULY_v/archive/2011/01/10/6127953.aspx1 f# ]4 K" A# ?. M8 Z( G9 o
% _6 H3 X* n0 U! @) P蒙特卡罗方法(Monte Carlo method),又称随机抽样或统计模拟方法,是一种以概率统计理论为指导
2 C7 I( X& X3 b" h$ i的一类非常重要的数值计算方法。此方法使用随机数(或更常见的伪随机数)来解决很多计算问题的方1 r7 N7 p! K+ C& r
法。
2 t/ F$ c; I+ P " {9 C$ R! M8 H/ |* A. P" k! Z
由于传统的经验方法由于不能逼近真实的物理过程,很难得到满意的结果,而蒙特卡罗方法由于能够真
& T& X+ S. ~8 Q7 m* S实地模拟实际物理过程,故解决问题与实际非常符合,可以得到很圆满的结果。# J2 |2 E; C/ ?
蒙特卡罗方法的基本原理及思想如下:
; R0 R0 j H( D" F+ x当所求解问题是某种随机事件出现的概率,或者是某个随机变量的期望值时,通过某种“实验”的方法
l9 c* n/ F/ _& M2 J. n! ?,以这种事件出现的频率估计这一随机事件的概率,或者得到这个随机变量的某些数字特征,并将其作
& Y: Q' f5 j8 ~# F' e L为问题的解。( i. s3 R3 L9 E3 I" X7 K
; t4 @3 J) p* ^, m( A. C# U+ t7 I, H有一个例子可以使你比较直观地了解蒙特卡洛方法:: h9 N D; l2 z2 ]0 Z1 z; ?
假设我们要计算一个不规则图形的面积,那么图形的不规则程度和分析性计算(比如,积分)的复杂程7 z- x# C1 z* M% {) R/ ~0 q
度是成正比的。蒙特卡洛方法是怎么计算的呢?假想你有一袋豆子,把豆子均匀地朝这个图形上撒,然- K3 @; g: p" j. x# ~# a7 p
后数这个图形之中有多少颗豆子,这个豆子的数目就是图形的面积。当你的豆子越小,撒的越多的时候# Q" _8 T9 D" I5 {/ k+ _+ G+ b
,结果就越精确。9 A# ]; C6 Q4 `! h) o+ N
在这里我们要假定豆子都在一个平面上,相互之间没有重叠。
4 A# o4 o F5 R6 {! Z
) E4 i! l5 c4 S, B% D1 W; w蒙特卡罗方法通过抓住事物运动的几何数量和几何特征,利用数学方法来加以模拟,即进行一种数字模) w! Q# }2 e( ]8 {
拟实验。它是以一个概率模型为基础,按照这个模型所描绘的过程,通过模拟实验的结果,作为问题的: l- J" ^8 Z& ~$ K' y+ X
近似解。
& a. C$ [7 C# i2 Q6 ^( [: }+ z ) r4 ]- S$ ^6 i8 n1 M
蒙特卡罗方法与一般计算方法有很大区别,一般计算方法对于解决多维或因素复杂的问题非常困难,而! f- K/ ~. b6 r& R- Y
蒙特卡罗方法对于解决这方面的问题却比较简单。其特点如下:
" i `+ G6 V f h) Q0 l B6 NI、 直接追踪粒子,物理思路清晰,易于理解。 / B5 R; b9 T) H& U3 @) |/ V) a" n
II、 采用随机抽样的方法,较真切的模拟粒子输运的过程,反映了统计涨落的规律。
5 H. z$ [3 e8 @+ rIII、不受系统多维、多因素等复杂性的限制,是解决复杂系统粒子输运问题的好方法。) N" k4 G& @9 d+ i- j, i) M K0 j. h
等等。. V" R% B! \) e: w* }
此算法,日后还会在本BLOG 内详细阐述。
+ ^6 u) i7 A( j7 m
$ k3 v! z6 M9 F0 {' z( g, B4 C! W; \9 N$ [
二、数据拟合、参数估计、插值等数据处理算法
G" t2 y. ~0 I8 Z a: c我们通常会遇到大量的数据需要处理, 而处理数据的关键就在于这些算法,通常使用Matlab作为工具。
, ^: Z M, f* q. V3 g/ Q数据拟合在数学建模比赛中中有应用,与图形处理有关的问题很多与拟合有关系,一个例子就是98年数+ f5 h. `% j. `% r$ `. X: C& K
学建模美国赛A题,生物组织切片的三维插值处理,94年A题逢山开路,山体海拔高度的插值计算,还有
/ _1 b' [) E1 P8 g/ ]7 \吵的沸沸扬扬可能会考的“非典”问题也要用到数据拟合算法,观察数据的走向进行处理。; @" m! T( n- K
: W5 r4 H, S% f2 k/ M$ G
此类问题在 MATLAB 中有很多现成的函数可以调用,熟悉MATLAB,这些方法都能游刃有余的用好。
: Z/ ^' H n- R2 [( F7 z
- Q9 h2 o4 A k0 Q8 E! T- f6 v- L! n. F. l+ `* N* I; P: S8 ~
三、线性规划、整数规划、多元规划、二次规划等规划类问题
' h2 c- ~4 E5 G; R! t数学建模竞赛中很多问题都和数学规划有关,可以说不少的模型都可以归结为一组不等式作为约束条件) G/ m8 }+ p9 T; i5 Y' p
、几个函数表达式作为目标函数的问题,遇到这类问题,求解就是关键了,比如98年B题,用很多不等式) i. O5 b. t/ F6 _
完全可以把问题刻画清楚,因此列举出规划后用 Lindo 、 Lingo 等软件来进行解决比较方便,所以还
* A( M2 ^# K. [4 h9 H需要熟悉这两个软件。
* Y9 o0 q5 [6 R
+ _ A" I8 p1 r9 A0 Z8 D" h6 o1 [ F9 o$ G5 m
四、图论算法
# s U' n- U5 o& U: l0 k) N这类问题算法有很多,
. {$ v& y0 {& _) D包括: Dijkstra 、 Floyd 、 Prim 、 Bellman-Ford ,最大流,二分匹配等问题。2 ?. y) a# f' S$ r+ `. |
- R' b _& X) E& v5 n% A关于此类图论算法,可参考Introduction to Algorithms--算法导论,关于图算法的第22章-第26章。2 H$ ^: F2 s* \6 s+ y3 l& Q
同时,本BLOG内经典算法研究系列,对Dijkstra算法有所简单描述,- K* l3 ^% F7 u
-----------
5 p1 P9 y0 {: Y, D0 Y6 u' k8 N经典算法研究系列:二、Dijkstra 算法初探$ N- S+ X. c' H! j* g4 [
http://blog.csdn.net/v_JULY_v/archive/2010/12/24/6096981.aspx
! Q8 {: C" a# P. |8 c9 n( e! a更多,请关注本BLOG 日后更新的博文。
( H" ?! ]7 A% p6 d, }$ c. [ * l; x/ A& |0 w4 u8 ~
3 Q& d. y3 }1 ^" v: q五、动态规划、回溯搜索、分治算法、分支定界等计算机算法3 p( W$ X9 M$ u0 ]6 S) ], V" {
在数学建模竞赛中,如:92 年B题用分枝定界法, 97年B题是典型的动态规划问题,
3 k, P' l& {7 X6 J) ^1 Q此外 98 年 B 题体现了分治算法。
4 O: ^& @( _) w1 i
2 @( e+ ]$ D+ h3 i这方面问题和 ACM 程序设计竞赛中的问题类似,
5 \/ H; }6 p9 K* I推荐看一下算法导论,与《计算机算法设计与分析》(电子工业出版社)等与计算机算法有关的书。6 i6 r9 e- d+ V& z$ S/ z
e& K. n4 C J9 a) @' W% k; l8 d3 c5 E/ Q9 I
六、最优化理论的三大经典算法:模拟退火法、神经网络、遗传算法 6 \: g. l. o- |: n: @. e. B
这十几年来最优化理论有了飞速发展,模拟退火法、神经网络、遗传算法这三类算法发展很快。
+ C7 ^; l" D0 Y3 a: y在数学建模竞赛中:比如97年A题的模拟退火算法,00年B题的神经网络分类算法,01年B题这种难题也可
! ]1 u0 c" x4 A/ [" s5 B以使用神经网络,还有美国竞赛89年A题也和 BP 算法有关系,当时是86年刚提出BP算法,89年就考了,& J! @: i, i( @
说明赛题可能是当今前沿科技的抽象体现。
* i" i% }* E$ J/ i3 V03 年 B 题伽马刀问题也是目前研究的课题,目前算法最佳的是遗传算法。( y$ s4 f1 X& K( ^
1 ]8 C7 B" G& F5 H( P4 n另,本人对人工智能非常感兴趣,遗传算法已在本BLOG内有所阐述,敬请参见。
9 l6 r/ Q @% i! |! D2 _8 O! y----------' D* k8 ?( U3 C
经典算法研究系列:七、深入浅出遗传算法,透析GA本质. l( J( j8 x' B
http://blog.csdn.net/v_JULY_v/archive/2011/01/12/6132775.aspx
& J4 y" \. U9 P+ }0 t( y+ x
$ g' ?' ]( o' f$ F }: r' ?, `其它俩大算法,模拟退火法,与神经网络,也定会在本BLOG内日后的博文更新中,详细阐述。
9 D, W3 N+ F9 I3 i: q: z 0 e- j/ K6 u+ ?9 Y
# G& c! y& c- g0 K. F' P- ]/ h七、网格算法和穷举法
- D. ?/ B2 f: O% M! C( @# y+ o网格算法和穷举法一样,只是网格法是连续问题的穷举。
( a4 z3 |% G2 ~: J% F* l' q比如要求在 N 个变量情况下的最优化问题,那么对这些变量可取的空间进行采点,
. l6 q! `0 P5 k8 i p! d/ Y1 w比如在 [ a; b ] 区间内取 M +1 个点,就是 a; a +( b ? a ) =M; a +2 ¢ ( b ? a ) =M ; …;b
* @: r1 z3 D. l( a" J那么这样循环就需要进行 ( M + 1) N 次运算,所以计算量很大。
0 e7 M+ ]; {3 X; v6 }, k
* \4 K u* x0 C' Z' J在数学建模竞赛中:比如 97 年 A 题、 99 年 B 题都可以用网格法搜索,这种方法最好在运算速度较
* R/ Z3 X" |" ^0 ]& T" m快的计算机中进行,还有要用高级语言来做,最好不要用 MATLAB 做网格,否则会算很久。( a" }3 z% |0 c' @5 u
穷举法大家都熟悉,自不用多说了。 4 {- N2 D4 A" j% h7 Q* `1 b) M
9 x% L' _- G5 U0 n8 j1 c3 U( P% K2 J8 ]
- f8 |9 b6 l. f- b# X, z1 B
八、一些连续离散化方法9 X) T' v' b7 m$ O9 @8 D' u1 |
大部分物理问题的编程解决,都和这种方法有一定的联系。物理问题是反映我们生活在一个连续的世界
# I( H9 `; v3 n中,计算机只能处理离散的量,所以需要对连续量进行离散处理。
7 }) ~; W7 P, Y/ \- i: G
! J4 o5 L0 S; j% l这种方法应用很广,而且和上面的很多算法有关。
: q$ S9 G& j: R0 r1 q事实上,网格算法、蒙特卡罗算法、模拟退火都用了这个思想。 $ R6 f3 i2 u6 M" H! C
8 x2 Z8 Z; I, C. h; L, F
% j) s8 O. k& f! ^, j$ h% m
九、数值分析算法
5 G& E& v3 O; I/ A+ z, L/ F1 ]数值分析(numerical analysis),是数学的一个分支,主要研究连续数学(区别于离散数学)问题的
]) L: e! W: D% J) h; N算法。
' p2 V# E% K. Z1 N/ J如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比 如方程组求解、矩阵运算、8 G" _2 ]- s! }* U: G/ E6 H& b
函数积分等算法就需要额外编写库函数进行调用。8 a1 K' t2 S; M6 Q5 A0 N
这类算法是针对高级语言而专门设的,如果你用的是 MATLAB 、 Mathematica ,大可不必准备,
6 ?7 o' b- t' ?" \因为像数值分析中有很多函数一般的数学软件是具备的。+ S, A7 y0 r1 m) U5 S# N3 k) m% J2 y& F
5 @* |( L6 d0 U2 a# V' ]2 o$ a2 L$ y3 Y& S
十、图象处理算法- C# a9 Q7 K: d: D2 Y7 Q/ f! ` c0 e* z6 O
在数学建模竞赛中:比如01 年 A 题中需要你会读 BMP 图象、美国赛 98 年 A 题需要你知道三维插值
h( `" s0 V0 k* z* x计算, 03 年 B 题要求更高,不但需要编程计算还要进行处理,而数模论文中也有很多图片需要展示,
' j0 T: l6 ^9 }9 |! X因此图象处理就是关键。做好这类问题,重要的是把MATLAB 学好,特别是图象处理的部分。; D% i1 o. X4 ^: ?
9 |6 w7 E- E6 Y1 O" ~0 l& A- K8 j
此数学建模十大算法的程序源码打包,请于此处下载:. Y7 `- t' S/ K- f# p9 K
http://download.csdn.net/source/30073361 S0 }$ Q( ^) i0 l5 w+ z6 Z! O# E' S
* L, o3 i1 [, ?0 q4 V& t9 u! `9 u
本人对算法,尤其感兴趣,且日渐愈浓, H% Q$ [1 z4 z) w1 p1 t
日后,更多的、好的、经典实用算法将会在本BLOG内有所详细而细致入微的阐述与深入研究。
& I2 h/ u# \3 q) @9 o% W+ F& ?4 L/ f完。9 l8 ^& B9 A5 v! }, _
- O: Z# |- {' \2 N( t1 P# c
% v$ {/ B; \3 h/ r作者声明:
/ \& W% W+ X6 h本人July对本博客所有任何文章、内容和资料享有版权,
' A F% d" n8 x转载请注明作者本人July及出处。谢谢。二零一一年一月二十九日。 |
zan
|