Seite 1 von 1

Total least square in python?

Verfasst: Sonntag 17. Juli 2011, 01:02
von p90
Hi,

folgendes Problem:
Ich habe messe y in Abhängigkeit von x.
Sowohl x als auch y sind fehlerbehaftet.
Ich weiß, das es einen linearen Zusammenhang geben sollte.
Nur wie fitte ich das ganze?
Das was ich bisher gelesen habe sagt, das ich einen Total least square fit (Achtung! Total nicht vergessen, least square bietet scipy ja an, total least square AFAIK nicht) durchführen muss.
Leider habe ich keine Ahnung wie das auszusehen hat und die einzigen Hinweise führen mich immer wieder zum selben Fortran Code.
Oder in python Code:

Code: Alles auswählen

x = [...,...,...]
y = [...,...,...] 
xerr = [...,...,...]
yerr = [...,...,...]
modelfunc = lamda x: m*x+b
Gesucht ist wie gesagt m und b, beides sind Skalare btw.
Hoffe jemand kann mir helfen.

bis später


p90

Re: Total least square in python?

Verfasst: Sonntag 17. Juli 2011, 05:31
von gkuhl
Hi,

wenn ich mich richtig erinnere, ist ein "Total Least Square" Problem am Ende nur ein Eigenwertproblem. Dafür gibt es bei SciPy dann Funktionen.

Grüße
Gerrit

Re: Total least square in python?

Verfasst: Sonntag 17. Juli 2011, 23:11
von CM
Oder wäre scipys ORD-Paket etwas für Dich?

Gruß,
Christian

Re: Total least square in python?

Verfasst: Dienstag 19. Juli 2011, 01:47
von p90
@CM
Genau das ist das was ich gesucht habe!
Danke!