matlab拟合韦布尔分布,Matlab 三参数Weibull分布拟合求解

function  a_b_c = wblthree(x)

% f(x) = b*a^(-b)*(x-c)^(b-1)*exp(-((x-c)/a)^b)

% a ------  尺度参数

% b ------  形状参数

% c ------  位置参数

disp('样本区间及最大值与最小值之比:')

x_range = [min(x) max(x) max(x)/min(x)]

alpha = [0.05]; %置信水平

c = linspace(0,min(x)-1,1000)';

Len_c = length(c);

for i = 1 : Len_c

[a_b(i,:),pci{i}] = wblfit(x-c(i),alpha);

lnL(i,1) = - wbllike([a_b(i,:)],x-c(i));

if a_b(i,2) <= 1

break;

end

end

c = c(1:i);

figure('name','wblthree--参数特征')

[max_lnL,position_max] = max(lnL);

subplot(311)

plot(c,lnL,'r',c(position_max),max_lnL,'rs')

title('c - lnL')

text(c(position_max),max_lnL,num2str(max_lnL));

subplot(312)

plot(c,a_b(:,1),'b&