Std::uniform_real_distribution包含范围

std::uniform_real_distribution inclusive range

本文关键字:包含 范围 distribution real Std uniform      更新时间:2023-10-16

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

如何得到[-1,1]范围内的均匀实分布?

实际上它可能无关紧要,但逻辑上我试图选择一个包含范围内的值

如果从整数开始,这更容易考虑。如果你通过[- 1,1]你会期望得到-1, 0。因为你想包含1,你会传递[-1,(1+1)),或[-1,2)。现在你得到-1, 0, 1

你想做同样的事情,但使用双精度:

借用这个答案:

#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;
}

(查看此处运行的代码)

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

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

下面是generate_cannonical的问题示例,其他real_发行版也出现类似的问题。