Python|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

  1. 从一个数字范围内产生随机数字
random.randrange(start, stop[, step])

从 range(start, stop, step) 返回一个随机选择的元素。
这相当于 choice(range(start, stop, step)),但实际上并没有构建一个 range 对象。
  1. 返回随机整数
random.randint(a, b)

相当于 randrange(a, b+1)
结果 N 满足: a <= N <= b
  1. 从非空序列 seq 返回一个随机元素。 如果 seq 为空,则引发 IndexError 异常。
random.choice(seq)

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

每一次运行会从列表中随机获得一个数字。
  1. 将序列 x 随机打乱
andom.shuffle(x[, random])

可选参数 random 是一个无参数函数,在 [0.0, 1.0) 中返回随机浮点数;默认情况下,这是函数 random()
import randomlst = [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)

  1. 返回从总体序列或集合中选择的唯一元素的 k 长度列表。 用于无重复的随机抽样。
random.sample(population, k, *, counts=None)

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

  1. 返回一个随机浮点数 N
random.uniform(a, b)

取决于等式 a + (b-a) * random() 中的浮点舍入,终点 b 可以包括或不包括在该范围内。
结果 N 满足:当 a <= ba <= N <= b ,当 b < ab <= N <= a
更多方法可查阅官方文档。
3. 不可预测之美
3.1 随机彩色点
Python|Python 随机(random)模块的不可预测之美
文章图片

解题思路: 可结合 turtle 模块绘制,随机小海龟出现的位置就可以了
import random import turtlecolors = ["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 即可得到 π 的近似值。
Python|Python 随机(random)模块的不可预测之美
文章图片

求解思路
  • 利用随机函数产生横坐标的值 x 和纵坐标的值 y(这两个值都应在0~1)
  • 判断由这两个随机数构成的点是否位于1/4圆的区域内(阴影部分),若该点位于阴影区域内则进行计数。
  • 不断产生新的点,由于随机函数生成的点坐标有一定的均匀性,当生成的点足够多时,就可得到阴影内和阴影外点的近似均匀分布。
  • 最后用在阴影内的点的数量除以总的点数,即可得到近似的阴影面积,也就得到了一个的1/4的近似值。
import randomi, 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 += 1print("PI=%f\n", 4 * s / n)

输出结果:
输入点的数量:9000000 PI= 3.141477777777778

输入的点数量越多,得到的 PI 的近似值就会越精确。
4 . 总结 【Python|Python 随机(random)模块的不可预测之美】随机数可以完美地模拟真实世界里的各种概率或随机事件。python 的随机数生成除了可以使用 random 模块外,还可以使用 numpy 库中所提供的方法。

    推荐阅读