蒙特卡罗方法是一种通过随机变量的数字模拟和统计分析来求取数学物理、工程技术问题近似解的数值方法,基本步骤:随机变量的抽样试验。按基本随机变量(输入随机变量)的已知概率分布进行随机抽样(数字模拟)。样本反应求解。对每个抽取的样本,按问题的性质采用确定性的控制数学、物理方程求取样本反应。计算反应量的统计量估计。对所有样本反应,按所求解答的类型分别求取输出随机变量的均值、方差或概率分布。
当求解确定性问题时,首先,要根据所提出的问题构造一个简单、适用的概率模型,使问题的解对应于该模型中随机变量的某些数字特征(如概率、数学期望、方差等);然后,在高速运行的计算机上生成随机数,并对随机数进行统计分析试验;最后,利用试验所获结果求出统计特征的估计值作为问题的近似解。总结以上思想,可以得出利用蒙特卡罗方法求解确定性问题的基本步骤为:
(1)根据所要求解的实际问题来构造概型,并使概型的某些统计特征恰好相当于所要求的问题的解。
(2)根据所建立的概率模型,设计、使用一些加速收敛的方法,以求加速收敛并提高计算精度。
(3)给出在计算机上产生概型中各种不同分布随机变量的方法。
(4)统计处理模拟结果,给出问题的近似解并做解的精度估计。
蒙特卡罗方法虽然可以求解许多确定性工程技术问题,但其独到之处还应该在于求解随机性问题。用蒙特卡罗方法求解随机性问题时,一般首先,根据问题的物理性质建立随机模型;然后,再根据模型中各个随机变量的分布,在计算机上产生随机数,进行大量的统计试验,以取得所求问题的大量试验值;最后,根据这些试验结果求它的统计特征量,从而获得所求问题的解。由此可见,用蒙特卡罗方法求解随机问题的步骤与求解确定性问题的步骤基本一致。
总之,蒙特卡罗方法的理论基础是概率论中的大数定律。设在N次独立试验中,n为事件A出现的次数,而P(A)为事件A在每次试验中出现的概率,贝努利大数定律指出,对于任意ε>0,当N→∞时,事件A出现的频率的概率收敛于事件的概率。
地下水系统随机模拟与管理
当随机变量满足独立分布时,若随机变量序列ξ1,ξ2,…,ξN的分布相同,ξi具有有限的数学期望E(ξi)=a,i=1,2,…,N,则根据柯欠莫哥洛夫大数定律,对于任意的ε>0,当N→∞时,变量ξi将以概率1收敛于期望值a。
在蒙特卡罗方法中,采用简单抽样方法进行随机变量的数字模拟,因此其所抽取的子样为具有同分布性质的独立随机变量,当抽取的样本个数足够大时,样本均值将以概率1收敛于分布均值,而事件A出现的频率则以概率收敛于事件A出现的概率,这样就保证了蒙特卡罗方法的概率收敛性。
根据所求解问题性质的不同,其基本随机变量可能属于不同的概率分布,为了产生不同分布类型的随机变量的抽样值(随机数),一般需先产生一个在[0,1]上均匀分布的随机变量的抽样值,然后按照给定的概率分布类型将其转化为所需随机变量的抽样值。因此,均匀分布随机变量随机数的生成是蒙特卡罗方法实现的基础。利用数值法产生的均匀随机变量的抽样值称之为伪随机数,这是因为数值方法的基础是某一数学递推公式,按这类递推公式产生的抽样与[0,1]均匀分布中的抽样在统计性质上不可能完全相同。
数学递推公式的一般形式是:
式中:f(xn,xn-1,…,xn-k)——某一给定的函数形式。根据这一函数式,当给定一组初值,x0,x-1,…,x-k后,便可依次求出x1,x2,…,xm…最常用的(0,1)均匀分布随机数生成的递推公式有:
乘同余法。用以产生(0,1)均匀分布随机数的递推公式为:
式中:λ,M和x0——预先给定的常数。
式(2、4)的意义是指以M除以λxi-1后得到的余数记为xi。由于是余数。
如此所得的随机数序列r1,r2,…,ri为具有(0,1)均匀分布的随机数。
由式(2、4)不难看出,不同的xi最多只能有M个,相应地不同的随机数ri也最多只能有M个。所以当产生的随机数ri个数多于M个时,就会出现循环数,这样,便再不能看成是随机数。为了使所产生的随机数能经得住数理统计中的独立性和均匀性检验,需要合理选择随机数生成参数x0,λ及M。
灰枣为什么叫灰枣 灰枣为什么叫灰枣图
1、灰枣在成熟变红之前,通体发灰所以得名灰枣。这层灰就像是挂了一层霜,但其实和柿子表面的灰差不多,大部分都是糖分,不影响食用。虽然它看似个头很小,而且还很脏,却是最好的枣。灰枣也称为新郑灰枣,果实呈长倒卵形,果皮为橙红色(0)人阅读时间:2024-10-05为什么有的艾草不苦有的苦 为什么有的艾草没香气
1、有的艾草不苦有的苦是因为品种不同。甜艾草指的是艾草,是一种多年生草本植物,别名有冰台、香艾、蕲艾、艾蒿、艾绒等;而苦艾草指的是苦艾,学名中亚苦蒿,多年生丛生植物,虽然都是菊科蒿属的植物,但两者不是同一种东西。2、草主(0)人阅读时间:2024-10-05腊梅种子是巴豆吗(腊梅种子是巴豆吗能吃吗)
1、腊梅籽不是巴豆。2、本质不同:巴豆:大戟科巴豆属植物巴豆树的干燥成熟果实。腊梅籽:植物腊梅(蜡梅)的果实。3、形态特征不同:巴豆形态特征:蒴果椭圆形,长约2厘米,直径1.4-2厘米,疏生短星状毛或近无毛;种子椭圆形,(0)人阅读时间:2024-10-05职场正能量励志短句(职场正能量短句霸气)
1、人生如歌,何必去过份的计较与留恋那些过往的东西,何必把那些已成定局的怨恨背在肩上,让心灵去承受这种不必要的负载,让精神再去遭受这种无谓的折磨!学会忘却吧!2、束缚于零星的事实而不敢想象的人缺乏创造力。3、青年人的眼睛(0)人阅读时间:2024-10-05灰枣为什么叫灰枣 灰枣为什么叫灰枣图
1、灰枣在成熟变红之前,通体发灰所以得名灰枣。这层灰就像是挂了一层霜,但其实和柿子表面的灰差不多,大部分都是糖分,不影响食用。虽然它看似个头很小,而且还很脏,却是最好的枣。灰枣也称为新郑灰枣,果实呈长倒卵形,果皮为橙红色..2024-10-05为什么有的艾草不苦有的苦 为什么有的艾草没香气
1、有的艾草不苦有的苦是因为品种不同。甜艾草指的是艾草,是一种多年生草本植物,别名有冰台、香艾、蕲艾、艾蒿、艾绒等;而苦艾草指的是苦艾,学名中亚苦蒿,多年生丛生植物,虽然都是菊科蒿属的植物,但两者不是同一种东西。2、草主..2024-10-05腊梅种子是巴豆吗(腊梅种子是巴豆吗能吃吗)
1、腊梅籽不是巴豆。2、本质不同:巴豆:大戟科巴豆属植物巴豆树的干燥成熟果实。腊梅籽:植物腊梅(蜡梅)的果实。3、形态特征不同:巴豆形态特征:蒴果椭圆形,长约2厘米,直径1.4-2厘米,疏生短星状毛或近无毛;种子椭圆形,..2024-10-05职场正能量励志短句(职场正能量短句霸气)
1、人生如歌,何必去过份的计较与留恋那些过往的东西,何必把那些已成定局的怨恨背在肩上,让心灵去承受这种不必要的负载,让精神再去遭受这种无谓的折磨!学会忘却吧!2、束缚于零星的事实而不敢想象的人缺乏创造力。3、青年人的眼睛..2024-10-05