The Modifikovaný zákon o normálním napájení (MLP) function is a three parameter function that can be used to model data that have characteristics of a normální distribuce protokolu a a mocenský zákon chování. Používá se k modelování funkční formy Funkce počáteční hmotnosti (MMF). Na rozdíl od ostatních funkčních forem MMF je MLP jedinou funkcí bez podmínek spojení.
Funkční forma distribuce MLP Uzavřená forma funkce hustoty pravděpodobnosti MLP je následující: F ( m ) = α 2 exp ( α μ 0 + α 2 σ 0 2 2 ) m − ( 1 + α ) erfc ( 1 2 ( α σ 0 − ln ( m ) − μ 0 σ 0 ) ) , m ∈ [ 0 , ∞ ) { displaystyle { begin {aligned} f (m) = { frac { alpha} {2}} exp left ( alpha mu _ {0} + { frac { alpha ^ {2} sigma _ {0} ^ {2}} {2}} vpravo) m ^ {- (1+ alpha)} { text {erfc}} left ({ frac {1} { sqrt {2} }} left ( alpha sigma _ {0} - { frac { ln (m) - mu _ {0}} { sigma _ {0}}} right) right), m v [0, infty) end {zarovnáno}}} kde α = δ y { displaystyle { begin {zarovnáno} alpha = { frac { delta} { gamma}} konec {zarovnáno}}} je asymptotický index mocninového práva distribuce. Tady μ 0 { displaystyle mu _ {0}} a σ 0 2 { displaystyle sigma _ {0} ^ {2}} jsou průměr a rozptyl podkladové lognormální distribuce, ze které je odvozena MLP.
Matematické vlastnosti distribuce MLP Následuje několik matematických vlastností distribuce MLP:
Kumulativní distribuce MLP kumulativní distribuční funkce ( F ( m ) = ∫ − ∞ m F ( t ) d t { displaystyle F (m) = int _ {- infty} ^ {m} f (t) , dt} ) darováno:
F ( m ) = 1 2 erfc ( − ln ( m ) − μ 0 2 σ 0 ) − 1 2 exp ( α μ 0 + α 2 σ 0 2 2 ) m − α erfc ( α σ 0 2 ( α σ 0 − ln ( m ) − μ 0 2 σ 0 ) ) { displaystyle { begin {aligned} F (m) = { frac {1} {2}} { text {erfc}} left (- { frac { ln (m) - mu _ {0 }} {{ sqrt {2}} sigma _ {0}}} vpravo) - { frac {1} {2}} exp left ( alpha mu _ {0} + { frac { alpha ^ {2} sigma _ {0} ^ {2}} {2}} vpravo) m ^ {- alpha} { text {erfc}} vlevo ({ frac { alpha sigma _ {0}} { sqrt {2}}} left ( alpha sigma _ {0} - { frac { ln (m) - mu _ {0}} {{ sqrt {2}} sigma _ {0}}} doprava) doprava) konec {zarovnáno}}} Můžeme to vidět jako m → 0 , { displaystyle m až 0,} že F ( m ) → 1 2 erfc ( − ln ( m − μ 0 ) 2 σ 0 ) , { displaystyle textstyle F (m) až { frac {1} {2}} operatorname {erfc} left (- { frac { ln (m- mu _ {0})} {{ sqrt {2}} sigma _ {0}}} vpravo),} což je kumulativní distribuční funkce pro lognormální distribuci s parametry μ 0 a σ 0 .
Zlá, odchylka, hrubé momenty The očekávaná hodnota z M { displaystyle M} k dává k { displaystyle k} th surový okamžik z M { displaystyle M} ,
⟨ M k ⟩ = ∫ 0 ∞ m k F ( m ) d m { displaystyle { begin {zarovnáno} langle M ^ {k} rangle = int _ {0} ^ { infty} m ^ {k} f (m) mathrm {d} m end {zarovnáno} }} To existuje tehdy a jen tehdy, když α> k { displaystyle k} , v takovém případě se stane:
⟨ M k ⟩ = α α − k exp ( σ 0 2 k 2 2 + μ 0 k ) , α > k { displaystyle { begin {aligned} langle M ^ {k} rangle = { frac { alpha} { alpha -k}} exp left ({ frac { sigma _ {0} ^ { 2} k ^ {2}} {2}} + mu _ {0} k right), alpha> k end {zarovnáno}}} který je k { displaystyle k} th hrubý moment lognormálního rozdělení s parametry μ0 a σ0 zmenšenα ⁄α- k { displaystyle k} v limitu α → ∞. To dává průměr a rozptyl distribuce MLP:
⟨ M ⟩ = α α − 1 exp ( σ 0 2 2 + μ 0 ) , α > 1 { displaystyle { begin {zarovnáno} langle M rangle = { frac { alpha} { alpha -1}} exp left ({ frac { sigma _ {0} ^ {2}} { 2}} + mu _ {0} right), alpha> 1 end {zarovnáno}}} ⟨ M 2 ⟩ = α α − 2 exp ( 2 ( σ 0 2 + μ 0 ) ) , α > 2 { displaystyle { begin {aligned} langle M ^ {2} rangle = { frac { alpha} { alpha -2}} exp left (2 left ( sigma _ {0} ^ { 2} + mu _ {0} right) right), alpha> 2 end {zarovnáno}}} Var ( M { displaystyle M} ) = ⟨ M { displaystyle M} 2 ⟩-(⟨ M { displaystyle M} ⟩)2 = α exp (σ0 2 + 2 μ0 ) (exp (σ0 2 ) / α-2 - α / (α-2)2 ), α> 2
Režim Řešení rovnice F ′ ( m ) { displaystyle f '(m)} = 0 (rovnice sklonu k nule v bodě maxim) pro m { displaystyle m} udává režim distribuce MLP.
F ′ ( m ) = 0 ⇔ K. erfc ( u ) = exp ( − u 2 ) , { displaystyle f '(m) = 0 Leftrightarrow K operatorname {erfc} (u) = exp (-u ^ {2}),} kde u = 1 2 ( α σ 0 − ln m − μ 0 σ 0 ) { displaystyle textstyle u = { frac {1} { sqrt {2}}} left ( alpha sigma _ {0} - { frac { ln m- mu _ {0}} { sigma _ {0}}} vpravo)} a K. = σ 0 ( α + 1 ) π 2 . { displaystyle K = sigma _ {0} ( alpha +1) { tfrac { sqrt { pi}} {2}}.}
K vyřešení této transcendentální rovnice jsou zapotřebí numerické metody. Je však třeba poznamenat, že pokud K. { displaystyle K} ≈1 pak u = 0 nám dává režim m { displaystyle m} * :
m ∗ = exp ( μ 0 + α σ 0 2 ) { displaystyle m ^ {*} = exp ( mu _ {0} + alpha sigma _ {0} ^ {2})} Náhodné variace Lognormální náhodná varieta je:
L ( μ , σ ) = exp ( μ + σ N ( 0 , 1 ) ) { displaystyle { begin {zarovnáno} L ( mu, sigma) = exp ( mu + sigma N (0,1)) end {zarovnáno}}} kde N ( 0 , 1 ) { displaystyle N (0,1)} je standardní normální náhodná variace. Exponenciální náhodná variace je:
E ( δ ) = − δ − 1 ln ( R ( 0 , 1 ) ) { displaystyle { begin {zarovnáno} E ( delta) = - delta ^ {- 1} ln (R (0,1)) end {zarovnáno}}} kde R (0,1) je jednotná náhodná variace v intervalu [0,1]. Pomocí těchto dvou můžeme odvodit náhodnou variaci pro distribuci MLP na:
M ( μ 0 , σ 0 , α ) = exp ( μ 0 + σ 0 N ( 0 , 1 ) − α − 1 ln ( R ( 0 , 1 ) ) ) { displaystyle { begin {sladěno} M ( mu _ {0}, sigma _ {0}, alpha) = exp ( mu _ {0} + sigma _ {0} N (0,1 ) - alpha ^ {- 1} ln (R (0,1))) end {zarovnáno}}} Reference Basu, Shantanu; Gil, M; Auddy, Sayatan (1. dubna 2015). „Distribuce MLP: upravený lognormální model zákonu moci pro hvězdnou počáteční hmotnostní funkci“ . MNRAS . 449 (3): 2413–2420. arXiv :1503.00023 . Bibcode :2015MNRAS.449.2413B . doi :10.1093 / mnras / stv445 .