Python 随机（random）模块的不可预测之美

1 . 概念

1.1 真、伪随机数

大部分的计算机语言都会提供 API 生成随机数，此类 API 称为随机数生成器。

计算机可以用随机数模拟现实世界中的各种随机概率问题，没有随机生成器的编程语言不是“好语言”。

什么是真随机数？

现实世界中的随机数：比如掷钱币、骰子、转轮、使用电子元件的噪音、核裂变等等。

计算机通过硬件技术摸拟现实世界中这种物理现象所生成的随机数，我们称其为真随机数。 这样的随机数生成器叫做物理性随机数生成器。生成真随机数对计算机的硬件技术要求较高。

真正随机数的特点：不可预测。

如在掷硬币时，你无法真正预测到下一次硬币的面向。

什么是伪随机数？

由算法摸拟生成的随机数称其为伪随机数。计算机编程语言中所生成的随机数基本上都是伪随机数。

伪随机数的特点：既然是由算法模拟的，虽然在一个较短的周期内是无法预测的，在一个较长的周期内的随机数具有可预测性。

1.2 随机数种子

生成伪随机数时，需要设置随机种子，种子作用就是在随机数的生成算法里注入一个动态变化量。

比如说使用系统的当前时间做随机种子，随机算法就可以在时间变化的基础上生成随机性更大的随机数。但是，如果不是在毫秒级别下生成随机数，同一时间点下所生成的大量随机数就有可能出现相等的情况。

选择种子时，可以考虑综合多维度的变化值进行运算。如在 UNIX 系统中，将系统时间、连入WIFI、甚至按下的键盘次数都量化为了seed。

参考指标越多，伪随机数就越接近真正的随机生成。

2. Python random 模块

random 模块实现了各种分布的伪随机数生成器。因为完全确定性，它不适用于所有目的，并且完全不适合加密目的。不应将此模块的伪随机生成器用于安全目的。 有关安全性或加密用途，可使用 Python 中的 secrets 模块。

使得之前需要导入 random 模块

import random

2.1 随机模块的方法

1. 初始化随机种子

random.seed(a=None, version=2)

• 如果 a 被省略或为 None ，则使用当前系统时间做随机种子。

• 如果操作系统提供随机源，则使用它们而不是系统时间。

• 如果 a 是 int 类型，则直接使用。

  当设置随机种子是一个常量，则每一次随机数是固定的。

  import random
  #设置随机种子是一个 int 常量
  random.seed(10)
  print(random.random())
  #设置随机种子是一个 int 常量
  random.seed(10)
  print(random.random())
  #设置随机种子是一个 int 常量
  random.seed(10)
  print(random.random())
  

  输出结果：

  0.5714025946899135
  0.5714025946899135
  0.5714025946899135
  

2. 从一个数字范围内产生随机数字

 random.randrange(start, stop[, step])

从 range(start, stop, step) 返回一个随机选择的元素。

这相当于 choice(range(start, stop, step))，但实际上并没有构建一个 range 对象。

3. 返回随机整数

   random.randint(a, b)

相当于 randrange(a, b+1)

结果 N 满足： a <= N <= b

4. 从非空序列 seq 返回一个随机元素。 如果 seq 为空，则引发 IndexError 异常。

random.choice(seq)

import random

lst = [5, 3, 90, 12, 4, 6]
r = random.choice(lst)
print(r)

每一次运行会从列表中随机获得一个数字。

5. 将序列 x 随机打乱

andom.shuffle(x[, random])

可选参数 random 是一个无参数函数，在 [0.0, 1.0) 中返回随机浮点数；默认情况下，这是函数 random()

import random

lst = [5.0, 3.0, 90.0, 12.0, 4.0, 6.0]
#使用 random.random 函数
random.shuffle(lst, random.random)
print(lst)
#输出结果
[3.0, 90.0, 6.0, 12.0, 5.0, 4.0]
#----------------------------------
def my_random():
    return float(random.randint(0, 1))

lst = [5.0, 3.0, 90.0, 12.0, 4.0, 6.0]
#使用用户自定义函数
random.shuffle(lst, my_random)
print(lst)

7. 返回从总体序列或集合中选择的唯一元素的 k 长度列表。 用于无重复的随机抽样。

random.sample(population, k, *, counts=None)

8. 返回 [0.0, 1.0) 范围内的下一个随机浮点数。

random.random()

9. 返回一个随机浮点数 N

random.uniform(a, b)

取决于等式 a + (b-a) * random() 中的浮点舍入，终点 b 可以包括或不包括在该范围内。

结果 N 满足：当 a <= b 时 a <= N <= b ，当 b < a 时 b <= N <= a 。

更多方法可查阅官方文档。

3. 不可预测之美

3.1 随机彩色点

解题思路： 可结合 turtle 模块绘制，随机小海龟出现的位置就可以了

import random
import turtle

colors = ["red", "blue", "green", "gray", "orange"]
for i in range(100):
    turtle.penup()
    x = random.randint(-300, 300)
    y = random.randint(-300, 300)

    turtle.goto(x, y)
    turtle.pendown()
    turtle.dot(20, colors[i % 5])

turtle.done()

3.2 求 π 的值

概率法又称为蒙特卡罗法，是一种非常重要的数值计算方法。

该方法是以概率和统计理论方法为基础的一种计算方法。将所求解的问题同一定的概率模型相联系，用计算机实现统计模拟或抽样，以获得问题的近似解。

假设有一个半径为 1 的圆，如图所示，则图中阴影部分（1/4圆）的面积就等于值的1/4。通过概率法计算出阴影部分的面积，也就得到了π 值的 1/4，将阴影部分面积乘以 4 即可得到 π 的近似值。

求解思路

• 利用随机函数产生横坐标的值 x 和纵坐标的值 y（这两个值都应在0~1）

• 判断由这两个随机数构成的点是否位于1/4圆的区域内（阴影部分），若该点位于阴影区域内则进行计数。

• 不断产生新的点，由于随机函数生成的点坐标有一定的均匀性，当生成的点足够多时，就可得到阴影内和阴影外点的近似均匀分布。

• 最后用在阴影内的点的数量除以总的点数，即可得到近似的阴影面积，也就得到了一个的1/4的近似值。

  import random
  
  i, n, s = 0, 0, 0
  x, y = 0.0, 0.0
  n = int(input("输入点的数量："))
  random.seed()
  for i in range(n):
      x = random.random()
      y = random.random()
      if (x * x + y * y) <= 1:
          s += 1
  
print("PI=%f\n", 4 * s / n)

输出结果：

输入点的数量：9000000
PI= 3.141477777777778

输入的点数量越多，得到的 PI 的近似值就会越精确。

4 . 总结

随机数可以完美地模拟真实世界里的各种概率或随机事件。python 的随机数生成除了可以使用 random 模块外，还可以使用 numpy 库中所提供的方法。