์ด ์์ฑ๊ธฐ๋ ๋ฌด์์ ํ๋์
์ด ๋๊ตฌ๋ ๋ก๊ทธ์ ๊ท๋ถํฌ(log-normal distribution)๋ฅผ ๋ฐ๋ฅด๋ ์์ฌ๋์ ๋ชฉ๋ก์ ๋ง๋ค์ด ์ค๋๋ค. ์ด๋ค ํ๋ฅ ๋ณ์ X์ ์์ฐ๋ก๊ทธ \(\ln(X)\)๊ฐ ์ ๊ท๋ถํฌ๋ฅผ ๋ฐ๋ฅผ ๋, X๋ฅผ ๋ก๊ทธ์ ๊ท๋ถํฌ๋ผ๊ณ ๋ถ๋ฆ ๋๋ค. ์ด ์์ฑ๊ธฐ๋ ๋จผ์ ๊ณ ์ ์ ์ธ Box-Muller ๋ณํ์ผ๋ก ํ์ค์ ๊ท ๋์๋ฅผ ๋ง๋ ๋ค, ์ด๋ฅผ ์ํ๋ ์ ๊ท๋ถํฌ \(N(\mu, \sigma^2)\)์ ๋ง๊ฒ ๋ณํํ๊ณ , ๋ง์ง๋ง์ผ๋ก ์ง์ํจ์๋ฅผ ์ทจํด ํญ์ ์์์ธ ๋ก๊ทธ์ ๊ท ๊ฐ์ ์ป์ต๋๋ค.
์ค์: ฮผ์ ฯ๊ฐ ์๋ฏธํ๋ ๊ฒ
๋งค๊ฐ๋ณ์ \(\mu\)์ \(\sigma\)๋ \(\ln(X)\), ์ฆ ๋ฐํ์ด ๋๋ ์ ๊ท๋ถํฌ์ ํ๊ท ๊ณผ ํ์คํธ์ฐจ์ ๋๋ค. X ์์ฒด์ ํ๊ท ๊ณผ ํ์คํธ์ฐจ๊ฐ ์๋๋ผ๋ ์ ์ ์ฃผ์ํ์ธ์. X์ ์ค์ ํ๊ท ์ \(\exp(\mu + \sigma^2/2)\), ์ค์๊ฐ์ \(\exp(\mu)\), ๋ถ์ฐ์ \((\exp(\sigma^2) - 1)\cdot\exp(2\mu + \sigma^2)\)๋ก ๊ณ์ฐ๋ฉ๋๋ค. ์ด ์ด๋ก ๊ฐ๋ค์ ์์ฑ๋ ํ๋ณธ๊ณผ ํจ๊ป ํ์๋๋ฏ๋ก, ์ถ๋ ฅ ๊ฒฐ๊ณผ๊ฐ ์ฌ๋ฐ๋ฅธ์ง ๊ฒ์ฐํ๋ ๋ฐ ํ์ฉํ ์ ์์ต๋๋ค.
$$\begin{aligned} \text{Mean} &= \exp\!\left(\mu + \tfrac{1}{2}\sigma^{2}\right) \\ \text{Median} &= \exp\!\left(\mu\right) \\ \text{Var} &= \left(e^{\sigma^{2}} - 1\right) e^{\,2\mu + \sigma^{2}} \end{aligned}$$์ฌ์ฉ ๋ฐฉ๋ฒ
\(\mu\)(์์์ ์ค์), \(\sigma\)(0 ์ด์์ ๊ฐ), ๊ทธ๋ฆฌ๊ณ ์์ฑํ ๋์์ ๊ฐ์(1~1000๊ฐ)๋ฅผ ์ ๋ ฅํ์ธ์. ํ์ํ ์ ํจ์ซ์ ์๋ฆฟ์๋ ์ ํํฉ๋๋ค. ์คํํ ๋๋ง๋ค ์๋ก์ด ๊ท ๋ฑ ๋์๋ฅผ ๋ฝ๊ธฐ ๋๋ฌธ์ ๊ฐ์ ๋งค๋ฒ ๋ฌ๋ผ์ง๋๋ค. Box-Muller ๋ฐฉ์์ ๊ท ๋ฑ ๋์ ํ ์์์ ์ ๊ท ๋์๋ฅผ ๋ ๊ฐ์ฉ ๋ง๋ค์ด ๋ด๋ฏ๋ก, ๊ฐ์๋ฅผ ํ์๋ก ์ง์ ํ๋ฉด ๋จ๋ ํ ๊ฐ์ ๋์๋ ๊ทธ๋ฅ ๋ฒ๋ฆฝ๋๋ค.
๊ณต์ ์ค๋ช
(0,1] ๊ตฌ๊ฐ์ ๊ท ๋ฑ ๋์ \(U_1, U_2\)์ ๋ํด \(Z = \sqrt{-2\ln U_1}\,\cos(2\pi U_2)\)๋ฅผ ๊ณ์ฐํ๋ฉด \(Z \sim N(0,1)\)์ ์ป์ต๋๋ค. ๊ทธ๋ค์ \(Y = \mu + \sigma\, Z\)๋ \(N(\mu, \sigma^2)\)๋ฅผ ๋ฐ๋ฅด๊ณ , \(X = \exp(Y)\)๋ ๋ก๊ทธ์ ๊ท๋ถํฌ๊ฐ ๋ฉ๋๋ค. \(\ln(0)\)์ ๊ณ์ฐํ๋ ์ผ์ ๋ง๊ธฐ ์ํด \(U_1\)์ ์์ฃผ ์์ ์ก์ค๋ก (epsilon) ๊ฐ์ผ๋ก ํํ์ ๋ก๋๋ค.
$$\begin{gathered} X = \exp\!\left(\mu + \sigma\, Z\right) \\[1.5em] \text{where}\quad \left\{ \begin{aligned} Z &= \sqrt{-2\ln U_1}\,\cos(2\pi U_2) \\ U_1, U_2 &\sim \text{Uniform}(0,1) \end{aligned} \right. \end{gathered}$$
๊ณ์ฐ ์์
\(\mu = 1\), \(\sigma = 2\)๋ก ๋๊ณ \(U_1 = 0.5\), \(U_2 = 0.25\)๋ฅผ ์ฌ์ฉํด ๋ด ์๋ค. ๊ทธ๋ฌ๋ฉด \(R = \sqrt{-2\ln 0.5} = 1.17741\)์ ๋๋ค. \(Z_1 = R\cos(\pi/2) = 0\)์ด๊ณ \(Z_2 = R\sin(\pi/2) = 1.17741\)์ด ๋ฉ๋๋ค. ๋ฐ๋ผ์
$$X_1 = \exp(1) = 2.71828, \quad X_2 = \exp(1 + 2\cdot 1.17741) = \exp(3.35482) \approx 28.64$$์ด๋ก ์ ํ๊ท \(\exp(3) = 20.0855\)์ ์ค์๊ฐ \(\exp(1) = 2.71828\)๊ณผ๋ ์ ๋ค์ด๋ง์ต๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
์ ์คํํ ๋๋ง๋ค ๊ฐ์ด ๋ค๋ฅธ๊ฐ์? Math.random()์ ์ฌ์ฉํ๋ ํ๋ฅ ์ ์์ฑ๊ธฐ์ด๊ธฐ ๋๋ฌธ์ ๋๋ค. ์๋(seed)๋ฅผ ๊ณ ์ ํ์ง ์์ผ๋ฉด ์คํํ ๋๋ง๋ค ๊ฒฐ๊ณผ๊ฐ ๋ฌ๋ผ์ง๋๋ค.
ฯ = 0์ด๋ฉด ์ด๋ป๊ฒ ๋๋์? ๋ถํฌ๊ฐ ํดํ(degenerate)๋์ด ๋ชจ๋ ๊ฐ์ด \(\exp(\mu)\)๋ก ๋์ผํด์ง๋๋ค.
๊ฐ์ด ์์๊ฐ ๋ ์๋ ์๋์? ์๋๋๋ค. ๋ก๊ทธ์ ๊ท๋ถํฌ์ ์ ์์ญ์ \((0, \infty)\)์ด๋ฏ๋ก ๋ชจ๋ ์ถ๋ ฅ๊ฐ์ ํญ์ ์์์ ๋๋ค.