三参数威布尔分布线性相关法拟合程序

function para = wbl3corrfit(x) % f(x) = b*a^(-b)*(x-c)^(b-1)*exp(-((x-c)/a)^b) % a ------ 尺度参数 % b ------ 形状参数 % c ------ 位置参数 % para = [a,b,c]

function  para = wbl3corrfit（x）
% f（x） = b*a^（-b）*（x-c）^（b-1）*exp（-（（x-c）/a）^b）
% a ------  尺度参数
% b ------  形状参数
% c ------  位置参数
% para = [abc]

if ~isvector（x）
   error（‘wbl3corrfit函数的输入必须为向量‘）;
end

if find（length（x）==size（x））==1
    x = x‘;
end

step = 0.01;
xSort = sort（x）;
R_xSort = 1-（1:length（xSort））/（length（xSort）+1）;
miu=0:s