C++生成正态分布随机数代码

最新推荐文章于 2021-05-24 16:36:12 发布

原创最新推荐文章于 2021-05-24 16:36:12 发布 · 2.5k 阅读

12 ·

CC 4.0 BY-SA版权

文章标签：

#C++ #随机数 #正态分布 #高斯分布

数学课程实验专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一个使用C++实现的高斯噪声生成算法，该算法利用Box-Muller变换从均匀分布转换到高斯分布。通过调整mu和sigma参数，可以生成具有不同均值和标准差的高斯随机数。


#include <cstdlib>
#include <cmath>
#include <limits>
double generateGaussianNoise(double mu, double sigma)
{
    const double epsilon = std::numeric_limits<double>::min();
    const double two_pi = 2.0*3.14159265358979323846;

    static double z0, z1;
    static bool generate;
    generate = !generate;

    if (!generate){
       return z1 * sigma + mu;
    }
    double u1, u2;
    do
     {
    	srand( (unsigned)time( NULL ) );
       u1 = rand() * (1.0 / RAND_MAX);
       u2 = rand() * (1.0 / RAND_MAX);
     }
    while ( u1 <= epsilon );

    z0 = sqrt(-2.0 * log(u1)) * cos(two_pi * u2);
    z1 = sqrt(-2.0 * log(u1)) * sin(two_pi * u2);
    return z0 * sigma + mu;
}