蒙特卡罗方法教学课件第三章由巳知分布的随机抽样.pptVIP

由于1－ξ1可用ξ1代替，该抽样方法可简化为： 对于ξ2＞ξ1的情况，可取Xf＝ξ1，因此 与β分布的推论相同。＞≤如下形式的分布称为乘分布：其中H(x)为非负函数，f1(x)为任意分布密度函数。 令M为H(x)的上界，乘抽样方法如下：抽样效率为：乘抽样方法≤＞例17.倒数分布抽样 倒数分布密度函数为： 其直接抽样方法为： 下面采用乘抽样方法，考虑如下分布族： 其中i=1，2，…，该分布的直接抽样方法为： 利用这一分布族，将倒数分布f(x)表示成: 其中， 乘法分布的抽样方法如下: 该分布的抽样效率为：＞≤例18.麦克斯韦(Maxwell)分布抽样 麦克斯韦分布密度函数的一般形式为： 使用乘抽样方法，令 该分布的直接抽样方法为： 此时 则麦克斯韦分布的抽样方法为: 该分布的抽样效率为：＞≤在实际问题中，经常会遇到如下形式的分布：其中Hn(x)为非负函数，fn(x)为任意分布密度函数，n=1，2，…。不失一般性，只考虑n=2的情况： 将f(x)改写成如下的加分布形式：乘加抽样方法其中乘加抽样方法为：该方法的抽样效率为：≤这种方法需要知道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相当于ξ≤，因此使用挑选抽样方法时，要注意以下两点：选取h(x)时要使得h(x)容易抽样且M的值要尽量小。因为M小能提高抽样效率。抽样效率是指在挑选抽样方法中进行挑选时被选中的概率。按此定义，该方法的抽样效率E为：所以，M越小，抽样效率越高。当f(x)在［0，1］上定义时，取h(x)=1，Xh=ξ，此时挑选