рдпрд╣ рдЬрдирд░реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ
рдпрд╣ рдЯреВрд▓ рдРрд╕реЗ рд╕реНрдпреВрдбреЛ-рд░реИрдВрдбрдо рдирдВрдмрд░реЛрдВ рдХреА рдПрдХ рд╕реВрдЪреА рддреИрдпрд╛рд░ рдХрд░рддрд╛ рд╣реИ рдЬреЛ рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рд╡рд┐рддрд░рдг (log-normal distribution) рдХрд╛ рдкрд╛рд▓рди рдХрд░рддреЗ рд╣реИрдВред рдХреЛрдИ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рдЪрд░ X рддрдм рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рдХрд╣рд▓рд╛рддрд╛ рд╣реИ рдЬрдм рдЙрд╕рдХрд╛ рдкреНрд░рд╛рдХреГрддрд┐рдХ рд▓рдШреБрдЧрдгрдХ, рдпрд╛рдиреА \(\ln(X)\), рдиреЙрд░реНрдорд▓ рд╡рд┐рддрд░рдг рдореЗрдВ рдЖрддрд╛ рд╣реИред рдЬрдирд░реЗрдЯрд░ рдкрд╣рд▓реЗ рдкреНрд░рд╕рд┐рджреНрдз Box-Muller рдЯреНрд░рд╛рдВрд╕рдлрд╝реЙрд░реНрдо рд╕реЗ рдорд╛рдирдХ-рдиреЙрд░реНрдорд▓ рдорд╛рди рдмрдирд╛рддрд╛ рд╣реИ, рдЙрдиреНрд╣реЗрдВ рд▓рдХреНрд╖рд┐рдд рдиреЙрд░реНрдорд▓ \(N(\mu, \sigma^2)\) рддрдХ рд╕реНрдХреЗрд▓ рдХрд░рддрд╛ рд╣реИ, рдФрд░ рдлрд┐рд░ рдЙрдирдХрд╛ exponential (рдШрд╛рддрд╛рдВрдХ) рд▓реЗрдХрд░ рд╣рдореЗрд╢рд╛ рдзрдирд╛рддреНрдордХ рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рдорд╛рди рдкреНрд░рд╛рдкреНрдд рдХрд░рддрд╛ рд╣реИред
рдЬрд╝рд░реВрд░реА рдмрд╛рдд: mu рдФрд░ sigma рдХрд╛ рдЕрд╕рд▓реА рдорддрд▓рдм
рдкреИрд░рд╛рдореАрдЯрд░ \(\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\) (рд╢реВрдиреНрдп рдпрд╛ рдзрдирд╛рддреНрдордХ), рдФрд░ рдЖрдк рдХрд┐рддрдиреЗ рдирдВрдмрд░ рдЪрд╛рд╣рддреЗ рд╣реИрдВ (1 рд╕реЗ 1000 рддрдХ) рджрд░реНрдЬ рдХрд░реЗрдВред рддрдп рдХрд░реЗрдВ рдХрд┐ рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХ (significant digits) рджрд┐рдЦрд╛рдиреЗ рд╣реИрдВред рд╣рд░ рдмрд╛рд░ рдЪрд▓рд╛рдиреЗ рдкрд░ рдирдП рдпреВрдирд┐рдлрд╝реЙрд░реНрдо рд░реИрдВрдбрдо рдирдВрдмрд░ рдирд┐рдХрд╛рд▓реЗ рдЬрд╛рддреЗ рд╣реИрдВ, рдЗрд╕рд▓рд┐рдП рдорд╛рди рд╣рд░ рдмрд╛рд░ рдмрджрд▓рддреЗ рд╣реИрдВред рдЪреВрдБрдХрд┐ Box-Muller рд╡рд┐рдзрд┐ рд╣рд░ рдпреВрдирд┐рдлрд╝реЙрд░реНрдо рдЬреЛрдбрд╝реА рд╕реЗ рджреЛ рдиреЙрд░реНрдорд▓ рдорд╛рди рджреЗрддреА рд╣реИ, рдЕрдЧрд░ рдЖрдк рд╡рд┐рд╖рдо (odd) рд╕рдВрдЦреНрдпрд╛ рдорд╛рдБрдЧрддреЗ рд╣реИрдВ рддреЛ рдПрдХ рдЕрддрд┐рд░рд┐рдХреНрдд рдорд╛рди рдХреЛ рдмрд╕ рдЫреЛрдбрд╝ рджрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред
рд╕реВрддреНрд░ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛
(0,1] рдореЗрдВ рдпреВрдирд┐рдлрд╝реЙрд░реНрдо рдорд╛рди U1, U2 рдХреЗ рд▓рд┐рдП: \(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)\) рд╕реЗ рдмрдЪрдиреЗ рдХреЗ рд▓рд┐рдП рд╣рдо U1 рдХреЛ рдПрдХ рдмрд╣реБрдд рдЫреЛрдЯреЗ 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\cdot\cos(\pi/2) = 0\) рдФрд░ \(Z_2 = R\cdot\sin(\pi/2) = 1.17741\)ред рдЗрд╕рд▓рд┐рдП \(X_1 = \exp(1) = 2.71828\) рдФрд░ \(X_2 = \exp(1 + 2\cdot 1.17741) = \exp(3.35482) \approx 28.64\)ред рд╕реИрджреНрдзрд╛рдВрддрд┐рдХ рдорд╛рдзреНрдп \(\exp(3) = 20.0855\) рдФрд░ рдордзреНрдпрд┐рдХрд╛ \(\exp(1) = 2.71828\) рдЗрди рдорд╛рдиреЛрдВ рдХреЗ рдЕрдиреБрд░реВрдк рд╣реИрдВред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рд╕рд╡рд╛рд▓
рд╣рд░ рдмрд╛рд░ рдЪрд▓рд╛рдиреЗ рдкрд░ рдореЗрд░реЗ рдорд╛рди рдЕрд▓рдЧ рдХреНрдпреЛрдВ рдЖрддреЗ рд╣реИрдВ? рдпрд╣ Math.random() рдХрд╛ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд░рдиреЗ рд╡рд╛рд▓рд╛ рдПрдХ рдпрд╛рджреГрдЪреНрдЫрд┐рдХ (stochastic) рдЬрдирд░реЗрдЯрд░ рд╣реИ; рдмрд┐рдирд╛ рдХрд┐рд╕реА рдирд┐рд╢реНрдЪрд┐рдд seed рдХреЗ рд╣рд░ рдмрд╛рд░ рдкрд░рд┐рдгрд╛рдо рдЕрд▓рдЧ рд╣реЛрддреЗ рд╣реИрдВред
рдЕрдЧрд░ sigma = 0 рд╣реЛ рддреЛ? рддрдм рд╡рд┐рддрд░рдг рдЕрдкрднреНрд░рд╖реНрдЯ (degenerate) рд╣реЛ рдЬрд╛рддрд╛ рд╣реИ рдФрд░ рд╣рд░ рдорд╛рди \(\exp(\mu)\) рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛрддрд╛ рд╣реИред
рдХреНрдпрд╛ рдХреЛрдИ рдорд╛рди рдЛрдгрд╛рддреНрдордХ рд╣реЛ рд╕рдХрддрд╛ рд╣реИ? рдирд╣реАрдВ тАФ рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рдХрд╛ рдкрд░рд┐рд╕рд░ \((0, \infty)\) рд╣реИ, рдЗрд╕рд▓рд┐рдП рд╕рднреА рдЖрдЙрдЯрдкреБрдЯ рд╣рдореЗрд╢рд╛ рдзрдирд╛рддреНрдордХ рд╣реЛрддреЗ рд╣реИрдВред