- 1、本文档共41页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Monte Carlo 模拟 内容提纲 1.引言 2.Monte Carlo模拟基本思想 3.随机数生成函数 4.应用实例举例 5.排队论模拟 6.Monte Carlo模拟求解规划问题 Monte Carlo方法的基本思想 一些人进行了实验,其结果列于下表 : 基本思想 由上面的例子可以看出,当所求问题的解是某个事件的概率,或者是某个随机变量的数学期望,或者是与之有关的量时,通过某种试验的方法,得出该事件发生的频率,再通过它得到问题的解。这就是蒙特卡罗方法的基本思想。 非常见分布的随机数的产生 逆变换方法 Acceptance-Rejection 方法最早由 Von Neumann提出,现在已经广泛应用于各种随机数的生成。 为要模拟服从给定分布的随机变量,用生成一个易于实现的不可约遍历链 作为随机样本,使其平稳分布为给定分布的方法,称为马氏链蒙特卡罗方法. [2] 离散型随机变量(以离散分布为例): 小实例一:投掷硬币 小实例二:投掷骰子 小实例五:生日问题 试验五源程序 小实例五的理论值计算 与Monte Carlo方法相似,但理论基础不同的方法—“拟蒙特卡罗方法”(Quasi-Monte Carlo方法)—近年来也获得迅速发展。这种方法的基本思想是“用确定性的超均匀分布序列(数学上称为Low Discrepancy Sequences)代替Monte Carlo方法中的随机数序列。对某些问题该方法的实际速度一般可比Monte Carlo方法提高数百倍。 排队问题随机模拟 排队论主要研究随机服务系统的工作过程。 在排队系统中,服务对象的到达时间和服务时间都是随机的。排队论通过对随机服务现象的统计研究,找出反映这些随机现象平均特性的规律指标,如排队的长度、等待的时间及服务利用率等。 [1] 系统的假设: (1) 顾客源是无穷的; (2) 排队的长度没有限制; (3)到达系统的顾客按先后顺序依次进入服务,“先到先服务”。 在某商店有一个售货员,顾客陆续来到,售货员逐个地接待顾客.当到来的顾客较多时,一部分顾客便须排队等待,被接待后的顾客便离开商店.设: 1.顾客到来间隔时间服从参数为0.1的指数分布. 2.对顾客的服务时间服从[4,15]上的均匀分布. 3.排队按先到先服务规则,队长无限制. 假定一个工作日为8小时,时间以分钟为单位。 [1]模拟一个工作日内完成服务的个数及顾客平均等待时间t. [2]模拟100个工作日,求出平均每日完成服务的个数及每日顾客的平 均等待时间。 单服务员的排队模型模拟 w:总等待时间; ci:第i个顾客的到达时刻; bi:第i个顾客开始服务时刻; ei:第i个顾客服务结束时刻; xi:第i-1个顾客与第i个顾客之间到达的间隔时间; yi:对第i个顾客的服务时间。 符号说明 c1 b1 c3 c4 c5 c2 e1 b2 e2 b3 e3 b4 e4 b5 ci=ci-1+ xi ei=bi+yi bi=max(ci,ei-1) t 思路分析 初始化:令i=1,ei-1=0,w=0 产生间隔时间随机数xi~参数为0.1的指数分布 ci=xi , bi=xi 产生服务时间随机数yi~[4,15]的均匀分布 ei=bi+yi 累计等待时间:w=w+bi-ci 准备下一次服务:i=i+1 产生间隔时间随机数xi~参数为0.1的指数分布 ci=ci-1+ xi 确定开始服务时间:bi=max(ci,ei-1) bi>480? Y N i=i-1,t=w/i 输出结果:完成服务个数:m=i 平均等待时间:t 停止 流程框图 用蒙特卡洛法解非线性规划问题 基本假设 试验点的第j个分量xj服从[aj ,bj]内的均匀分布. 符号假设 求解过程 先产生一个随机数作为初始试验点,以后则将上一个试验点的第j个分量随机产生,其它分量不变而产生一新的试验点.这样,每产生一个新试验点只需一个新的随机数分量.当KMAXK或PMAXP时停止迭代. 框 图 初始化:给定MAXK,MAXP;k=0,p=0,Q:大整数 xj=aj+R(bj-aj) j=1,2,…,n j=0 j=j+1,p=p+1 PMAXP? Y N xj=aj+R(bj-aj) gi(X)≥0? i=1,2…n Y N jn? f(X)≥Q? Y N X*=X,Q=f(X) k=k+1 kMAXK? Y N 输出X,Q,停止 Y N 在Matlab软件包中编程,共需三个M-文件:randlp.m, mylp.m, lpconst.m.主程序为randlp.m. %
您可能关注的文档
- MSA第三版讲义.ppt
- 毕业设计规范性.doc
- 毕业设计——回转轴零件的数控车床加工.doc
- 毕业设计-火力发电厂锅炉给水处理设计.doc
- 2017版高考历史大二轮总复习与增分策略板块二中国近现代史第7讲改革开放前的现代中国(1949~1978年).ppt
- 毕业设计——基于NX6.0的过滤板工艺设计及仿真加工.doc
- MSA测量系统(第四版).ppt
- 毕业设计——皮带输送机滚筒轴的加工工艺.doc
- MRO行业基本情况.ppt
- 毕业设计——塑料水杯的注塑模设计.doc
- 新的一年工作展望.docx
- 医生年终个人工作的述职报告(3篇).docx
- 2023年消防设施操作员之消防设备中级技能考前冲刺练习题附答案详解.docx
- 2022-2023年环境影响评价工程师之环评技术导则与标准通关练习题包括详细解答.docx
- 2023年中级注册安全工程师之安全生产管理考前冲刺检测卷和答案.docx
- 2023年中级银行从业资格之中级银行管理考前冲刺测试卷提供答案解析.docx
- 2023年公共营养师之二级营养师通关模拟考试试卷附带答案.docx
- 证券分析师之发布证券研究报告业务考前冲刺模拟题库.docx
- 2022-2023年二级建造师之二建建设工程法规及相关知识综合提升测试卷附答案.docx
- 2023年二级建造师之二建机电工程实务通关模拟考试试卷提供答案解析.docx
文档评论(0)