二项分布的表达式就不赘述了
与泊松分布类似,关键是要退出二项分布的递推公式:
P{X=i+1} = ( (n-i)/(i+1) )*( p/(1-p) )*P{X=i}
随机变量的产生步骤:
1. 生成随机数U
2. c=p/(1-p) , i=0, pr=(1-p)^n , F=pr (pr的值由二项分布表达式令i=0得到的,即初始项)
3. if U<F 则X=i,停止搜索
4.pr=( c*(n-i)/(i+1) )*pr, F=F+pr, i=i+1
5.返回步骤3
平均搜索长度为1+np
当p很大时可以优化
优化: 求与其互补的二项分布(n, 1-p), 然后用n减去得到的随机变量(想想为什么,成功的次数与不成功的次数加起来是n)
此外,也可以像泊松分布一样从均值附近开始搜索而不是从头来,即从Int(np)附近开始