我们利用具有复共轭的实三次代数整数上伯努利映射的混沌真轨道,开发了一个伪随机比特生成器。虽然该生成器具有较高的计算成本,但它精确地模拟了能够生成理想随机二进制序列的伯努利映射。特别地,我们阐明了一种种子选择方法,该方法可以毫无偏见地选择初始点(即种子),并可以避免由其派生的伪随机二进制序列的后部分中的重叠。此外,为了评估生成器的内存使用情况,我们给出了真实轨道上点表示增长的上限。我们还报告了大量测试结果,表明生成的伪随机序列具有良好的统计特性,并且比目前最流行的生成器之一Mersenne Twister MT19937具有优势。

1
M。
第页。
维坦伊
,
Kolmogorov复杂性及其应用简介
,第2版。(
施普林格
,
纽约
,
1997
).
2
H。
杉田
,
蒙特卡罗方法、随机数和伪随机数
(
日本数学学会
,
东京
,
2011
).
三。
D.E.博士。
克努特
,
计算机编程的艺术
,第3版。(
出版商
,
雷丁,马萨诸塞州
,
1998
),第2卷,第3章。
4
D.H.博士。
莱默
, “
大规模计算单元中的数学方法
,“in
第二届大型数字计算机械研讨会论文集
(
哈佛大学出版社
,
1951
),第页。
141
146
.
5
南威尔士。
戈洛姆
,
移位寄存器序列
(
爱琴公园出版社
,
加利福尼亚州拉古纳山
,
1982
).
6
M。
松本
T。
西村
, “
Mersenne Twister:一种623维均匀伪随机数发生器
,”
ACM事务处理。模型。计算。模拟。
8
,
30
(
1998
).
7
第页。
比林斯利
,
遍历理论与信息
(
威利
,
纽约
,
1965
).
8
E.公司。
亚特兰蒂·杰克逊
,
非线性动力学的观点
(
剑桥大学出版社
,
剑桥
,
1991
),第1卷,第4章。
9
F、。
拉努
, “
离散平面面积保护映射的数值研究
,”
阿童木。天体物理学。
31
,
289
301
(
1974
).
10
下午。
粘合剂
R.V.公司。
延森
, “
用有限状态机模拟混沌行为
,”
物理学。版次A
34
,
4460
4463
(
1986
).
11
C、。
贝克
G.公司。
勒普斯托夫
, “
相空间离散化对动力系统长期行为的影响
,”
物理D
25
,
173
180
(
1987
).
12
C、。
格雷博吉
,
S.M.公司。
哈梅尔
,
J.A.公司。
约克
、和
T。
索尔
, “
混沌动力学中物理轨迹的阴影:遏制和细化
,”
物理学。修订稿。
65
,
1527
1530
(
1990
).
13
D.J.D.博士。
赚取
美国。
屈里曼
, “
哈密顿映射的精确数值研究:无舍入误差迭代
,”
物理D
56
,
1
22
(
1992
).
14
美国。
道森
,
C、。
格雷博吉
,
T。
索尔
、和
J.A.公司。
约克
, “
Lyapunov指数在零左右波动时阴影的障碍
,”
物理学。修订稿。
73
,
1927
1930
(
1994
).
15
第页。
钻石
,
第页。
克洛登
,
答:。
波克罗夫斯基
、和
答:。
弗拉迪米洛夫
, “
一类混沌动力系统和单吸引中心随机映射数值模拟中的坍塌效应
,”
物理D
86
,
559
571
(
1995
).
16
T。
索尔
,
C、。
格雷博吉
、和
J.A.公司。
约克
, “
数值混沌解保持有效多久?
,”
物理学。修订稿。
79
,
59
62
(
1997
).
17
G.公司。
J.A.公司。
约克
, “
一维映射中混沌的崩溃
,”
物理D
136
,
18
30
(
2000
).
18
答:。
斋藤
, “
修正Bernoulli映射的计算方面
,”
掠夺。西奥。物理学。供应商。
161
,
328
331
(
2006
).
19
出于类似的原因,也没有提出使用帐篷映射混沌轨道的伪随机数生成器0,1由提供M(M)T型(x个)=1负极|2x个负极1|或者面包师的转变[0,1)2由提供
M(M)b条(x个,)=2x个,2如果x个0,1/22x个负极1,+12如果x个1/2,1
尽管这些地图,连同M(M)B类,提供(字面上)混沌映射的教科书示例。
20
答:。
斋藤
答:。
山口
, “
利用Bernoulli映射的混沌真轨道生成伪随机数
,”
混乱
26
,
063122
(
2016
).
21
S.M.公司。
乌拉姆
J。
冯·诺依曼
, “
随机过程与确定性过程的结合
,”
牛市。美国数学。Soc公司。
53
,
1120
(
1947
).
22
T.Y.公司。
J.A.公司。
约克
, “
上的遍历地图[0,1]和非线性伪随机数生成器
,”
非线性分析。
2
,
473
481
(
1978
).
23
美国。
上好佳
H。
井上
, “
伪随机数发生器与混沌
,”
事务处理。IECE标准。
E65型
,
534
541
(
1982
).
24
美国。
冗长的
H。
Trotter公司
, “
一些代数数的连分式
,”
J.Reine Angew。数学。
255
,
112
134
(
1972
);附录,J.Reine Angew。数学。267, 219–220 (1974).
25
F、。
维瓦尔第
, “
不可约多项式上的动力学
,”
非线性
5
,
941
960
(
1992
).
26
J.H。
洛文斯坦
,
G.公司。
波贾斯帕拉
、和
F、。
维瓦尔第
, “
kicked-oscillator模型中的粘性轨道
,”
动态。系统。
20
,
413
451
(
2005
).
27
美国。
秋山
, “
beta展开的有限性和周期性-数论和动力学开放问题
,”
CIRM条例
1
,
9
(
2009
).
28
M。
古角
,
美国。
伊藤
,
答:。
斋藤
,
J。
塔穆拉
、和
美国。
靖国神社
, “
一种新的多维慢连分式算法和阶梯曲面
,”
实验数学。
23
,
390
410
(
2014
).
29
答:。
斋藤
美国。
伊藤
, “
利用三次非理性计算真实混沌轨道
,”
物理D
268
,
100
105
(
2014
).
30
答:。
斋藤
,
美国。
靖国神社
,
J。
塔穆拉
、和
美国。
伊藤
, “
用代数数模拟任意维分段线性和线性分式映射的真轨道
,”
混乱
25
,
063103
(
2015
).
31
E.公司。
见鬼去吧
,
代数数理论讲座
(
施普林格
,
纽约
,
1981
),第5章。
32
É.
博雷尔
, “
希弗尔苏尔(Sur les chiffres décimaux de)2et潜水员面临的概率问题
,”
C.R.学院。科学。巴黎
230
,
591
593
(
1950
).
33
B。
阿达姆切夫斯基
年。
比若
, “
关于代数数的复杂性I.整数基的展开
,”
安。数学。
165
,
547
565
(
2007
).
34
G.公司。
马沙里亚
,Marsaglia随机数CDROM,包括随机性测试的Diehard电池组(1995),可在https://web.archive.org/web/20160125103112/http://stat.fsu.edu/pub/diehard/。
35
答:。
鲁欣
,
J。
索托
,
J。
Nechvatal公司
,
M。
史密斯
,
E.公司。
巴克
,
美国。
,
M。
莱文森
,
M。
万格尔
,
D。
银行
,
答:。
赫克特
,
J。
德雷
、和
美国。
Vo公司
《密码应用随机和伪随机数生成器的统计测试套件》,NIST特别出版物800-22第1a版(2010年)。
36
第页。
欧库耶
R。
西马德
, “
TestU01:随机数生成器经验测试的C库
,”
ACM事务处理。数学。柔和。
33
,
22
(
2007
).
37
F、。
帕内东
,
第页。
L'Ecuyer公司
、和
M。
松本
, “
基于模2线性递归的改进长周期生成器
,”
ACM事务处理。数学。柔和。
32
,
1
16
(
2006
).
38
美国。
骚扰
, “
F类2-Mersenne Twister伪随机数生成器的线性关系
,”
数学。计算。模拟。
100
,
103
113
(
2014
).
39
H。
尼德雷特
, “
多项式的因式分解与有限域上的一些线性代数问题
,”
线性代数应用。
192
,
301
328
(
1993
).
40
H。
尼德雷特
, “
产生伪随机数的多重递归矩阵方法
,”
有限域应用。
1
,
30
(
1995
).
41
对于初始化,我们执行了与程序中采用的相同的程序,该程序位于http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/MT2002/CODES/mt19937ar。c(c).
您当前无权访问此内容。