会计学蒙特卡罗(MonteCarlo)方法(fāngfǎ)随机(suíjī)模拟—MonteCarlo方法1.蒙特卡罗方法(fāngfǎ)概述蒙特卡罗方法的基本(jīběn)思想例1.蒲丰氏问题(wèntí)一些(yīxiē)人进行了实验，其结果列于下表：例2.射击问题（打靶(dǎbǎ)游戏）现假设该运动员进行了Ｎ次射击，每次射击的弹着点依次为r1，r2，…，rN，则Ｎ次得分g(r1)，g(r2)，…，g(rN)。其算术平均值为代表了该运动员的成绩。换言之，为积分(jīfēn)<g>的估计值，或近似值。在该例中，用Ｎ次试验所得成绩的算术平均值作为数学期望<g>的估计值（积分(jīfēn)近似值）。基本(jīběn)思想蒙特卡罗方法是用随机试验的方法计算积分，即将所要计算的积分看作服从某种分布密度函数f(r)的随机变量ｇ(r)的数学期望通过某种试验，得到Ｎ个观察值r1，r2，…，rN（用概率(gàilǜ)语言来说，从分布密度函数f(r)中抽取Ｎ个子样r1，r2，…，rN，），将相应的Ｎ个随机变量的值g(r1)，g(r2)，…，g(rN)的算术平均值作为积分的估计值（近似值）。为了得到具有一定精确度的近似解，所需试验的次数是很多的，通过人工(réngōng)方法作大量的试验相当困难，甚至是不可能的。电子计算机的出现，使得人们可以通过电子计算机来模拟随机试验过程，把巨大数目的随机试验交由计算机完成，使得蒙特卡罗方法得以广泛地应用，在现代化的科学技术中发挥应有的作用。计算机模拟试验(shìyàn)过程例１．蒲丰氏问题(wèntí)如何产生任意的（x,θ）？x在［0，a］上任意取值，表示x在［0，a］上是均匀分布的，其分布密度函数为：θ的分布密度函数为：因此，产生任意的（x,θ）的过程(guòchéng)就变成了由f1(x)抽样x及由f2(θ)抽样θ的过程(guòchéng)了。由此得到：其中ξ1，ξ2均为（0,1）上均匀分布的随机变量。每次投针试验，实际上变成在计算机上从两个均匀分布的随机变量中抽样得到（x,θ），然后定义描述针与平行线相交(xiāngjiāo)状况的随机变量s(x,θ)，为如果投针Ｎ次，则是针与平行线相交(xiāngjiāo)概率Ｐ的估计值。理论上，于是有例２．射击(shèjī)问题蒙特卡罗方法(fāngfǎ)的收敛性，误差收敛性误差(wùchā)当N充分大时，有如下的近似(jìnsì)式其中α称为显著性，1－α称为置信水平。这表明，不等式近似(jìnsì)地以概率1－α成立，且误差收敛速度的阶为。蒙特卡罗方法的误差ε定义为误差的两点说明：第一，蒙特卡罗方法的误差为概率误差，这与其他数值计算方法是有区别的；第二，误差中的均方差σ是未知的，必须(bìxū)使用其估计值来代替，在计算所求量的同时，可计算出；减小方差的各种(ɡèzhǒnɡ)技巧效率(xiàolǜ)蒙特卡罗方法(fāngfǎ)的特点能够比较逼真地描述具有随机性质的事物的特点及物理实验过程从这个意义上讲，蒙特卡罗方法可以部分代替物理实验，甚至可以得到物理实验难以得到的结果。用蒙特卡罗方法解决(jiějué)实际问题，可以直接从实际问题本身出发，而不从方程或数学表达式出发。它有直观、形象的特点。在计算s维空间中的任一区域Ds上的积分时，无论区域Ds的形状多么特殊，只要能给出描述Ds的几何特征的条件(tiáojiàn)，就可以从Ds中均匀产生N个点，得到积分的近似值。其中Ds为区域Ds的体积。优点3：收敛速度(sùdù)与问题的维数无关优点4:具有(jùyǒu)同时计算多个未知量的能力优点5:误差容易(róngyì)确定优点6:程序结构简单(jiǎndān)，易于实现缺点(quēdiǎn)1:收敛速度慢缺点(quēdiǎn)2:误差具有概率性蒙特卡罗方法(fāngfǎ)的主要应用范围