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 で、exp の代わりに power(10,...) を試しましたが、それは問題ではないようです。

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項も含まれています。