C++生成随机数：高斯/正态分布（gaussian/normal distribution）

最新推荐文章于 2025-02-01 00:15:00 发布

yxpandjay

最新推荐文章于 2025-02-01 00:15:00 发布

阅读量3.6k

点赞数 2

CC 4.0 BY-SA版权

分类专栏： C++

本文链接：https://blog.youkuaiyun.com/yxpandjay/article/details/80607311

C++ 专栏收录该内容

21 篇文章

订阅专栏

本文介绍了一种由Marsaglia和Bray提出的生成高斯分布(正态分布)随机数序列的经典方法，并提供了C++实现代码。该方法能够生成期望为0.0，方差为1.0的随机数序列，通过简单的数学变换还可以得到任意期望和方差的高斯分布随机数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

高斯分布也称为正态分布（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;
}