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