! r6 Y) f" Y. }6 ^) s- ]数学建模比赛是本科生和研究生阶段最重要的比赛之一,包括全国大学生数学建模竞赛(俗称“国赛”)和美国大学生数学建模竞赛(俗称“美赛”)。在这些比赛中取得好成绩,不仅有助于保研、有助于找工作,更重要的是形成科学的思维模式。下面列举了十大算法,在数学建模竞赛中有着无比广泛而重要的应用。4 Z8 r( q' w) Y; }8 p
$ f- Q9 s G* c+ d: E' c/ H 01 - J* q- X" R& D& {
& }2 ^, {( j. p2 z* Y7 o
蒙特卡罗算法 ! B1 l$ R$ w4 h3 Q B$ P
1946 年,美国拉斯阿莫斯国家实验室的三位科学家 JohnvonNeumann, Stan Ulam和 Nick Metropolis 共同发明了蒙特卡罗方法。
# { b6 g9 M* D. _蒙特卡罗方法(Monte Carlo method),又称随机抽样或统计模拟方法,是一种以概率统计理论为指导的一类非常重要的数值计算方法。此方法使用随机数(或更常见的伪随机数)来解决很多计算问题的方法。
3 h! r7 \% e$ K" B( z6 f% R" E$ H: n C
由于传统的经验方法由于不能逼近真实的物理过程,很难得到满意的结果,而蒙特卡罗方法由于能够真实地模拟实际物理过程,故解决问题与实际非常符合,可以得到很圆满的结果。
: `4 F% O6 e# I7 }3 r- t
& d# l) U" v8 j; A9 v蒙特卡罗方法的基本原理及思想如下:
) o! S: j1 j, f# o
1 R3 T5 o% N% E' Y3 c当所求解问题是某种随机事件出现的概率,或者是某个随机变量的期望值时,通过某种“实验”的方法,以这种事件出现的频率估计这一随机事件的概率,或者得到这个随机变量的某些数字特征,并将其作为问题的解。
. a" o0 d7 P) w# `8 T) n) R4 Q6 i: T" E& P9 I! n
举个栗子,直观了解蒙特卡洛方法:
1 v# G6 y5 Q9 Z3 I& B# ]" l
4 V; _7 A5 `" u! P ?* I" w0 Q; D假设我们要计算一个不规则图形的面积,那么图形的不规则程度和分析性计算(比如:积分)的复杂程度是成正比的。蒙特卡洛方法是怎么计算的呢?假想你有一袋豆子,把豆子均匀地朝这个图形上撒,然后数这个图形之中有多少颗豆子,这个豆子的数目就是图形的面积。当你的豆子越小,撒的越多的时候,结果就越精确。在这里我们要假定豆子都在一个平面上,相互之间没有重叠。
" g# p% }' v& P7 v9 Q
7 V: U' P7 x. ^% L i7 ?4 K9 f) v$ c+ C7 a$ \
蒙特卡罗方法通过抓住事物运动的几何数量和几何特征,利用数学方法来加以模拟,即进行一种数字模拟实验。它是以一个概率模型为基础,按照这个模型所描绘的过程,通过模拟实验的结果,作为问题的近似解。
% _: w$ ?; I9 g9 u4 a: u" j
/ Z# s# R+ O! o* I蒙特卡罗方法与一般计算方法有很大区别,一般计算方法对于解决多维或因素复杂的问题非常困难,而蒙特卡罗方法对于解决这方面的问题却比较简单。其特点如下:
/ l# E+ `+ I1 X+ f7 a
* Q1 G% @5 }& w- ]a、直接追踪粒子,物理思路清晰,易于理解;
' x5 m4 A8 K4 p+ e! {2 v( q5 i' R3 v5 J! F4 x, m( c7 a
b、采用随机抽样的方法,较真切的模拟粒子输运的过程,反映了统计涨落的规律;. y. f# `2 z; F" G/ W
# e0 M" {% t2 ^# j0 h mc、不受系统多维、多因素等复杂性的限制,是解决复杂系统粒子输运问题的好方法等等
3 l6 x6 ^" T6 P
, d t3 S. z ^* T# [4 m% m$ A 02
# J2 c2 V' R2 t8 H% g' l$ v, D* M9 p& y
数据拟合、参数估计、插值等数据处理算法 ' u1 L k( c( v# l# U1 m3 L1 O
. B. B' Q3 r. s* z$ k. n3 W" P& x O" {# Z
我们通常会遇到大量的数据需要处理,而处理数据的关键就在于这些算法,通常使用 Matlab 作为工具。/ |* j5 Z( S8 N1 Y) a
. d( `, A. k% g) R# h& b) @9 [
数据拟合在数学建模比赛中中有应用,与图形处理有关的问题很多与拟合有关系,一个例子就是 98 年数学建模美国赛 A 题,生物组织切片的三维插值处理,94 年 A 题逢山开路,山体海拔高度的插值计算,还有吵的沸沸扬扬可能会考的“非典”问题也要用到数据拟合算法,观察数据的走向进行处理。( S5 l( H/ T. H+ g# ~$ C6 K! e
% ]" _ t+ `$ g$ L* D
. E3 g v* G' g" }' ]! R. u% ]
此类问题在 MATLAB 中有很多现成的函数可以调用,熟悉 MATLAB,这些方法都能游刃有余的用好。
2 s. ]% A) z) y" M3 k7 y
% F6 H# f* _" `( |7 x. u9 Q) l' } 03
2 e! P0 j! o$ }+ ?# r7 D n$ ]
8 N: Y( h. j) p线性规划、整数规划、多元规划、二次规划等规划类问题 3 [& N1 o7 B+ [/ F3 ^7 _
: l+ R6 g( |4 i$ n
数学建模竞赛中很多问题都和数学规划有关,可以说不少的模型都可以归结为一组不等式作为约束条件、几个函数表达式作为目标函数的问题。
7 V1 K2 H9 V6 z- @! ?
9 T0 b$ F% A) H; p# u& i遇到这类问题,求解就是关键了,比如 98 年 B 题,用很多不等式完全可以把问题刻画清楚,因此列举出规划后用 Lindo、Lingo 等软件来进行解决比较方便,所以还需要熟悉这两个软件。0 z& b/ j5 q8 z9 y4 e$ i
- y' m. @/ D+ l- _$ m8 y 04
$ o' J8 P: \- F
) H3 c9 u U8 x8 n1 P* w9 I+ b- g9 Z图论算法
& x6 r6 _( l, `' \$ ?9 ^) K
/ C9 {+ w# x$ t; H% h这类问题算法有很多,包括:Dijkstra、Floyd、Prim、Bellman-Ford,最大流,二分匹配等问题。
: V2 {# l/ Z- c9 N- {2 |6 i! {- D! W, Y2 A- @$ C/ C3 e/ a
关于此类图论算法,可参考 IntroductiontoAlgorithms--算法导论,关于图算法的第22章-第26章。
8 j! q6 a1 @7 T: t
3 S6 r; f4 r* e# t8 z3 x
) Q2 q i& Q; Z- F/ ~- _' N) [( t' ]/ M% r" g8 I6 X4 B5 J2 E
0 @* V0 r( D5 L* j8 C: C1 G1 u
05
" L9 [& m" S+ X" G9 n! q8 P. R d$ a3 T) f8 }
动态规划、回溯搜索、分治算法、分支定界等计算机算法 F. O( f' a' @ \8 g$ C: K5 H- x
6 w0 F, e: V2 o! C5 {" u
在数学建模竞赛中,如:92 年 B 题用分枝定界法,97年 B 题是典型的动态规划问题,此外 98 年 B 题体现了分治算法。+ u* t5 K! s- o+ J
/ y$ y f- w2 N! ~
6 ?) r2 J& P! r0 c1 U
这方面问题和 ACM 程序设计竞赛中的问题类似,推荐看一下算法导论,与《计算机算法设计与分析》(电子工业出版社)等与计算机算法有关的书。
5 U2 O4 y H" ?* z, }" K* A$ y9 q8 F6 y$ `( e" d2 J, b* J
06 - R$ w0 C0 t6 x' }
/ s- o, C9 o! z- T+ m最优化理论的三大经典算法:模拟退火法、神经网络、遗传算法
\4 {2 K0 L- ~6 v9 L4 Y1 h) T# s' L
这十几年来最优化理论有了飞速发展,模拟退火法、神经网络、遗传算法这三类算法发展很快。8 E6 E* \- D( c; r6 Y7 }. m* I' ?
8 p. r8 ^( h5 a6 X( d$ o4 s在数学建模竞赛中:比如 97 年 A 题的模拟退火算法,00 年 B 题的神经网络分类算法,01 年 B 题这种难题也可以使用神经网络。3 Y+ O2 s) y; K4 |0 z$ F" R c
- u& e4 H9 p" A# s
还有美国竞赛 89 年 A 题也和 BP 算法有关系,当时是 86 年刚提出 BP 算法,89 年就考了,说明赛题可能是当今前沿科技的抽象体现。6 }: T. U9 B) r, s
8 ~2 Z5 G! x! @
03 年 B 题伽马刀问题也是目前研究的课题,目前算法最佳的是遗传算法。8 Y. L# O4 ]: |$ ?
% z, o: {, z- g% o6 U
% ~- V" e9 n- h) a. i2 `+ i
4 D6 f: G% n* f v$ Y/ v; R
) Z2 x! d1 I# c2 Z! i0 ]
07
3 A' N- r- Z+ I' i/ U1 a$ r: Y/ T# R* ^* h6 E# T$ T+ o% Z, Q# j& u
网格算法和穷举法 + l4 }1 H) M5 f1 H M
* u% m' ]6 p9 |$ D
网格算法和穷举法一样,只是网格法是连续问题的穷举。1 K1 y" v+ s$ w6 |3 }* w# y& M
+ \ e `/ y! _6 |0 w; _比如要求在 N 个变量情况下的最优化问题,那么对这些变量可取的空间进行采点,比如在 [a;b] 区间内取 M+1 个点,那么这样循环就需要进行 (M+1)N 次运算,所以计算量很大。( `. v1 b% G! Y$ X6 k$ Y
1 s" a. |3 `8 n7 g# g( J# P在数学建模竞赛中:比如 97 年 A 题、99 年 B 题都可以用网格法搜索,这种方法最好在运算速度较快的计算机中进行,还有要用高级语言来做,最好不要用MATLAB 做网格,否则会算很久。. d, g1 S6 p% E( G
$ B( |! E+ c0 b% y' E5 g! o4 S) n% n3 Q ^6 ?
7 E3 Y& c8 x& Y) X: k
穷举法大家都熟悉,自不用多说了。
3 J2 I L" _; c7 g" c4 P1 A7 |: @9 u& x7 I e& H
08
* |- I( D! d% k! A6 f# i7 j0 e: h; a# N8 s) `8 E! j
一些连续离散化方法
1 ~0 S3 s( M5 m. a" H f8 V7 d% d4 n) \/ S; g
大部分物理问题的编程解决,都和这种方法有一定的联系。物理问题是反映我们生活在一个连续的世界中,计算机只能处理离散的量,所以需要对连续量进行离散处理。: u1 `' x1 Y3 c
" K L2 t& k* D4 j# k5 U6 o
这种方法应用很广,而且和上面的很多算法有关。事实上,网格算法、蒙特卡罗算法、模拟退火都用了这个思想。/ c& V8 l+ e3 v$ U8 {
9 k% y' f! K1 h" W) A, W7 [: ? 09 5 x/ O$ L( G5 F7 c0 W: x
) h: h- C8 N! N: b2 U: R
数值分析算法
) z: X- r2 b8 ], L
4 ^1 H& c! Q* k3 w9 z数值分析(numericalanalysis),是数学的一个分支,主要研究连续数学(区别于离散数学)问题的算法。
3 L! t' x* Y$ J3 {1 c+ Z) C, q
2 Q+ g0 Y% T v" n如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比如方程组求解、矩阵运算、函数积分等算法就需要额外编写库函数进行调用。3 f) U( g" k/ p4 V
) |, h \& n* J" \这类算法是针对高级语言而专门设的,如果你用的是 MATLAB、Mathematica,大可不必准备,因为像数值分析中有很多函数一般的数学软件是具备的。; @+ z/ h. v$ V" d
0 q2 S& l- d# a$ R 10
0 n6 {* C7 Y& g$ o8 f* C
( W) N1 w8 ]7 m. U. N图象处理算法
! i" n2 w1 D+ z5 {5 w2 A
" O- {/ O8 a% O& o) a在数学建模竞赛中:比如 01 年 A 题中需要你会读 BMP 图象、美国赛 98 年 A 题需要你知道三维插值计算,03 年 B 题要求更高,不但需要编程计算还要进行处理,而数模论文中也有很多图片需要展示,因此图象处理就是关键。做好这类问题,重要的是把 MATLAB 学好,特别是图象处理的部分。
- X0 ~4 F; T _/ }) g1 m1 L
8 m6 R/ o# A3 b3 }7 s, e- Z# s8 j' H
$ x2 C* g+ j; f3 ?
' ]0 O( d$ Y% ^5 I0 u/ k1 U% z" ] |