在C ++中是否有使用三角分布生成随机数的函数?

时间:2019-04-15 00:42:22

标签: c++ simulation distribution

我正在寻找c ++中的函数(我给它一个函数(min,mode,max),然后它返回由三角分布生成的随机数)。如果有实现此功能的代码,那就太好了。

2 个答案:

答案 0 :(得分:1)

std::piecewise_linear_distribution可用于建模三角形分布。

这是一个基于链接的cppreference页面上的示例代码的示例,该示例生成一个三角形分布,该三角形分布生成一个0到30之间的数字,一个峰值为20:

#include <random>
#include <iostream>
#include <iomanip>
#include <array>
#include <map>

std::piecewise_linear_distribution<double> triangular_distribution(double min, double peak, double max)
{
    std::array<double, 3> i{min, peak, max};
    std::array<double, 3> w{0, 1, 0};
    return std::piecewise_linear_distribution<double>{i.begin(), i.end(), w.begin()};
}

int main() {
    std::random_device rd;
    // create a mersenne twister PRNG seeded from some implementation-defined random source
    std::mt19937 gen(rd());

    // create a triangular distribution with a minimum of 0, a peak at 20, and a maximum of 30
    auto dist = triangular_distribution(0, 20, 30);

    std::map<int, int> hist;

    // use our distribution to generate 10,000 random numbers
    // (truncated to integers for the sake of output; the generated numbers are actually real numbers)
    for (int i = 0; i < 10000; ++i) {
        double num = dist(gen);
        ++hist[num];
    }

    // print out a nice histogram of the numbers generated
    for(auto p : hist) {
        std::cout << std::setw(2) << std::setfill('0') << p.first << ' '
            << std::string(p.second/10,'*') << '\n';
    }
}

可能的输出:

00 **
01 *****
02 ******
03 ************
04 **************
05 ******************
06 **********************
07 *************************
08 **************************
09 *********************************
10 ************************************
11 **************************************
12 *************************************
13 ********************************************
14 **************************************************
15 **************************************************
16 *******************************************************
17 *******************************************************
18 ************************************************************
19 *****************************************************************
20 **************************************************************
21 *******************************************************
22 ************************************************
23 *******************************************
24 ***************************************
25 ******************************
26 **************************
27 ****************
28 ***********
29 ***

答案 1 :(得分:0)

您可以通过汇总两个均值U₁ U2 简单三角形分布 T 中进行采样。用C术语表示为:

 float U1 = (float)rand() / RAND_MAX;
 float U2 = (float)rand() / RAND_MAX;
 float T  = U1 + U2;

相应缩放:结果在0到2之间,且模式为1。