(完整word版)蒙特卡洛方法及其在风险评估中的应用

1、123123123蒙特卡洛方法及其应用1风险评估及蒙特卡洛方法概述1 1蒙特卡洛方法。蒙特卡洛方法，又称随机模拟方法或统计模拟方法，是在20世纪40年代随着电子计算机的发明而提出的。 它是以统计抽样理论为基础，利用随机数，经过对随机变量已有数据的统计进行抽样实验或随机模拟，以求得统计量的某个数字特征并将其作为待解决问题的数值 解。蒙特卡洛模拟方法的基本原理是：假定随机变量X1、X2、X3Xn、Y，其中X1、X2、X3Xn的概率分布已知， 且X1、X2、X3Xn、Y有函数关系：Y=F （ X1、X 2、X3Xn）, 希望求得随机变量 Y的近似分布情况及数字特征。 通过抽取符合其概率分布的随机数列

2、X1、X2、X3Xn带入其函数关系式计算获得 Y的值。当模拟的次数足够多的时候，我们就可 以得到与实际情况相近的函数 Y的概率分布和数字特征。蒙特卡洛法的特点是预测结果给出了预测值的最大值，最小值和最可能值，给出了预测值的区间范围及分布规律。1 . 2风险评估概述。风险表现为损损益的不确定性，说明风险产生的结果可能带来损失、获利或是无损失也无获利，属于广义风险。正是因为未来的不确定性使得每一个项目都存在风险。对于一个公司而言，各种投资项目通常会具有不同程度的风险，这些风险对于一个公司的影响不可小视，小到一个项目投资资本的按时回收，大到公司的总风险、公司正常运营。因此，对于风险的测量以及控制是非

3、常重要的一个环节。风险评估就是量化测评某一事件或事物带来的影响的可能程度。根据“经济人”假设，收益最大化是投资者的主要追求目标，面对不可避免的风险时， 降低风险，防止或减少损失，以实现预期最佳是投资的目标。当评价风险大小时， 常有两种评价方式：定性分析与定量分析法。 定性分析一般是根据 风险度或风险大小等指标对风险因素进行优先级排序，为进一步分析或处理风险提供参考。这种方法适用于对比不同项目的风险程度，但这种方法最大的缺陷是在于，在多个项目中风险最小者也有可能亏损。 而定量分析法则是将一些风险指标量化得到一系列的量化指标。通过这些简单易懂的指标，才能使公司的经营者、投资者对于项目分风险有正确的

4、评估与判断，采取有针对性的措施，最终做出有利于公司的决策。2蒙特卡洛方法在风险评估中的运用2.1方法简介在定量分析法下，选取一个合适的量化指标是非常重要的。对于一般的项目投资而言，项目投资回报是否能按时收回，项目是否能够为公司带来利润是决策者需要考察的问题，也就是风险。在这种情况下，这一投资未来的收益（净现值或内部收益率）以及其相对于预期 的偏离程度常常被用作衡量风险的指标。针对一个投资项目，影响未来收益的因素很多，例如，随着时间的推移，需要追加投 资数额可能会发生变化；在实业中，随着生产规模的扩大，可能出现规模经济或者规模不经 济，使得成本有所下降或上升；由于受到资金量限制， 追加投资的量会

5、受到项目回报的影响，若项目已实现的收益率达到某一标准后才继续投资，否则就退出市场。2.1.1蒙特卡洛模拟的一般步骤蒙特卡洛模拟的一般步骤如下：1. 选取随机变量，即对净现值最敏感的变量。2. 确定随机变量的概率分布3. 为各随机变量抽取随机数4. 将抽得的随机数转化为各输入变量的抽样值5. 将抽样值构成一组项目评价基础数据6. 根据基础数据计算出一种随机状况下的评价指标值7. 重复上述过程，进行反复多次模拟，得出多组评价指标值8. 整理模拟结果所得评价指标的期望值、方差、标准差、概率分布及累计概率分布，绘制累计概率图，同时，检验模拟次数是否满足预定的精度要求根据上述结果，分析各随机变量对项目收

6、益的影响。2.1.2蒙特卡洛模拟结果的分析与应用根据所得的概率分布以及概率分布图，我们可以获得关于项目未来收益的一系列评价指标，例如未来现金流的净现值的期望值、方差、标准差、净现值可能的区间以及概率。但是，公司的决策者根据公司总体情况，需要综合考虑很多其他的因素，故在了解该项目风险指标的基础上，可以根据公司现金流的需求状况、公司整体运营情况来决策。具体而言，首先，公司投资的回报需要用于弥补公司除成本外的各项费用开支，因此， 仅仅要求项目未来收益的现值为正还不能够使得公司盈利，决策者需要在了解总成本的基础上确定一个收支相抵的净现值额，再结合模拟的结果进行决策。其次，对于一个公司而言， 公司可能同

7、时有数个项目在运营中，决策者就需要考虑整个公司所有项目之间的平衡。例如，公司的某一其他项目在未来的一时间点需要一笔现金投入， 这笔现金投入来源于我们目前考察项目的资金回收。为了保证公司资金链的流畅，就需要了解项目资金回收的情况。2.2模型改进2.2.1输入变量关联性改进在项目评估中，可能有多个风险敏感变量会对目标变量造成影响，尽管蒙特卡洛方法可以设置多个风险敏感变量，但是传统的蒙特卡洛方法不考虑变量之间的关系，那么对于部分案例，我们就无法观察到风险敏感变量之间的关系。关联性改进就是通过研究风险敏感变量之间的关系，试图将变量之间的关系嵌入模型，使得模型更加完善。 最典型的例子就是规模效应。规模效

8、应是指销售量或者产量与单位可变成本之间的关系，可以分为规模经济、 规模不变以及规模不经济。规模经济就是说随着销售量或者产量的增加，单位可变成本是呈现递减的趋势；同样的，规模不经济就是单位成本随着销售量或者产量的增加而递增。2.2.2伪随机数列的改进在软件Matlab中，命令rand ()可以用来产生 0到1之间服从均匀分布随机数列，然 而这种随机数是根据一定的算法，如逆同余法、乘同余法、线性同余法等产生服从均匀分布的随机数。但上述各方法均存在一定的不足，如高维不均匀性和长周期相关性现象，会导致仿真收敛速度慢及结果波动大等一系列问题。基于上述原因，传统蒙特卡洛方法往往会造成“空隙和簇”的现象，造

9、成对采样空间的搜索不充分。为了获得分布更加均匀的数列，可以采用分布更加均匀的拟随机数列，可以使用精选的确定的样本点。而且由于拟随即序列的收敛速度要高于伪随机序列，它可以用较少的样本数就可以达到相对高的精度。3案例分析3.1案例某饮料企业现准备开发一种新型果汁饮料的投资项目，其初始投资额为200万元该项目一旦投入运营后，第一年产品的销量是一个服从均值为 200万件而标准差为60万件的正 态分布，根据这种产品的生命周期规律， 第二年销量将在第一年的基础上增长 30%,而第三 年销量将在第二年基础上增长 -20%三年内每年还需投入固定成本 100万元。新产品单位可变成本为服从2到4的均匀分配。商品零

10、售价格为服从期望为4,方差为2的正态分布。试分析此项目的风险。1 考虑将项目投资后三年内的现金流的净现值作为评估风险的依据，其中，净现值的计算公式如下：P = /+（/? -C 严卩+i）J其中：B-年现金流入，C-年现金流出，i-贴现率，n-项目寿命周期，I-项目初始投资额。2对项目的已知数据进行简单整理，考虑金钱的时间价值，设年贴现率为10%此例题只有一个变量，数据整理如下表：固定敌值的输几参数初始投熒额（百万元）2初始治堡均值（百万件）2初始销量标准差（百万杵）0.5销售第二年増长率磁销售第二年増长率-20%年固定成年（百万兀）1年贴现率10%变动数值的输入養数初始谓量（百万件）正态外苻

11、单位可变成本（元）2到4均勻分布价格（元）2至恪正态分布3依据案例，设销售量为x。x服从均值为2，标准差为0.6的正态分布；零售价格为m,服从期望为 4,方差为2的正态分布；单位可变成本为n,服从 2到4的均匀分布。写出目标函数：g=-l+ g 1（x）+ g 2（x）+ g 3（x）第一年：g1(x)=(m*x-1- n*x)*(1i) 1i(1i)第二年：g2(x)=(m*x_1- n*x)*(1i)21i(1i)2第三年：g3(x)=(mx-1- nx)*(1 i)：3 1i(1i)34.按照构建的模型，使用 Matlab进行编程、计算模拟、绘图。首先可以统计出所有模拟中的最值，期望和方

12、差，如下表：特定净现值一百万欢模押if现值均值百万兀）2.42一百万跑模拟角现值标准差陌万兀）441- 27-百万欢模拟现值最大值百万元）185- 32百万欢模拟淨现值最小值f百万兀）-153.91其次画出概率分布图和累计概率图:累计概率国由上图可以看出，在该饮料企业当前的运营情况和经济环境下，此项目投资的值大部分都落在（0,1000000 ）区间内,均值为2.42百万元，净现值大于 0的概率为53%净现值约为10,000,000的概率约为33%收益净现值达到 20,000,000的概率约为18%这表明投资项目的可行性比较高，项目投资经济上基本是安全的。而且，此案例中所设定的年贴现率为10%这

13、是一个非常高的费率，而往往年贴现率要远远低于这一水平，因此，这个项目是 个比较具有投资价值的项目。3.2案例改进 3.2.1输入变量关联性改进对于饮料项目而言，随着销量的增加，对于原材料的需求会增大，这使得企业在采购原材料时能够降低成本，同时在进行生产时所耗费的各种费用分摊到单个产品后的单位成本就 会降低，这就是规模效应。在上述案例中，当销量小于800,000时，单位可变成本为 3元，然而当销量大于800,000 后，销量每增加200,000,单位成本会降低 0.2元。假设单位可变成本为n,销量为x，则随着销量的增加，他们之间的关系可以表达为：n=3-（x-0.8）/ 0.3）*0.2）。则数

14、据表格变更为下表:jKi十概1FES呵常-ft. Fl亏列I固定数值的弔声入参数初始投资额（右万元2初始销量均值百万件12初始销量标准差（白万件）0.6销售第二年增长率30%销售第三年増长率-20%年固定成本（百万兀1年贴现率10%变动数值的辛市入参数销量（:百万件）正态分布单位可变成本（兀）随着销星降低的増加价格C元）2到&正态分布再进行模拟，统计出所有模拟中的最值，期望和方差，如下表:1特定净观值一百7T次模拟净现值均值百万兀11.24一百万次模拟淨现值方差百万兀】432. 70一百厅次複艇班植最夫值百万兀）198, 0D一百万次模拟净现值最小值（百万兀）-9&.玮其次画出概率分布图和累计

15、概率图:由概率分布图和累计概率图中我们可以看出，此项目的净现值主要落在(-5,000,000 ,40,000,000 )区间内，均值为 14.24百万元，净现值大于 0的概率为68% 接近70%可以 说项目投资的整体风险较小，适宜投资。但是，对于饮料生产企业，前期需要大量的资金投 入用于采购生产线、开拓市场，但是在后期，尤其是产销量出现大幅增长以后呈现规模效 应，单位可变成本下降，同时市场占有率不断提高会使得后期的资金投入与产出比增加， 收益也会有所增加。3.2.2伪随机数列的改进在软件Matlab中，使用命令rand ()可以产生 0到1之间的随机数，同时 matlab中 还有很多其他产生随机数的命令，例如利用命令unidrnd ( N)可以产生均匀分布(离散)随机数，利用 unifrnd(A,B) 可以产生A,B上均匀分布的随机数。这两个命令都可以使得 所产生的随机数列分布更加均匀。4总结