3.1由已知分布的随机抽样分析.ppt

蒙特卡罗方法及应用 第三章 由已知分布的随机抽样 3.1 随机抽样及其特点 3.2 直接抽样方法 3.2.1 离散型分布的直接抽样方法 3.2.1 离散型分布的直接抽样方法 3.2.1 离散型分布的直接抽样方法 3.2.1 几个典型的例子 3.2.1 几个典型的例子 3.2.1 几个典型的例子 3.2.1 几个典型的例子 3.2.2 连续型分布的直接抽样方法 3.2.2 几个典型的例子 3.2.2 几个典型的例子 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 3.3 挑选抽样方法 例7. 圆内均匀分布抽样 令圆半径为R0，点到圆心的距离为r，则r的分布密度函数为 分布函数为 容易知道，该分布的直接抽样方法是 由于开方运算在计算机上很费时间，该方法不是好方法。下面使用挑选抽样方法：取 则抽样框图为 显然，没有必要舍弃ξ1＞ξ2的情况，此时，只需取 就可以了，即 3.4 替换抽样方法 例7. 散射方位角余弦分布的抽样 直接抽样方法为： 令φ=2θ，则θ在［0,π］上均匀分布，作变换 其中0≤ρ≤1，0≤ρ≤π，则 (x,y) 表示上半个单位圆内的点。如果 (x,y) 在上半个单位圆内均匀分布，则θ在［0,π］上均匀分布，由于 因此抽样sinφ和cosφ的问题就变成在上半个单位圆内均匀抽样 (x,y) 的问题。 为获得上半个单位圆内 的均匀点，采用挑选法，在 上半个单位圆的外切矩形内 均匀投点（如图）。 舍弃圆外的点，余下的就是所要求的点。 抽样方法为： * * 随机数序列是由单位均匀分布的总体中抽取的简单子样，属于一种特殊的由已知分布的随机抽样问题。 由己知分布的随机抽样指的是由己知分布的总体中抽取简单子样。 本章所叙述的由任意已知分布中抽取简单子样，是在假设随机数为已知量的前提下，使用严格的数学方法，从已知分布中随机抽样产生。 只要所采用的随机数满足独立性和均匀性要求，那么由上述方法产生的简单子样严格服从具有相同的总体分布并相互独立。 对于任意给定的累积分布函数F(x)，要从中抽取随机变量序列X1,X2,…,Xn ，则可以产生一[0,1]间均匀分布的随机数序列， ?1, ?2,…, ?n ，按下式来确定Xi值 也就是说:当抽取?j后，判断所要取得各点值中，所对应的分布函数值大于?j，则其最小值/下界为所取值。 可以证明随机变量序列X1,X2,…,XN具有相同分布。 直接抽样法适应于离散型的随机变量 设离散型随机变量X的可能取值为x1, x2, …, xN, 其概率为 累积分布函数： 0 x1 xN-1 xN x2 xk-1 xk 1 F(x) ? pk 0 x1 xN-1 xN p1 p2 pN x2 pk xk-1 xk 方法： 计算yk = yk-1 + pk，k = 2,3,…,N, y1 = p1 产生在[0,1]区间上均匀分布的随机数? ； 求满足yk-1?≤yk 的k值; 随机变量的第k个取值即为欲抽取的值。 0 x1 xN-1 xN x2 xk-1 xk 1 F(x)=y ? pk 0 x1 xN-1 xN p1 p2 pN x2 pk xk-1 xk 证明： 0 x1 xN-1 xN x2 xk-1 xk 1 F(x) ? pk 0 x1 xN-1 xN p1 p2 pN x2 pk xk-1 xk 即:所产生的随机数的概率密度函数/pdf为pk p3=0.2 B3+c3 p2=0.3 B2+c2 p1=0.5 B1+c1 X? 例1、某核素衰变末态的随机抽样 设某核素X有三种衰变方式，其分支比如下 随机选取每次衰变的衰变方式（衰变道）?直接抽样法 ? = U[0,1] 例2、中子与核的反应类型确定 假设中子与核的反应类型有如下几种:弹性散射，非弹性散射，辐射俘获，裂变，相应的反应截面分别为σel ,σin,σf ,σa。则发生每一种反应类型的概率依次为 ： 其中反应总截面σt＝ σel＋ σin＋ σf＋ σa。 反应类型的确定方法为：产生一个随机数? 例3、碰撞目标核种类的确定 中子或光子在介质中发生碰撞时，如介质是由多种元素组成，需要确定到底与那种元素发生碰撞。假定入射粒子与介质中每种核的发生反应的宏观总截面分别为Σ1, Σ2, … , Σn，则中子或光子