一.功能
生成jqdl分布的随机数。
二、方法介绍
jqdl分布的概率密度函数如下
() )是
f(x )=(FRAC ) x ) (sigma^{2} ) e ^ {-x ^ {2}/2 (sigma ^ {2} ) ) x0
() ) ) )。
JDL分布的平均值为((sigma ) sqrt ) Frac{pi}{2} ) ),方差为() ) left )2-frac { pi } {2} (right ) )
首先,用逆变换法生成参数 (贝塔=2) )的指数分布的随机变量((y ) ),其概率密度函数为() ) f ) y ) ) frac {1} {2} e ^ {- frac { y } 具体方法如下。
均匀分布随机数(u ),即) (u ) simu (0,1 ) );
计算((y=-2(ln ) u ) );
(x=) sigma ) sqrt{y} )。
三.使用说明
是一种用c语言实现生成jqdl分布随机数的方法。
/**************************
sigma ---jqdl分布的参数sigma
seed ---随机数种子
* * * * * * * * * * * * * * *
#include 'math.h '
#include 'uniform.c '
双阵列(双sigma,长int * s ) )。
{
u=统一(0.0、1.0、s );
x=-2.0*log(u;
x=sigma*sqrt(x;
返回(x;
}
统一. c文件请参阅均匀分布的随机数
标签: frac、sqrt、jqdl、分布、随机数、sigma
资料来源: https://www.cn blogs.com/Liam-Ji/p/11631644.html