一、功能
产生贝努利-高斯分布的随机数。
二、方法简介
贝努利-高斯分布的随机变量(x)是贝努利分布的随机变量(y)与高斯分布的随机变量(z)的乘积,即(x=y*x)。因此,贝努利-高斯分布的随机数可视为:每当贝努利序列中有1出现时,打开高斯随机数发生器,并用其输出代替1。贝努利-高斯分布的均值为(pmu),方差为(p),其中(p)是贝努利分布的参数,(mu)是高斯分布的均值。
在地震勘探信号处理中,常用贝努利-高斯序列描述地下主要层状结构的反射作用。
产生贝努利-高斯分布随机变量(x)的具体算法如下:
产生贝努利分布的随机数(y),即(y sim BN(p));
产生高斯分布的随机数(z),即(z sim N(mu,sigma));
计算(x=y*x)。
三、使用说明
是用C语言实现产生贝努利-高斯分布随机数的方法如下:
/************************************
p ---贝努利分布参数p
mean ---高斯分布的均值mu
sigma ---高斯分布的均方差sigma
s ---随机数种子
************************************/
#include "gauss.c"
double bg(double p, double mean, double sigma, long int *s)
{
double u;
double x;
u = uniform(0.0, 1.0, s);
if(u < p)
x = gauss(mean, sigma, s);
else
x = 0.0
return(x);
}
gauss.c文件参见正态分布的随机数