Hallo Peter,
Kann mir bitte jemand sagen, warum der Fit nicht läuft?
Fit ist fit, ein ähnliches Problem wurde neulich behandelt, also
In[44]:= Clear[sol]
sol = FindFit[
SetPrecision[data, 20],
\[Nu] PDF[NormalDistribution[\[Mu], \[Sigma]], x],
{\[Mu], \[Sigma], \[Nu]}, x]
Out[45]= {\[Mu] -> 3602.2150727800663179, \[Sigma] ->
8.2168928566957955527, \[Nu] -> 1.7970073213172530926*10^6}
man muss einen Normierungsfaktor zulassen. Der Ansatz ist trotzdem
unpassend, weil bei der Normalverteilung eine grössere Standardabweichung
unweigerlich einen kleineren Maximalwert nach sich zieht. Der Fit ist
daher grotesk schlecht; mit einem polynomialen Ansatz ist man schon ohne
grosses Probieren besser bedient
In[74]:= Clear[sol]
sol = FindFit[SetPrecision[data, 10],
Subscript[a, 0] + Subscript[a, 1] x + Subscript[a, 2 ] x^2 +
Subscript[a, 3] x^3, {Subscript[a, 0], Subscript[a, 1], Subscript[
a, 2], Subscript[a, 3]}, x]
Out[75]= {Subscript[a, 0] -> -117.9655198,
Subscript[a, 1] -> -0.4307277147, Subscript[a, 2] -> 0.0003611036317,
Subscript[a, 3] -> -2.530333882*10^-8}
Show[{ListPlot[data, DisplayFunction -> Identity],
Plot[(Subscript[a, 0] + Subscript[a, 1] x + Subscript[a, 2 ] x^2 +
Subscript[a, 3] x^3) /. sol, {x, 0, 13500},
DisplayFunction -> Identity]},
DisplayFunction -> $DisplayFunction]
Gruss
Udo.
|