蒙特卡罗模拟方法课件解析.ppt

蒙特卡罗模拟方法 经济管理学院 任课教师 ：莫展宏 蒙特卡罗模拟方法 一、蒙特卡罗方法概述 二、蒙特卡罗方法模型 三、蒙特卡罗方法的优缺点及其适用范围 四、相关案例分析及软件操作 五、问题及相关答案 Monte Carlo方法的发展历史 早在17世纪，人们就知道用事件发生的“频率”来决定事件的“概率”。从方法特征的角度来说可以一直追溯到1777年法国科学家蒲丰（Buffon）提出的一种计算圆周率的方法——随机投针试验，即著名的蒲丰投针问题。 1777年，古稀之年的蒲丰在家中请来好些客人玩投针游戏（针长是线距之半），他事先没有给客人讲与π有关的事。客人们虽然不知道主人的用意，但是都参加了游戏。他们共投针2212次，其中704次相交。蒲丰说，2212/704=3.142，这就是π值。这着实让人们惊喜不已。 投针步骤 这一方法的步骤是： 1） 取一张白纸，在上面画上许多条间距为d的平行线。 2） 取一根长度为l（ld） 的针，随机地向画有平行直线的纸上掷n次，观察针与直线相交的次数，记为m 3）计算针与直线相交的概率． 例．蒲丰氏问题 设针投到地面上的位置可以用一组参数（x,θ）来描述，x为针中心的坐标，θ为针与平行线的夹角，如图所示。 任意投针，就是意味着x与θ都是任意取的，但x的范围限于［0，a］，夹角θ的范围限于［0，π］。在此情况下，针与平行线相交的数学条件是 一些人进行了实验，其结果列于下表 ： 20世纪四十年代，由于电子计算机的出现，利用电子计算机可以实现大量的随机抽样的试验，使得用随机试验方法解决实际问题才有了可能。 其中作为当时的代表性工作便是在第二次世界大战期间，为解决原子弹研制工作中，裂变物质的中子随机扩散问题，美国数学家冯.诺伊曼（Von Neumann）和乌拉姆（Ulam）等提出蒙特卡罗模拟方法。 由于当时工作是保密的，就给这种方法起了一个代号叫蒙特卡罗，即摩纳哥的一个赌城的名字。用赌城的名字作为随机模拟的名称，既反映了该方法的部分内涵，又易记忆，因而很快就得到人们的普遍接受。 蒙特卡罗方法的基本思想 蒙特卡罗方法又称计算机随机模拟方法。它是以概率统计理论为基础的一种方法。 由蒲丰试验可以看出，当所求问题的解是某个事件的概率，或者是某个随机变量的数学期望，或者是与概率、数学期望有关的量时，通过某种试验的方法，得出该事件发生的频率，或者该随机变量若干个具体观察值的算术平均值，通过它得到问题的解。这就是蒙特卡罗方法的基本思想。 因此，可以通俗地说，蒙特卡罗方法是用随机试验的方法计算积分，即将所要计算的积分看作服从某种分布密度函数f(r)的随机变量ｇ(r)的数学期望 通过某种试验，得到Ｎ个观察值r1，r2，…，rN（用概率语言来说，从分布密度函数f(r)中抽取Ｎ个子样r1，r2，…，rN，），将相应的Ｎ个随机变量的值g(r1)，g(r2)，…，g(rN)的算术平均值 作为积分的估计值（近似值）。 计算机模拟试验过程 计算机模拟试验过程，就是将试验过程（如投针问题）化为数学问题，在计算机上实现。 模拟程序 l=1; d=2; m=0; n=10000 for k=1:n; x=unifrnd(0,d/2); y=unifrnd(0,pi); if x0.5*1*sin(y) m=m+1 else end end p=m/n pi_m=1/p 例子 某投资项目每年所得盈 利额A由投资额P、劳动 生产率L、和原料及能 源价格Q三个因素。 模型建立的两点说明 Monte Carlo方法在求解一个问题是，总是需要根据问题的要求构造一个用于求解的概率统计模型，常见的模型把问题的解化为一个随机变量 的某个参数 的估计问题。 要估计的参数 通常设定为 的数学期望（亦平均值，即 ）。按统计学惯例， 可用 的样本 的平均值来估计，即 这时就必须采用主观概率,即由专家做出主观估计得到的概率。 另一方面,在对估测目标的资料与数据不足的情况下,不可能得知风险变量的真实分布时,根据当时或以前所收集到的类似信息和历史资料,通过专家分析或利用德尔菲法还是能够比较准确地估计上述各风险因素并用各种概率分布进行描述的。 Crystal ball软件对各种概率分布进行拟合以选取最合适的分布。 抽样次数与结果精度 解的均值与方差的计算公式： 随机数 随机数的定义 用Monte Car