2010-08-02 20 views
1

Je suis en train de jouer avec des codes pour résoudre ce en utilisant un solveur Levenberg-Marquardt mais je suis intrigué par ce qui est donné ici http://www.itl.nist.gov/div898/strd/nls/data/LINKS/DATA/Misra1a.datNIST de référence de montage non linéaire Misra1a

Lorsque je branche les valeurs certifiées et Regardez les valeurs prédites, elles sont très différentes des valeurs Y réelles ... sûrement je fais quelque chose de mal ici - quoi?

 
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 

Je pensais que peut-être la base était de 10 et de puissance essayé (10, ...) au lieu de exp, mais cela ne semble pas être le problème.

Répondre

0

Vos X et Y les colonnes sont permutées (vraiment, aller chercher près à leur table à nouveau):

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

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

Soit dit en passant, leur modèle comprend aussi un terme + e final.

+0

Oui, ... l'ai trouvé et était sur le point de rayer ma question – CrunchedNumber

+0

Le terme d'erreur e est juste le terme à minimiser, ce n'est pas un paramètre. – CrunchedNumber

+0

Ah! J'ai réussi à oublier plusieurs années de cours de statistiques ... – sarnold