The O'zgartirilgan Lognormal kuch qonuni (MLP) funktsiyasi - bu uchta parametr funktsiyasidir, u a ning xususiyatlariga ega bo'lgan ma'lumotlarni modellashtirish uchun ishlatilishi mumkin normal taqsimot va a kuch qonuni xulq-atvor. Ning funktsional shaklini modellashtirish uchun ishlatilgan Dastlabki ommaviy funktsiya (XVF). XVFning boshqa funktsional shakllaridan farqli o'laroq, MLP qo'shilish shartlari bo'lmagan yagona funktsiyadir.
MLP taqsimotining funktsional shakli
MLP ning zichlik funktsiyasining yopiq shakli quyidagicha: f ( m ) = a 2 tugatish ( a m 0 + a 2 σ 0 2 2 ) m − ( 1 + a ) erfc ( 1 2 ( a σ 0 − ln ( m ) − 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}} o'ng) m ^ {- (1+ alfa)} { text {erfc}} left ({ frac {1} { sqrt {2} }} chap ( alfa sigma _ {0} - { frac { ln (m) - mu _ {0}} { sigma _ {0}}} o'ng) o'ng), m ichida [0, infty) end {hizalangan}}} qayerda a = δ γ { displaystyle { begin {aligned} alpha = { frac { delta} { gamma}} end {aligned}}} taqsimotning asimptotik kuch-qonun ko'rsatkichidir. Bu yerda m 0 { displaystyle mu _ {0}} va σ 0 2 { displaystyle sigma _ {0} ^ {2}} MLP kelib chiqadigan asosiy lognormal taqsimotning o'rtacha va dispersiyasi.
MLP taqsimotining matematik xususiyatlari
Quyida MLP taqsimotining bir nechta matematik xususiyatlari keltirilgan:
Kümülatif taqsimot MLP kümülatif taqsimlash funktsiyasi ( F ( m ) = ∫ − ∞ m f ( t ) d t { displaystyle F (m) = int _ {- infty} ^ {m} f (t) , dt} ) tomonidan berilgan:
F ( m ) = 1 2 erfc ( − ln ( m ) − m 0 2 σ 0 ) − 1 2 tugatish ( a m 0 + a 2 σ 0 2 2 ) m − a erfc ( a σ 0 2 ( a σ 0 − ln ( m ) − m 0 2 σ 0 ) ) { displaystyle { begin {aligned} F (m) = { frac {1} {2}} { text {erfc}} left (- { frac { ln (m) - mu _ {0) }} {{ sqrt {2}} sigma _ {0}}} o'ng) - { frac {1} {2}} exp left ( alpha mu _ {0} + { frac { alfa ^ {2} sigma _ {0} ^ {2}} {2}} o'ng) m ^ {- alpha} { text {erfc}} chap ({ frac { alpha sigma _ {0}} { sqrt {2}}} chap ( alfa sigma _ {0} - { frac { ln (m) - mu _ {0}} {{ sqrt {2}} sigma _ {0}}} right) right) end {aligned}}} Buni biz shunday ko'rishimiz mumkin m → 0 , { displaystyle m dan 0 gacha,} bu F ( m ) → 1 2 erfc ( − ln ( m − m 0 ) 2 σ 0 ) , { displaystyle textstyle F (m) to { frac {1} {2}} operatorname {erfc} left (- { frac { ln (m- mu _ {0})} {{ sqrt {2}} sigma _ {0}}} o'ng),} parametrlarga ega lognormal taqsimot uchun kümülatif taqsimlash funktsiyasi m 0 va σ 0 .
O'rtacha, o'zgaruvchanlik, xom lahzalar The kutish qiymati ning M { displaystyle M} k beradi k { displaystyle k} th xom lahza ning M { displaystyle M} ,
⟨ M k ⟩ = ∫ 0 ∞ m k f ( m ) d m { displaystyle { begin {aligned} langle M ^ {k} rangle = int _ {0} ^ { infty} m ^ {k} f (m) mathrm {d} m end {aligned} }} Bu faqat a> bo'lsa, mavjud bo'ladi k { displaystyle k} , bu holda:
⟨ M k ⟩ = a a − k tugatish ( σ 0 2 k 2 2 + m 0 k ) , a > 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 {hizalangan}}} qaysi k { displaystyle k} th m parametrlari bilan lognormal taqsimotning xom momenti0 va σ0 miqyosia ⁄a- k { displaystyle k} a → ∞ chegarasida. Bu MLP taqsimotining o'rtacha va farqlanishini beradi:
⟨ M ⟩ = a a − 1 tugatish ( σ 0 2 2 + m 0 ) , a > 1 { displaystyle { begin {aligned} langle M rangle = { frac { alpha} { alpha -1}} exp left ({ frac { sigma _ {0} ^ {2}} { 2}} + mu _ {0} right), alpha> 1 end {hizalangan}}} ⟨ M 2 ⟩ = a a − 2 tugatish ( 2 ( σ 0 2 + m 0 ) ) , a > 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 {hizalangan}}} Var ( M { displaystyle M} ) = ⟨ M { displaystyle M} 2 ⟩-(⟨ M { displaystyle M} ⟩)2 = a exp (b0 2 + 2m0 ) (exp (σ.)0 2 ) / a-2 - a / (a-2)2 ), a> 2
Rejim Tenglamaning echimi f ′ ( m ) { displaystyle f '(m)} = 0 (qiyalikni maksimal darajadagi nuqtaga nolga tenglashtirish) uchun m { displaystyle m} MLP tarqatish rejimini beradi.
f ′ ( m ) = 0 ⇔ K erfc ( siz ) = tugatish ( − siz 2 ) , { displaystyle f '(m) = 0 Leftrightarrow K operator nomi {erfc} (u) = exp (-u ^ {2}),} qayerda siz = 1 2 ( a σ 0 − ln m − m 0 σ 0 ) { displaystyle textstyle u = { frac {1} { sqrt {2}}} chap ( alfa sigma _ {0} - { frac { ln m- mu _ {0}} { sigma _ {0}}} o'ng)} va K = σ 0 ( a + 1 ) π 2 . { displaystyle K = sigma _ {0} ( alfa +1) { tfrac { sqrt { pi}} {2}}.}
Ushbu transandantal tenglamani echish uchun raqamli usullar talab qilinadi. Biroq, agar shunday bo'lsa K { displaystyle K} -1, keyin u = 0 bizga rejimni beradi m { displaystyle m} * :
m ∗ = tugatish ( m 0 + a σ 0 2 ) { displaystyle m ^ {*} = exp ( mu _ {0} + alpha sigma _ {0} ^ {2})} Tasodifiy o'zgaruvchanlik Lognormal tasodifiy o'zgarish quyidagicha:
L ( m , σ ) = tugatish ( m + σ N ( 0 , 1 ) ) { displaystyle { begin {aligned} L ( mu, sigma) = exp ( mu + sigma N (0,1)) end {aligned}}} qayerda N ( 0 , 1 ) { displaystyle N (0,1)} standart odatiy tasodifiy o'zgaruvchidir. Eksponensial tasodifiy miqdor quyidagicha:
E ( δ ) = − δ − 1 ln ( R ( 0 , 1 ) ) { displaystyle { begin {aligned} E ( delta) = - delta ^ {- 1} ln (R (0,1)) end {aligned}}} bu erda R (0,1) - [0,1] oralig'idagi bir xil tasodifiy o'zgaruvchidir. Ushbu ikkitadan foydalanib, MLP taqsimotining tasodifiy o'zgarishini quyidagicha olishimiz mumkin:
M ( m 0 , σ 0 , a ) = tugatish ( m 0 + σ 0 N ( 0 , 1 ) − a − 1 ln ( R ( 0 , 1 ) ) ) { displaystyle { begin {aligned} M ( mu _ {0}, sigma _ {0}, alpha) = exp ( mu _ {0} + sigma _ {0} N (0,1) ) - alpha ^ {- 1} ln (R (0,1))) end {hizalangan}}} Adabiyotlar
Basu, Shantanu; Gil, M; Auddy, Sayatan (2015 yil 1-aprel). "MLP taqsimoti: yulduzlarning boshlang'ich massasi funktsiyasi uchun modifikatsiyalangan kuch-qonun modeli" . MNRAS . 449 (3): 2413–2420. arXiv :1503.00023 . Bibcode :2015MNRAS.449.2413B . doi :10.1093 / mnras / stv445 .