合适的参数威布尔分布

37视图(30天)

显示旧的评论

Kalle 2023年5月4日8点

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1957859-fit-parameter-to-weibull-distribution

编辑: Kalle2023年5月4日八点二十分

你好,我的伤口为什么我定制3-parametric威布尔分布函数不能符合这一特定数据但2-parametric威布尔分布函数。

                          % %输入数据
                         
                          = (5.05002673467000 4.03148539861043 5.39815092086792 5.80624039967855 6.40216302871704 5.76023737589518 4.58155600229899 5.00490283966065 6.05405394236247 3.56419595082601 4.50625197092692 6.45612986882528 5.02723805109660 5.93072589238485 5.20658095677694 6.97059233983358 6.17937358220418 4.19038120905558 5.45702393849691 4.59787972768148 2.83733654022217 5.06135972340902 4.45322751998902 3.77954316139221 - 6.04727951685588 5.98482131958008 - 4.78870089848836 5.37800931930542 - 5.78453779220581);
                         
                          = ';
                         
                          % %威布尔的设置
                         
                          custompdf = @ (x, a, b, c) (x > c) * (b / a)。*(((得到)/ a) ^ (b - 1))。* exp(-((得到)/ a)。^ b);
                         
                          customcdf = @ (x, a, b, c) (x > c)。* 1-exp(-((得到)/ a)。^ b);
                         
                          选择= statset (“麦克斯特”1 e5,“MaxFunEvals”1 e5,“FunValCheck”,“关闭”);
                         
                          params =大中型企业(,“pdf”custompdf,“开始”,(意思是(A)性病(A)分钟(A)],“选项”选择,下界的[0 0负),“UpperBound”,正正敏(A)]);
                         
                          一个= params (1,1);
                         
                          b = params (1、2);
                         
                          c = params (1、3);
                         
                          x = [c + eps (c): 0.1: max () * 4.5);
                         
                          % %分位数
                         
                          p = wblfit(一个);
                         
                          [nlogl, pcov] = wbllike (p);
                         
                          [quantile_Wbl2, q95lo q95up] = wblinv (0.9991666, p (1), p (2), pcov);
                         
                          Q =分位数(customcdf (x, a, b, c), 0.9991666);
                         
                          quantile_Wbl3 = fsolve (@ (x) customcdf (x, a, b, c) - 0.9991666, quantile_Wbl2);
                         
                          % %数据
                         
                          情节(x, customcdf (x, a, b, c),“线宽”2);
                         
                          持有在
                         
                          情节(x, wblcdf (x, p(1),(2)页));
                         
                          持有在
                         
                          cdfplot(一个);
                         
                          传奇(“WBL-3-Parameter”,“WBL-2-Parameter”,“daten”)
                         
                          持有从
                         
                          % % Anpassungstest
                         
                          % kolmogorov-smirnov Anpassungstest
                         
                          [h_wbl3, p_wbl3] =键糟(A, (customcdf (A, A, b, c)], 0.05, 0);
                         
                          [h_wbl2, p_wbl2] =键糟(A, (wblcdf p (p (1), (2))), 0.05, 0);