梅森旋转随机数生成实例

最新推荐文章于 2023-05-23 20:14:50 发布
最新推荐文章于 2023-05-23 20:14:50 发布
阅读量4.9k 收藏 10
点赞数 1
分类专栏： c/c++ 文章标签： c++ 梅森旋转算法 随机数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/chinabhlt/article/details/45535565
版权
梅森旋转算法是1997年由松本真和西村拓士开发的伪随机数生成器，以其在有限二进制字段上的矩阵线性递归特性，提供高质量的随机数。该算法以梅森素数为周期，并有两个常见变体：MT19937（32位）和MT19937-64（64位）。本文将介绍如何用C++实现梅森旋转算法。
摘要由CSDN通过智能技术生成

梅森旋转演算法（Mersenne twister）是一个伪随机数发生算法。由松本真和西村拓士[1]在1997年开发，基于有限二进制字段上的矩阵线性递归F_{2}。可以快速产生高质量的伪随机数， 修正了古典随机数发生算法的很多缺陷。

Mersenne Twister这个名字来自周期长度取自梅森素数的这样一个事实。这个算法通常使用两个相近的变体，不同之处在于使用了不同的梅森素数。一个更新的和更常用的是MT19937, 32位字长。 还有一个变种是64位版的MT19937-64。 对于一个k位的长度，Mersenne Twister会在[0,2^k-1]的区间之间生成离散型均匀分布的随机数。

优点：

最为广泛使用Mersenne Twister的一种变体是MT19937，可以产生32位整数序列。具有以下的优点:

  1. 周期非常长，达到219937−1。尽管如此长的周期并不必然意味着高质量的伪随机数，但短周期（

最低0.47元/天 解锁文章
jlaij
关注
专栏目录
梅森旋转产生随机数c语言实现,C++生成随机数的实现代码
weixin_31062993的博客
05-23 687
C风格C为随机数提供的工具是rand、srand和RAND_MAX，定义在中。srand为rand设置种子，如果不设置，相当于调用过srand(1)。rand产生伪随机数，其范围为0到RAND_MAX，RAND_MAX至少是32767，在MSVC和GCC中这个值都是32767。伪随机数看似随机，实则是有规律可循的，对于相同的种子值，rand产生的序列完全相同，也就是说无论你给srand一个什么数字...
梅森旋转算法
06-04
梅森旋转算法是R,Python,Ruby,IDL,Free Pascal,PHP,Maple,Matlab,GMP和GSL的默认伪随机数产生器。从C++11开始，C++也可以使用这种算法。在Boost C++,Glib和NAG数值库中，作为插件提供。 在SPSS中，梅森选旋转算法是两个PRNG中的一个：另一个是产生器仅仅为保证旧程序的兼容性，梅森旋转被描述为”更加可靠“。梅森旋转在SAS中同样是PRNG中的一个，另一个产生器是旧时的且已经被弃用。
梅森旋转法产生随机数
shuimu12345678的专栏
06-26 3154
直接贴代码了： #include #include typedef unsigned long uint32; #define N (624) // length of state vector #define M (397) #define K (0x9908B0DFU) // a magic cons
随机数：梅森旋转算法
weixin_46596227的博客
03-24 1936
两个日本学者研究提出的伪随机数算法，在这个算法当中用到了梅森素数，所以称为梅森旋转算法。 简单介绍一下梅森素数，梅森素数的意思是形如：2n−12^n-12n−1的素数。利用梅森素数的性质可以设计出周期长度为梅森素数长度的随机数周期。比如目前Python、C++11等语言当中用的随机数计算包都是用的这种算法。目前常用的版本周期是：219937−12^{19937}-1219937−1 梅森旋转算法的周期非常非常长，但是仍不是安全的随机数算法，仍然有可能会被黑客破解。只不过和其他算法相比，被破解的概率以及难度增
【转】梅森旋转（Mersenne twister）及其他随机数生成算法
WangYouJin321的博客
10-31 785
例如(32,7,5,3,2,1,0)是指本原多项式x^32+x^7+x^5+x^3+x^2+x+1，把它转化为最大周期LFSR就是在LFSR的第32，7，5，2，1位抽头。Mersenne Twister有以下优点：随机性好，在计算机上容易实现，占用内存较少(mt19937的C程式码执行仅需624个字的工作区域)，与其它已使用的伪随机数发生器相比，产生随机数的速度快、周期长，可达到2^19937-1，且具有623维均匀分布的性质，对于一般的应用来说，足够大了，序列关联比较小，能通过很多随机性测试。
Mersenne Twister（梅森旋转算法）
最新发布
斯坦福-MIT-Harvard-CMU
05-23 1266
Mersenne Twister引擎是C++标准库中最常用的随机数引擎之一，具有良好的统计性质和长周期，因此在众多应用中得到了广泛的应用。高质量的随机数：Mersenne Twister引擎产生的随机数序列在统计上具有良好的性质，如均匀性、独立性和周期性。快速和高效：Mersenne Twister引擎的算法设计使其具有良好的性能，生成随机数的速度较快。在上面的示例中，生成的随机数似乎都是相同的值，这是因为没有添加产生随机数的分布函数。它是一种高质量的随机数引擎，具有良好的统计性质和较长的周期。
实例008——产生随机数.zip
03-22
本实例008——产生随机数，旨在探讨如何在不同的编程语言中生成随机数，以及随机数生成器的基本原理和使用技巧。 首先，我们要理解随机数的定义：它是一个在一定范围内无法预测的数值。在计算机科学中，虽然数字是...
random_num_gen.zip_random_randomnum_random（num）_随机数 verilog_随机数生
09-24
随机数生成通常基于确定性的算法，如线性同余法、梅森旋转算法等。这些算法虽然本质上是确定的，但通过精心选择参数，可以使得产生的数列在统计上看起来像随机序列。在硬件实现中，随机数生成器通常使用混沌电路或者...
随机数生成与蒙特卡洛模拟：在数值计算中的应用
在计算机科学和数学领域，随机数生成和蒙特卡洛模拟是非常重要的技术。随机数生成指的是通过特定的算法和方法产生一系列看似无规律的数字。而蒙特卡洛模拟则是一种基于随机数生成的计算方法，通过大量的随机模拟来...
梅森旋转随机算法
破楼兰
01-27 1万+
梅森旋转随机算法，C++和Java代码。 C++： #ifndef _MersenneTwister_H_ #define _MersenneTwister_H_ #include #include #include typedef int32_t MS_INT; class MersenneTwister { public: void rseed(MS_INT seed){
伪随机数生成——梅森旋转（Mersenne Twister/MT）算法笔记
热门推荐
tick_tock97的博客
11-28 3万+
前言 最近在看吴军博士的《数学之美》一书，把很多之前没注意到，没用到，甚至不知道怎么用的数学知识和实际问题联系了起来，感觉打开了新世界的大门一样。这本书很多知识点还有技术都是点到为止，并没有深入，所谓师傅领进门，修行在个人吧。所以从本篇开始，博主将对数学之美一书中的一些提到的东西做个总结。不对的地方希望各位看官大佬们多多指正。 在本书第16章《信息指纹及其应用》一文中，介绍到了现在常用的
伪随机数生成算法-梅森旋转（Mersenne Twister/MT）算法介绍
liinux-Talk is cheap,show me the code.
10-21 1万+
今天主要是来研究梅森旋转算法，它是用来产生伪随机数的，实际上产生伪随机数的方法有很多种，比如线性同余法， 平方取中法等等。但是这些方法产生的随机数质量往往不是很高，而今天介绍的梅森旋转算法可以产生高质量的伪随 机数，并且效率高效，弥补了传统伪随机数生成器的不足。梅森旋转算法的最长周期取自一个梅森素数， 由此命名为梅森旋转算法。常见的两种为基于32位的MT19937-32和基于64位的MT19...
生成n套数位加减乘除_伪随机数生成：梅森旋转（Mersenne Twister/MT）算法
weixin_39616348的博客
12-03 250
(给算法爱好者加星标，修炼编程内功)作者：tik_tokc97https://blog.csdn.net/tick_tock97/article/details/78657851前言最近在看吴军博士的《数学之美》一书，把很多之前没注意到，没用到，甚至不知道怎么用的数学知识和实际问题联系了起来，感觉打开了新世界的大门一样。这本书很多知识点还有技术都是点到为止，并没有深入，所谓师傅领进门，修...
MT19937（梅森旋转算法）
shshss64的博客
10-14 9079
MT19937算法
JAVA梅森旋转随机算法_梅森旋转算法（Mersenne twister）生成伪随机数
weixin_35006433的博客
02-16 422
import java.util.Random;/*** MT19937的Java实现*/public class MTRandom extends Random {// Constants used in the original C implementationprivate final static int UPPER_MASK = 0x80000000;private final stat...
用C语言写一个生成梅森旋转随机数
03-31
梅森旋转算法是一种常见的伪随机数生成算法，其原理是通过一系列位运算和循环移位操作来生成随机数。下面是一个用C语言实现的梅森旋转随机数生成器： ```c #include <stdint.h> #define N 624 #define M 397 #define MATRIX_A 0x9908b0dfUL #define UPPER_MASK 0x80000000UL #define LOWER_MASK 0x7fffffffUL static uint32_t mt[N]; static int mti = N + 1; void init_genrand(uint32_t seed) { mt[0] = seed; for (mti = 1; mti < N; mti++) { mt[mti] = 1812433253UL * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti; } } void generate_numbers() { int i; uint32_t y; static uint32_t mag01[2] = {0x0UL, MATRIX_A}; for (i = 0; i < N; i++) { y = (mt[i] & UPPER_MASK) | (mt[(i+1) % N] & LOWER_MASK); mt[i] = mt[(i+M) % N] ^ (y >> 1) ^ mag01[y & 0x1UL]; } mti = 0; } uint32_t extract_number() { if (mti >= N) { generate_numbers(); } uint32_t y = mt[mti++]; y ^= (y >> 11); y ^= (y << 7) & 0x9d2c5680UL; y ^= (y << 15) & 0xefc60000UL; y ^= (y >> 18); return y; } ``` 这个实现中，我们使用了一个624个元素的数组 `mt` 来存储当前的状态，用 `mti` 来记录当前状态的位置。在 `init_genrand` 函数中，我们初始化了这个数组，以 `seed` 为种子生成初始状态。然后，在 `generate_numbers` 函数中，我们使用了梅森旋转算法来更新状态数组。最后，在 `extract_number` 函数中，我们根据当前状态生成一个随机数，并更新状态位置。 要使用这个随机数生成器，可以按照以下方式调用： ```c init_genrand(12345); for (int i = 0; i < 10; i++) { uint32_t rand_num = extract_number(); printf("%u\n", rand_num); } ``` 这个例子中，我们以 `12345` 为种子初始化了随机数生成器，并生成了10个随机数。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值