C++生成随机数:高斯/正态分布(gaussian/normal distribution)随机数函数
  admin 发表于 2015-03-16 17:27  分类:技术资料  4,703 次阅读  0条评论

高斯分布也称为正态分布(normal distribution)

常用的成熟的生成高斯分布随机数序列的方法由Marsaglia和Bray在1964年提出,C++版本如下:

#include <stdlib.h>

#include <math.h>

double gaussrand()

{

   static double V1, V2, S;

   static int phase = 0;

   double X;

   

   if ( phase == 0 ) {

       do {

           double U1 = (double)rand() / RAND_MAX;

           double U2 = (double)rand() / RAND_MAX;

           

           V1 = 2 * U1 – 1;

           V2 = 2 * U2 – 1;

           S = V1 * V1 + V2 * V2;

       } while(S >= 1 || S == 0);

       

       X = V1 * sqrt(-2 * log(S) / S);

   } else

       X = V2 * sqrt(-2 * log(S) / S);

       

   phase = 1 – phase;

   return X;

}

  这样生成的高斯分布随机数序列的期望为0.0,方差为1.0。若指定期望为E,方差为V,则只需增加:

X = X * V + E;

文章来自:http://www.cnblogs.com/yeahgis/archive/2012/07/13/2590485.html

本文固定链接: http://www.youduoshao.com/2015-03-16/201503162397.html

填写您的称呼和邮箱即可发布评论

快捷键:Ctrl+Enter