蒙特卡罗模拟方法课件[77p][338mb]

蒙特卡罗模拟方法报告人：杨林吴颖科目：项目风险管理任课教师：尹志军蒙特卡罗模拟方法.一、蒙特卡罗方法概述.二、蒙特卡罗方法模型.三、蒙特卡罗方法的优缺点及其适用范围.四、相关案例分析及软件操作.五、问题及相关答案MonteCarlo方法的发展历史.早在17世纪，人们就知道用事件发生的“频率”来决定事件的“概率”。从方法特征的角度来说可以一直追溯到18世纪后半叶的蒲丰（Buffon）随机投针试验，即著名的蒲丰问题。1707-1788.1777年，古稀之年的蒲丰在家中请来好些客人玩投针游戏（针长是线距之半），他事先没有给客人讲与π有关的事。客人们虽然不知道主人的用意，但是都参加了游戏。他们共投针2212次，其中704次相交。蒲丰说，2212/704=3.142，这就是π值。这着实让人们惊喜不已。例．蒲丰氏问题.设针投到地面上的位置可以用一组参数（x,θ）来描述，x为针中心的坐标，θ为针与平行线的夹角，如图所示。.任意投针，就是意味着x与θ都是任意取的，但x的范围限于［0，a］，夹角θ的范围限于［0，π］。在此情况下，针与平行线相交的数学条件是针在平行线间的位置xlsin1,当xlsins(x,)0,其他1NsNs(xi,i)Ni1Ps(x,)f(x)f()dxd12dlsindx2l00aa2l2laPasN.一些人进行了实验，其结果列于下表：实验者年份投计次数π的实验值沃尔弗(Wolf)185050003.1596斯密思(Smith)185532043.1553福克斯(Fox)189411203.1419拉查里尼190134083.1415929(Lazzarini)20世纪四十年代，由于电子计算机的出现，利用电子计算机可以实现大量的随机抽样的试验，使得用随机试验方法解决实际问题才有了可能。其中作为当时的代表性工作便是在第二次世界大战期间，为解决原子弹研制工作中，裂变物质的中子随机扩散问题，美国数学家冯.诺伊曼（VonNeumann）和乌拉姆（Ulam）等提出蒙特卡罗模拟方法。由于当时工作是保密的，就给这种方法起了一个代号叫蒙特卡罗，即摩纳哥的一个赌城的名字。用赌城的名字作为随机模拟的名称，既反映了该方法的部分内涵，又易记忆，因而很快就得到人们的普遍接受。蒙特卡罗方法的基本思想.蒙特卡罗方法又称计算机随机模拟方法。它是以概率统计理论为基础的一种方法。.由蒲丰试验可以看出，当所求问题的解是某个事件的概率，或者是某个随机变量的数学期望，或者是与概率、数学期望有关的量时，通过某种试验的方法，得出该事件发生的频率，或者该随机变量若干个具体观察值的算术平均值，通过它得到问题的解。这就是蒙特卡罗方法的基本思想。因此，可以通俗地说，蒙特卡罗方法是用随机试验的方法计算积分，即将所要计算的积分看作服从某种分布密度函数f(r)的随机变量ｇ(r)的数学期望gg(r)f(r)dr0通过某种试验，得到Ｎ个观察值r1，r2，…，rN（用概率语言来说，从分布密度函数f(r)中抽取Ｎ个子样r1，r2，…，rN，），将相应的Ｎ个随机变量的值g(r1)，g(r2)，…，g(rN)的算术平均值1NgNg(ri)Ni1作为积分的估计值（近似值）。计算机模拟试验过程计算机模拟试验过程，就是将试验过程（如投针问题）化为数学问题，在计算机上实现。模拟程序.l=1;.d=2;.m=0;.n=10000.fork=1:n;.x=unifrnd(0,d/2);.y=unifrnd(0,pi);.ifx<0.5*1*sin(y).m=m+1.else.end.end.p=m/n.pi_m=1/p①建立概率统计模型②收集模型中风险变量的数据，确定风⑤根据随机数在各风险因数的分布函数险变量的概率分布中随机抽样，代入第一步中建立的数学模型③根据风险分析的精度要求，确定模拟次数N④建立对随机变量的抽样方法，产生随机数。⑥N个样本值⑦统计分析，估计均值， 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 差例子1N某投资项目每年所得盈AaPbL2cQ2d利额A由投资额P、劳动生产率L、和原料及能收集P,L,Q数据，确定分布函源价格Q三个因素。数f(P),f(L),f(Q)抽取P,L,Q一模拟次数N；根据分组随机布函数，产生随机数数，带根据历史数据，预测未来。入模型产生N个A值统计分析，估计均值，标准差X模型建立的两点说明.MonteCarlo方法在求解一个问题是，总是需要根据问题的要求构造一个用于求解的概率统计模型，常见的模型把问题的解化为一个随机变量X的某个参数的估计问题。.要估计的参数通常设定为的数学期望（亦平均值，即EX()）。按X(XXX...)统计学惯例，可用的样本1,2,n1nXX的平均值来估计，即knk1收集模型中风险变量的数据，确定风险因数的分布函数.这时就必须采用主观概率,即由专家做出主观估计得到的概率。.另一方面,在对估测目标的 资料 新概念英语资料下载李居明饿命改运学pdf成本会计期末资料社会工作导论资料工程结算所需资料清单 与数据不足的情况下,不可能得知风险变量的真实分布时,根据当时或以前所收集到的类似信息和历史资料,通过专家分析或利用德尔菲法还是能够比较准确地估计上述各风险因素并用各种概率分布进行描述的。.Crystalball软件对各种概率分布进行拟合以选取最合适的分布。抽样次数与结果精度1.解的均值与方差的计算公式：E(),()XVarX2nx2x是随机变量X的方差，而称Var()X为估计量方差。通常蒙特卡罗模拟中的样本量n很大，由统计学的中心极限定理知X渐进正态分布，即：xn1Xx1t2limp(x)e2dtxxn2X从而p()1xn式中α位小概率，1-α称为置信度：是标准正态分布中与α对应的临界值，可有统计分布表查得。X得统计学上称为与置信水平α对应的置信区间：由p()1xnXxx//nXn我们就把x/n记做是误差得到人们习惯的结果误差表示：X对于指定的误差ε,模拟所需抽样次数n可由导出：2xn随机数.随机数的定义用MonteCarlo方法模拟某过程时，需要产生各种概率分布的随机变量。最简单、最基本、最重要的随机变量是在［0，1］上均匀分布的随机变量。由该分布抽取的简单子样称为随机数序列，其中每一个体称为随机数。随机数属于一种特殊的由已知分布的随机抽样问题。随机数是随机抽样的基本工具。［0，1］上均匀分布（单位均匀分布），其分布密度函数为：1,0x1fx()0,其他分布函数为：0,x0F(x)x,0x1特征：独立性、均匀性1,x1随机数的产生方法.随机数表.物理方法.计算机方法随机数表.随机数表是由0,1,2,…,9十个数字组成，每个数字以0.1的概率出现，数字之间相互独立。.方法：如果要得到n位有效数字的随机数，只需将表中每n个相邻的随机数字合并在一起，且在最高位的前边加上小数点即可。例如：某随机数表第一行数字为7634258910…，要想得到三位有效数字的随机数依次为：0.763，0.425，0.891物理方法.基本原理：利用某些物理现象，在计算机上增加些特殊设备，可以在计算机上直接产生随机数。.缺点：无法重复实现费用昂贵计算机方法.在计算机上产生随机数最实用、最常见的方法是数学方法，即用如下递推公式：nn1T()产生随机数序列，对于给定的初始值n，确定n1，n=1,2…存在的问题：1，不满足相互独立的要求2，不可避免的出现重复问题所以成为伪随机数问题的解决：1.选取好的递推公式2.不是本质问题产生伪随机数的乘同余方法.乘同余方法是由Lehmer在1951年提出来的，它的一般形式是：对于任一初始值x1，伪随机数序列由下面递推公式确定：xii1a.x(modM)xi12,i1,2,Ma为乘子，X为种子（初值）；M成为模数。上式表示x是ax.iii1被M整除后的余数，叫做x1与对模M的同余。•利用乘同余法产生伪随机数的步骤如下：（1）取种子x1、乘子、和模数M；（2）由式（1）获得一系列，x2...；（3）由式（2）得到一系列，…。这就是所要产生的伪随机数2的序列1乘同余方法在计算机上的使用.为了便于在计算机上使用，通常取：Ｍ=2s.其中s为计算机中二进制数的最大可能有效位数.x1=奇数.a=52k+1.其中k为使52k+1在计算机上所能容纳的最大整数，即a为计算机上所能容纳的5的最大奇次幂。一般地，s=32时，a=513；s=48，a=515等。伪随机数序列的最大容量λ(M)=2s-2。.乘同余方法是使用的最多、最广的方法，在计算机上被广泛地使用。用MATLAB产生随机数.语言：连续均匀分布的函数表达式为R=unifrnd(A,B).演示：forn=1:100;k=unifrnd(0,1)end随机抽样及其特点.由巳知分布的随机抽样指的是由己知分布的总体中抽取简单子样。随机数序列是由单位均匀分布的总体中抽取的简单子样，属于一种特殊的由已知分布的随机抽样问题。下表所叙述的由任意已知分布中抽取简单子样，是在假设随机数为已知量的前提下，使用严格的数学方法产生的。.直接抽样方法.对于任意给定的分布函数F(x)，直接抽样方法如下：Xninft,n1,2,,NF(t)n.其中，ξ1，ξ2，…，ξN为随机数序列。为方便起见，将上式简化为：XFinftF(t).若不加特殊说明，今后将总用这种类似的简化形式表示，ξ总表示随机数。离散型分布的直接抽样方法.对于任意离散型分布：F(x)Pixix.其中x1，x2，…为离散型分布函数的跳跃点，P1，P2，…为相应的概率，根据前述直接抽样法，有离散型分布的直接抽样方法如下：I－1I当XFxI,PiPii＝1i＝1.该结果表明，为了实现由任意离散型分布的随机抽样，直接抽样方法是非常理想的。例1.二项分布的抽样.二项分布为离散型分布，其概率函数为：P(xn)PCnPn(1P)NnnN.其中，P为概率。对该分布的直接抽样方法如下：n－1n当XFn,PiPii＝0i＝0例2.掷骰子点数的抽样.掷骰子点数X=n的概率为：1P(Xn)6.选取随机数ξ，如n1n66.则XFnXF[6]1.在等概率的情况下，可使用如下更简单的方法：.其中［］表示取整数。连续型分布的直接抽样方法.对于连续型分布，如果分布函数F(x)的反函数F－1(x)存在，则直接抽样方法是：1XFF()例3.在［a，b］上均匀分布的抽样.在［a，b］上均匀分布的分布函数为：0当xaxaF(x)当axbba1当xb.则XFa(ba)由任意已知分布中抽取简单子样的方法还包括，挑选抽样方法，复合抽样方法，复合挑选抽样方法，替换抽样方法。圆内均匀分布抽样要用到挑选抽样方法，指数分布函数抽样要用到复合抽样方法，正态分布的抽样和β分布的抽样要用到替换抽样方法等。每种方法各有其优缺点和使用范围。常用概率分布的抽样公式分布名称抽样公式注[a，b]均匀分布abar指数分布lnr12正态分布ri6i1caabacar,0rbaa，b，c为三角分布三角分布cabbabc1r,r1的参数baf[()]()abar12fmrbr11as分布a()barm1rs2r，s为函数参数三角分布三角形概率分布是一种应用较广连续型概率分布,它是一种3点估计:特别适用于对那些风险变量缺乏历史统计资料和数据,但可以经过咨询专家意见,得出各参数变量的最乐观值(a),最可能出现的中间值(b)以及最悲观值(m),这3个估计值(a,b,m)构成一个三角形分布。实际上，Matlab软件为我们提供了一种简单快捷的产生各种常用分布随机数的方法。其功能和特点：（1）界面友好，编程效率高。（2）功能强大，可扩展性强。（3）强大的数值计算功能和符号计算功能。（4）图形功能灵活方便。Matlab常用的随机数产生函数函数名调用形式函数注释betarndR=betarnd(A,B)分布随机数产生函数binorndR=binornd(N,P,MM,NN)二项分布随机数产生函数chi2rndR=chi2rnd(v)卡方分布随机数产生函数frndR=frnd(v1,v2)F分布随机数产生函数georndR=geornd(p)几何分布随机数产生函数hygerndR=hygernd(M,K,N)超几何分布随机数产生函数mvnrndR=mvnrnd(mu,sigma,cases)多元正态分布随机数产生函数normrndR=normrnd(mu,sigma)正态分布随机数产生函数trndR=trnd(v)t分布随机数产生函数有了这些随机产生函数，就可以直接产生满足分布F(x)的随机数了，而无需通过先求出连续均匀分布的随机数，再通过抽样公式得出所求分布函数的随机抽样。演示：forn=1:100;k=betarnd(0.1,100)end蒙特卡罗方法的特点优点缺点①能够比较逼真地描述具①收敛速度慢。有随机性质的事物的特②误差具有概率性。点及物理实验过程。③进行模拟的前提是②受几何条件限制小。各输入变量是相互独③收敛速度与问题的维数立的。无关。④误差容易确定。⑤程序结构简单，易于实现。①能够比较逼真地描述具有随机性质的事物的特点及物理实验过程.从这个意义上讲，蒙特卡罗方法可以部分代替物理实验，甚至可以得到物理实验难以得到的结果。用蒙特卡罗方法解决实际问题，可以直接从实际问题本身出发，而不从方程或数学表达式出发。它有直观、形象的特点。②受几何条件限制小.在计算s维空间中4的任一区域Ds上的2积分，无论区域Ds的形状多么特殊，-4-224只要能给出描述Ds的几何特征的条件，-2-4就可以从Ds中均匀产生N个点③收敛速度与问题的维数无关.由误差定义可知，在给定置信水平情况下，蒙特卡罗方法的收敛速度为O(N1/2)，与问题本身的维数无关。维数的变化，只引起抽样时间及估计量计算时间的变化，不影响误差。也就是说，使用蒙特卡罗方法时，抽取的子样总数N与维数s无关。维数的增加，除了增加相应的计算量外，不影响问题的误差。这一特点，决定了蒙特卡罗方法对多维问题的适应性。程序结构简单，易于实现.在计算机上进行蒙特卡罗方法计算时，程序结构简单，分块性强，易于实现。①收敛速度慢.如前所述，蒙特卡罗方法的收敛为O(N1/2)，一般不容易得到精确度较高的近似结果。对于维数少（三维以下）的问题，不如其他方法好。②误差具有概率性.由于蒙特卡罗方法的误差是在一定置信水平下估计的，所以它的误差具有概率性，而不是一般意义下的误差。蒙特卡罗方法的主要应用范围.蒙特卡罗方法所特有的优点，使得它的应用范围越来越广。它的主要应用范围包括：粒子输运问题，统计物理，典型数学问题，真空技术，激光技术以及医学，生物，探矿等方面，特别适用于在计算机上对大型项目、新产品项目和其他含有大量不确定因素的复杂决策系统进行风险模拟分析。随着科学技术的发展，其应用范围将更加广泛。第五节项目风险案例分析现以成都某房地产开发公司对一综合开发用地进行投资开发为例，用基于蒙特卡罗模拟方法为原理的EXCEL插件——CrystalBall工具对该开发项目进行风险决策分析。一、项目概况和基本数据的确定该项目位于成都市锦江区，占地面积47亩；该房地产公司根据市场状况调查，结合该地块的规划说明，在做了充分的 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 之后，确定了两套主要的投资方案。甲方案：该地块主要以小高层电梯住宅开发为主，辅以车库和部分商业配套设施，开发期共三年。甲方案预测出的的主要经济技术指标见表5-1。表5-1甲方案的主要经济技术指标建设经营期序号项目合计200520062007一现金流入45306018064272421销售收入4530601806427242二现金流出413531627712329127471开发建设投资2658316277850218042营业税金及附加25140100315123土地增值税22920022924所得税9964028257139三净现金流量（税后）3953-16277573514495累计净现金流量（税后）-16277-105423953四现值系数（i=10%)10.9090.826五净现值（税后）915-16277521411979累计净现值（税后）-16277-11064915乙方案：将该地块开发为商业类地产为主，外设露天停车场，配以部分小户型电梯公寓，开发期仍为三年。乙方案预测出的的主要经济技术指标见表5-2。表5-2乙方案的主要经济技术指标建设经营期序号项目合计200520062007一现金流入54660032082218401销售收入5466003208221840二现金流出492151762819391121961开发建设投资30626176281095520432营业税金及附加30340182212123土地增值税41900041904所得税11365066144750三净现金流量（税后）5445-17628134299644累计净现金流量（税后）-17628-41995445四现值系数（i=10%)10.9090.826五净现值（税后）2550-17628122087970累计净现值（税后）-17628-54202550根据该表5-1第五项，我们可以得出甲方案的财务净现值NPV=915万元，同样根据该表5-2第五项，我们可以得出乙方案的财务净现值NPV=2550万元。通过对两种方案动态财务指标的比较，我们可以很明确的断定采用乙方案将是开发商最佳的选择。但不容忽略的一点是，以商业类开发为主的乙方案，在销售期间，销售面积和销售价格具有较大的不确定性；而以住宅类开发为主的甲方案在对未来的销售面积和销售价格方面将有更大的把握度。仅从这点上我们就可以判断乙方案的风险大于甲方案。为了做出精准的判断，需要在此基础之上进行更精准的风险分析。二、采用蒙特卡罗方法进行风险决策分析（一）、识别项目风险在投资开发项目时，实际情况千差万别，重要的风险变量也各不相同，这就需要分析人员根据项目的具体情况，运用适当的风险辨识的方法从影响投资的众多因素中找出关键的风险变量。本案例采用“德尔菲法”确定影响该项目的7个主要风险变量：住宅销售收入（P1*S1）、商业销售收入（P2*S2）、土地费用（K1）、前期费用（K2）、开发建设费用（K3）、营销费用（K4）、其他费用（K5）。（二）、确定每个风险变量的概率分布同样采用“德尔菲法”估计出以上7个风险变量概率分布和其分布函数中的具体参数,如下表所示：表5-3甲方案风险变量概率分布第一年分布参数住宅类销售收入三角分布无销售收入商业类销售收入三角分布无销售收入土地费用均匀分布a：11182b:12105前期费用正态分布均值：911方差：50开发建设费用三角分布a：3112b：3374m：3276营销费用三角分布a：235b：329m：313其他费用正态分布均值：249方差：15第二年分布参数住宅类销售收入三角分布a：13710b：18762m：14432商业类销售收入三角分布a：759b：1036m：1012土地费用均匀分布无支出前期费用正态分布均值：727方差：30开发建设费用三角分布a：6027b：6813m：6551营销费用三角分布a：251b：326m：313其他费用正态分布均值：911方差：55第三年住宅类销售收入三角分布a：21569b：28515m：22704商业类销售收入三角分布a：1304b：1739m：1656土地费用均匀分布无支出前期费用正态分布无支出开发建设费用三角分布a：1085b：1136m：1092营销费用三角分布a：334b：443m：418其他费用正态分布均值：294方差：20表5-4乙方案风险变量概率分布第一年分布参数住宅类销售收入三角分布无销售收入商业类销售收入三角分布无销售收入土地费用均匀分布a：11182b:12105前期费用正态分布均值：1249方差：80开发建设费用三角分布a：4007b：4555m：4218营销费用三角分布a：258b：413m：368其他费用正态分布均值：265方差：30第二年分布参数住宅类销售收入三角分布a：3996b：5328m：4440商业类销售收入三角分布a：14190b：28948m：28380土地费用均匀分布无支出前期费用正态分布均值：1003方差：90开发建设费用三角分布a：7760b：9110m：8435营销费用三角分布a：472b：565m：491其他费用正态分布均值：1025方差：100第三年住宅类销售收入三角分布a：1080b：1440m：1200商业类销售收入三角分布a：10526b：21053m：20640土地费用均匀分布无支出前期费用正态分布无支出开发建设费用三角分布a：1397b：1518m：1405营销费用三角分布a：350b：442m：368其他费用正态分布均值：269方差：30三、定义模型并确定模拟次数定义财务净现值NPV的模型为：nn25NPV(P*S)(/,,)llPFittcklPFitt(/,,)ct1l1t1l11(/,,)PFit其中，ct，i为基准折现率，n为1i项目的生命周期。c为了确保模拟结果与实际分布最大限度的接近一致，我们取95%的置信度，拟进行10000次的模拟实验。进行10000次的模拟，得出甲、乙方案的NPV的统计数据。表5-5甲方案的 评价 LEC评价法下载LEC评价法下载评价量规免费下载学院评价表文档下载学院评价表文档下载 指标统计值统计值NPV模拟次数10000均值672.24中值604.66标准差1052.27方差1107271.23偏差0.3347峰度2.72Coeff.ofVariability1.57最小值-1833.45最大值4448.76标准误差1052表5-6乙方案的评价指标统计值统计值NPV模拟次数10000均值432.59中值617.6标准差2157.44方差4654568.25偏差-0.3882峰度2.66Coeff.ofVariability4.99最小值-7334.47最大值5529.92标准误差21.57(四）、分析决策1、通过表5-5甲方案的财务净现值统计值和表5-6乙方案的财务净现值统计值，我们可以出，两个方案的NPV期望值均大于零，但甲方案的值大于乙方案。2、进一步对各方案的风险度进行比较，甲方案NPV的标准差为1052.27，而乙的标准差为2157.44，说明乙方案的偏离程度较大；并且甲方案NPV介于[min:-1833.45,max:4448.76]之间,乙方案NPV在[min:-7334.47,max:5529.92]之间，再次说明乙方案NPV的风险度大于甲方案。3、利用EXCEAL可以很容易评价指标具体的概率分布，如表5-7，表5-7甲乙方案风险概率分布甲方案的概率分布乙方案的概率分布概率分布NPV概率分布NPV0℅-1955.550℅-7322.82929710℅-635.3310℅-2546.58849120℅-260.6820℅-1446.0021328.24℅030℅-649.928374430℅52.1339.33℅040℅342.1640℅37.7284323150℅62350℅648.925504960℅913.2760℅1242.51551870℅1214.7670℅1810.41007580℅1585.5480℅2404.75315290℅2098.3990℅3149.852139100℅4534.23100℅5477.691348因此，应该采用甲方案。4、总结通过上面的分析，利用蒙特卡罗方法模拟分析得出的结果与使用传统的分析技术得出的结果相比，不仅能够分析风险因素对整个项目预期收益的影响程度，而且还能科学地估计出风险发生的概率大小，并且这样的估计是建立在充分考虑了多个风险变量共同影响、共同作用的基础之上，能够为风险决策者提供有实用价值的决策依据。因此有助于我们对多套投资方案进行筛选比较。CrystalBall软件简介CrystalBall软件是由美国Decisioneering公司开发的,为Excel电子表格提供的功能强大的加载宏。它充分利用微软视窗环境,提供了含有易学易用的图形包的高级模拟技术的独特组合。该软件包主要有计算机仿真模拟功能、时间序列数据生成预测和OptQuest功能,使其可以在运行结果中自动搜索仿真模型的最优解。CrystalBall软件的使用步骤.①定义随机的输入单元格:加载CrystalBall到Excel中,并且建立一个工作表,将投资预测的相关变量输入电子表格中;.②定义随机单元格的概率分布:利用软件的DefineAssumption功能为相应变量设定概率分布,利用DefineDecision定义决策变量;.③定义预测的输出单元格:利用DefineForecast功能定义输出变量的单元格;.④设定运行参数:在RunPreference功能中定义模拟次数、敏感度分析等参数;.⑤运行仿真:点击Run进行模拟运算,分析模拟结果。问题1、蒙特卡罗方法的基本思想是什么？2、用蒙特卡罗模型解决实际问题的基本步骤是什么？3、蒙特卡罗方法的优缺点各有哪些？4、由蒙特卡罗方法的误差公式我们可以推断出其有那些优缺点？5蒙特卡罗模拟与随机抽样统计分析有什么区别？Theanswer.1、当所求问题的解是某个事件的概率，或者是某个随机变量的数学期望，或者是与概率、数学期望有关的量时，通过某种试验的方法，得出该事件发生的频率，或者该随机变量若干个具体观察值的算术平均值，通过它得到问题的解。这就是蒙特卡罗方法的基本思想。.2、（1）建立数学模型（2）收集模型中风险变量的数据，确定风险因数的分布函数。（3）确定模拟次数、产生随机数。（4）由产生的随机数在各风险变量的分布函数中随机抽样，带入模型求出目标变量的一个样本值。（5）重复第4步N次，产生N个样本值，对得到的N个样本值进行统计分析。.3、优点①能够比较逼真地描述具有随机性质的事物的特点及物理实验过程。②受几何条件限制小。③收敛速度与问题的维数无关。④误差容易确定。⑤程序结构简单，易于实现。.缺点①收敛速度慢。②误差具有概率性。③进行模拟的前提是各输入变量是相互独立的。、通常，蒙特卡罗方法的误差定义为.4εN.上式中与置信度α是一一对应的，根据问题的要求确定出置信水平后，查标准正态分布表，就可以确定出。蒙特卡罗方法的误差为概率误差由此性质我们可以得知蒙特卡罗的优点：①收敛速度与问题的维数无关。②误差容易确定。以及缺点：误差具有概率性。.5、蒙特卡罗模拟是一种simulation，通过建立模型，产生相应分布的随机数(实际是伪随机变量)，来模拟实际存在的过程，并且分析相关的结果。.比如，首先观察客观世界，某一个过程有几个步骤，每个步骤所用的时间符合哪种分布(可以是估计)，数学特征是什么(可能通过测量)，然后建立起相关的模型。.通过一定的算法，产生符合标准正态分布的伪随机变量，然后对变量进行运算，来产生其他分布的变量，比如泊松分布，指数分布等。.对于模型中的每个步骤，通过产生的随机变量进行模拟，来得到整个过程所用的时间的数学特征，比如平均值，方差等。根据这些数学特征，就可以对以后的发展做出预测。.而随机抽样统计分析是对实际数据的抽样分析。通过抽样的数学特征，来估计样本的数学特征。也就是通过局部认识总体。.一句话概括，蒙特卡罗模拟对计算机产生的，而不是对实际存在的数据统计，而随机抽样分析是对实际存在的数据分析。思考题.对蒙特卡洛方法中风险变量相关性的探讨？主要文献1陈立文.项目投资风险分析理论与方法.北京：机械工业出版社，20042尹志军.工程项目投资风险度量与模拟及收益优化模型研究：[学位论文]，天津：河北工业大学，20023王卓甫.工程项目风险管理.北京：中国水利水电出版社，20024石海均.房地产投资分析.大连：大连理工大学出版社，19945闫雪晶.蒙特卡罗模拟方法房地产投资风险分析的应用：[学位论文]，成都：西南财经大学，20066杨衡.蒙特卡罗模拟优化与风险决策分析的应用研究：[学位论文]，天津：天津大学，20047刘军.蒙特卡罗方法及应用：课件，哈佛大学，20058林谦.蒙特卡罗方法及应用：课件，清华大学，2007请老师和各位同学批评指正！谢谢！