贝努利分布的随机数

一、功能

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

二、方法简介

贝努利分布的概率密度函数为

[f(x)=left{egin{matrix} p, &x = 1 \ 1-p, & x = 0 end{matrix} ight. ]

通常用(BN(p))表示。贝努利分布的均值为(p)，方差为(p(1-p))。

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

1. 产生均匀分布的随机数(u)，即(u sim U(0,1))；
2. 如果(xleqslant 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文件参见均匀分布的随机数