Borrar filtros
Borrar filtros

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 visualizaciones (últimos 30 días)
Hi
  3 comentarios
Mayank Agrawal
Mayank Agrawal el 6 de Mzo. de 2020
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

Iniciar sesión para comentar.

Respuesta aceptada

Torsten
Torsten el 26 de Nov. de 2015
Editada: Torsten el 26 de Nov. de 2015
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 comentarios
kader
kader el 27 de Nov. de 2015
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 el 27 de Nov. de 2015
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.

Iniciar sesión para comentar.

Más respuestas (1)

Torsten
Torsten el 26 de Nov. de 2015
Editada: Torsten el 26 de Nov. de 2015
  2 comentarios

Iniciar sesión para comentar.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by