DMUG-Archiv 2009

Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

Re: Fit läuft nicht

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.


Verweise:
Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

DMUG DMUG-Archiv, http://www.mathematica.ch/archiv.html