also die funktion ist

f = a*x^c

und gefittet werden a und c.

Ingesamt sind es 10 einzelne Datenpunkte...

Hallo,
hat einer von euch schon mal mit dem (Scipy) Befehl leastsq() gearbeitet?
Das Fitten eines Leastsquare-Fit laeuft ohne Probleme, aber wuerde mir auch gerne den Fehler/Ungenauigkeit des Fits ausgeben lassen.
Weiss jemand wie das geht?

Danke

Ich habe das jetzt durch rausschreiben in eine neue Datei geloest, aber ich weiss nicht wie ich ihm sage, dass er die letze Zeile loeschen soll, vorallem da alle Datein unterschiedlich lang sind...




fin = open( 'test.dat', "r" )
data_list = fin.readlines()
fin.close()

del data_list[0:2]

fout ...

Hallo,

ich habe das Problem das ich in Python eine ganze Menge an Dateien plotten moechte (mit matplotlib). Die ganzen ascii-Daten fangen aber alle mit 3 Zeilen text an und hoehren mit unterstrichen auf also:


dateinamte.txt
xachse yachse zachse
------------------------------
1 2 3
1 2 3 ...

stimmt das war die falsche Datei copy/paste :oops:

neuer Versuch



import numpy, asciidata, pylab, matplotlib.pyplot
import matplotlib
from pylab import *


matplotlib.rc('text', usetex = True)
gaussian = lambda x: 3*exp(-(30-x)**2/20.)

X, data = loadtxt('test.dat', unpack=True)


plot(data ...

from pylab import *

gaussian = lambda x: 3*exp(-(30-x)**2/20.)

data = gaussian(arange(100))

plot(data)

X = arange(data.size)
x = sum(X*data)/sum(data)
width = sqrt(abs(sum((X-x)**2*data)/sum(data)))

max = data.max()

fit = lambda t : max*exp(-(t-x)**2/(2*width**2))

plot(fit(X))

show()


Das ...

Aber das plottet doch das ganze nur und macht keinen gauss Fit oder?

8.08305e+03 5.18096e-01
8.08310e+03 5.65463e-01
8.08315e+03 6.09921e-01
8.08320e+03 5.58614e-01
8.08325e+03 5.42059e-01
8.08330e+03 5.22219e-01
8.08335e+03 5.47670e-01
8.08340e+03 5.46657e-01
8.08345e+03 4.82227e-01
8.08350e+03 5.09462e-01
8.08355e+03 5.36048e-01
8.08360e+03 4.99334e-01 ...

die Datein sind 10 000 Zeilen lang.
Mein Problem ist das ich so gut wie keine Ahnung von python habe.
Ich dachte man mueste das
data = gaussian(arange(100)) mit sowas wie


f = open('test1.dat', 'r')
data = f.readlines()
f.close oder data1 = asciidata.open('test.dat')




ersetzen, damit er ...

ja das weiss ich. Aber ich habe doch eine Gauss Verteilung. Da bin ich mir ziemlich sicher.
Mein Problem ist einfach, dass Programm so unzuschreiben, dass es meine Datei benutzt und diese Gauss fittet.

Ich bin mir nicht sicher was du mit Struktur meinst.
Ich habe eine Tabelle und das sind Messungen von Spektrallinien. Und ich brauche am Ende von denen die Halbwertsbreite und die 'Flaeche'.
Jedenfalls kann man die Linien gut mit Gauss kurven approximieren.

Es sind einfach zwei Spalten in einer Ascii Tabelle. X und Y. Und dann jeweils durch Leerzeichen getrennt.

Hallo, ich bin gerade dabei Python zu lernen und moechte nun Datenpunkte in einer Ascii Tabelle mit einer Gauss function fitten. Im Scipy cookbook habe ich folgende Anweisung gefunden.

from pylab import *

gaussian = lambda x: 3*exp(-(30-x)**2/20.)

data = gaussian(arange(100))

plot(data)

X ...

super danke !

Mit numpy funktioniert es.

Hallo ,
ich moechte eine for loop prgramieren und dazu eine Schrittweite einsetzen die keine Integer ist. Also
for x in range(6000,8000,0.15):

Bekomme aber immer ein Fehler.
Gibt es dafuer vielleicht eine Loesung?

(Achso: der Fehler liegt wirklich(!) an der 0.15, wenn ich 1 einsetze ...

Hallo,

ich bin neu hier und auch neu bei Python und habe ein Problem, bei dem ich nicht weiterkomme. Ich habe eine Ascii-Tabelle und möchte diese neu binnen.
Dazu habe ich bei http://www.scipy.org/Cookbook/Rebinning
ein Cookbook gefunden, verstehe es leider aber nicht :roll: .


def rebin( a ...