三参数威布尔分布拟合程序

function a_b_c = wbl3fit(x) % f(x) = b*a^(-b)*(x-c)^(b-1)*exp(-((x-c)/a)^b) % a ------ 尺度参数 % b ------ 形状参数 % c ------ 位置参数

function  a_b_c = wbl3fit（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（0min（x）1000）‘;
Len_c = length（c）-1;
for i = 1 : Len_c
    [a_b（i:）pci{i}] = wblfit（x-c（i）alpha）;
    lnL（i1） = - wbllike（[a_b（i:）]x-c（i））;
    if a_b（i2） <= 1
        break;
    end
end
c = c（1:i）;
figure（‘name‘‘wblthree--参数特征‘）
[max_lnLposition_max] = max（lnL）;
subplot（311）
plot（clnL‘r‘c（position_max）max_lnL‘rs‘）
title（‘c - lnL‘）
text（c（position_max）max_lnLnum2str（max_lnL））;

subplot（312）
plot（ca_b（:1）‘b‘c（position_max）a_b（position_max1）‘bs‘）
text（c（position_max）median（a_b（:1））num2str（a_b（positi