第十一讲 Monte-Carlo模拟方法在风险分析中的应用.ppt

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,本讲目标：,了解,Monte-Carlo,模拟的基本思路和步骤、计算机的工作模拟流程；,熟悉,crystal,ball,软件的使用。,第十一讲,Monte-Carlo,模拟方法在风险分析中的应用,第一节 蒙特卡罗法概述,蒙特卡罗方法也称为随机模拟、统计表试验方法，是一种依据统计理论，利用计算机来研究风险发生概率或风险损失的数值计算方法。在目前的工程项目风险分析中，是一种应用广泛、相对较精确的方法。,蒙特卡罗方法源于第二次世界大战期间，为解决原子弹研制工作中，裂变物质的中子随机扩散问题，美国数学家冯,.,诺伊曼（,Von Neumann,）和乌拉姆（,Ulam,）等提出蒙特卡罗模拟方法。由于当时工作是保密的，就给这种方法起了一个代号叫蒙特卡罗，即摩纳哥的一个赌城的名字。用赌城的名字作为随机模拟的名称，既反映了该方法的部分内涵，又易记忆，因而很快就得到人们的普遍接受。,该方法的应用十分广泛，如：项目融资、投标报价、投资和房地产等。,蒙特卡罗方法的基本思想是：,(1),将目标变量用一数学模型表示，该数学模型可被称为模拟模型，模型中尽可能地包含影响该目标变量的主要风险变量。,(2),每个风险变量的风险结果及其相应的概率值均可用一具体的概率分布来描述，然后利用抽样技术来产生随机数，再根据这一随机数在各风险变量的概率分布中随机取一值。,(3),当各风险变量的取值确定后，目标变量就可根据所建立的模拟模型计算得出。这样重复,N,次，便可得到,N,组目标变量值，通过产生随机数得出目标变量具体值的过程便是蒙特卡罗模拟过程。,例如，某投资项目每年所得盈利额,A,由投资额,P,、劳动生产率,L,和原料及能源价格,Q,三个因素确定，表达式为：,式中,a,，,b,，,c,，,d,均为常数；,P,、,L,、,Q,均为随机变量，其分布函数分别为,f(P),、,f(L),、,f(Q),。根据其分布函数每次随机抽取一组,P,、,L,、,Q,，通过上式便可得到,A,，重复取,N,次（模拟次数），便可得到,N,个,A,值，然后对,N,个样本值进行统计分析，得到分布曲线，并检验其概率分布，估计其均值和标准差。,蒙特卡罗模拟方法的优点：,（,1,）模拟算法简单，过程灵活；,（,2,）可模拟分析多元风险因素变化对结果的影响；,（,3,）模拟成本低，并可方便地补充更新数据。,蒙特卡罗模拟方法的局限性：,（,1,）蒙特卡罗方法要求的数据信息较多。,（,2,）进行模拟的前提是各输入变量是相互独立的。,（,3,）对一些复杂问题，要想达到较高的模拟精度需要进行较多的模拟次数。,蒙特卡罗模拟方法的步骤,(1),明确要解决的问题的实质，确定分析计算的目标，把握与分析计算相关的环境、条件等基本问题；将被分析的项目目标变量用一系列风险变量以一具体的数学模型表示出来。,(2),对第一步建立的数学模型中的风险变量进行风险识别和分析，收集风险变量的相关数据，对其加工分析。在风险分析基础上对各风险变量进行风险结果及相应的概率分析，确定风险因素的分布函数。,(3),根据风险分析的精度要求，确定模拟次数、产生随机数，并根据随机数在各风险变量的概率分布中随机取一值。将各风险变量的取值代入第一步中建立的数学模型，求得项目目标变量一具体值，即得到一个随机事件的样本值。,(4),重复第三步的实验操作，取得,N,个目标变量值。一般说来，重复的试验次数越多，目标值的分布越接近实际情况。,(5),对得到的,N,个样本值进行统计分析，得到分布曲线，并检验其概率分布，估计其均值和标准差。,第二节 模拟次数的确定,在蒙特卡罗方法中，模拟次数,N,是首先要确定的参数。模拟次数,N,和蒙特卡罗方法的计算精度直接相关，下面给出了在给定置信度和误差条件下，模拟次数,N,的计算公式。,抽样次数与结果精度,如果均值与方差的计算公式：,是随机变量,X,的方差，而称 为估计量方差。通常蒙特卡罗模拟中的样本量,n,很大，由统计学的中心极限定理知 渐进正态分布，即：,从而,式中,为小概率，,1-,称为置信度：是标准正态分布中与,对应的临界值，可由统计分布表查得。,则：,由,得统计学上称为与置信水平,1-,对应的置信区间：,我们就把 记做是误差,得到人们习惯的结果误差表示：,对于指定的误差,模拟所需抽样次数,n,可由 导出：,3,、随机变量的抽样,随机变量的抽样指的就是由已知分布的总体中产生简单子样。本文随机变量的抽样采用了直接抽样法，即把已取得的伪随机数,r,，通过给定的抽样模型，使伪随机数,r,变为给定分布的随机变量的样本。下表列出了常用概率分布的抽样公式。,常用概率分布的抽样公式,分布名称,抽样公式,注,a,，,b,均匀分布,指数分布,正态分布,三角分布,a,，,b,，,c,为三角分布的参数,a,c,为范围，均值为,b,。,有了这些随机产生函数，就可以直接产生满足分布,F(x),的随机数了，而无需通过先求出连续均匀分布的随机数，再通过抽样公式得出所求分布的随机数。下面来通过一个实例来加深对蒙特卡罗模拟方法的理解。,4,、项目风险案例分析,现以成都某房地产开发公司对一综合开发用地进行投资开发为例，用基于蒙特卡罗模拟方法为原理的,EXCEL,插件,Crystal Ball,工具对该开发项目进行风险决策分析。,(1),项目概况和基本数据的确定,该项目位于成都市锦江区，占地面积,47,亩；该房地产公司根据市场状况调查，结合该地块的规划说明，在做了充分的方案设计之后，确定了两套主要的投资方案。,甲方案：该地块主要以小高层电梯住宅开发为主，辅以车库和部分商业配套设施，开发期共三年。甲方案预测出的的主要经济技术指标见表,1,。,表,1,甲方案的主要经济技术指标,序号,项目,合计,建设经营期,2005,2006,2007,一,现金流入,45306,0,18064,27242,1,销售收入,45306,0,18064,27242,二,现金流出,41353,16277,12329,12747,1,开发建设投资,26583,16277,8502,1804,2,营业税金及附加,2514,0,1003,1512,3,土地增值税,2292,0,0,2292,4,所得税,9964,0,2825,7139,三,净现金流量（税后）,3953,-16277,5735,14495,累计净现金流量（税后）,-16277,-10542,3953,四,现值系数（,i=10%),1,0.909,0.826,五,净现金流量折现值（税后）,915,-16277,5214,11979,累计净现金流量折现值（税后）,-16277,-11064,915,乙方案：将该地块开发为商业类地产为主，外设露天停车场，配以部分小户型电梯公寓，开发期仍为三年。乙方案预测出的的主要经济技术指标见表,2,。,表,2,乙方案的主要经济技术指标,序号,项目,合计,建设经营期,2005,2006,2007,一,现金流入,54660,0,32082,21840,1,销售收入,54660,0,32082,21840,二,现金流出,49215,17628,19391,12196,1,开发建设投资,30626,17628,10955,2043,2,营业税金及附加,3034,0,1822,1212,3,土地增值税,4190,0,0,4190,4,所得税,11365,0,6614,4750,三,净现金流量（税后）,5445,-17628,13429,9644,累计净现金流量（税后）,-17628,-4199,5445,四,现值系数（,i=10%),1,0.909,0.826,五,净现金流量折现值（税后）,2550,-17628,12208,7970,累计净现净现金流量折现值（税后）,-17628,-5420,2550,根据该表,1,第五项，我们可以得出甲方案的财务净现值,NPV=915,万元，同样根据该表,2,第五项，我们可以得出乙方案的财务净现值,NPV=2550,万元。通过对两种方案动态财务指标的比较，我们可以很明确的断定采用乙方案将是开发商最佳的选择。但不容忽略的一点是，以商业类开发为主的乙方案，在销售期间，销售面积和销售价格具有较大的不确定性；而以住宅类开发为主的甲方案在对未来的销售面积和销售价格方面将有更大的把握度。仅从这点上我们就可以判断乙方案的风险大于甲方案。为了做出精准的判断，需要在此基础之上进行更精准的风险分析。,（,2),采用蒙特卡罗方法进行风险决策分析,识别项目风险,在投资开发项目时，实际情况千差万别，重要的风险变量也各不相同，这就需要分析人员根据项目的具体情况，运用适当的风险辨识的方法从影响投资的众多因素中找出关键的风险变量。本案例采用“德尔菲法”确定影响该项目的,7,个主要风险变量：住宅销售收入（,P1,S1,）、商业销售收入（,P2,S2,）、土地费用（,K1,）、前期费用（,K2,）、开发建设费用（,K3,）、营销费用（,K4,）、其他费用（,K5,）。,确定每个风险变量的概率分布,同样采用“德尔菲法”估计出以上,7,个风险变量概率分布和其分布函数中的具体参数，如下表所示：,表,3,甲方案风险变量概率分布,第一年,分布,参数,住宅类销售收入,三角分布,无销售收入,商业类销售收入,三角分布,无销售收入,土地费用,均匀分布,a,：,11182 b:12105,前期费用,正态分布,均值：,911,标准差：,50,开发建设费用,三角分布,a,：,3112 b,：,3374 m,：,3276,营销费用,三角分布,a,：,1235 b,：,1329 m,：,1313,其他费用,正态分布,均值：,249,标准差：,15,第二年,分布,参数,住宅类销售收入,三角分布,a,：,13710 b,：,18762 m,：,14432,商业类销售收入,三角分布,a,：,759 b,：,1036 m,：,1012,土地费用,均匀分布,无支出,前期费用,正态分布,均值：,727,标准差：,30,开发建设费用,三角分布,a,：,6027 b,：,6813 m,：,6551,营销费用,三角分布,a,：,2510 b,：,3260 m,：,3130,其他费用,正态分布,均值：,2911,标准差：,55,第三年,住宅类销售收入,三角分布,a,：,21569 b,：,28515 m,：,22704,商业类销售收入,三角分布,a,：,1304 b,：,1739 m,：,1656,土地费用,均匀分布,无支出,前期费用,正态分布,无支出,开发建设费用,三角分布,a,：,1085 b,：,1136 m,：,1092,营销费用,三角分布,a,：,3340 b,：,4430 m,：,4180,其他费用,正态分布,均值：,1294,标准差：,20,表,4,乙方案风险变量概率分布,第一年,分布,参数,住宅类销售收入,三角分布,无销售收入,商业类销售收入,三角分布,无销售收入,土地费用,均匀分布,a,：,11182 b:12105,前期费用,正态分布,均值：,1249,方差：,80,开发建设费用,三角分布,a,：,4007 b,：,4555 m,：,4218,营销费用,三角分布,a,：,2580 b,：,4130 m,：,3680,其他费用,正态分布,均值：,265,方差：,30,第二年,分布,参数,住宅类销售收入,三角分布,a,：,3996 b,：,5328 m,：,4440,商业类销售收入,三角分布,a,：,14190 b,：,28948 m,：,28380,土地费用,均匀分布,无支出,前期费用,正态分布,均值：,1003,方差：,90,开发建设费用,三角分布,a,：,7760 b,：,9110 m,：,8435,营销费用,三角分布,a,：,4720 b,：,5650 m,：,4910,其他费用,正态分布,均值：,825,方差：,100,第三年,住宅类销售收入,三角分布,a,：,1080 b,：,1440 m,：,1200,商业类销售收入,三角分布,a,：,10526 b,：,21053 m,：,20640,土地费用,均匀分布,无支出,前期费用,正态分布,无支出,开发建设费用,三角分布,a,：,1397 b,：,1518 m,：,1405,营销费用,三角分布,a,：,3500 b,：,4420 m,：,3680,其他费用,正态分布,均值：,869,方差：,30,定义模型并确定模拟次数,定义财务净现值,NPV,的模型为：,其中，,i,为基准折现率，,n,为项目的生命周期。,为了确保模拟结果与实际分布最大限度的接近一致，我们取,95%,的置信度，拟进行,10000,次的模拟实验。进行,10000,次的模拟，得出甲、乙方案的,NPV,的统计数据。,表,5,甲方案的评价指标统计值,统计值,NPV,模拟次数,10000,均值,1870.87,中值,1764.68,标准差,1669.47,方差,2787120.65,偏差,0.3640,峰度,2.78,Coeff.of Variability,0.8923,最小值,-2107.73,最大值,7909.98,标准误差,16.69,表,6,乙方案的评价指标统计值,统计值,NPV,模拟次数,10000,均值,1269.00,中值,1587.07,标准差,3744.51,方差,14021330.40,偏差,-0.3999,峰度,2.68,Coeff.of Variability,2.95,最小值,-12762.89,最大值,9593.81,标准误差,37.45,分析决策,A,、通过表,5,甲方案的财务净现值统计值和表,6,乙方案的财务净现值统计值，我们可以看出，两个方案的,NPV,期望值均大于零，但甲方案的值大于乙方案。,B,、进一步对各方案的风险度进行比较，甲方案,NPV,的标准差为,1669.47,，而乙的标准差为,3744.51,，说明乙方案的偏离程度较大；并且甲方案,NPV,介于,min-2107.73,，,max:7909.98,之间,乙方案,NPV,在,min-12762.89,，,max:9593.81,之间，再次说明乙方案,NPV,的风险度大于甲方案。,C,、利用,EXCEAL,可以很容易评价指标具体的概率分布，如表,7,，,表,7,甲乙方案风险概率分布,甲方案的概率分布,乙方案的概率分布,概率分布,NPV,概率分布,NPV,0,-2107.73,0,-12762.89,10,-237.19,10,-3897.51,17,0,20,-2000.29,20,359.32,30,-604.08,30,837.27,35,0,40,1306.85,40,509.65,50,1764.68,50,1587.07,60,2220.22,60,2574.99,70,2711.55,70,3609.49,80,3285.19,80,4694.82,90,4124.39,90,5947.03,100,7909.98,100,9593.81,因此，应该采用甲方案。,(4),总结,通过上面的分析，利用蒙特卡罗方法模拟分析得出的结果与使用传统的分析技术得出的结果相比，不仅能够分析风险因素对整个项目预期收益的影响程度，而且还能科学地估计出风险发生的概率大小，并且这样的估计是建立在充分考虑了多个风险变量共同影响、共同作用的基础之上，能够为风险决策者提供有实用价值的决策依据。因此有助于我们对多套投资方案进行筛选比较。,作业,以净现值模型为例说明蒙特卡洛方法的模拟过程。,