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，并尝试使用 power(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 术语。