Hallo,
nach der gestrigen Hilfe habe ich ein kleines Programm
gebastelt, aber ich mache Syntax-Fehler denke ich!
Wer kann mir einen Tip zur Hilfe geben?
ClearAll["Global`*"]
<< Statistics`NonlinearFit`
d[1] = 1;
d[2] = 2;
d[3] = 1;
startv = {3216, 675, 1280, -55.69, -2.8, -5.86, -16.6, -4.07, 5.7};
data = {{0, 1, 0, 1, 674.699582`}, {0, 2, 0, 0, 1332.3452`}, {0, 2, 0, 2,
1355.156`}, {0, 1, 1, 1, 1951.535`}, {1, 0, 0, 0, 3216.889192`}, {2, 0,
0, 0, 6322.3868`}, {1, 1, 0, 1, 3874.971832`}, {0, 0, 1, 0,
1280.9125`}, {0, 0, 2, 0, 2550.1`}, {0, 3, 0, 1, 1990.212`}, {0, 4, 0,
0, 2642.773`}, {0, 3, 0, 3, 2041.5759`}}
gv[v_, l_] :=
Sum[omega[i]*(v[i] + d[i]/2), {i, 1, 3}] +
Sum[Sum[x[i, j]*(v[i] + d[i]/2)*(v[j] + d[i]/2), {i, 1, j}], {j, 1, 3}] +
g22 l^2
NonlinearRegress[data,
gv[v_, l_], {v[1], v[2], v[3], l}, {omega[1], omega[2], omega[3], x[1, 1],
x[2, 2], x[3, 3], x[1, 2], x[1, 3], x[2, 3], g22}, InitialGauss -> startv,
Errors -> 1]
Als Fehlermeldung bekomme ich, the model is not numerical at.....
x[1,3] ist durch den data Datensatz nicht bestimmt, wird zwar in
gv[v_,l_] aufgestellt aber weil v[1] und v[2] = null sind,
ist x[1,3], nicht definiert - ich weiss nicht, ob ich es unter den
parametern von NonlinearRegress ueberhaupt auffuehren darf.
b) eine kleine Frage noch zu dem Messwert (letzte Position in jedem data
Vektor) dafuer habe ich auch einen Fehler, wie kann ich diesen in
data eingeben, um ihn fuer die Fehlerrechnung zu benutzen?
Gruesse, Michael
--------------------------------------------------------------------------
Michael Jung Room 1002
Justus Liebig Universitaet Giessen Phone ++49 (0)641 99 34529
Institut fuer Physikalische Chemie Fax ++49 (0)641 99 34579
Heinrich-Buff-Ring 58
35392 Giessen
Germany
email: michael.jung@XXXXXXX.de
|