伪随机数生成算法

写在前面

伪随机数生成算法在计算机科学领域应用广泛,比如枪击游戏里子弹命中扰动、数据科学里对样本进行随机采样、密码设计、仿真领域等等,背后都会用到伪随机数生成算法。

骰子

说随机,那什么是随机呢?随机意味着不可预测,没有任何规律。谈随机数,一定是在序列当中,单拿出一个数谈随机是没有意义的。给一个数字序列,如果能在其中发现规律可以预测或以一定概率(大于“猜”的概率)预测接下来的数,那么这个序列就不是随机的。

在20世纪早期科学工作中就开始需要使用随机数,为了获取随机数,研究人员通过物理方式采集了成千上万的随机数,并发布给他人使用,比如RAND公司在1955年发布的《A Million Random Digits with 100,000 Normal Deviates》(百万乱数表)——亚马逊美国现在还有卖~。但是,通过物理方式采集“真”随机数并不高效,实时获取需要附加额外的随机数发生装置,而且获取速度缓慢、序列不可复现,如果将采集到随机数全保存下来则需要占用额外的存储空间,而且数量终究是有限的,于是大家开始寻求生成“伪”随机数的数学方法。伪随机数,顾名思义,即看起来是随机的但实际上不是,在不知其背后生成方式的情况下,生成的序列看上去毫无规律可言。

本文源自个人兴趣通过查阅参考文献整理所得,再加上个人的理解,大部分图片来自WIKI。

统计学检验

如何判断一个序列是否够随机呢?伪随机数生成算法多种多样,总要分出个孰好孰差,如何对各自的随机性进行定量评估呢?主要有两类方式,其出发点都是试图定量评估序列中是否隐含某种规律或模式

  • 实证检验。给定一个随机序列而不告知其背后的生成方式,尝试对观测到的分布进行拟合,以预测后面的序列,或者查看其中是否具有某些统计规律,比如查看分布是否均匀、连续值的相关性、某些数出现位置的间隔分布是否有规律等等。具体有 χ 2 \chi ^2 χ2检验KS检验、Frequency test、Serial test等等。

  • 理论检验。直接分析生成器的理论性质(已知生成方式),生成器通常需要配置一些参数,不同的参数会影响生成序列的质量,比如考察不同参数对随机序列周期的影响。

可在下一小节对理论检验一窥一二,但本文的重点不在于此,就不详细展开了,详细内容可见参考资料。

线性同余法

lin­ear con­gru­en­tial generator(LCG)线性同余法是最早最知名的伪随机数生成算法之一,曾被广泛应用,后逐渐被更优秀的算法替代,其通过如下递推关系定义:

X n + 1 = ( a X n + c )   m o d   m X_{n+1} = (aX_n + c)\ mod \ m Xn+1=(aXn+c) mod m

其中, X X X为伪随机序列,

  • m m m m > 0 m > 0 m>0,模数,显然其也是生成序列的最大周期
  • a a a 0 &lt; a &lt; m 0 &lt; a &lt; m 0<a<m,乘数
  • c c c 0 ≤ c &lt; m 0 \leq c &lt; m 0c<m,增量
  • X 0 X_0 X0 0 ≤ X 0 &lt; m 0 \leq X_0 &lt; m 0X0
  • 14
    点赞
  • 52
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
学术研究中,随机数生成算法是用于产生看起来像随机数的数字流的算法。由于有限状态机不能产生真正的随机数,现有的随机数生成算法通过生成具有较大重复周期的数字流来降低重复率,从而使得生成的数字看起来像是真正的随机数。这些算法通常以较小的数字作为输入,然后生成巨大甚至无限的数字流。这个输入被称为种子。随机数生成器被认为是安全的,如果没有计算机程序能够区分它们生成的数字流和真正的随机流。为了确定随机数生成器的质量,统计检验被用来检验生成的数字流是否满足随机性的要求。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* [关于真随机数随机数](https://blog.csdn.net/eyeder/article/details/48732185)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *2* *3* [随机函数 密码学_密码学与理论2:什么是随机](https://blog.csdn.net/danpu0978/article/details/107276441)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值