这是X[1 0.962756100000000 0.940400000000000 0.810055900000000 0.638733700000000 0.564200000000000 0.495344000000000 0.456200000000000 0.374300000000000 0.324000000000000 0.288600000000000 0.247672300000000]
这是Y[0 7.64216000000000 84.3694500000000 730.794300000000 943.654000000000 1203.54400000000 1489.30400000000 1814.87200000000 2541.68100000000 3399.59500000000 4385.21900000000 5469.34100000000]
想拟合这个函数,
y=a*[(((x-b)/(1-b))^(-1/m))-1]^(1-m)
本来用的是这几句代码,但是报错
fun = @(p,x)p(1).*((((x-p(2))./(1-p(2)))^(-1/p(3)))-1)^(1-p(3));
x0 = [100;0;0]; %% 初始值
p = lsqcurvefit(fun,x0,x,y); |
|