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