Không phải là một câu trả lời đầy đủ, xin lỗi, nhưng một vài ý tưởng (để mong nhận xét). Lưu ý rằng những gì bạn có là một sản phẩm của IID biến ngẫu nhiên, nơi K là một biến ngẫu nhiên (rv) với một bản phân phối Poisson với tham số λ . Điều đó có thể được sử dụng cho một "kiểm tra vệ sinh" khác, một mô phỏng (sử dụng số mũ của tỷ lệ 1):K+ 1Kλ
set.seed(7*11*13)
N <- 1000000
prods <- rep(0, N)
ks <- rpois(N, 1)+1
for (i in 1:N) {
k <- ks[i]
prods[i] <- prod( rexp(k, 1))
}
qqplot( qexp(ppoints(N)), prods)
Kết quả qqplot
(không được hiển thị ở đây) nằm xa một đường thẳng, do đó, điều này không có vẻ là một hàm mũ của tỷ lệ 1. Giá trị trung bình là đúng, phương sai lớn, có đuôi bên phải dài hơn nhiều so với số mũ. Những gì có thể được thực hiện trên lý thuyết? Biến đổi Mellin https://en.wikipedia.org/wiki/Mellin_transform được điều chỉnh phù hợp với các sản phẩm của các biến ngẫu nhiên độc lập. Tôi sẽ chỉ tính cho mũ với tốc độ 1. Mellin biến đổi của sau đó là
M 1 ( s ) = E V s 0 = ∫ ∞ 0 x s e - xV0
nên biến đổi Mellin của sản phẩm của k + 1 iid lũy thừa là
M k + 1 ( s ) = Γ ( s + 1 ) k + 1
Vì K có phân phối poisson với tham số λ , biến đổi Mellin của sản phẩm ngẫu nhiên có số nhân ngẫu nhiên K + 1 , là
M ( s ) = E M K
M1( s ) = EVS0= ∫∞0xSe- xdx = Γ ( s + 1 )
k + 1Mk + 1( s ) = Γ ( s + 1 )k + 1
KλK+ 1
nhưng tôi không thể tìm thấy một nghịch đảo của biến đổi này. Nhưng lưu ý rằng nếu
Xlà biến ngẫu nhiên không âm với biến đổi Mellin
MX(t), thì xác định
Y=logXchúng ta thấy rằng
KY(t)=EetM( s ) = EMK+ 1( s ) = EΓ ( s + 1 )K+ 1= Γ ( s + 1 ) e- λΣk = 0∞λkk !Γ ( s + 1 )k= e- λΓ ( s + 1 ) eλ Γ ( s + 1 )
XMX( t )Y= nhật kýX
vì vậy biến đổi Mellin của
Xlà hàm tạo mô men của logarit
Ycủa nó. Vì vậy, bằng cách sử dụng chúng ta có thể xấp xỉ phân phối
Xvới các phương pháp xấp xỉ yên ngựa,
Làm thế nào gần đúng yên ngựa hoạt động? và tìm kiếm trang web này.
KY( t ) = Eet Y= Eet logX= Eeđăng nhập( Xt)= EXt= MX( y)
XYX