+ L |6 U$ a* E( X( \/ D, g4 ^有一个例子可以使你比较直观地了解蒙特卡洛方法: $ a/ q& I- P* {5 V假设我们要计算一个不规则图形的面积,那么图形的不规则程度和分析性计算(比如,积分)的复杂程& }7 m$ |% J' s( R! a* H/ b% T
度是成正比的。蒙特卡洛方法是怎么计算的呢?假想你有一袋豆子,把豆子均匀地朝这个图形上撒,然+ f6 D+ V* D& b! Z
后数这个图形之中有多少颗豆子,这个豆子的数目就是图形的面积。当你的豆子越小,撒的越多的时候 ; y8 P( R2 G9 l, T( Y! u5 f,结果就越精确。$ l4 Q: y, z6 y/ r ^# i; q- |
在这里我们要假定豆子都在一个平面上,相互之间没有重叠。 , b. E `( j; f+ t! F/ X0 G% W6 ^/ h9 X5 _! N/ K
蒙特卡罗方法通过抓住事物运动的几何数量和几何特征,利用数学方法来加以模拟,即进行一种数字模- M9 u- L% X9 I" ~" b Z) R1 o3 }' v
拟实验。它是以一个概率模型为基础,按照这个模型所描绘的过程,通过模拟实验的结果,作为问题的+ k, W5 `4 u4 J& x* A0 C
近似解。 ) n! C6 K0 E* H$ U$ s5 |4 w ; I& s- L% D2 R; m; [; _, X6 z蒙特卡罗方法与一般计算方法有很大区别,一般计算方法对于解决多维或因素复杂的问题非常困难,而 2 P# A1 g6 L1 o6 I4 T4 g3 X- a) M蒙特卡罗方法对于解决这方面的问题却比较简单。其特点如下: / a, e( u7 B; M, M% ?I、 直接追踪粒子,物理思路清晰,易于理解。 / D/ B1 C, f6 o7 X( H# G7 i7 sII、 采用随机抽样的方法,较真切的模拟粒子输运的过程,反映了统计涨落的规律。, g% g2 c) N U6 i
III、不受系统多维、多因素等复杂性的限制,是解决复杂系统粒子输运问题的好方法。$ Z- d4 X2 @+ H6 b
等等。 % h- T* n7 ~1 S# G7 c) F7 i4 E" q此算法,日后还会在本BLOG 内详细阐述。 % L! t3 m. k, h9 @$ j h2 {9 Z8 S! m
" k, ^6 `) F! Y+ u( F& G! v
二、数据拟合、参数估计、插值等数据处理算法 : u0 r+ Y- T& p% q我们通常会遇到大量的数据需要处理, 而处理数据的关键就在于这些算法,通常使用Matlab作为工具。- l$ E4 _' o( p
数据拟合在数学建模比赛中中有应用,与图形处理有关的问题很多与拟合有关系,一个例子就是98年数 & x0 [" X6 ^% X* c6 ]/ r. Q学建模美国赛A题,生物组织切片的三维插值处理,94年A题逢山开路,山体海拔高度的插值计算,还有 9 m7 m+ }; ?$ k4 K) A' o吵的沸沸扬扬可能会考的“非典”问题也要用到数据拟合算法,观察数据的走向进行处理。 ) P- B' k8 e6 a! v* b 4 `: H7 _ S, q N4 h! a
此类问题在 MATLAB 中有很多现成的函数可以调用,熟悉MATLAB,这些方法都能游刃有余的用好。 1 k& o$ H: g: F) n5 G 3 M, V" K- E& j $ t+ K6 c; u9 J! j1 S' b" x, r三、线性规划、整数规划、多元规划、二次规划等规划类问题 9 j7 i6 M. R9 N- p4 A! D/ A数学建模竞赛中很多问题都和数学规划有关,可以说不少的模型都可以归结为一组不等式作为约束条件 / z, b) ?" {; p: D、几个函数表达式作为目标函数的问题,遇到这类问题,求解就是关键了,比如98年B题,用很多不等式8 w; F, }. s! R7 ?3 {0 L
完全可以把问题刻画清楚,因此列举出规划后用 Lindo 、 Lingo 等软件来进行解决比较方便,所以还+ Q% R" _7 D& n8 e1 {) \" Y' O
需要熟悉这两个软件。; W9 h6 v$ C9 p) J6 {% n9 W! Q% z# u
0 Z6 \- F; H+ D H8 }& e8 P# v
. J, ~) |$ @5 ^' z9 E7 c: V
四、图论算法& e6 I+ V# }0 g/ T: y
这类问题算法有很多, 9 s3 a4 X. w y7 q包括: Dijkstra 、 Floyd 、 Prim 、 Bellman-Ford ,最大流,二分匹配等问题。9 E% L @# p" w5 D+ j4 U
6 p7 ]" w3 D* F
关于此类图论算法,可参考Introduction to Algorithms--算法导论,关于图算法的第22章-第26章。 + @2 f' O1 O! H同时,本BLOG内经典算法研究系列,对Dijkstra算法有所简单描述, u0 J1 ~ \$ i/ W. b2 b
-----------/ ?- s: \7 z6 A _
经典算法研究系列:二、Dijkstra 算法初探 % k- M! A: a: {1 Q/ M* ^4 f- Lhttp://blog.csdn.net/v_JULY_v/archive/2010/12/24/6096981.aspx ! R5 k4 ]" t% k, ^# ]2 G更多,请关注本BLOG 日后更新的博文。+ _( M: z/ [1 u, g# n
' K# x6 P6 p% A
" G" C% l0 W- ]0 H
五、动态规划、回溯搜索、分治算法、分支定界等计算机算法; l6 ?: [# ~. `0 S( Q9 B
在数学建模竞赛中,如:92 年B题用分枝定界法, 97年B题是典型的动态规划问题, ; ?! V, B( Z: Z" a' J此外 98 年 B 题体现了分治算法。 3 K7 u; }' _+ S4 p+ h' K( R$ x: ^" R6 y4 ]# r
这方面问题和 ACM 程序设计竞赛中的问题类似, % b! h6 |$ R/ R4 X$ H2 w8 W& Z推荐看一下算法导论,与《计算机算法设计与分析》(电子工业出版社)等与计算机算法有关的书。, H& x8 V& h3 |
) L: G, U3 ~, v1 f! X( N' C- r
. l5 L% G5 Y3 d5 V+ N六、最优化理论的三大经典算法:模拟退火法、神经网络、遗传算法 ; T, Y* c7 Q6 J! P! ~, V/ [
这十几年来最优化理论有了飞速发展,模拟退火法、神经网络、遗传算法这三类算法发展很快。, T% R* ^" D. r9 e
在数学建模竞赛中:比如97年A题的模拟退火算法,00年B题的神经网络分类算法,01年B题这种难题也可; @6 q5 Y( [4 m/ R1 M
以使用神经网络,还有美国竞赛89年A题也和 BP 算法有关系,当时是86年刚提出BP算法,89年就考了,( b" I' O) t$ ~& ~$ s. c8 Y
说明赛题可能是当今前沿科技的抽象体现。 ) ` `9 ^# T, C+ {7 k
03 年 B 题伽马刀问题也是目前研究的课题,目前算法最佳的是遗传算法。 1 _" G1 @5 K' k6 H9 B% R : S$ ~7 j% H( C# f
另,本人对人工智能非常感兴趣,遗传算法已在本BLOG内有所阐述,敬请参见。 . ?) K7 y( D6 }0 C---------- ^4 ?6 {+ H9 E3 e+ }' L经典算法研究系列:七、深入浅出遗传算法,透析GA本质& F4 p h! I E! ?; W1 M http://blog.csdn.net/v_JULY_v/archive/2011/01/12/6132775.aspx. v- F% N: W7 o" O" Z/ \/ P7 F
R3 m) a* W( O2 b/ ]% `$ C
其它俩大算法,模拟退火法,与神经网络,也定会在本BLOG内日后的博文更新中,详细阐述。7 V0 r: T, ]9 s) b% h( _
( r' V/ A; x7 `& j
0 }6 Y8 B0 t6 k+ T9 \" k3 n# o七、网格算法和穷举法. X4 [, ^* O' n
网格算法和穷举法一样,只是网格法是连续问题的穷举。) m) d( x# b$ S2 c0 f) O
比如要求在 N 个变量情况下的最优化问题,那么对这些变量可取的空间进行采点,) c) n4 o$ [# z! g2 C
比如在 [ a; b ] 区间内取 M +1 个点,就是 a; a +( b ? a ) =M; a +2 ¢ ( b ? a ) =M ; …;b ; M( b2 ]+ p3 Y( p那么这样循环就需要进行 ( M + 1) N 次运算,所以计算量很大。# Q! i; O( S$ N; l9 Y% Z
. H8 X. v" n) B+ Z! Y+ d3 M在数学建模竞赛中:比如 97 年 A 题、 99 年 B 题都可以用网格法搜索,这种方法最好在运算速度较( q* N! I8 b$ F- t: @( i) d6 `- _
快的计算机中进行,还有要用高级语言来做,最好不要用 MATLAB 做网格,否则会算很久。( Q" V% f/ H: l
穷举法大家都熟悉,自不用多说了。 % @1 c( L0 f# [1 a" B " _* g5 ^+ X) w5 D5 r! f1 o % y" n, V) U. F& l八、一些连续离散化方法5 |; Q5 B3 U0 k. |! r
大部分物理问题的编程解决,都和这种方法有一定的联系。物理问题是反映我们生活在一个连续的世界 7 @. V$ k. L) k) V |中,计算机只能处理离散的量,所以需要对连续量进行离散处理。# M/ ~9 `% x$ d
& z- `. }' x! @
这种方法应用很广,而且和上面的很多算法有关。3 D" C( V8 z* e. F; s2 o
事实上,网格算法、蒙特卡罗算法、模拟退火都用了这个思想。 2 k+ b# L- n/ k3 d + X! Z6 z9 f- ^: b1 |6 w + B1 @( o3 S1 T U+ m! b$ n F九、数值分析算法 + V' o! N8 X6 p, j* O, P0 [6 M数值分析(numerical analysis),是数学的一个分支,主要研究连续数学(区别于离散数学)问题的8 [8 v7 k: m* d! ~8 N' y
算法。 0 c) ~, z5 J% v如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比 如方程组求解、矩阵运算、8 U+ u8 {5 |' E- Y: |
函数积分等算法就需要额外编写库函数进行调用。% d5 ~. ^7 g! g# F3 ~
这类算法是针对高级语言而专门设的,如果你用的是 MATLAB 、 Mathematica ,大可不必准备,, p) q# K/ C& `) z
因为像数值分析中有很多函数一般的数学软件是具备的。8 B8 y6 \+ V: \) V3 C- N9 q
3 `# W3 K1 R, Q( E% U( ^
2 X/ ` ^) D1 w( b十、图象处理算法 0 c& \$ K s9 `' N6 X; `- F7 B在数学建模竞赛中:比如01 年 A 题中需要你会读 BMP 图象、美国赛 98 年 A 题需要你知道三维插值 u; b, i1 G. q3 U) m) R计算, 03 年 B 题要求更高,不但需要编程计算还要进行处理,而数模论文中也有很多图片需要展示,8 N2 R* f: X# p) L- L
因此图象处理就是关键。做好这类问题,重要的是把MATLAB 学好,特别是图象处理的部分。; W" @0 R. f2 H1 F+ Q: j$ R& J& {