第三章 蒙特卡罗仿真.ppt

第二章蒙特卡罗模拟及随机数产生一 蒙特卡罗模拟二 随机数的产生 蒙特卡罗 MonteCarlo 方法 或称计算机随机模拟方法 是一种基于 随机数 的计算方法 这一方法源于美国在第二次世界大战进研制原子弹的 曼哈顿计划 该计划的主持人之一 数学家冯 诺伊曼用驰名世界的赌城 摩纳哥的MonteCarlo 来命名这种方法 为它蒙上了一层神秘色彩 MonteCarlo方法的基本思想很早以前就被人们所发现和利用 早在17世纪 人们就知道用事件发生的 频率 来决定事件的 概率 19世纪人们用投针试验的方法来决定圆周率 考虑平面上的一个边长为1的正方形及其内部的一个形状不规则的 图形 如何求出这个 图形 的面积呢 MonteCarlo方法是这样一种 随机化 的方法 向该正方形 随机地 投掷N个点落于 图形 内 则该 图形 的面积近似为M N 随机数的产生 EXCEL中关于模拟的函数的简介 1 Rand 产生 01 上的均匀分布的随机数 ab 上均匀分布的随机数 b a rand a 2 Norminv rand mean standard dev 产生均值为mean 标准差为standard dev的正态分布随机变量 x1 x2 20 60 60 20 20 X1 X2 20 60 60 40 40 60 60 5 9 模拟技术在经营管理系统的应用 风险分析问题存储订货问题排队问题 ATM机设置 其他问题 风险分析 风险分析是在不确定条件下预测某项决策后果的过程 例子 某电脑公司开发新型打印机 初步的市场调查和财务分析给出了如下的相关数据 零售价 249美元 台第一年管理费用 400000美元第一年广告费用 600000美元 预测市场需求量为 15000台劳动力成本 45美元零件费用 90美元以上数据若均为确定数值则无风险 利润 售价 劳动力成本 零件费用 销售量 广告费用 管理费用 销售量 劳动成本 零件成本均为外部因素 是随机变量 各因素的分布情况如下 直接劳动力成本的概率分布 离散分布 零件成本 服从均匀分布 80 100元 80 90 100 需求量 服从正态分布平均值 15000标准差 4500 15000 库存问题 简单描述 公司销售某种产品 每月初订货一次 若该月货物有剩余 则每单位货物支付一定的存储费 但若不够销售则要支付不足部分的赔偿费 每月的需求量为一随机变量 预测某一存货水平下的盈利水平 某公司销售产品 售价125元 进货价格为75元 单位的存储费用为15元 单位的缺货费用为30元 需求量为服从均值为100 标准差为20的正态分布 设其补充水准为100 模拟其盈利情况 排队系统有某种服务台一个 服务某种业务 来此接受服务的顾客实行先到先服务的规则 服务台忙则排入一个队列 顾客皆为单个相互独立的到达 相继到达的时间间隔是一个随机变量 服务台对每位顾客的服务时间也是一个随机变量 对此排队服务系统的问题是想要了解服务台业务忙闲情况和排队等候情况 某银行 设有一台ATM机 顾客到达为 05 上的均匀分布 交易时间服从均值为2 标准差为0 5的正态分布 对此系统进行模拟 分析顾客必须等待的概率 及最长等待时间等问题 一风险投资机会 成功失败概率都是0 5 每投资1元 成功得到1 6元回报 原投的资本金仍归还给你 若失败 则损失1元 投资次数 金额不限 你为了不把钱输光 采用如下策略 总是拿你所持的一半去投资 假设钱是无限可分的 一直可以投资下去 不会破产 假设开始有1000000元 假设初始有a元钱 一次投资后有如下可能 1 成功 a 0 5 a 0 5 a 0 5 1 6 1 8a2 失败 0 5a一次投资后的资本期望值 E 1 8a 0 5 0 5a 0 5 1 15a两次投资后E 1 152a三次投资后E 1 153an次投资后E 1 15na10000次后E 1 1510000 1000000约等于10606元 相似性 模型研究基础相似性是模型概念和模型建立的基础 这里所说的模型与对象 系统 之间的相似性是广义的 它可以是外形上 行为上或结构上的 相似性的概念适用于非常广泛的一类物质对象 包括物理的 生物的 工程的 非工程的 当一个对象相似于另一个对象 它们之间就具有了原型和模型间的关系 如果对象A是对象B的模型 则A相似于B 可写为A B 按照相似原理 相似性是相互的 因此B A同样是存在和正确的 相似性 模型研究基础当研究对象的某一特殊问题时 常常引入简化模型的概念 简化模型仅保留对象原型的本质特征 而忽略一切非本质的特征 这一要求使得在研究复杂的对象时易于实现仿真的目的 由于简化模型忽略了对象