正态分布的随机数发生器 in C#

转载 2006年06月08日 17:04:00

来自:http://dev.csdn.net/author/Solstice/d99caee900e4400780d3177ac64402d2.html

Box 和 Muller 在 1958 年给出了由均匀分布的随机变量生成正态分布的随机变量的算法。设 U1, U2 是区间 (0, 1) 上均匀分布的随机变量,且相互独立。

 

主要参考《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;
      }
    }
  }

 

相关文章推荐

C#产生正态分布、泊松分布、指数分布、负指数分布随机数(原创)

http://blog.sina.com.cn/s/blog_76c31b8e0100qskf.html 在编程过程中,由于数据仿真模拟的需要,我们经常需要产生一些随机数,在C#中,产生一般...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

正态分布的随机数发生器 in C#

 主要参考《Numerical Recipes in C++ 2/e》p.292~p.294 和《Simulation Modeling and Analysis 3/e》p.465~p.466。Bo...
  • holym
  • holym
  • 2006年05月16日 13:30
  • 1423

手把手教你写C#窗体程序,画出泊松分布、正态分布曲线 函数图像

想要源码的童鞋请在下方留言~~以泊松分布,产生50个随机数,画函数图像为例进行详细描述: 1. 首先打开你的vs2010,建立项目,选择窗体程序,点击确定 将需要的控件从工具箱中拖到窗体,如图所示,...
  • huplion
  • huplion
  • 2017年04月09日 08:50
  • 948

C语言产生满足正态分布的随机数

C语言中可以通过rand函数生成满足均匀分布的随机数,但是生成满足正太分布的随机数就没有那么简单了,下面对常用的几种方法进行总结并用C++编程实现。 方法一:由均匀分布的随机数来产生   ...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

[C#] 查标准正态分布表

C#里面要计算正态分布是一件比较麻烦的事情,一般是通过查表来实现的。
  • wwwsq
  • wwwsq
  • 2015年11月05日 16:53
  • 3998

C#Random得随机数求均值、方差、正态分布

最近在做中小学试卷分析系统,其中数据的分析让自己很头疼,整个系统采用B/S架构。在分析试卷难度梯度的时候需要用到正态分布,自己做了一些,也查阅了一些资料,终于掌握了将一组数据分析检验,最后生成正态分布...

C#产生正态分布、泊松分布、指数分布、负指数分布随机数(原创)

http://blog.sina.com.cn/s/blog_76c31b8e0100qskf.html 在编程过程中,由于数据仿真模拟的需要,我们经常需要产生一些随机数,在C#中,产生一般...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:正态分布的随机数发生器 in C#
举报原因:
原因补充:

(最多只允许输入30个字)