我注意到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
如何获得范围而不是范围中的最小数字?
也许更优雅一点,因为你真的想从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()
的第一个参数)
您需要一个步骤参数:
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
你可以做到
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