<iframe align="top" marginwidth="0" marginheight="0" src="http://www.zealware.com/46860.html" frameborder="0" width="468" scrolling="no" height="60"></iframe>
主要参考《Numerical Recipes in C++ 2/e》p.292~p.294 和《Simulation Modeling and Analysis 3/e》p.465~p.466。
Box 和 Muller 在 1958 年给出了由均匀分布的随机变量生成正态分布的随机变量的算法。设 U1, U2 是区间 (0, 1) 上均匀分布的随机变量,且相互独立。令
X1 = sqrt(-2*log(U1)) * cos(2*PI*U2);
X2 = sqrt(-2*log(U1)) * sin(2*PI*U2);
那么 X1, X2 服从 N(0,1) 分布,且相互独立。等于说我们用两个独立的 U(0,1) 随机数得到了两个独立的 N(0,1)随机数。
Marsaglia 和 Bray 在 1964 年提出了一种改进算法,避免使用三角函数。以下的实现代码用的就是这种改进算法。
//
// Gaussian Random Number Generator class
// ref. ``Numerical Recipes in C++ 2/e'', p.293 ~ p.294
//
public class GaussianRNG
{
int iset;
double gset;
Random r1, r2;
public GaussianRNG()
{
r1 = new Random(unchecked((int)DateTime.Now.Ticks));
r2 = new Random(~unchecked((int)DateTime.Now.Ticks));
iset = 0;
}
public double Next()
{
double fac, rsq, v1, v2;
if (iset == 0) {
do {
v1 = 2.0 * r1.NextDouble() - 1.0;
v2 = 2.0 * r2.NextDouble() - 1.0;
rsq = v1*v1 + v2*v2;
} while (rsq >= 1.0 || rsq == 0.0);
fac = Math.Sqrt(-2.0*Math.Log(rsq)/rsq);
gset = v1*fac;
iset = 1;
return v2*fac;
} else {
iset = 0;
return gset;
}
}
}
Trackback: http://tb.blog.csdn.net/TrackBack.aspx?PostId=178666
分享到:
相关推荐
使用 C 语言产生正态分布的随机数。给出了全部的源代码。产生的随机数保存在文本文件中,并给出了 MATLAB 直方图验证的源代码。
根据平均值和标准差生成符合正态分布的随机数
实现正态分布的伪随机数发生器
用c语言编写的产生正态随机数的源代码: 算法一产生12个(0,1)平均分布的随机函数,用大数定理可以模拟出正态分布。 算法二用到了数学中的雅可比变换,直接生成正态分布,但此算法在计算很大规模的数时 会出现...
用labview实现的正态分布的随机数并生成对应的直方图(Achieved with labview normally distributed random numbers and generate the corresponding histogram)
30种随机数发生器,计算统计学课上作业,包含5种单独的随机数发生器和两两组合的25种随机数发生器共计30种随机数发生器的matlab代码。
最近上通信建模这门课时范平志老师布置的这个作业,我做好后传上来方便后人参考,包括了三个C语言程序,分别产生服从正态分布、瑞利分布、泊松分布的随机数。程序是是用的C语言编写,备有大量注释,浅显易懂,且全部...
用蒙特卡洛方法获取标准正态分布的伪随机数
用该程序可输出一个(0,1)区间内的服从高斯分布的随机数,精度为15位字长。
使用现代优化理论中的混合同余法,用C语言产生服从标准正态分布的随机数。
正态散布随机数生成算法.pdf
利用boxmuller算法生成符合正态分布的随机数。输入方差、期望、数目后运行即可。
实现正态分布随机数,可以控制平均值和标准差
主要介绍了C#利用Random得随机数求均值、方差、正态分布的方法,涉及C#数学运算及概率统计的相关技巧,需要的朋友可以参考下
用来产生指数,泊松,正态分布的随机数 有详细注释 这些函数的数学原理请参看浙大4版的概率统计教材
monte carlo随机数发生器 软件MATLAB 服从正态分布
国密标准 《GM/T 0105-2021 软件随机数发生器设计指南》 正式版
随机数发生器的C语言实现,包括平方取中法随机数发生器、积取中法随机数发生器、常数乘子随机数发生器、斐波那契随机数发生器。