随机抽样培训资料课件.pptVIP

这种方法需要知道P1的值(P2=1－P1），这对有些分布是很困难的。下面的方法可以不用计算P1 ： 对于任意小于1的正数P1 ，令P2=1－P1 ； 则采用复合挑选抽样方法，有： 第六十二页，共一百一十页。 当取 时，抽样效率最高 这时，乘加抽样方法为： ≤ 第六十三页，共一百一十页。 由于 可知第一种方法比第二种方法的抽样效率高。 第六十四页，共一百一十页。 例19. 光子散射后能量分布的抽样 令光子散射前后的能量分别为 和 （以 m0c2 为单位，m0为电子静止质量，c 为光速）， ， 则 x 的分布密度函数为： 该分布即为光子散射能量分布，它是由著名的Klin－Nishina 公式确定的。其中 K(α) 为归一因子： 第六十五页，共一百一十页。 把光子散射能量分布改写成如下形式： 在［1, 1+2α］上定义如下函数: 第六十六页，共一百一十页。 则有 使用乘加抽样方法： 第六十七页，共一百一十页。 光子散射能量分布的抽样方法为： 该方法的抽样效率为： ＞ ＞ ＞ ≤ ≤ ≤ 第六十八页，共一百一十页。 乘减分布的形式为： 其中H1(x) 、H2(x)为非负函数，f1(x)、f2(x) 为任意分布密度函数。 与减抽样方法类似，乘减分布的抽样方法也分为两种。 乘减抽样方法 第六十九页，共一百一十页。 （1）将 f (x) 表示为 令H1(x)的上界为M1， 的下界为m，使用乘抽 样方法得到如下乘减抽样方法： 第七十页，共一百一十页。 （2）将 f (x) 表示为 令H2(x)的上界为M2，使用乘抽样方法，得到另一种乘减抽样方法： 第七十一页，共一百一十页。 例20. 裂变中子谱分布抽样 裂变中子谱分布的一般形式为： 其中A，B，C，Emin，Emax 均为与元素有关的量。令 其中λ为归一因子，γ为任意参数。 第七十二页，共一百一十页。 相应的 H1(E)，H2(E) 为： 于是裂变中子谱分布可以表示成乘减分布形式: 容易确定 H1(E) 的上界为： 为提高抽样效率，应取γ使得M1 达到最小，此时 第七十三页，共一百一十页。 取 m＝0，令 则裂变中子谱分布的抽样方法为: 抽样效率 ＞ ≤ 第七十四页，共一百一十页。 对称分布的一般形式为： 其中 f1(x) 为任意分布密度函数，满足偶函数对称条件，H(x) 为任意奇函数，即对任意x满足： 对称分布的抽样方法如下：取η=2ξ－1 对称抽样方法 ≤ 第七十五页，共一百一十页。 证明： 因为η=2ξ－1，η≤x 相当于ξ≤ ，因此 第七十六页，共一百一十页。 例21. 质心系各向同性散射角余弦分布抽样 在质心系各向同性散射的假设下，为得到实验室系散射角余弦，需首先抽样确定质心条散射角余弦： 再利用下面转换公式: 得到实验室系散射角余弦μL。其中A为碰撞核质量，θC、θL 分别为质心系和实验室系散射角。 第七十七页，共一百一十页。 为避免开方运算，可以使用对称分布抽样。 根据转换公式可得： 依照质心系散射各向同性的假定，可得到实验室系散射角余弦μL 的分布如下： 该密度函数中的第一项为偶函数，第二项为奇函数，因而是对称分布。其中 第七十八页，共一百一十页。 复合挑选抽样方法 考虑另一种形式的复合分布如下： 其中0≤H(x,y)≤M，f2(x/y)表示与参数y有关的条件分布密度函数，F1(y)表示分布函数。抽样方法如下： 第三十页，共一百一十页。 证明： 抽样效率为：E=1/M 第三十一页，共一百一十页。 为了实现某个复杂的随机变量 y 的抽样，将其表示成若干个简单的随机变量 x1，x2，…，xn 的函数 得到 x1，x2，…，xn 的抽样后，即可确定 y 的抽样，这种方法叫作替换法抽样。即 替换抽样方法 第三十二页，共一百一十页。 例11. 散射方位角余弦分布的抽样 散射方位角φ在［0,2π］上均匀分布，则其正弦和余弦sinφ和cosφ服从如下分布： 直接抽样方法为： 第三十三页，共一