How can I Calculate the PDF and CDF of a product of two i.i.d exponentially distributed random variables with mean a and b respectively

22 ビュー (過去 30 日間)
Hi
  3 件のコメント
Torsten
Torsten 2018 年 6 月 7 日
The links from below are still valid.
Best wishes
Torsten.
Mayank Agrawal
Mayank Agrawal 2020 年 3 月 6 日
hello sir
I have a similar question but a little advanced one where I want to plot pdf and cdf of Z = (x*y)/(y+b) where 'b' is a positive real constant and x & y have same exponential distributions.
Can you please help me out how to plot the pdf and cdf of 'z' in MATLAB?
CDF of Z: P(Z<=z) = 1 - (e^(-x/lambda1)/lambda2)*(sqrt(4*b*z*lambda2/lambda1))*K1(sqrt(4*b*z/(lambda1*lambda2)))
so we differentiated this cdf w.r.t. z using chain rule and found the pdf which had 3 terms in addition.
we are not able to plot the pdf and cdf of z in MATLAB. please do help us out.
Thanks

サインインしてコメントする。

採用された回答

Torsten
Torsten 2015 年 11 月 26 日
編集済み: Torsten 2015 年 11 月 26 日
You get
F(x)=1-2*sqrt((lambda1)*(lambda2)*x)*K1(2*sqrt((lambda1)*(lambda2)*x))
f(x)=2*(lambda1)*(lambda2)*K0(2*sqrt((lambda1)*(lambda2)*x))
with
lambda1 = 1/m
lambda2 = 1/n
K0, K1 : Modified Bessel functions of the second kind
Best wishes
Torsten.
  2 件のコメント
kader
kader 2015 年 11 月 27 日
Dear Torsten, I really grateful to you. Thanks a lot for this answer. In MatLab, the function is like as follows. For K1, it is besselK(1,x). Right? Though I got the expected result from the mentioned changes in MatLab. So, for K0 it will be besselK(0,x)? F(x)=1-2*sqrt((lambda1)*(lambda2)*x)*besselk(1,2*sqrt((lambda1)*(lambda2)*x))
Torsten
Torsten 2015 年 11 月 27 日
Yes, in MATLAB notation it's
f=@(x)2*(lambda1)*(lambda2)*besselk(0,2*sqrt((lambda1)*(lambda2)*x))
F=@(x)1-2*sqrt((lambda1)*(lambda2)*x)*besselk(1,2*sqrt((lambda1)*(lambda2)*x))
Best wishes
Torsten.

サインインしてコメントする。

その他の回答 (1 件)

Torsten
Torsten 2015 年 11 月 26 日
編集済み: Torsten 2015 年 11 月 26 日
  2 件のコメント
kader
kader 2015 年 11 月 26 日
But that is for the same rate parameter for two random variable!!!

サインインしてコメントする。

タグ

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by