numerical - 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


XY Y 估计  
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术语。