如何在Julia中生成一系列随机浮点数?

问题描述 投票:1回答:3

我注意到rand(x)是一个整数的x给了我一个随机浮点数组。我想知道如何在一定范围内生成一组随机浮点型变量。我尝试使用如下范围:

rand(.4:.6, 5, 5)

我得到:

 0.4  0.4  0.4  0.4  0.4
 0.4  0.4  0.4  0.4  0.4
 0.4  0.4  0.4  0.4  0.4
 0.4  0.4  0.4  0.4  0.4
 0.4  0.4  0.4  0.4  0.4

如何获得范围而不是范围中的最小数字?

random julia
3个回答
8
投票

也许更优雅一点,因为你真的想从Uniform发行版中采样,你可以使用Distribution包:

julia> using Distributions
julia> rand(Uniform(0.4,0.6),5,5)
5×5 Array{Float64,2}:
 0.547602  0.513855  0.414453  0.511282  0.550517
 0.575946  0.520085  0.564056  0.478139  0.48139
 0.409698  0.596125  0.477438  0.53572   0.445147
 0.567152  0.585673  0.53824   0.597792  0.594287
 0.549916  0.56659   0.502528  0.550121  0.554276

然后,同样的方法适用于来自其他众所周知或用户定义的分布的采样(只需将分布作为rand()的第一个参数)


3
投票

您需要一个步骤参数:

rand(.4:.1:.6, 5, 5)

.1将为您的范围提供一个步骤,这对于浮点数是必需的,而不是递增1所必需的。问题是无论隐式精度如何,它都将为1。如果您需要比以下更精确的增量:

rand(.4:.0001:.6, 5, 5)

这将为您提供类似于以下内容的结果:

 0.4587  0.557   0.586   0.4541  0.4686
 0.4545  0.4789  0.4921  0.4451  0.4212
 0.4373  0.5056  0.4229  0.5167  0.5504
 0.5494  0.4068  0.5316  0.4378  0.5495
 0.4368  0.4384  0.5265  0.5995  0.5231

0
投票

你可以做到

julia> map(x->0.4+x*(0.6-0.4),rand(5,5))
5×5 Array{Float64,2}:
 0.455445  0.475007  0.518734  0.463064  0.400925
 0.509436  0.527338  0.566976  0.482812  0.501817
 0.405967  0.563425  0.574607  0.502343  0.483075
 0.50317   0.482894  0.54584   0.594157  0.528844
 0.50418   0.515788  0.5554    0.580199  0.505396

一般规则是

julia> map(  x -> start + x * (stop - start), rand(5,5)  )

其中start为0.4,stop为0.6

您甚至可以通过这样的方式生成六面骰子,其中x范围从1到7,即1 <x <7,因为x恰好是1.0或7.0的概率为零

julia> map(x->Integer(floor(1+x*(7-1))),rand(5,5))
5×5 Array{Int64,2}:
 2  6  6  3  2
 3  1  3  1  6
 5  4  6  1  5
 3  6  5  5  3
 3  4  3  5  4

或者你可以使用

julia> rand(1:6,5,5)
5×5 Array{Int64,2}:
 3  6  3  5  5
 2  1  3  3  3
 1  5  4  1  5
 5  5  5  5  1
 3  2  1  5  6
© www.soinside.com 2019 - 2024. All rights reserved.