在Julia中获取[-1,1]中的随机浮点值

Question

是

p = rand(-1.:eps():1., 100000)

在[-1,1]中获取随机Float值的好方法？

一个常见的建议似乎是2.*兰特（100000）-1 相反，但是（a）这永远不会返回1，因为rand的范围是[0，1），（b）这跳过了很多值：为了参数的缘故，假设eps（）==0.1，然后rand从（0.1，0.2，0.3，…，0.9）返回，在这个计算之后，您将从（-0.8，-0.6，-0.4，…，0.8）得到结果，因此结果在该范围内不再是均匀随机的。

（注：就性能而言，我在顶部的版本似乎比后一个慢4倍。）

在给定范围内获得统一随机浮点数的一般推荐方法是什么？




Answer 1

使用Distributions.jl包在（-1,1）之间创建一个统一分布，并使用rand对其进行采样。

julia> using Distributions

julia> rand(Uniform(-1, 1), 10000)
10000-element Vector{Float64}:
  0.2497721424626267
  ...
 -0.27818099962886844

如果您不需要向量，只需要一个标量数，您可以这样称呼它（感谢@DNF指出这一点）：

julia> rand(Uniform(-1,1))
-0.02748614119728021

您还可以对不同形状的矩阵/向量进行采样：

julia> rand(Uniform(-1, 1), 3, 3)
3×3 Matrix{Float64}:
 -0.290787  -0.521785    0.255328
  0.621928  -0.775802   -0.0569048
  0.987687   0.0298955  -0.100009

查看分发文档。这里是jl。