金母鸡量化教学场：浅析蒙特卡洛模拟（二）

　　本篇文章继续带大家浅析蒙特卡洛模拟。

　　

　　由于设计、建造和测试这三个要素都是呈标准正态分布，我们可以根据上面表格中的各自的均值和标准差数据大致画出这三个要素工期的概率分布图如下面的样子：

　　

　　我们要用蒙特卡洛模拟来定量分析整个项目的工期进度风险。于是我们用计算机来模拟项目的实施，我们的思路是： 第一步：随机选取每个WBS要素的工期值作为输入（因为每个要素的工期不是恒定的，本身就是一个估计的分布区间）； 第二步：然后把三个WBS要素的值相加得到整个项目的工期值，这样就完成了一次模拟； 第三步：重复第一二步，然后就这样一次一次的模拟，需要模拟成千上万次最终得到成千上万个整个项目总工期的数值； 第四步：再对这些海量模拟次数得到海量总工期数值进行统计分析，得出其最终的项目总工期估计的概率分布。

　　我们先做第一步。第一步需要我们先产生这些每个要素的随机工期值。Excel里面有个函数可以生成呈正态分布的随机数，就是NORMINV。我们的设计要素的第一个随机工期取值的公式就是这么写的：=ROUND(NORMINV(RAND(),$E$3,$F$3),0)，如下图所示：

　　

　　解释一下这个公式：ROUND(NORMINV(RAND(),$E$3,$F$3),0)，RAND() 是生成0到1之间的随机数，NORMINV(RAND(),$E$3,$F$3) 是生成呈均值为E3（图中为14）、标准差为F3（图中为2）的正态分布的随机数，ROUND 是四舍五入的意思，这样回车我们就生成了设计这个要素的第一个随机工期值17。同理我们把这个公式值往下拉，复制400次（我们此例中模拟400次），就得到了400个呈正态分布的随机工期值；然后建造和测试的随机工期值也是照葫芦画瓢，这样我们就得到了这3个要素的400次模拟的随机值，再每次的3个要素的随机值相加得到总工期的模拟值，如下图所示：

　　

　　此时前三步就做完了得到了总工期的一组数据（400个）。现在开始做第四步对这组数据做统计分析和作图。

　　4.1 先把总工期这一列（图中E列）400个值拷贝一份，粘贴数值到另外一列（注意粘贴的时候选择“选择性粘贴”然后选“值”，因为随机数随时变动，这儿需要把值固定下来），用MAX和MIN函数计算出这一列400个值的最大值为76，和最小值为45，作为分组依据，然后在旁边 I 列依次升序排列42-78（前后多取几个数值图像更完整）这部分数值作为分组数据，如下图所示：

　　

　　4.2 然后计算每个分组数据在总工期这组数据中出现的概率，这儿需要用到函数FREQUENCY，计算概率的公式为：=FREQUENCY(H8:H407,I8:I40)/400，意思是统计每个分组数据在总工期这一组数据中出现的次数，再除以模拟总次数400就得到这个分组数据出现的概率。再计算一个累积概率值，累积概率值就是前面的所有单个概率值加起来的概率，比如算分组数据46的累积概率值就是把46以下的数值的概率值全部加起来，这样我们就得到关于分组数据在总工期这组数据中出现的概率和累积概率的数据，如下表格所示：

　　

　　4.3 通过对上面表格的数据，以分组数据为X轴，出现概率和累积概率的值为Y轴于是可以做出下面关于总工期的概率分布图：

　　

　　红色柱状图是整个项目估计刚好多少天完工的概率数据，比如图中60天对应的概率大约是11%，表示整个项目刚好60天完工的概率是11%；蓝线就是我们PMBOK上定量风险分析得到的那张S曲线图，也是我们最终蒙特卡洛模拟需要得到的最终的模拟输出：总工期的概率分布图。通过这个S曲线，我们可以预测整个项目在多少天内完工的概率。比如我们需要预测整个项目在56天完工的概率，通过S曲线了解到56天对应的累积概率是34%左右，也就是整个项目56天内完工的概率是34%，那么56天内不能完工的概率就是1-34%=66%，这就是风险。如果觉得风险太高无法接受，那么我们可以把工期适当规划长一些，比如60天，这样查询S曲线可以得到60天内整个项目完工的概率是70%，这样就只有剩下30%的不能按时完工的风险，项目在进度方面的风险就大大降低了。

　　

　　好了，至此关于蒙特卡洛模拟的基本概念和操作流程就说完了，希望能通过本文对大家学习和理解蒙特卡洛模拟有所助益。

　　举报/反馈