当前位置: 首页 > news >正文

蒙特卡洛随机模拟

蒙特卡洛随机模拟

简介

蒙特卡洛模拟是在计算机上模拟项目实施了成千上万次,每次输入都随机选择输入值。由于每个输入很多时候本身就是一个估计区间,因此计算机模型会随机选取每个输入的该区间内的任意值,通过大量成千上万甚至百万次的模拟次数,最终得出一个累计概率分布图,这个就是蒙特卡洛模拟。

模拟举例

我们以定量风险分析项目总成本为例来简要介绍一下蒙特卡洛模拟。你做为项目经理,项目主要有三大阶段需要花费成本,分别是计划,编码,测试。假设这三块内容预估费用的概率是正态分布(单位:万RMB)。分别对三要素得到如下最乐观,最可能和最悲观的成本数值。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-yRhWrJsw-1677461323969)(C:\Users\why\AppData\Roaming\Typora\typora-user-images\image-20230116195552464.png)]

现在我们需要用蒙特卡洛模拟来以这三个要素的成本分布为输入,来模拟得到整个项目的成本概率分布图。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-LQxish6J-1677461323970)(C:\Users\why\AppData\Roaming\Typora\typora-user-images\image-20230116195624223.png)]

我们要用蒙特卡洛模拟来定量分析整个项目的成本风险。于是我们用电脑来模拟项目的实施。

第一步:随机抽取三要素(计划、编码、测试)的成本值作为输入;

(对于第一步中,产生每个要素的随机成本值,Excel里面有个函数可以生成呈正态分布的随机数,就是NORMINV。)

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-nOSab5an-1677461323971)(C:\Users\why\AppData\Roaming\Typora\typora-user-images\image-20230116195801827.png)]

第二步:然后把三个要素的值求和Σ得到整个项目的成本值,这样就完成了一次模拟;

第三步:对第一二步进行重复套娃,然后就这样一次一次的模拟,需要模拟成千上万次最终得到成千上万个整个项目成本的数值(我们这个例子中套娃500次);

第三步中对三要素重复500次生成随机值,就得到了500组三要素的随机值,并将它们的数值相加得到500组随机总成本。如下图:

img

第四步:再对这些海量成本值进行统计分析,得出最终的项目总成本概率分布。

在这里插入图片描述

通过对上面表格的数据可以计算出每个随机总成本出现的概率(和累计概率)。将分组数据为X轴,出现概率和累积概率的值为Y轴,可以得到总成本的概率分布图:

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-0nWPiSfR-1677461323975)(C:\Users\why\AppData\Roaming\Typora\typora-user-images\image-20230116195942101.png)]

柱状图是整个项目估计刚好花多少成本完工的概率数据,比如图中花28万对应的概率大约是14%,表示整个项目刚好花28万完工的概率是14%。

图中的曲线是每个柱状图(每一次总成本)概率的累计值得到的曲线,是总成本的累计概率分布图。比如还是以28万这个点为例,X轴28万,曲线在Y轴上对应的是70%的概率,这就代表着总成本在28万的话,有70%的概率能够完成项目。

这个曲线就是模特卡洛模拟所得到的交付物,最终我们可以根据这个模拟曲线,来判断花多少成本,能够完成项目的概率是多大,这样就将风险量化了。

http://www.lryc.cn/news/22425.html

相关文章:

  • Android从屏幕刷新到View的绘制(三)之Handler异步消息与同步屏障
  • 最新版axios@1.3.x取消请求-AbortController-初体验-番茄出品
  • Git的简述
  • webpack实战,手写loader和plugin
  • STM32CubeMX按键模块化 点灯
  • C#专栏目录(长期更新)
  • BurpSuite配置抓取HTTPS数据包
  • 图片转base64格式返回给前端,前端如何展示?
  • C++入门知识【超详解】
  • 零基础、非计算机系学Python该如何上手?
  • 关于 vue3 模板引用
  • Redis | 安装Redis和启动Redis服务
  • 博客要考虑的最佳WordPress主题
  • C 学习笔记 —— 函数指针
  • FastDDS-3. DDS层
  • 9.2 IGMPv2
  • 巨头混战,抢着“兜底”自动驾驶安全
  • RightCapital 第一轮面试题
  • Python曲线肘部点检测-膝部点自动检测
  • 【算法题】最大矩形面积,单调栈解法
  • 活动策划|深度分析年货节活动该如何策划!
  • Idea启动遇到 Web server failed to start. Port 8080 was already in use. 报错
  • Python3中zip()函数知识点总结
  • 过滤器,监听器,拦截器的原理与在Servlet和Spring的应用
  • minio spring boot 秒传、分片上传、断点续传文件实现
  • MTK平台使用Omnipeek分析空口协议讲解
  • string和自动推断类型
  • 【软件测试】从功能到自动化测试,测试人的进阶之路细节,这些必不可少......
  • C语言青蛙跳台阶【图文详解】
  • 笔记(五)——list容器的基础理论知识