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&