一様乱数とは? わかりやすく解説

Weblio 辞書 > 同じ種類の言葉 > 人文 > 高等数学 > 乱数 > 一様乱数の意味・解説 

一様乱数

読み方いちようらんすう
【英】:uniform random numbers

概要

すべての値の出現確率等し乱数のこと. 離散型の一様乱数の場合は,とりうる値の集合としてふつう\{0,1,2,\cdots,m-1\} \,(m-1 \,計算機表現できる最大自然数に近い数), あるいは, この中から等間隔抽出した数の集合想定する. 離散型の一様乱数をm \,割ったものを近似的に区間[0,1)上の連続一様分布にしたがう乱数見なして, 標準一様乱数という.

詳説

 確率変数実現値見なしうる数の列のことを乱数 (または乱数列)(random numbersまたはrandom number sequence)という. 確率変数の従う分布として一様分布想定する場合には, 対応する乱数のことを一様乱数 (uniform random numbers)という. 例えば, サイコロ振って出る目(数)の系列は, 典型的な一様乱数である. この系列には次の2つ性質がある. 1) 系列長ければ各数の相対出現頻度ほぼ等しい. 2) 次に出る数を確実に予測することは不可能である.

 大量乱数使用する実験では, サイコロ振って乱数作るのは現実的でないので, 簡単なアルゴリズム生成される乱数もどきの数(擬似乱数)で代用するのがふつうである. そして, これを単に乱数と呼ぶことが多い. この意味での乱数は, 上記の2) の性質持たないが, 1) の性質近似的に満たしているものと考えられている. このような乱数の生成法数多く提案されているが, 現在比較的よく使われているものを以下にあげる.

線形合同法

 1948年頃にレーマー(Lehmer)によって提案され, その後多数人々によって研究され方法であり, 線形漸化式



X_n=aX_{n-1}+c \ \ \ \ \ (\mbox{mod}\ \ m) \,


使って非負整数\{X_n\} \,生成する. ここで, a \,およびm \,正整数であり, c \,非負整数である. 特にc=0 \,場合には, 乗算合同法と呼ぶ. パラメタ選び方に関して多く研究結果があるが, 現在比較良いとされているものをいくつか表1に示す.

表1線形合同法使われるパラメタの例
a \, c \, m \, X_0 \, \mu_2 \, \mu_3 \, \mu_4 \, \mu_5 \, \mu_6 \, 周期
1\,664\,525 \, * \, 2^{32} \, 任意 16.1 10.6 8.0 6.0 5.0 2^{32} \,
1\,566\,083\,941 \, 0 2^{32} \, 奇数 14.8 9.7 7.5 5.6 4.2 2^{30} \,
48\,828\,125 \, 0 2^{32} \, 奇数 14.8 8.8 7.4 5.7 4.9 2^{30} \,
2\,100\,005\,341 \, 0 2^{31}-1 \, 正整数 15.4 10.2 7.7 6.0 5.0 2^{31}-2 \,
397\,204\,094 \, 0 2^{31}-1 \, 正整数 14.8 9.7 7.4 6.0 5.0 2^{31}-2 \,
314\,159\,369 \, 0 2^{31}-1 \, 正整数 15.2 9.9 7.6 5.9 5.1 2^{31}-2 \,

c \,の列の * \,印は, 任意の奇数使用してよいことを表す.


 線形合同法によって生成される乱数欠点として, 多次元結晶構造と言われているものがある. これは, k \,次元超立方体内にランダムに点を配置する目的で, 点の座標(x_n, x_{n+1},\cdots, x_{n+k-1}), n=1,2,\cdots, \,定めたとすると, これらの点はすべて比較少数等間隔並んだ超平面の上規則的にのってしまい, ランダムにならないという性質である. m \,超平面枚数の上界との関係を表2に示す. また, 表1\mu_k \,は, k \,次元の点配置作ったとき, \mu_k \,ビット精度ではほぼ一様な配置になることを意味している.


表2超平面枚数の上
m \, k=3 \, 4 5 6 7 8 9 10
2^{24} \, 465 141 72 47 36 30 26 23
2^{32} \, 2\,953 \, 566 220 120 80 60 48 41
2^{36} \, 7442 \, 1133 \, 383 191 119 85 66 54
2^{48} \, 119086 \, 9065 \, 2021 \, 766 391 240 167 126


M系列法

 ガロア体GF(2)上の任意の原始多項式



x^p+c_1x^{p-1}+c_2x^{p-2}+\cdots + c_p


選び, その係数係数とする漸化式



a_n = c_1a_{n-1}+c_2a_{n-2}+\cdots + c_pa_{n-p}
\ \ \ \ \ (\mbox{mod}\ \ 2)


によって生成される系列\{a_n\} \,考える. この系列M系列(M-sequence)あるいはシフトレジスタ系列, 極大多項式系列などと呼ばれ, 硬貨投げて表が出れば1, 裏が出れば0として得られる系列類似の性質を持つことが知られている. 表3実用的な原始多項式の例を示す.


表3原始多項式の例

x^{521}+x^{32}+1 \,
x^{607}+x^{273}+1 \,
x^{1279}+x^{418}+1 \,
x^{521}+x^{455}+x^{437}+x^{350}+1 \,
x^{607}+x^{461}+x^{307}+x^{167}+1 \,


 これは1ビット系列なので, 多数ビット系列\{X_n\} \,作るためには, 漸化式



X_n = c_1X_{n-1} \oplus c_2X_{n-2} \oplus \cdots \oplus c_pX_{n-p}


用いる. ここで, 記号\oplus \,2進法でのけた上りなしの足し算(ビットごとの排他的論理和)を表す. これによって生成される系列は, トーズワース(Tausworthe)系列あるいはGFSR(Generalized Feedback Shift Register)系列呼ばれている.

 この漸化式を使う場合初期値設定法については, [1] を参照するとよい. これを使って例え32ビット整数系列\{X_n\} \,生成すると, 合同法のような多次元結晶構造生じことはなく, 32ビット精度[p/32] \,次元まで一様に分布する. また, この系列周期2^p-1 \,であり, 自己相関関数の値は位相差2^p/32 \,以下ならばほぼ0に等しい.

 M系列使ったもう少し複雑な乱数生成法として, 最近提案されメルセンヌ・ツイスター(Mersenne Twister) [3] がある. これは上記のものに比べて, 同じ記憶容量で, はるかに長い周期と高い次元一様性達成できるという特徴有する. 詳細

http://www.math.keio.ac.jp/

参照するとよい.




参考文献

[1] 伏見正則, 『乱数』(UP応用数学選書12), 東京大学出版会, 1989.

[2] D. E. Knuth, The Art of Computer Programming, Vol.2: Seminumerical Algorithms, 2nd ed., Addison-Wesley, 1981. 渋谷政昭訳, 『準数値算法/乱数』, サイエンス社, 1981.

[3] M. Matsumoto and T. Nishimura, "Mersenne Twister: A 623-Dimensionally Equidistributed Uniform Pseudo-Random Number Generator," ACM Transactions on Modeling and Computer Simulation, 8 (1998), 3-30.





一様乱数と同じ種類の言葉


英和和英テキスト翻訳>> Weblio翻訳
英語⇒日本語日本語⇒英語
  

辞書ショートカット

すべての辞書の索引

「一様乱数」の関連用語

一様乱数のお隣キーワード
検索ランキング

   

英語⇒日本語
日本語⇒英語
   



一様乱数のページの著作権
Weblio 辞書 情報提供元は 参加元一覧 にて確認できます。

   
日本オペレーションズ・リサーチ学会日本オペレーションズ・リサーチ学会
Copyright (C) 2024 (社)日本オペレーションズ・リサーチ学会 All rights reserved.

©2024 GRAS Group, Inc.RSS