《计算物理(本科)》[第11章].pptVIP

. 取九个随机数，按式计算第一个中子的（X1，Y1 ，Z1） ，判断此坐标点是否在铀块内，若是在内部则Nin增加1，接着同样方法计算第二个中子； 重复上步N次； 计算 f= Nin / N； 调节常数M和s使 f=1，得到临界时的Mc和s的值。 （6）模拟计算： . 第十一章 蒙特卡罗方法 §11.1 蒙特卡罗方法概述 蒙特卡罗(Monte Carlo)方法：利用随机数统计地去计算和模拟给定的问题。 说明： 1、也称统计模拟法、随机抽样法或统计试验法。 2、Monte Carlo方法的命名：世界上著名的赌城摩洛哥的Monte Carlo。 3、方法：先构造一个与物理问题等价的随机过程，当完成大量的随机试验后，结果由多次事件的平均值给出。 . 5、求解确定性物理问题：如微分方程、定积分、线性方程等。同其它数值计算方法相比是速度慢。 6、求解复杂物理问题：如果物理学问题的严格算法不知道，或非常复杂，则蒙特卡罗方法有意想不到的成功。特别在分子运动学、输运现象、布朗运动、放射性衰变等问题中，由于本身有一定的统计规律性，这种方法很奏效。 7、计算机模拟：蒙特卡罗方法广泛应用于“计算机模拟”，在计算机上模拟真实过程。 4、随机数的抽样：现都用计算机程序来产生，一般不用物理方法抽样。计算机产生的是伪随机数。 . §11.2 蒙特卡罗方法的原理 用蒙特卡罗方法寻找某个未知量x时，利用计算机产生的随机变量ξ，得到的期望值E(ξ) x= E(ξ) 为了估算ξ的平均值，构造随机变量ξ的若干独立的实验数序列{ξi|i=1，2，3，……n}，得 . § 11.3 伪（赝）随机变量的抽样 实际上，大多数伪随机变量不满足[0，1]均匀分布，而是具有符合分布密度函数为f(x)的分布。 抽取符合f(x) 随机数的步骤如下： （2）从上面随机数总体中抽取一个简单子样，使它满足分布密度函数f(x) ξ ：ξ 1，ξ2，ξ3 …… （1）在[0，1]之间抽取均匀分布的伪随机数序列 γ： γ1，γ2，γ3 …… . 一、离散型分布随机变量的直接抽样法 随机变量ξ的概率表 其中xi为离散型随机变量ξ的跳跃点,Pi为相应概率。 变换方法： 第一步 构造一个矢量 . 第二步 产生一个[0，1]间随机数γ，找到区间 使γ恰好落在这个区间内。 第三步 ξP取对应的值xj，其表达式为 . 例1.γ光子与物质相互作用的抽样问题。 物理过程：γ光子与物质作用有康普顿效应、光电效应和电子对效应三种类型 。其中光电效应和电子对效应为光子吸收过程。 设三种过程的碰撞截面分别σs 、σe和σp，则总截面σT = σs +σe+σp 。 根据给定的[0，1]之间均匀分布的随机数γ ，求应产生那种效应？ . （3）产生[0,1]间均匀随机数γ： 若γ＜ σs /σT，则发生康普顿散射； 若σs/σT≤γ＜(σs+σe)/σT ， 则发生光电过程； 若γ≥(σs + σe ) /σT ，则发生电子对产生过程。 （2）构造一个矢量 解：（1）随机变量的概率表 . 二、连续型分布随机变量的直接抽样法 一个γ对应一个ξ。如果积分式能求出解析形式的反函数，则得到变换公式 ξF=F -1(γ) 11.12 f(x) x a ξ 1-γ γ b 对连续型随机变量ξ的值，可直接解连续分布方程 其中f(x)为在（a，b）区间的归一化分布密度函数。 . 例1.均匀分布：在区间（a，b）内找出符合均匀分布的随机变量 。 解法2：线性变换，令 解法1： . 例2．指数分布：指数分布的一般形式为 f(x)=λe-λx x ≥0其中λ0 ，求符合此分布的随机变量ξ 求反函数 ξ=-1/λln(1-F(ξ)) 其中0≤F ≤ 1，由上式得