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.
Oui, ... l'ai trouvé et était sur le point de rayer ma question –
CrunchedNumber
Le terme d'erreur e est juste le terme à minimiser, ce n'est pas un paramètre. – CrunchedNumber
Ah! J'ai réussi à oublier plusieurs années de cours de statistiques ... – sarnold