我将如何在具有高斯分布的两个值之间产生随机数
How would I produce random numbers between two values with a Gaussian distribution
刚接触 Python,正在做书中的一些练习。我需要生成 200 到 600 之间的 800 个随机数,具有高斯分布。我已经走到这一步了:
x = pylab.zeros(800,float)
for x in range (0,800):
y = random.gauss(550,30)
但是,只要所有800都符合高斯分布,这不是会产生任何数字吗?我需要 200 到 600 之间的值。
高斯分布没有界限,但您可以使您不太可能在您的范围之外进行采样。例如,您可以对平均值为 400 且标准差为 200/3 的数字进行采样,这意味着超出范围 [200, 600] 将超出 3 个标准差。
mean = 400
stdev = 200/3 # 99.73% chance the sample will fall in your desired range
values = [random.gauss(mean, stdev) for _ in range(800)]
如果你想要一个有界的伪高斯分布,你可以这样做
values = []
while len(values) < 800:
sample = random.gauss(mean, stdev)
if sample >= 200 and sample < 600:
values.append(sample)
因此,如果您对超出所需范围的值进行采样,则将其丢弃并重新采样。
我认为最好的方法是生成高斯分布,然后 re-scale 所有值,以便 max=600
和 min=200
.
刚接触 Python,正在做书中的一些练习。我需要生成 200 到 600 之间的 800 个随机数,具有高斯分布。我已经走到这一步了:
x = pylab.zeros(800,float)
for x in range (0,800):
y = random.gauss(550,30)
但是,只要所有800都符合高斯分布,这不是会产生任何数字吗?我需要 200 到 600 之间的值。
高斯分布没有界限,但您可以使您不太可能在您的范围之外进行采样。例如,您可以对平均值为 400 且标准差为 200/3 的数字进行采样,这意味着超出范围 [200, 600] 将超出 3 个标准差。
mean = 400
stdev = 200/3 # 99.73% chance the sample will fall in your desired range
values = [random.gauss(mean, stdev) for _ in range(800)]
如果你想要一个有界的伪高斯分布,你可以这样做
values = []
while len(values) < 800:
sample = random.gauss(mean, stdev)
if sample >= 200 and sample < 600:
values.append(sample)
因此,如果您对超出所需范围的值进行采样,则将其丢弃并重新采样。
我认为最好的方法是生成高斯分布,然后 re-scale 所有值,以便 max=600
和 min=200
.