非一様乱数とは - OR事典 Weblio辞書

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

OR事典

非一様乱数

読み方：ひいちようらんすう
【英】：non-uniform random numbers

概要

一様分布以外の分布をする乱数の総称. 一様乱数に適当な変換を施して作るのがふつうである. 変換法として, 逆関数法, 棄却法, 合成法, 別名法, 一様乱数の比を利用する方法, 等の種々の分布に適用可能な方法と, 個々の分布に特有の方法が多数提案されている.

詳説

一様分布以外の確率分布に従う乱数を非一様乱数という. これらの乱数は, 一様乱数に適当な変換を施して作るのがふつうである. 変換方法には, 種々の分布に対して適用可能な一般的なものと, 個々の分布の特徴を利用するものとがある. これらの方法を網羅的に扱っているのが [3]である. ここでは, 一般的な変換法を二つ紹介した後, 実用上大事な非一様乱数生成法を述べる. 以下では, (整数型の)一様乱数を正規化して得られる区間[0,1)上の(実数型の)一様乱数を $U, U_1, U_2\,$ 等で表し, これらは互いに独立であるものと仮定する. また, 発生させたい乱数を $X\,$ とし, その分布関数を $F(x)\,$ , それが連続分布ならば, その密度関数を $f(x)\,$ で表す.

[逆関数法]

$F(x)\,$ の逆関数 $F^{-1}(x)\,$ を使って $X=F^{-1}(U)\,$ とするものである. 例えば, 指数分布, ワイブル分布, ロジスティック分布(スケール・パラメータはいずれも1とする)なら, それぞれ $X=-\log U\,$ , $X=(-\log U)^\alpha\,$ , $X=\log(U/(1-U))\,$ とすればよい. 逆関数が解析的に求まらない場合には, 近似式等を用いることになる.

[二者択一法]

ウォーカー(A. J. Walker) [4] によって提案されたもので, 別名法(alias method)ともいう.

取りうる値の個数 $(n)\,$ が有限個の離散分布に従う乱数を, $n\,$ に無関係なO(1)の速度で生成できる, 簡単で効率的な方法である. ただし, 初期設定にO $(n)\,$ の時間とO $(n)\,$ のメモリを必要とする. 例えばポアソン分布のように, 取りうる値の個数が無限の場合には, 分布の裾を適当に打ち切って適用すればよい.

[正規分布]

逆関数法を適用する場合には, 逆関数 $F^{-1}(x)\,$ が解析的には求まらないので, 例えば次の山内二郎の近似式を使う.

$\begin{array}{rll} F^{-1}(x)& = &\left\{ \begin{array}{ll} -w, & \ \ \ \ \ 0 < x < 0.5 \\ +w, & \ \ \ \ \ 0.5 \leq x < 1.0 \end{array}\right.\\ w & = & \{ z(2.0611786-\frac{5.7262204}{z+11.640595})\}^{1/2}\\ z & = & -\log\{4x(1-x)\} \end{array} \,$

ボックス・ミュラー(Box-Muller)法は, 変換公式

$X_1=\sqrt{-2\log U_1}\cos(2\pi U_2), \,$

$X_2=\sqrt{-2\log U_1}\sin(2\pi U_2) \,$

を使って, 2個の一様乱数 $U_1,U_2\,$ から互いに独立な標準正規乱数 $X_1,X_2\,$ を作るものである.

平均ベクトルが $\boldsymbol{\mu}\,$ , 分散共分散行列が $V\,$ の $n\,$ 次元正規乱数ベクトル ${\mathbf{ Y}}\,$ を生成するためには, 互いに独立な $n\,$ 個の標準正規乱数からなるベクトル ${\mathbf{ X}}\,$ に変換

$\mathbf{Y} = \boldsymbol{\mu} + A \mathbf{X} \,$

を施せばよい. ただし, ベクトルはすべて列ベクトルとし, $A\,$ は $V=AA^\top\,$ ( $A^\top\,$ は $A\,$ の転置行列)を満たす下三角行列であり, コレスキー(Cholesky)分解によって求める.

[アーラン分布]

フェイズが $k\,$ (正整数)のアーラン分布の密度関数は

$f(x)=\lambda \mathrm{e}^{-\lambda x}(\lambda x)^{k-1}/(k-1)! \ \ \ \ \ \ \ (x\geq 0) \,$

である. これは, パラメータ(平均値の逆数)が $\lambda\,$ の指数分布に従う $k\,$ 個の確率変数の和の分布であるから,

$X=-\log(U_1\cdots U_k)/\lambda \,$

とすればよい.

[ポアソン分布]

平均値が $\lambda\,$ のポアソン乱数は, 平均が $1/\lambda\,$ の指数乱数との関係を利用して作ることができる. すなわち, $U_1U_2\cdots U_n > {\mbox{e}}^{-\lambda}\,$ が成り立つ最大の整数 $n\,$ を $X\,$ とすればよい. 1個のポアソン乱数を発生するのに必要な一様乱数の個数の平均値は $\lambda+1\,$ であるから, $\lambda\,$ が大変に大きいときには, この方法は効率が悪いので, 二者択一法を使うほうがよい.

参考文献

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

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

[3] L. Devroye, Non-Uniform Random Variate Generation, Springer-Verlag, 1986.

[4] A. J. Walker, "An Efficient Method for Generating Discrete Random Variables with General Distributions," ACM Transactions on Mathematical Software, 3 (1977), 253-256.

「OR事典」の他の用語

シミュレーション：

負相関変量法連続型シミュレーション連続型シミュレーション言語重点サンプリング離散型シミュレーション離散型シミュレーション言語非一様乱数

索引トップ用語の索引ランキングカテゴリー

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

>>同じ種類の言葉 >>高等数学に関連する言葉

>> 「非一様乱数」を含む用語の索引

非一様乱数のページへのリンク

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

1		井上和子
2		沓名環希
3		井上和子 (女優)
4		井上和子_(女優)
5		イノウエカズコ
6		行列のできる法律相談所
7		同居人カップルの殺人推理旅行
8		メガゾーン２３　PARTⅡ　秘密く・だ・さ・い
9		TPP
10		井上和子 (言語学者)

11		捗る
12		リーガル・ハイ
13		出不精
14		福利厚生
15		アセスメント
16		奇跡の人
17		自由民主党幹事長
18		友引
19		ダイバーシティ
20		蛍/少年

21		言葉
22		功罪
23		未満
24		タカラダニ
25		石松千恵美
26		ハクビシン
27		先負
28		土曜ワイド劇場
29		ご指導ご鞭撻
30		注意喚起


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

非一様乱数とは？

OR事典

非一様乱数

概要

詳説