= (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)]);
x = [c + eps (c): 0.1: max () * 4.5);
[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)页));
传奇(“WBL-3-Parameter”,“WBL-2-Parameter”,“daten”)
[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);