" }$ c) |& i1 C$ s! D* l实地模拟实际物理过程,故解决问题与实际非常符合,可以得到很圆满的结果。 ' C- ]$ R0 u, d W ; F4 o Y* n( n: x蒙特卡罗方法的基本原理及思想如下: & r+ j2 n |& W% z当所求解问题是某种随机事件出现的概率,或者是某个随机变量的期望值时,通过某种“实验”的方法$ I* @# M8 b, s! \
f* R: k8 L* z% d1 G* O. M,以这种事件出现的频率估计这一随机事件的概率,或者得到这个随机变量的某些数字特征,并将其作' E {$ G( x0 ]" P# S$ ?5 E0 [
% @* i U. B9 }7 q5 c+ ?) I$ x" {
为问题的解。( N8 T# K% \; L6 \- u
6 V6 |; o) `" a, e& V9 ^ 3 g1 P- d$ G0 [# V; c: o- k/ C4 l5 t1 q
有一个例子可以使你比较直观地了解蒙特卡洛方法: - A9 I; Z; e3 R' h% L8 P- O) T' D假设我们要计算一个不规则图形的面积,那么图形的不规则程度和分析性计算(比如,积分)的复杂程 9 G1 {' u/ j( ]( X$ C! R ~ : b! g1 i8 r& {+ w度是成正比的。蒙特卡洛方法是怎么计算的呢?假想你有一袋豆子,把豆子均匀地朝这个图形上撒,然 : y1 b' @) {, _2 N- J2 j) L4 _; c7 {8 Q5 C' b9 _
后数这个图形之中有多少颗豆子,这个豆子的数目就是图形的面积。当你的豆子越小,撒的越多的时候 3 G6 n2 U0 g! i$ [! s) t 4 R% ~; ^' Q% E* u+ q0 j,结果就越精确。 4 M4 h0 ], K2 c4 u在这里我们要假定豆子都在一个平面上,相互之间没有重叠。. y/ }$ a3 ~$ @9 k g
$ G, |3 l& J+ K
; s0 C2 |7 l0 S
蒙特卡罗方法通过抓住事物运动的几何数量和几何特征,利用数学方法来加以模拟,即进行一种数字模 $ ]- A1 E1 F) a4 ^. E & e" T4 k0 O; K" J5 o拟实验。它是以一个概率模型为基础,按照这个模型所描绘的过程,通过模拟实验的结果,作为问题的 " n2 b* H7 [# a1 O7 F& j9 V1 g: s: Y
近似解。: R1 y5 Q! M$ e6 E, s8 h- a
! c" M8 y- a- s# b
/ B* O6 N& o0 F! _
+ o3 S' f. ~) \' X' I7 t T
蒙特卡罗方法与一般计算方法有很大区别,一般计算方法对于解决多维或因素复杂的问题非常困难,而6 L" [9 v$ p4 X4 r
4 {* \0 M, _( N$ Z/ u: c$ v8 T
蒙特卡罗方法对于解决这方面的问题却比较简单。其特点如下: 2 o9 o( Z0 o {( m1 ]0 a) ZI、 直接追踪粒子,物理思路清晰,易于理解。 3 Q6 S( |/ S2 ]
II、 采用随机抽样的方法,较真切的模拟粒子输运的过程,反映了统计涨落的规律。' y6 x: u4 e" E7 C
III、不受系统多维、多因素等复杂性的限制,是解决复杂系统粒子输运问题的好方法。3 Q3 t a5 e) O
等等。 9 r$ q' T, ?6 o* G 8 w5 U5 u2 O. d此算法,日后还会在本BLOG 内详细阐述。 , {3 S! K6 q4 J* f 2 N2 {5 @& b! E8 t6 e- n; k) A% t - W) u% w" v( Z% Z$ E) ~1 R/ o # r7 O* j' H- w0 p , V) i! K/ t. N二、数据拟合、参数估计、插值等数据处理算法 ; F: p3 s( W! I( r) |我们通常会遇到大量的数据需要处理, 而处理数据的关键就在于这些算法,通常使用Matlab作为工具。 ' B4 j5 W5 W" ~7 O- t) Q+ q' S( [" l$ @# c
数据拟合在数学建模比赛中中有应用,与图形处理有关的问题很多与拟合有关系,一个例子就是98年数 1 ~( l3 N1 e& O# K* ?/ v% d + l1 s+ \9 P2 i U5 |0 M0 f% P# W* j学建模美国赛A题,生物组织切片的三维插值处理,94年A题逢山开路,山体海拔高度的插值计算,还有 0 i( V. g( D' L+ }2 k# [; ]' Q3 }0 u; X+ r5 i ]
吵的沸沸扬扬可能会考的“非典”问题也要用到数据拟合算法,观察数据的走向进行处理。 " z0 w" {' ?3 H. q! L, p9 }; U6 u7 q" _
: W: w Q3 ]/ u: r+ ]3 f8 p$ S
+ J/ I" v) {4 a! b5 C3 S
此类问题在 MATLAB 中有很多现成的函数可以调用,熟悉MATLAB,这些方法都能游刃有余的用好。5 ~* Z4 \; A i. X0 F
" X$ S2 \1 M3 W+ j7 _
0 }& V4 }3 @8 R
2 g9 r- o( K" u 4 C: U3 x5 j& x$ _' ?三、线性规划、整数规划、多元规划、二次规划等规划类问题# p/ J7 s% H* a$ }. c& p
数学建模竞赛中很多问题都和数学规划有关,可以说不少的模型都可以归结为一组不等式作为约束条件; V. S: |, ] [) q
6 u7 }. }6 \" ], O" A0 z$ C7 Y
、几个函数表达式作为目标函数的问题,遇到这类问题,求解就是关键了,比如98年B题,用很多不等式( b1 m" Q+ W0 M J: A ~2 [/ x
% {6 c& b* n# q! o- B
完全可以把问题刻画清楚,因此列举出规划后用 Lindo 、 Lingo 等软件来进行解决比较方便,所以还3 e' V) r1 S4 L: V5 z. O( O
* ?7 K0 F* \" n( @1 h9 V1 I
需要熟悉这两个软件。 6 [8 i$ r2 f& g+ _( O- ~, R! n4 `$ }6 P1 P y. V# `( M& G
; d6 }* ~& X" `" t: d$ o) k - T- j$ w% x Z4 ]9 \ , C$ e: z/ M) a. y6 l6 I2 G- P四、图论算法5 q2 i0 X% H+ U6 R+ e9 o
这类问题算法有很多,; p, M5 ]; I8 ]6 _& b6 m( o2 G
包括: Dijkstra 、 Floyd 、 Prim 、 Bellman-Ford ,最大流,二分匹配等问题。& y6 n' f8 x4 @7 p, H1 j
$ e8 g' v6 s2 T) W 5 B. a4 g- Z. @8 ~4 B% ^- b8 E1 C u 4 v6 d- {0 ^. a$ T) V6 d关于此类图论算法,可参考Introduction to Algorithms--算法导论,关于图算法的第22章-第26章。 ! {7 z) V" q* J$ L同时,本BLOG内经典算法研究系列,对Dijkstra算法有所简单描述, , x, q) Q$ N* s4 Z-----------* G4 `. o& p4 D/ J* E8 V& Z4 F
经典算法研究系列:二、Dijkstra 算法初探 , X; ]9 i' W) k, Ihttp://blog.csdn.net/v_JULY_v/archive/2010/12/24/6096981.aspx 4 u! `: G8 g4 F7 ]0 E( [ 1 u9 R6 u7 p! ]2 I更多,请关注本BLOG 日后更新的博文。 5 P' m- c, O5 d0 O$ C+ |9 Q! a 2 {' K* A- u4 E. W: |* ^$ V5 V* g) o) u) Z3 N% {1 D
8 X7 ~6 D2 R5 |# r, }; ?4 W( I
' M T2 c4 v9 M% ^ h2 ~
五、动态规划、回溯搜索、分治算法、分支定界等计算机算法8 k* K' e+ o( L$ J# h
在数学建模竞赛中,如:92 年B题用分枝定界法, 97年B题是典型的动态规划问题, ) m) w2 Z; }% G此外 98 年 B 题体现了分治算法。$ y P: v& N; E/ `9 l. \4 H
* l2 j3 c. C7 w* Z& @( y7 d& S( x : p; l: ~) d7 W3 w$ L' P. O这方面问题和 ACM 程序设计竞赛中的问题类似, 6 y: c* L/ F+ I2 I推荐看一下算法导论,与《计算机算法设计与分析》(电子工业出版社)等与计算机算法有关的书。3 ^) A5 N3 |3 H
1 B, ]% j- F8 C2 Q4 w P ; k1 E9 m) N; Q2 H2 @( `6 g, W: I0 Y4 m- o% o
# C% D* G6 d9 F六、最优化理论的三大经典算法:模拟退火法、神经网络、遗传算法 & r$ k- x, z5 t这十几年来最优化理论有了飞速发展,模拟退火法、神经网络、遗传算法这三类算法发展很快。 9 @# C* I3 z; X ; \" O1 X4 P6 M4 F: a0 P在数学建模竞赛中:比如97年A题的模拟退火算法,00年B题的神经网络分类算法,01年B题这种难题也可 - Y4 T- C4 v, V3 i 2 @, g# S8 C" O以使用神经网络,还有美国竞赛89年A题也和 BP 算法有关系,当时是86年刚提出BP算法,89年就考了, % i8 P2 c) I+ ^, u/ ?0 [' y# J' L# V4 E! ^6 D+ g. B1 \
说明赛题可能是当今前沿科技的抽象体现。 . Q+ N/ G* b# J( b5 i* b5 L6 h03 年 B 题伽马刀问题也是目前研究的课题,目前算法最佳的是遗传算法。 9 b+ S3 F9 z& m+ H; p ) w+ b. e' @2 H% e$ t$ T9 k1 H* }9 Y9 W& ?
/ Z# H1 u# l/ ? F7 ?, _ {
另,本人对人工智能非常感兴趣,遗传算法已在本BLOG内有所阐述,敬请参见。$ U% Z6 ?# a, w$ n: a2 h
---------- - H4 h/ N( L H$ ^( [经典算法研究系列:七、深入浅出遗传算法,透析GA本质9 d8 W; Z/ l5 [2 V/ R7 F6 W
http://blog.csdn.net/v_JULY_v/archive/2011/01/12/6132775.aspx * y# J- r _& v+ R2 H7 d. j8 c 3 Q" f5 P5 m* N4 ^! h6 }3 ] % q' U- k+ f: G4 a. ~ ( ~1 Z( h4 R0 P/ G0 J其它俩大算法,模拟退火法,与神经网络,也定会在本BLOG内日后的博文更新中,详细阐述。) s- t4 h& P+ B( D
! c- g: L. Z! A! x
6 A1 E4 N: R) G* O" N
* T; |3 h2 c9 W) H7 [) z) o- H. f o9 s! o& P- N t
七、网格算法和穷举法 D! N2 ?" v+ [' Y3 C; P& c$ e* N
网格算法和穷举法一样,只是网格法是连续问题的穷举。) J2 o6 Y0 Y5 t9 ~: v
比如要求在 N 个变量情况下的最优化问题,那么对这些变量可取的空间进行采点,/ T2 ~; a1 p# I/ G% t
比如在 [ a; b ] 区间内取 M +1 个点,就是 a; a +( b ? a ) =M; a +2 ¢ ( b ? a ) =M ; …;b 4 F5 F- u% i/ @3 T4 {- h4 ?3 t 6 t7 T M/ J7 R+ }, h! e那么这样循环就需要进行 ( M + 1) N 次运算,所以计算量很大。 . T1 g2 H" G4 x- I+ F8 N# j5 z* H6 c4 `& H
% o" `; K9 A- J4 f' S& f
在数学建模竞赛中:比如 97 年 A 题、 99 年 B 题都可以用网格法搜索,这种方法最好在运算速度较: q$ R* q& u, a H) N: }8 I4 H2 @
: D% w9 T. \& ]: Z4 X7 |6 d快的计算机中进行,还有要用高级语言来做,最好不要用 MATLAB 做网格,否则会算很久。" T# ]$ W9 x8 U: k% g* V
- |5 ^5 v4 g& C2 A, T- ^% c穷举法大家都熟悉,自不用多说了。 * z& d5 x: Q8 a! X; e + o& W% s* H% _8 Q' \) P* a. s) Z9 h1 X. |4 ~3 }0 R, x
6 m2 F" |4 B" \9 [" S* C: [3 M
! Q0 e! d8 c# o( w$ ?八、一些连续离散化方法- `3 e6 G5 O+ F6 e" N6 s
大部分物理问题的编程解决,都和这种方法有一定的联系。物理问题是反映我们生活在一个连续的世界& Z! `$ i3 t1 z9 n- ~; ^4 A# r0 D
+ K' T& [( J! W% S! _1 }! Z" A
中,计算机只能处理离散的量,所以需要对连续量进行离散处理。 - E$ X0 y6 P( m% k0 | - V; L5 w9 x' K& o9 ]3 x' H/ S 7 ?3 ~( d. r1 U8 ^& Q这种方法应用很广,而且和上面的很多算法有关。 ( X, S, v3 R/ \- }* @" T事实上,网格算法、蒙特卡罗算法、模拟退火都用了这个思想。 + @' E" i* e# ^$ Y _& c
0 ~5 O2 {5 o4 i, F5 q% u( H: e3 v
( v, G; {. j, p/ O2 y8 ^; Q- E$ I+ K7 D) p2 E
# t) A8 U" ~9 E2 W九、数值分析算法4 F3 ?6 h. I- g/ C% z
数值分析(numerical analysis),是数学的一个分支,主要研究连续数学(区别于离散数学)问题的 . y, D. i, i, _+ q7 D: v/ ?4 g/ t+ C9 M* v' e$ L
算法。* f: p A9 b9 J; O
6 S7 @+ [3 ^1 D4 t如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比 如方程组求解、矩阵运算、 $ T# a3 g6 N! X U% x # J6 u" b) ]6 B5 K7 G函数积分等算法就需要额外编写库函数进行调用。 9 U6 Y) C* k2 b3 g$ j $ x/ j5 E' u$ `+ c2 Q2 k这类算法是针对高级语言而专门设的,如果你用的是 MATLAB 、 Mathematica ,大可不必准备,3 _" |. J- o, g
因为像数值分析中有很多函数一般的数学软件是具备的。1 q- A3 V; }8 ?; _1 j7 W
% r& ]+ C0 X" s3 n
: A$ s. z+ d' f3 V( e6 F% r* C% k
' W- l* v/ d& c1 d9 `% t8 [十、图象处理算法 # U9 ~) R: s! g' Y$ {' g在数学建模竞赛中:比如01 年 A 题中需要你会读 BMP 图象、美国赛 98 年 A 题需要你知道三维插值5 X, j5 o' h& m3 i) x4 M
7 ^& |8 _- O( Q0 l
计算, 03 年 B 题要求更高,不但需要编程计算还要进行处理,而数模论文中也有很多图片需要展示, # A% i: T% @- Q; M) [+ _+ ^( [' |1 }" q9 H0 t/ s G
因此图象处理就是关键。做好这类问题,重要的是把MATLAB 学好,特别是图象处理的部分。 3 T$ i' h0 N# ^. p( E4 v. M' w$ V! X7 v) ]