加入收藏

蒙特卡罗方法的核心思想 蒙特卡罗方法

2023-08-27 17:37:25 来源:互联网

1、蒙特卡罗方法(Monte Carlo method)是一种通过随机变量的数字模拟和统计分析来求取数学物理、工程技术问题近似解的数值方法,利用这种方法求解问题的过程可以归纳为下列三个基本步骤:(1)随机变量的抽样试验。


(资料图片)

2、按基本随机变量(输入随机变量)的已知概率分布进行随机抽样(数字模拟)。

3、(2)样本反应求解。

4、对每个抽取的样本,按问题的性质采用确定性的控制数学、物理方程求取样本反应。

5、(3)计算反应量的统计量估计。

6、对所有样本反应,按所求解答的类型分别求取输出随机变量的均值、方差或概率分布。

7、当求解确定性问题时,首先,要根据所提出的问题构造一个简单、适用的概率模型,使问题的解对应于该模型中随机变量的某些数字特征(如概率、数学期望、方差等);然后,在高速运行的计算机上生成随机数,并对随机数进行统计分析试验;最后,利用试验所获结果求出统计特征的估计值作为问题的近似解。

8、总结以上思想,可以得出利用蒙特卡罗方法求解确定性问题的基本步骤为:(1)根据所要求解的实际问题来构造概型,并使概型的某些统计特征恰好相当于所要求的问题的解。

9、(2)根据所建立的概率模型,设计、使用一些加速收敛的方法,以求加速收敛并提高计算精度。

10、(3)给出在计算机上产生概型中各种不同分布随机变量的方法。

11、(4)统计处理模拟结果,给出问题的近似解并做解的精度估计。

12、蒙特卡罗方法虽然可以求解许多确定性工程技术问题,但其独到之处还应该在于求解随机性问题。

13、用蒙特卡罗方法求解随机性问题时,一般首先,根据问题的物理性质建立随机模型;然后,再根据模型中各个随机变量的分布,在计算机上产生随机数,进行大量的统计试验,以取得所求问题的大量试验值;最后,根据这些试验结果求它的统计特征量,从而获得所求问题的解。

14、由此可见,用蒙特卡罗方法求解随机问题的步骤与求解确定性问题的步骤基本一致。

15、总之,蒙特卡罗方法的理论基础是概率论中的大数定律。

16、设在N次独立试验中,n为事件A出现的次数,而P(A)为事件A在每次试验中出现的概率,贝努利大数定律指出,对于任意ε>0,当 N→∞时,事件 A 出现的频率的概率收敛于事件的概率。

17、即地下水系统随机模拟与管理当随机变量满足独立分布时,若随机变量序列ξ1,ξ2,…,ξN的分布相同,ξi具有有限的数学期望E(ξi)=a,i=1,2,…,N,则根据柯欠莫哥洛夫大数定律,对于任意的ε>0,当N→∞时,变量ξi 将以概率1收敛于期望值 a,即地下水系统随机模拟与管理在蒙特卡罗方法中,采用简单抽样方法进行随机变量的数字模拟,因此其所抽取的子样为具有同分布性质的独立随机变量,当抽取的样本个数足够大时,样本均值将以概率1收敛于分布均值,而事件 A 出现的频率则以概率收敛于事件A 出现的概率,这样就保证了蒙特卡罗方法的概率收敛性。

18、2.1.1 均匀分布随机数的生成根据所求解问题性质的不同,其基本随机变量可能属于不同的概率分布,为了产生不同分布类型的随机变量的抽样值(随机数),一般需先产生一个在[0,1]上均匀分布的随机变量的抽样值,然后按照给定的概率分布类型将其转化为所需随机变量的抽样值。

19、因此,均匀分布随机变量随机数的生成是蒙特卡罗方法实现的基础。

20、利用数值法产生的均匀随机变量的抽样值称之为伪随机数,这是因为数值方法的基础是某一数学递推公式,按这类递推公式产生的抽样与[0,1]均匀分布中的抽样在统计性质上不可能完全相同。

21、数学递推公式的一般形式是:地下水系统随机模拟与管理式中:f(xn,xn-1,…,xn-k)——某一给定的函数形式。

22、根据这一函数式,当给定一组初值,x0,x-1,…,x-k后,便可依次求出x1,x2,…,xm…最常用的(0,1)均匀分布随机数生成的递推公式有:(1)乘同余法。

23、用以产生(0,1)均匀分布随机数的递推公式为:地下水系统随机模拟与管理式中:λ,M和x0——预先给定的常数。

24、式(2.4)的意义是指以 M 除以λxi-1后得到的余数记为 xi。

25、由于是余数,所以,即有:地下水系统随机模拟与管理如此所得的随机数序列r1,r2,…,ri为具有(0,1)均匀分布的随机数。

26、由式(2.4)不难看出,不同的xi最多只能有M个,相应地不同的随机数ri也最多只能有M个。

27、所以当产生的随机数ri个数多于M个时,就会出现循环数,这样,便再不能看成是随机数。

28、为了使所产生的随机数能经得住数理统计中的独立性和均匀性检验,需要合理选择随机数生成参数x0,λ及M。

29、表2.1所列为几个经过检验的参数,以供参考。

30、表2.1(2)混合同余法。

31、混合同余法的递推公式为:地下水系统随机模拟与管理通过适当地选取参数,可以改变伪随机数的统计性质。

32、其他有关伪随机数的生成技术读者可参阅文献[32,41]。

33、2.1.2 任意分布随机数的生成任意分布随机数的生成是以(0,1)均匀分布随机数为基础,通过适当的数学变换来形成。

34、可以证明有下列任意分布随机数生成公式。

35、(1)(a,b)上均匀分布随机数的生成公式为:地下水系统随机模拟与管理(2)具有指数分布概率密度f(x)=λe-λx(x≥0)的随机数生成公式为:地下水系统随机模拟与管理(3)正态分布N(0,1)随机数生成公式为:地下水系统随机模拟与管理(4)正态分布N(μ,σ)随机数生成公式为:将式(2.8)的xi代入式:地下水系统随机模拟与管理即可得 N(μ,σ)分布随机数上述各式中的ri 为(0,1)均匀分布随机数。

36、2.1.3 随机数的统计检验为了进一步了解所生成的随机数是否具有我们所需要的随机数特性,往往需要对所生成的随机数进行参数检验,均匀性检验和独立性检验。

37、参数检验主要是为了检验随机数的子样均值和理论均值的差异是否显著,(0,1)上均匀分布的随机变量R的期望值和方差分别为:地下水系统随机模拟与管理地下水系统随机模拟与管理设随机变数R共有n个观测值r1,r2,…,rn,则由中心极限定理得知:式中:地下水系统随机模拟与管理渐近服从标准正态分布 N(0,1),可以进行 U 检验。

38、当给定显著性水平后,即可根据正态分布表确定临界值,据此判断-r 与其期望值E(R)之差异是否显著,从而决定能否把 r1,r2,…,rn看做是(0,1)均匀分布随机变量 R 的n 个独立取值。

39、均匀性检验又称频率检验,它检验随机数的经验频率与理论频率的差异是否显著。

40、把(0,1)区间分成 k 等份,以(i=1,2,…,k)表示第 i 个小区间,如 rs 是(0,1)上均匀分布的随机变量 R 的一个取样值,则它落在任一小区间的概率 Pi均匀等于这些小区间的长度,故 n 个值落在任一个小区间的平均数为mi=nPi=n/k,设 n 个rs 值落入第i 个小区间有ni个,则统计量:地下水系统随机模拟与管理渐近地服从χ2(k-1)分布。

41、据此可进行显著性检验。

42、独立性检验主要是检验随机数r1,r2,…,中前后各数的统计相关性是否显著。

43、两个随机变数的相关系数反映它们之间的线性相关程度,若两个随机变数相互独立,则它们的相关系数ρK=0,故可通过相关系数来检验随机数的独立性。

44、设给定n个随机数r1,r2,…,rn,前后距离为k的样本相关系数的计算公式为:式中:地下水系统随机模拟与管理当独立性假设(ρ=0)成立时,则当 n 充分大(如 n>50+k)时,统计量 U=渐近地服从标准正态分布N(0,1),故可进行 U 检验。

本文到此分享完毕,希望对大家有所帮助。

关键词:

相关新闻

资讯

新疆呼图壁:海棠果获丰收 农户采摘忙
新疆呼图壁:海棠果获丰收 农户采摘忙

26日,在新疆昌吉回族自治州呼图壁县五工台镇龙王庙村......更多>

139亿元!湖北又一高速公路中标信息公布
139亿元!湖北又一高速公路中标信息公布

139亿元!湖北又一高速公路中标信息公布【来源:中国......更多>

俄国防部称在布良斯克州上空击落乌无人机 乌方暂无回应
俄国防部称在布良斯克州上空击落乌无人机 乌方暂无回应

8月27日电,俄罗斯国防部称,莫斯科时间8月26日21时45......更多>

阳光股份(000608.SZ)半年报净亏损2206.94万元,同比收窄
阳光股份(000608.SZ)半年报净亏损2206.94万元,同比收窄

格隆汇8月25日丨阳光股份(000608 SZ)披露2023年半年......更多>

蓝色协议错误代码102解决方法攻略
蓝色协议错误代码102解决方法攻略

相信有许多蓝色协议用户不知道蓝色协议错误代码102解......更多>

青海尖扎:乡村公路联通美景 串起群众幸福生活
青海尖扎:乡村公路联通美景 串起群众幸福生活

走进青海省黄南藏族自治州尖扎县昂拉乡德吉村,这座位......更多>

孕妈妈科学饮食与营养配餐金典(关于孕妈妈科学饮食与营养配餐金典简述)
孕妈妈科学饮食与营养配餐金典(关于孕妈妈科学饮食与营养配餐金典简述)

,你们好,今天0471房产来聊聊一篇妈妈科学饮食与营养......更多>

惠程科技2023年半年度董事会经营评述
惠程科技2023年半年度董事会经营评述

惠程科技2023年半年度董事会经营评述内容如下:  一......更多>

淡水鲍鱼价格行情 日本水产品价格暴跌 鲍鱼已降价30% 基本情况讲解
淡水鲍鱼价格行情 日本水产品价格暴跌 鲍鱼已降价30% 基本情况讲解

大家好,今日关于【淡水鲍鱼价格行情日本水产品价格暴......更多>

端午节的由来 端午节的习俗
端午节的由来 端午节的习俗

你们好,最近0471房产发现有诸多的对于端午节的由来,......更多>

关注

【短视频】甘肃紧急下拨5000万资金用于抗旱救灾
【短视频】甘肃紧急下拨5000万资金用于抗旱救灾
甘肃省财政厅、省农业农村厅、省水利厅三部门日前从农... 更多>
【短视频】甘肃紧急下拨5000万资金用于抗旱救灾
甘肃省财政厅、省农业农村厅、省水利厅三部门日前从农... 更多>
地铁五号线东延段完成综合监控系统电力大联调
近日,广州地铁五号线东延段(文冲-黄埔新港)综合监... 更多>
三预警齐发!南方水汽充沛 川渝陕豫等地将有强降雨
预计,8月24日20时至25日20时,四川东部、重庆西部和... 更多>
宁波银行(002142):盈利增长稳健 存贷扩张强劲
宁波银行(002142):盈利增长稳健存贷扩张强劲 更多>
格利特是什么材质
格利特是一种由纤维素纤维制成的材质。纤维素是由植物... 更多>
拒绝广告打扰!魅族Flyme 10系统,给你0广告的清爽体验!
在当今这个信息爆炸的时代,广告无孔不入,无处不在。... 更多>
山楂树瓶里怎么养水仙花
水仙花是一种美丽的花卉,它的花朵娇艳欲滴,散发出淡... 更多>
限时优惠!Apple MacBook Air 13.6英寸笔记本促销
力度空前。原售价8199元的AppleMacBookAir13 6英寸笔... 更多>
郑州发布地质灾害气象风险预警
郑报融媒把握正确舆论导向,关注百姓生活,侧重报道社... 更多>
新天地上市后首份半年报:净利润同比增长30.91%
新天地表示,随着募投项目的实施,产能未来有望进一步... 更多>
机器人:目前公司核应急机器人处于示范应用阶段
机器人在投资者互动平台表示,公司核应急机器人可应用... 更多>
法国拟加强管控关键原材料产业,外媒炒作“瞄准中国”
不过接受《环球时报》记者采访的专家认为,法国的相关... 更多>