蒙特卡罗方法在实验核物理中的应用(蒙特卡罗方法)

导读 今天菲菲来为大家解答以上的问题。蒙特卡罗方法在实验核物理中的应用,蒙特卡罗方法相信很多小伙伴还不知道,现在让我们一起来看看吧!1、蒙...

今天菲菲来为大家解答以上的问题。蒙特卡罗方法在实验核物理中的应用,蒙特卡罗方法相信很多小伙伴还不知道,现在让我们一起来看看吧!

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 检验。

本文就为大家分享到这里,希望小伙伴们会喜欢。

免责声明:本文由用户上传,如有侵权请联系删除!

猜你喜欢