随机数模拟方法

合集下载

随机模拟

随机模拟

随机模拟(蒙特卡罗算法)一 随机模拟法随机模拟法也叫蒙特卡罗法,它是用计算机模拟随机现象,通过大量仿真试验,进行分析推断,特别是对于一些复杂的随机变量,不能从数学上得到它的概率分布,而通过简单的随机模拟就可以得到近似的解答。

M onte Carlo 法也用于求解一些非随机问题,如重积分、非线性方程组求解、最优化问题等。

需要指出的是,Monte Carlo 计算量大,精度也不高,因而主要用于求那些解析方法或常规数学方法难解问题的低精度解,或用于对其他算法的验证。

蒙特卡罗方法的基本思想是:当所求解问题是某种随机事件出现的概率,或者是某个随机变量的期望值时,通过某种“实验”的方法,以这种事件出现的频率估计这一随机事件的概率,或者得到这个随机变量的某些数字特征,并将其作为问题的解。

在解决实际问题的时候应用蒙特·卡罗方法主要有两部分工作: 用蒙特卡罗方法模拟某一过程时,需要产生各种概率分布的随机变量。

用统计方法把模型的数字特征估计出来,从而得到实际问题的数值解。

使用蒙特卡罗方法进行分子模拟计算是按照以下步骤进行的:使用随机数发生器产生一个随机的分子构型。

对此分子构型的其中粒子坐标做无规则的改变,产生一个新的分子构型。

计算新的分子构型的能量。

比较新的分子构型于改变前的分子构型的能量变化,判断是否接受该构型。

若新的分子构型能量低于原分子构型的能量,则接受新的构型。

若新的分子构型能量高于原分子构型的能量,则计算玻尔茲曼常数,同时产生一个随机数。

若这个随机数大于所计算出的玻尔兹曼因子,则放弃这个构型,重新计算。

若这个随机数小于所计算出的玻尔兹曼因子,则接受这个构型,使用这个构型重复再做下一次迭代。

如此进行迭代计算,直至最后搜索出低于所给能量条件的分子构型结束。

二 随机模拟法应用实例考虑二重积分(,)AI f x y dxdy =⎰⎰,其中(,)0,(,)f x y x y A ≥∀∈根据几何意义,它是以(,)f x y 为曲面顶点,A 为底的柱体C 的体积。

随机模拟的方法和应用

随机模拟的方法和应用

随机模拟的方法和应用随机模拟是一种重要的数学方法,可以用来模拟各种现实世界中复杂的系统、行为和事件。

它的应用领域广泛,包括金融、统计学、天气预测、交通规划、工程设计等多个领域。

本文将简要介绍随机模拟的基础知识以及其在不同领域的应用。

1. 随机模拟的基础知识随机模拟的实质是通过计算机程序生成的一系列随机数,来模拟真实的随机过程。

因此,随机模拟的核心是随机数生成器。

随机数生成器需要生成能够代表真实随机事件的随机数,这需要考虑一些关键问题:如何确定随机数的分布、如何生成不相关的随机数、如何满足特定的统计性质等。

常用的随机数生成方法包括线性同余发生器、Marsaglia发生器、梅森旋转游程测试以及基于物理过程的随机数发生器。

这些方法在不同场合下各有优缺点,可以根据具体需求进行选择。

随机模拟的另一个基础是随机过程的建模。

随机过程是一组与时间有关的随机变量序列,用来描述某个系统、事件或行为的随机性质。

在进行随机模拟前,需要根据实际应用建立相应的随机过程模型,通常包括确定随机变量的分布、相关性结构以及参数等。

2. 随机模拟在金融中的应用在金融领域,随机模拟被广泛应用于风险管理、资产定价、投资组合优化等方面。

随机模拟可以通过模拟不断变化的金融市场来评估不同投资策略的风险水平和收益率。

其中,蒙特卡罗模拟是一种常用的方法,它通过生成随机数对股票价格进行模拟,以此来分析不同投资组合在不同市场情况下的表现。

此外,随机模拟还可以用来构建金融风险模型,包括VaR、CVaR等风险指标。

通过随机模拟的方法,可以不断地生成样本数据,并结合实际数据来计算风险指标,从而更加准确地评估金融投资风险。

3. 随机模拟在天气预测中的应用天气预测是一项非常重要的应用领域,也是随机模拟的重要应用之一。

天气系统具有复杂的非线性关系,因此难以建立确定性模型。

随机模拟通过计算机程序模拟大气系统、海洋系统等自然系统的复杂变化,提供了一种高效、准确的天气预测方法。

蒙特卡洛随机模拟

蒙特卡洛随机模拟

蒙特卡洛随机模拟蒙特卡洛模拟法简介蒙特卡洛(Monte Carlo)方法是一种应用随机数来进行计算机摸你的方法。

此方法对研究对象进行随机抽样,通过对样本值的观察统计,求得所研究系统的某些参数。

作为随机模拟方法,起源可追溯到18世纪下半叶蒲峰实验。

蒙特卡洛模拟法的应用领域蒙特卡洛模拟法的应用领域主要有:1.直接应用蒙特卡洛模拟:应用大规模的随机数列来模拟复杂系统,得到某些参数或重要指标。

2.蒙特卡洛积分:利用随机数列计算积分,维数越高,积分效率越高。

蒙特卡洛模拟法求解步骤应用此方法求解工程技术问题可以分为两类:确定性问题和随机性问题。

解题步骤如下:1.根据提出的问题构造一个简单、适用的概率模型或随机模型,使问题的解对应于该模型中随机变量的某些特征(如概率、均值和方差等),所构造的模型在主要特征参量方面要与实际问题或系统相一致2 .根据模型中各个随机变量的分布,在计算机上产生随机数,实现一次模拟过程所需的足够数量的随机数。

通常先产生均匀分布的随机数,然后生成服从某一分布的随机数,方可进行随机模拟试验。

3. 根据概率模型的特点和随机变量的分布特性,设计和选取合适的抽样方法,并对每个随机变量进行抽样(包括直接抽样、分层抽样、相关抽样、重要抽样等)。

4.按照所建立的模型进行仿真试验、计算,求出问题的随机解。

5. 统计分析模拟试验结果,给出问题的概率解以及解的精度估计。

在可靠性分析和设计中,用蒙特卡洛模拟法可以确定复杂随机变量的概率分布和数字特征,可以通过随机模拟估算系统和零件的可靠度,也可以模拟随机过程、寻求系统最优参数等。

一. 预备知识:1.随机数的产生提示:均匀分布(0, 1)U 的随机数可由C 语言或Matlab 自动产生,在此基础上可产生其他分布的随机数. 2.逆变换法:设随机变量U 服从(0,1)上的均匀分布,则)(1U F X -=的分布函数为)(x F . 步骤:(1) 产生)1,0(U 的随机数U ;(2) 计算)(1U F X -=, 则X 服从)(x F 分布. 问题:练习用此方法产生常见分布随机数.例如“指数分布,均匀分布),(b a U ”.还有其它哪种常见分布的随机数可用此方法方便产生? 3.产生离散分布随机数已知离散随机变量X 的概率分布:)2,1(,)( ===K P x X P k k ,产生随机变量X 的随机数可采用如下算法:a) 将区间[0.1]依次分为长度为 ,,21p p 的小区间 ,,21I I ;b) 产生[0,1]均匀分布随机数R ,若k I R ∈则令k x X =,重复(b),即得离散随机变量X 的随机数序列.问题:(1) 下表给出了离散分布X 的概率分布表,试产生100个随机数.X 的概率分布表(2) 用此方法给出100个二项分布(20, 0.1)B 的随机数及10个泊松分布P(1)的随机数. 4. 正态分布的抽样提示:设21,U U 是独立同分布的)1,0(U 变量,令)2sin()ln 2()2cos()ln 2(22/11222/111U U X U U X ππ-=-=则1X 与2X 独立 ,均服从标准正态分布. 步骤:(1) 由)1,0(U 独立抽取1122,U u U u ==(2) 用(*)式计算21,x x .用此方法可同时产生两个标准正态分布的随机数.问题: 有关随机数产生方法很多,查阅相关材料进行系统总结.二. 随机决策问题1.某小贩每天以一元的价格购进一种鲜花,卖出价为b 元/束,当天卖不出去的花全部损失,顾客一天内对花的需求量是随机变量, 服从泊松分布,(),0, 1, 2,,!kP X k ek k λλ-=== .其中常数λ由多日销售量的平均值来估计, 问小贩每天应购进多少束鲜花?(准则:期望收入S(u)最高) 问题:(1) 在给定 1.25, 50b λ==的值后, 画出目标函数S(u)连线散点图, 观察单调性,给出最优决策*u ;(2) 选取其他的λ,b ,再观察S(u)的单调性;(3) 用计算机模拟方法来求出最优决策*u .对固定的u ,例如,u=40,对随机变量X 模拟100次,每次模拟得到一个收入,求出100个收入的平均值,即得到在决策u=40情况下的可能收入;(4) 对所有的可能的u ,重复(3),从中找最大的,并与(1)的结果相比较. 3.一重定积分的蒙特卡罗算法问题描述:假设函数()f x 在[,]a b 内有界连续,且()0f x ≥,求解定积分()baI f x dx =⎰.为计算出其值,可构造概率模型如下:取一个边长分别为b a -和c 的矩形D ,使曲边梯形在矩形域之内,如图2,并在矩形内随机投点,假设随机点均匀地落在整个矩形之内,则落在图中灰色区域内的随机点数k 与投点总数N 之比k/N 就近似地等于曲线下方面积(即阴影面积)与矩形面积之比,从而得出近似积分()kI b a c N≈-.图2例 求211x e--⎰由于2x e -是非初等函数,我们很难求出其原函数,所以用牛顿-莱布尼茨公式无法求解,但可以运用蒙特卡罗方法求出其近似值.将上述方法推广到一般情况:假设函数()f x 在[a ,b]内有界连续,对于定积分()baI f x dx =⎰,为计算出其值,可构造如下概率模型:取一个边长分别为b a -和c d -的矩形D ,使曲线[,]a b 段的值在矩形域之内,如图3,并在矩形内随机投点,假设随机点均匀地落在整个矩形之内,则落在图中x 轴上下灰色区域内的随机点数m 与n 的差与投点总数p 之比(m-n)/P 就近似地等于曲线上下方面积之差(即阴影面积之差)与矩形面积之比,从而得出近似积分()()m nI b a c d P-≈--.图34. 二重积分的蒙特卡罗算法问题描述:实际计算中常常要遇到如(,)Df x y dxdy ⎰⎰的二重积分,发现被积函数的原函数往往很难求出,或者原函数根本就不是初等函数,对于这样的重积分,蒙特卡罗方法也有成熟的计算方法. 方法1: 步骤:1,取一个包含D 的矩形区域Ω:,a x b c y d ≤≤≤≤,面积()()A b a d c =--;2,(,), 1,2,,i i x y i n = ,为Ω上的均匀分布随机数列,不妨设(,),1,2,i i x y i n = ()为落在D 中的n 个随机数,则n 充分大时,有1(,)(,)ki i i DA f x y dxdy f x y n =≈∑⎰⎰.方法2: 对二重积分(,)AI f x y dxdy =⎰⎰,假设(,)f x y 为区域A 上的有界函数,且(,)0f x y ≥,几何意义对应的是以(,)f x y 为曲面顶, A 为底的曲顶柱体C 的体积.因此,用均匀随机数计算二重积分的蒙特卡罗方法基本思路为:假设曲顶柱体C 包含在己知体积为DV的几何体D 的内部,在D 内产生N 个均匀随机点,统计出在C 内部的随机点数目C N ,则DC V I N N=.例:计算(1Adxdy +⎰⎰,其中22{(,)|1}A x y x y =+≤.分析:该二重积分可以看作以1+曲顶柱体在一个边长为2的立方体内,用数学分析方法可计算出其精确值为π.。

随机模拟方法总结

随机模拟方法总结

随机模拟方法总结引言随机模拟方法是一种基于概率和统计的数值计算方法,通过模拟随机事件的方式,来求解实际问题。

随机模拟方法在各个领域中都有广泛的应用,特别是在金融、物理、计算机科学和工程等领域。

本文将总结随机模拟方法的基本原理和常用的应用场景。

基本原理随机模拟方法的基本原理是通过生成服从某种概率分布的随机数,并在该分布上进行采样,来模拟实际问题。

其基本步骤如下:1.确定概率分布:根据实际问题的特点和要求,选择合适的概率分布,如均匀分布、正态分布等。

2.生成随机数:利用确定的概率分布,生成服从该分布的随机数序列。

3.采样模拟:根据具体问题,对生成的随机数进行采样模拟,得到问题的解或近似解。

4.分析结果:对采样模拟得到的结果进行统计分析,评估其准确性和可靠性。

常用应用场景随机模拟方法在各个领域中都有广泛的应用,下面列举几个常见的应用场景:金融风险评估在金融领域,随机模拟方法常用于风险评估。

通过模拟随机的市场变动、利率变化等因素,来评估投资组合的风险水平。

这些模拟结果可以帮助投资者做出更加准确的决策,降低投资风险。

物理系统模拟在物理学领域,随机模拟方法广泛应用于物理系统的建模和模拟。

通过随机模拟方法可以模拟分子动力学、粒子运动等复杂的物理现象,进一步深入理解和预测实验中观察到的现象。

计算机网络性能评估随机模拟方法可以用于评估计算机网络的性能。

通过模拟网络中的随机事件,如消息传输延迟、丢包率等,可以评估网络的性能指标,从而优化网络架构和改进网络协议。

工程系统仿真在工程领域,随机模拟方法可用于工程系统的仿真和优化。

通过模拟随机因素对工程系统的影响,可以评估系统的可靠性和性能,并进行系统优化设计。

常用模拟算法实际应用中,常用的随机模拟算法包括:•蒙特卡洛方法:通过随机采样和统计学方法,进行数值计算和模拟,如求解积分、求解微分方程等。

•马尔可夫链蒙特卡洛方法:利用马尔可夫链的性质,进行随机抽样和模拟,如在复杂系统中进行参数估计和优化。

蒙特卡洛随机模拟方法

蒙特卡洛随机模拟方法

蒙特卡洛随机模拟方法一、概述蒙特卡洛随机模拟方法是一种基于随机数的数值计算方法,它通过随机抽样来模拟实验过程,从而得到实验结果的概率分布。

在金融、物理、工程等领域有着广泛的应用。

二、基本思想蒙特卡洛随机模拟方法的基本思想是通过大量的随机抽样来模拟实验过程,从而得到实验结果的概率分布。

其主要步骤包括:1. 确定问题和目标:确定需要解决的问题和目标,例如计算某个事件发生的概率或者某个变量的期望值。

2. 建立模型:建立与问题相关的数学模型,并将其转化为计算机程序。

3. 生成随机数:根据所选用的分布函数生成符合要求的随机数。

4. 进行模拟实验:利用生成的随机数进行多次重复实验,并记录每次实验结果。

5. 统计分析:对多次重复实验结果进行统计分析,得到所需结果。

三、常用应用1. 金融领域中对衍生品价格进行估值;2. 工程领域中对结构可靠性进行评估;3. 物理领域中对粒子运动进行模拟;4. 生物领域中对药物作用机制进行研究。

四、具体步骤1. 确定问题和目标:首先需要明确需要解决的问题和目标,例如计算某个事件发生的概率或者某个变量的期望值。

2. 建立模型:建立与问题相关的数学模型,并将其转化为计算机程序。

例如,如果需要计算某个事件发生的概率,可以采用蒙特卡洛方法生成符合要求的随机数,并根据随机数判断事件是否发生。

如果需要计算某个变量的期望值,可以通过多次重复实验得到该变量在不同条件下的取值,并根据统计学原理计算其期望值。

3. 生成随机数:根据所选用的分布函数生成符合要求的随机数。

常见的分布函数包括均匀分布、正态分布、指数分布等。

4. 进行模拟实验:利用生成的随机数进行多次重复实验,并记录每次实验结果。

通常情况下,需要进行大量重复实验才能得到准确可靠的结果。

5. 统计分析:对多次重复实验结果进行统计分析,得到所需结果。

常见的统计分析方法包括求和、平均值、方差等。

五、优缺点1. 优点:蒙特卡洛随机模拟方法具有灵活性、精度高、适用范围广等优点,可以处理各种复杂问题,并且可以通过增加样本容量来提高精度。

蒙特卡罗(Monte Carlo)方法简介

蒙特卡罗(Monte Carlo)方法简介

蒙特卡罗(Monte Carlo)方法简介蒙特卡罗(Monte Carlo)方法简介蒙特卡罗(Monte Carlo)方法,也称为计算机随机模拟方法,是一种基于"随机数"的计算方法。

一起源这一方法源于美国在第二次世界大战进研制原子弹的"曼哈顿计划"。

Monte Carlo方法创始人主要是这四位:Stanislaw Marcin Ulam, Enrico Fermi, John von Neumann(学计算机的肯定都认识这个牛人吧)和Nicholas Metropolis。

Stanislaw Marcin Ulam是波兰裔美籍数学家,早年是研究拓扑的,后因参与曼哈顿工程,兴趣遂转向应用数学,他首先提出用Monte Carlo方法解决计算数学中的一些问题,然后又将其应用到解决链式反应的理论中去,可以说是MC方法的奠基人;Enrico Fermi是个物理大牛,理论和实验同时都是大牛,这在物理界很少见,在“物理大牛的八卦”那篇文章里提到这个人很多次,对于这么牛的人只能是英年早逝了(别说我嘴损啊,上帝都嫉妒!);John von Neumann可以说是计算机界的牛顿吧,太牛了,结果和Fermi一样,被上帝嫉妒了;Nicholas Metropolis,希腊裔美籍数学家,物理学家,计算机科学家,这个人对Monte Carlo方法做的贡献相当大,正式由于他提出的一种什么算法(名字忘了),才使得Monte Carlo方法能够得到如此广泛的应用,这人现在还活着,与前几位牛人不同,Metropolis很专一,他一生主要的贡献就是Monte Carlo方法。

蒙特卡罗方法的名字来源于摩纳哥的一个城市蒙地卡罗,该城市以赌博业闻名,而蒙特•罗方法正是以概率为基础的方法。

与它对应的是确定性算法。

二解决问题的基本思路Monte Carlo方法的基本思想很早以前就被人们所发现和利用。

早在17世纪,人们就知道用事件发生的"频率"来决定事件的"概率"。

随机数生成公式

随机数生成公式

随机数生成公式随机数生成公式是一种计算机程序中常用的技术,可以生成随机的数字,用于模拟和实验等场景中。

本文将介绍几种常见的随机数生成公式及其应用场景。

一、线性同余法(Linear Congruential Method)线性同余法是一种简单而又高效的随机数生成方法,其公式为:Xn+1 = (aXn + c) mod m其中Xn为当前随机数,a、c、m为常数,mod为模运算符。

该公式的原理是通过不断迭代计算,每次得到一个新的随机数。

该方法的优点是计算速度快,缺点是会产生周期性重复的随机数序列。

该方法常用于模拟和实验场景中。

二、梅森旋转算法(Mersenne Twister)梅森旋转算法是一种广泛应用的随机数生成方法,其公式为:Xn+1 = Xn⊕(Xn >> u)其中Xn为当前随机数,⊕为异或运算符,>>为右移运算符,u为常数。

该公式的原理是通过对当前随机数进行位运算,得到一个新的随机数。

该方法的优点是生成的随机数序列较为均匀,缺点是计算速度较慢。

该方法常用于加密和安全场景中。

三、高斯分布随机数生成公式(Gaussian Distribution)高斯分布随机数生成公式是一种生成符合正态分布(高斯分布)的随机数的方法,其公式为:X = μ + σ * Z其中μ为均值,σ为标准差,Z为符合标准正态分布的随机数。

该公式的原理是通过对标准正态分布进行线性变换,得到符合正态分布的随机数。

该方法的优点是生成的随机数符合实际分布规律,缺点是计算量较大。

该方法常用于金融和统计场景中。

四、指数分布随机数生成公式(Exponential Distribution)指数分布随机数生成公式是一种生成符合指数分布的随机数的方法,其公式为:X = -ln(U) / λ其中U为符合均匀分布的随机数,ln为自然对数函数,λ为指数分布的参数。

该公式的原理是通过对均匀分布进行变换,得到符合指数分布的随机数。

随机现象模拟

随机现象模拟

一、 随机数的生成与随机现象的模拟从建模的角度看,现实对象都不可避免的具有一些不确定性,我们所涉及的变量都是随机变量,往往建立的是随机模型。

我们对随机模型的变量进行模拟,或者对随机系统进行模拟,例如渡口模型,超市收银台顾客到达时间模型,都需要产生服从某种分布的一系列随机数,亦或者服从某种分布的一系列样本值。

随机数在概率算法设计中扮演着十分重要的角色。

在现实计算机上无法产生真正的随机数,因此在概率算法中使用的随机数都是一定程度上随机的,即伪随机数。

在产生随机数的数学方法中,同余法是使用较为广泛的方法,其中以乘同余法和混合同余法应用最广,因为他们具有很好的统计性质。

产生(0,1)内均匀分布的随机数的递推公式:)(mod 1M X X n n λ=+M X r n n /=其中,λ为乘因子,M 是模数。

第一式称做以M 为模数的同余式,当给定了一个初始值0X ,计算出的n r r r ,,21即(0,1)上均匀分布的随机数序列。

例题:取1000,7,10===M X λ,有:807.0,401.0,343.0,049.0,007.054321=====r r r r r很多数学软件都可以产生随机数,但对于特殊分布或者要求数据量很大时就不太有效,仍然需要寻求一种简便、经济、可靠的并能在计算机上实现产生随机数的方法。

Matlab产生随机数Unifrnd unifrnd ( A,B,m,n) [A,B]上均匀分布(连续) 随机数Unidrnd unidrnd(N,m,n) 均匀分布(离散)随机数Exprnd exprnd(Lambda,m,n) 参数为Lambda的指数分布随机数Normrnd normrnd(MU,SIGMA,m,n) 参数为MU,SIGMA的正态分布随机数chi2rnd chi2rnd(N,m,n) 自由度为N的卡方分布随机数Trnd trnd(N,m,n) 自由度为N的t分布随机数Frnd frnd(N1, N2,m,n) 第一自由度为N1,第二自由度为N2的F分布随机数gamrnd gamrnd(A, B,m,n) 参数为A, B的分布随机数betarnd betarnd(A, B,m,n) 参数为A, B的分布随机数lognrnd lognrnd(MU, SIGMA,m,n) 参数为MU, SIGMA的对数正态分布随机数nbinrnd nbinrnd(R, P,m,n) 参数为R,P的负二项式分布随机数ncfrnd ncfrnd(N1, N2, delta,m,n) 参数为N1,N2,delta的非中心F分布随机数nctrnd nctrnd(N, delta,m,n) 参数为N,delta的非中心t分布随机数ncx2rnd ncx2rnd(N, delta,m,n) 参数为N,delta的非中心卡方分布随机数raylrnd raylrnd(B,m,n) 参数为B的瑞利分布随机数weibrnd weibrnd(A, B,m,n) 参数为A, B的韦伯分布随机数binornd binornd(N,P,m,n) 参数为N, p的二项分布随机数geornd geornd(P,m,n) 参数为p的几何分布随机数hygernd hygernd(M,K,N,m,n) 参数为M,K,N的超几何分布随机数Poissrnd poissrnd(Lambda,m,n) 参数为Lambda的泊松分布随机数(1)随机现象的模拟例:超市出口有若干个收款台,两项服务:收款、装袋。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

知识迁移
例1 天气预报说,在今后的三天中,每一 天下雨的概率均为40%,用随机模拟方法估 计这三天中恰有两天下雨的概率约是多少?
要点分析: (1)今后三天的天气状况是随机的, 共有四种可能结果,每个结果的出现 不是等可能的. (2)用数字1,2,3,4表示下雨,数 字5,6,7,8,9,0表示不下雨,体现 下雨的概率是40%.
思考5:一般地,如果一个古典概型的基 本事件总数为n,在没有试验条件的情况 下,你有什么办法进行m次实验,并得到 相应的试验结果?在几何概型中的变量 呢?
思考6:如果一次试验中各基本事件不都 是等可能发生,如何利用上述方法获得 的试验结果?
探究(二):随机模拟方法 思考1:对于古典概型,我们可以将随机 试验中所有基本事件进行编号,利用计 算器或计算机产生随机数,从而获得试 验结果.这种用计算器或计算机模拟试 验的方法,称为随机模拟方法或蒙特卡 罗方法(Monte Carlo).你认为这种方 法的最大优点是什么? 不需要对试验进行具体操作,可以广 泛应用到各个领域.
(3)用计算机产生三组随机数,代表 三天的天气状况. (4)产生30组随机数,相当于做30次 重复试验,以其中表示恰有两天下雨的 随机数的频率作为这三天中恰有两天下 雨的概率的近似值. Excel演示
(5)据有关概率原理可知,这三天中 恰有两天下雨的概率 P=3掷两粒骰子,计算出现点数之 和为7的概率,利用随机模拟方法试验 200次,计算出现点数之和为7的频率, 并分析两个结果的联系和差异. 模拟试验
思考2:随机数表中的数是0~9之间的随 机数,你有什么办法得到随机数表? 我们可以利用计算器产生随机数,其 操作方法见教材P130及计算器使用说 明书. 我们也可以利用计算机产生随机数, 用Excel演示:
思考3:若抛掷一枚均匀的骰子30次,如 果没有骰子,你有什么办法得到试验的 结果? 思考4:若抛掷一枚均匀的硬币50次,如 果没有硬币,你有什么办法得到试验的 结果?
y
1
模拟试验
1 x
-1
0
以直线x=1,x=-1,y=0,y=1为边界作矩形, 用随机模拟方法计算落在抛物区域内的均匀 随机点的频率,则所求区域的面积=频率×2.
小结作业 1.用计算机或计算器产生的随机数,是 依照确定的算法产生的数,具有周期性 (周期很长),这些数有类似随机数的 性质,但不是真正意义上的随机数,称 为伪随机数.
3.4
随机数的产生
问题提出
1.通过大量重复试验,反复计算事件 发生的频率,再由频率的稳定值估计概 率,是十分费时的.对于实践中大量非古 典概型的事件概率,又缺乏相关原理和 公式求解.因此,我们设想通过计算机模 拟试验解决这些矛盾.
探究1:随机数的产生
思考1:对于某个指定范围内的整数,每 次从中有放回随机取出的一个数都称为 随机数. 那么你有什么办法产生1~20之 间的随机数 . 抽签法
思考2:用随机模拟方法抛掷一枚均匀的 硬币100次,那么如何统计这100次试验 中“出现正面朝上”的频数和频率.演示. 思考3:把抛掷两枚均匀的硬币作为一次 试验,则一次试验中基本事件的总数为 多少?若把这些基本事件数字化,可以 怎样设置? 思考4:用随机模拟方法抛掷两枚均匀的 硬币100次,如何估计出现一次正面和一 次反面的概率? 演示.
2.随机模拟方法是通过将一次试验所有
等可能发生的结果数字化,由计算机或 计算器产生的随机数,来替代每次试验 的结果,其基本思想是用产生随机数的 频率估计事件发生的概率,这是一种简 单、实用的科研方法,在实践中有着广 泛的应用.
作业:
P134 A组: 5,6. B组: 1,2.
例3 在下图的正方形中随机撒一把豆 子,如何用随机模拟的方法估计圆周率 的值.
(1)圆面积︰正方形面 积=落在圆中的豆子数︰ 落在正方形中的豆子数. (2)设正方形的边长为2, 则 落在圆中的豆子数 p= ÷落在正方形中的豆子数 ×4.
例4 利用随机模拟方法计算由y=1和 y=x2 所围成的图形的面积.
思考5:假设你家订了一份报纸,送报人 可能在早上6:30~7:30之间把报纸送到 你家,你父亲离开家去上班的时间在早 上7:00~8:00之间,如果把“你父亲在 离开家之前能得到报纸”称为事件A,那 么事件A是哪种类型的事件? 随机事件
思考6:设X、Y为[0,1]上的均匀随机数, 6.5+X表示送报人到达你家的时间,7+ Y表示父亲离开家的时间,若事件A发生, 则X、Y应满足什么关系? 7+Y >6.5+X,即Y>X-0.5. 思考7:如何利用计算机做100次模拟试 验,计算事件A发生的频率,从而估计事 件A发生的概率? 模拟试验
相关文档
最新文档