NIST非线性拟合基准Misra1a

Question

我一直用一些代码来解决这个问题，使用Levenberg-Marquardt求解器 但我对这里给出的内容感到困惑 http://www.itl.nist.gov/div898/strd/nls/data/LINKS/DATA/Misra1a.dat

当我插入认证值并查看预测值时，它们看起来与实际Y值完全不同......我肯定在这里做错了什么？

Y= b1*(1-exp[-b2*X])
b1=238.94212918
b2=0.00055015643181


X       Y       Y-estimate  
10.07   77.6    1.32E+00
14.73   114.9   1.93E+00
17.94   141.1   2.35E+00
23.93   190.8   3.13E+00
29.61   239.9   3.86E+00
35.18   289     4.58E+00
40.02   332.8   5.20E+00
44.82   378.4   5.82E+00
50.76   434.8   6.58E+00
55.05   477.3   7.13E+00
61.01   536.8   7.89E+00
66.4     593.1   8.57E+00
75.47   689.1   9.72E+00
81.78   760     1.05E+01

我想也许基数是10并且尝试了力量（10，...）而不是exp但是这似乎不是问题。

Answer 1

您的X和Y列已被交换（实际上，再次在他们的桌子上看密切）：

238.94212918*(1-e(-0.00055015643181*77.6))
9.98626636447322174420

238.94212918*(1-e(-0.00055015643181*10.07))
1.32009728485679509663

顺便提一下，他们的模型还包括最终的+ e术语。