贝努利分布的随机数

一、功能

产生贝努利分布的随机数。

二、方法简介

贝努利分布的概率密度函数为
\[ f(x)=\left\{\begin{matrix} p, &x = 1 \\ 1-p, & x = 0 \end{matrix}\right. \]
通常用\(BN(p)\)表示。贝努利分布的均值为\(p\)，方差为\(p(1-p)\)。

产生贝努利分布随机变量\(x\)的具体算法如下：

1. 产生均匀分布的随机数\(u\)，即\(u \sim U(0,1)\)；
2. 如果\(x\leqslant p\)，那么\(x=1\)；否则\(x=0\)。

三、使用说明

是用C语言实现产生贝努利分布随机数的方法如下：

/************************************
    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文件参见均匀分布的随机数