2010-08-02 90 views
1

我一直在与一些代码玩弄这种使用列文伯格 - 马夸特求解 解决,但我通过什么这里给出的疑惑 http://www.itl.nist.gov/div898/strd/nls/data/LINKS/DATA/Misra1a.datNIST非线性拟合基准Misra1a

当我插上标准值和看看预测值,他们看起来与实际的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,但这似乎并不是问题。

回答

0

XY列交换(真的,去再看看在他们的桌子密切):

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

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

顺便说一下,他们的模型还包括最终+ e项。

+0

是的, ...发现它,并且正准备脱离我的问题 – CrunchedNumber 2010-08-02 10:57:53

+0

误差项e只是最小化的术语,它不是一个参数。 – CrunchedNumber 2010-08-02 10:59:35

+0

啊!我设法忘记了几年的统计课程...... – sarnold 2010-08-02 11:20:40