如何使用标准::normal_distribution

How to use std::normal_distribution

本文关键字:distribution normal 标准 何使用      更新时间:2023-10-16

我是一个经验丰富的C程序员,偶尔会被迫使用一点C++。

我需要从具有各种均值和方差的正态分布生成随机数。如果我有一个 C 函数来执行此操作,称为 normal(float mean, float var)那么我可以编写以下代码:

int i;
float sample;
for(i = 0;i < 1000;i++)
{
  sample = normal(mean[i],variance[i]);
  do_something_with_this_value(sample);
}

请注意,i 的每个值都有不同的均值和方差。

C 不包含一个叫做 normal 的函数,但C++包含,实际上它叫做 std::normal_distribution。不幸的是,我的C++不足以理解文档中的语法。谁能告诉我如何使用 std::normal_distribution 实现我的 C 代码的功能。

std::normal_distribution不是函数而是模板化类

你可以像这样使用它:

#include <random>
int main(int, char**)
{
    // random device class instance, source of 'true' randomness for initializing random seed
    std::random_device rd; 
    // Mersenne twister PRNG, initialized with seed from previous random device instance
    std::mt19937 gen(rd()); 
    
    int i;
    float sample;
    for(i = 0; i < 1000; ++i)
    {
        // instance of class std::normal_distribution with specific mean and stddev
        std::normal_distribution<float> d(mean[i], stddev[i]); 
        // get random number with normal distribution using gen as random source
        sample = d(gen); 
        // profit
        do_something_with_this_value(sample); 
    }
    return 0;
}