how to fit exponential distribution function on data?

Question

Mos_bad 2020-1-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/501971-how-to-fit-exponential-distribution-function-on-data

回答： Image Analyst 2020-1-26

The vector m follows the truncated exponential equation (F_M) and it is shown by solid black line in figure. I intend to fit an exponential distribution function to data and find the parameter lambda (1/mean). Even though I've used fitdist(x,distname), the fitted exp. dist. shown in dashed line which is way different from the data. here is the code:

M_min=4.5; M_max=8.0; 
m=M_min:0.0001:M_max;
a=4.56; b=1.0; 
alpha=a*log(10);beta=b*log(10);
nu=exp(alpha-beta*M_min);
F_M=(1-exp(-beta*(m-M_min))) / (1-(exp(-beta*(M_max-M_min)))); % CDF of Mag.
pd = fitdist(m','Exponential');
figure(1); plot(m,1-F_M,'-','linewidth',2); 
hold on; plot(m,1-cdf(pd,m),'--'); 
legend('data','fitted dist')

2 个评论
显示无隐藏无

Walter Roberson 2020-1-25

a and b are not defined in the third line.

Mos_bad 2020-1-25

just editted. thanks for pointing that out.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2020-1-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/501971-how-to-fit-exponential-distribution-function-on-data#answer_412031

在 MATLAB Online 中打开

You do not have an exponential distribution. (1 minus an exponential) is not an exponential.

On the other hand if you fit using the equation

a*exp(-b*x)+c

instead of

a*exp(-b*x)

then you get pretty much a perfect fit.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Mos_bad 2020-1-26

All I want to do is to devide vertical axis to 1000 intervals and pick a random value of magnitude (horizental axis) at each interval. Kind of Latin hypercube sampling.

请先登录，再进行评论。