Question

我正在尝试使用Boost根据使用C ++的beta分布生成随机数。我在网上看到很多根据random.hpp中的分布生成随机数的例子（例如this book）。但是，我无法将它们翻译为使用beta.hpp中的beta发行版。

感谢。

Answer 1

您首先要从范围（0,1）中均匀地绘制一个随机数。给定任何分布，然后您可以将该数字插入分布的“分位数函数”，结果就像从分布中抽取随机值一样。来自here：

从具有不跳转的cdf的任意分布生成随机数的一般方法是使用cdf的反函数：G（y）= F ^ { - 1}（y）。如果u（1），...，u（n）是来自均匀（0,1）分布的随机数，则G（u（1）），...，G（u（n））是随机的使用cdf F（x）进行分布的样本。

那么我们如何获得beta分布的分位数函数？ beta.hpp的文档是here。你应该可以使用这样的东西：

#include <boost/math/distributions.hpp>
using namespace boost::math;

double alpha, beta, randFromUnif; 
//parameters and the random value on (0,1) you drew

beta_distribution<> dist(alpha, beta);
double randFromDist = quantile(dist, randFromUnif);

Answer 2

根据boost的随机数库演示 Random_demo.cpp和Generating integers with different probabilities

您应该做的是使用“ variate_generator”类绑定您的随机数生成器和分布。

一个例子可能看起来像

#include <iostream>
#include "boost/random.hpp"

int main(int argc, char *argv[])
{
    int seed = 2018;
    typedef boost::random::mt19937 RandomNumberGenerator;
    typedef boost::random::beta_distribution<> BetaDistribution;
    typedef boost::variate_generator<RandomNumberGenerator&, BetaDistribution> 
    Generator;
    RandomNumberGenerator Rng(seed);
    BetaDistribution distribution(2,5);
    Generator getRandomNumber(Rng,distribution);
    for (int idx = 0 ; idx < 1000 ; ++idx) 
    {
        std::cout << getRandomNumber() << std::endl;
    }
return 0;
}

但是，在最新的文档enter link description here中，boost似乎建议将生成器直接传递到分发对象。以下代码的结果是相同的。

#include <iostream>
#include "boost/random.hpp"

int main(int argc, char *argv[])
{
    int seed = 2018;
    typedef boost::random::mt19937 RandomNumberGenerator;
    typedef boost::random::beta_distribution<> BetaDistribution;
    RandomNumberGenerator Rng(seed);
    BetaDistribution distribution(2,5);
    for (int idx = 0 ; idx < 1000 ; ++idx) 
    {
        std::cout << distribution(Rng) << std::endl;
    }
return 0;
}

使用Boost生成基于Beta分布的随机数

2 个答案: