MCP рдХреЗ рдорд╛рдзреНрдпрдо рд╕реЗ рдХрдиреЗрдХреНрдЯ рдХрд░реЗрдВ тЖТ

рдЧрдгрдирд╛ рджрд░реНрдЬ рдХрд░реЗрдВ

рд╕реВрддреНрд░ (рдлреЙрд░реНрдореВрд▓рд╛)

Show calculation steps (1)
  1. Distribution Statistics

    Distribution Statistics: рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рд╡рд┐рддрд░рдг рд░реИрдВрдбрдо рдирдВрдмрд░ рдЬрдирд░реЗрдЯрд░

    Theoretical mean, median, and variance of the log-normal distribution from mu and sigma.

рд╡рд┐рдЬреНрдЮрд╛рдкрди

рдкрд░рд┐рдгрд╛рдо

рддреИрдпрд╛рд░ рдХрд┐рдП рдЧрдП рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рд░реИрдВрдбрдо рдирдВрдмрд░
10
рдзрдирд╛рддреНрдордХ рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рдорд╛рди (Box-Muller)
0.0.04424369296
1.2.742452308
2.0.8228509476
3.10.73468553
4.8.181534841
5.0.7208715469
6.5.599260011
7.0.7510678292
8.29.95311855
9.3.336275718
Theoretical mean E[X] = exp(╬╝ + ╧Г┬▓/2) 20.085537
Theoretical median = exp(╬╝) 2.718282
рд╕реИрджреНрдзрд╛рдВрддрд┐рдХ рдкреНрд░рд╕рд░рдг Var[X] 21,623.037001

рдпрд╣ рдЬрдирд░реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ

рдпрд╣ рдЯреВрд▓ рдРрд╕реЗ рд╕реНрдпреВрдбреЛ-рд░реИрдВрдбрдо рдирдВрдмрд░реЛрдВ рдХреА рдПрдХ рд╕реВрдЪреА рддреИрдпрд╛рд░ рдХрд░рддрд╛ рд╣реИ рдЬреЛ рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рд╡рд┐рддрд░рдг (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}$$
рдмреЙрдХреНрд╕-рдореБрд▓рд░ рд░реВрдкрд╛рдВрддрд░рдг рдХрд╛ рдкреНрд░рд╡рд╛рд╣ рдЖрд░реЗрдЦ рдЬреЛ рджреЛ рд╕рдорд╛рди рд╕рдВрдЦреНрдпрд╛рдУрдВ рдХреЛ рд▓реЙрдЧ-рдиреЙрд░реНрдорд▓ рдорд╛рди рдореЗрдВ рдмрджрд▓рддрд╛ рд╣реИ
рджреЛ рд╕рдорд╛рди рдпрд╛рджреГрдЪреНрдЫрд┐рдХ рд╕рдВрдЦреНрдпрд╛рдПрдБ рдмреЙрдХреНрд╕-рдореБрд▓рд░ рд╕реЗ рдорд╛рдирдХ рд╕рд╛рдорд╛рдиреНрдп Z рдмрдирддреА рд╣реИрдВ, рдлрд┐рд░ рд╕реНрдХреЗрд▓ рдФрд░ рдШрд╛рддрд╛рдВрдХрд┐рдд рд╣реЛрдХрд░ X рдмрди рдЬрд╛рддреА рд╣реИрдВред

рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг

\(\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)\) рд╣реИ, рдЗрд╕рд▓рд┐рдП рд╕рднреА рдЖрдЙрдЯрдкреБрдЯ рд╣рдореЗрд╢рд╛ рдзрдирд╛рддреНрдордХ рд╣реЛрддреЗ рд╣реИрдВред

рдЕрдВрддрд┐рдо рдЕрдкрдбреЗрдЯ: