ボックス＝ミュラー法

ファイル:Box-Muller transform visualisation.svg

ボックス=ミュラー法のイメージ。単位正方形(u1, u2)内部の色の付いた点は円状に散布され、2次元ガウス分布となる。上と右の余白に点で示される曲線は変換後の確率密度関数である。図は有限の区間でのプロットであるが、実際には変換後の分布は無限に広がる。the SVG fileでは、カーソルを合わせた点と関係する点をハイライトする。

ボックス＝ミュラー法（ボックス＝ミュラーほう、英: Box–Muller's method）とは、一様分布に従う確率変数から標準ガウス分布に従う確率変数を生成させる手法。計算機シミュレーションにおいて、ガウス分布に従う擬似乱数の発生に応用される。統計学者ジョージ・ボックスとマーヴィン・マラー（ミュラー）によって考案された。

概要

確率変数 X 及び Y が互いに独立で、ともに（0, 1）上での一様分布に従うものとする。このとき、

[math] \begin{align} Z_1 &= \sqrt{-2 \log{X}} \cos{2\pi Y}, \\ Z_2 &= \sqrt{-2 \log{X}} \sin{2\pi Y} \end{align} [/math]

で定義される Z₁, Z₂ は、平均 0、分散 1 の標準ガウス分布N(0,1)に従う互いに独立な確率変数となる。一様分布に従う X 及び Y からガウス分布に従う Z₁, Z₂ を与えるこの変換をボックス＝ミュラー変換という。また、このガウス分布に従う確率変数を生成させる方法のことをボックス＝ミュラー法という。ボックス＝ミュラー法によって、比較的生成が容易な一様分布に従う乱数から、応用上、重要なガウス分布に従う乱数を生成させることができる。

発想

2次元の標準ガウス分布に従う (Z₁, Z₂) において、2変数が互いに独立であれば、結合確率密度関数

[math] f(z_1,z_2) = \frac{1}{2\pi}e^{-\frac{z_1^{\, 2} + z_2^{\, 2}}{2} } [/math]

は、円周上で定数値を与えることから、偏角

[math] \Theta = \arctan{\frac{Z_2}{Z_1}} [/math]

は (0, 2π) 上で、一様分布をなす。一方、2次元ベクトル (Z₁, Z₂) の大きさの2乗

[math] R^2= Z_1^{\, 2} + Z_2^{\, 2} [/math]

は自由度2のカイ二乗分布に従う。ここで、カイ二乗分布の性質から e^−R2/2 は、(0, 1) 上の一様分布となる。

これらのことから、逆に (0, 1) 上で一様分布する2つの独立な確率変数 X, Y により、

[math] \begin{align} \Theta &= 2 \pi Y ,\\ R^2 &= - 2\log{X} \end{align} [/math]

とすれば、

[math] \begin{align} Z_1 &= R \cos{\Theta} ,\\ Z_2 &= R \sin{\Theta} \end{align} [/math]

で定義される確率変数 Z₁, Z₂ は標準ガウス分布 N(0, 1) に従うこととなる。

証明

ボックス=ミュラー変換が標準ガウス分布を与えることは、特性関数を調べることで確認できる。実際、Z₁ については、その特性関数は、

[math] \Phi_{Z_1}(\xi)= \langle e^{i\xi Z_1} \rangle = \int_0^1 \int_0^1 e^{i \xi \sqrt{-2 \log{x}} \cos{2\pi y}}dxdy [/math]

であり、変数変換

[math] x = e^{ -\frac{r^2}{2} }, \, y = \frac{\theta}{2\pi}, \, dxdy=|J|drd\theta = \frac{r}{2\pi} e^{- \frac{r^2}{2} } drd\theta [/math]

によって、

[math] \Phi_{Z_1}(\xi) = \frac{1}{2\pi} \int_0^{\infty} \int_0^{2\pi} e^{i \xi r \cos{\theta}-\frac{r^2}{2} } r dr d\theta [/math]

となるが、さらに変数変換

[math] z=r \cos{\theta}, \, w=r \sin{\theta}, \, dzdw=|J|dr d \theta= r dr d \theta [/math]

を行えば、

[math] \begin{align} \Phi_{Z_1}(\xi) &= \frac{1}{2\pi} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} e^{i \xi z -\frac{z^2+w^2}{2} } dz dw \\ &=e^{-\frac{\xi^2}{2} } \frac{1}{2\pi} \int_{-\infty}^{\infty} e^{-\frac{(z-i \xi)^2}{2} } dz \int_{-\infty}^{\infty} e^{-\frac{w^2}{2}} dw \\ &=e^{-\frac{\xi^2}{2}} \end{align} [/math]

を得る。これは、標準ガウス分布 N(0, 1) の特性関数にほかならない。

参考文献

原論文

• G.E.P. Box and M.E. Muller ,“A note on the generation of random normal deviates,” Annals Math. Stat, Vol. 29, pp. 610-611 (1958) 　doi:10.1214/aoms/1177706645

参考書籍

• 四辻哲章 『計算機シミュレーションのための確率分布乱数生成法 』 プレアデス出版 (2010年),　ISBN 978-4903814353
• 添田喬, 太田光雄, 大松繁 『数理統計の基礎と応用』 日新出版 (2000年),　ISBN 978-4817301079

関連項目

• 乱数列
• ガウス分布