Convolution of two log normal distributions

Question

Joshua Woodard 2019-9-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/480815-convolution-of-two-log-normal-distributions

编辑： Bruno Luong 2019-9-18

Greetings. I am trying to do the convolutions of two lognormal distributions however, I am getting errors. I started to question my method but I cannot find a mistake in my script. Is there a better way? I seem to be getting unexpected numerical values. I would expect like any CDF to approach 1, but this one is not.

%% Convolution of two LogNormal Distributions.

clc
clear all
format long
muX = 9.7224; % For Ni Distribution (X-domain)
sigmaX = 0.3332; % For Ni Distribution (X-domain)
muY = 8.6878; % For Np Distribution (Y-domain)
sigmaY = 0.2454; % For Np Distribution (Y-domain)
t = inf; % Cycles input (would expect an answer of 1 with t = inf.)
fun = @(y,x) exp(-0.5.*((log(x) - muX).^2)./(sigmaX.^2))./(x.*sigmaX.*sqrt(2.*pi)) .* exp(-0.5.*((log(y) - muY).^2)./(sigmaY.^2))./(y.*sigmaY.*sqrt(2.*pi));
P = integral2(fun,-inf,t - 'x',-inf,inf,'RelTol',1e-12,'AbsTol',1e-12)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jeff Miller 2019-9-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/480815-convolution-of-two-log-normal-distributions#answer_392406

在 MATLAB Online 中打开

LognormalsConv.jpg

Convolutions are pretty easy to do in Cupid. For example, the following code gives the attached figure

muX = 9.7224; % For Ni Distribution (X-domain)
sigmaX = 0.3332; % For Ni Distribution (X-domain)
muY = 8.6878; % For Np Distribution (Y-domain)
sigmaY = 0.2454; % For Np Distribution (Y-domain)
conv = Convolution(Lognormal(muX,sigmaX),Lognormal(muY,sigmaY));
conv.PlotDens

This might be handy if you also want to try other distributions, try fitting data, etc.