std :: uniform_real_distribution包含范围

时间:2013-04-25 20:49:35

标签: c++ random c++11

C ++ 11 std :: uniform_real_distribution(-1,1)给出了[-1,1]范围内的数字。

如何在[-1,1]范围内获得统一的实际分布?

实际上它可能并不重要,但从逻辑上讲,我正在尝试选择包含范围内的值。

2 个答案:

答案 0 :(得分:18)

如果从查看整数开始,这会更容易思考。如果你传递[-1,1),你会得到-1, 0。由于您想要包含1,因此您将传递[-1,(1 + 1))或[-1,2]。现在你得到-1, 0, 1

你想做同样的事情,但是有双打:

借鉴this answer

#include <cfloat> // DBL_MAX
#include <cmath> // std::nextafter
#include <random>
#include <iostream>

int main()
{
  const double start = -1.0;
  const double stop = 1.0;

  std::random_device rd;
  std::mt19937 gen(rd());

  // Note: uniform_real_distribution does [start, stop),
  //   but we want to do [start, stop].
  //   Pass the next largest value instead.
  std::uniform_real_distribution<> dis(start, std::nextafter(stop, DBL_MAX));

  for (auto i = 0; i < 100; ++i)
  {
    std::cout << dis(gen) << "\n";
  }
  std::cout << std::endl;
}

(参见代码运行here

也就是说,在你想要的之后找到下一个最大的double值,然后将其作为结束值传递。

答案 1 :(得分:2)

不幸的是,浮点分布的实际实现不允许您如此精确。例如,uniform_real_distribution<float>应该在给定的半范围内产生值,但由于舍入问题,它实际上可能产生包含范围内的值。

Here's generate_cannonical问题的一个示例,其他real_distributions也会出现类似的问题。