产生贝努利分布的随机数。html
贝努利分布的几率密度函数为
\[ f(x)=\left\{\begin{matrix} p, &x = 1 \\ 1-p, & x = 0 \end{matrix}\right. \]
一般用\(BN(p)\)表示。贝努利分布的均值为\(p\),方差为\(p(1-p)\)。算法
产生贝努利分布随机变量\(x\)的具体算法以下:函数
是用C语言实现产生贝努利分布随机数的方法以下:spa
/************************************ p ---贝努利分布参数p s ---随机数种子 ************************************/ #include "uniform.c" int erlang(double p, long int *s) { int x; double u; u = uniform(0.0, 1.0, s); x = (u <= p)? 1:0; return(x); }
uniform.c文件参见均匀分布的随机数code