Erklären von Kalman-Filtern in Zustandsraummodellen

10

Welche Schritte sind bei der Verwendung von Kalman-Filtern in Zustandsraummodellen erforderlich?

Ich habe einige verschiedene Formulierungen gesehen, bin mir aber über die Details nicht sicher. Zum Beispiel beginnt Cowpertwait mit diesem Satz von Gleichungen:

yt=Ftθt+vt
θt=Gtθt1+wt

wobei und , sind unsere unbekannten Schätzungen und sind die beobachteten Werte.θ0N(m0,C0),vtN(0,Vt)wtN(0,Wt)θtyt

Cowpertwait definiert die beteiligten Verteilungen (vorherige, Wahrscheinlichkeits- bzw. hintere Verteilung):

θt|Dt1N(at,Rt)
yt|θtN(Ftθt,Vt)
θt|DtN(mt,Ct)

mit

at=Gtmt1,Rt=GtCt1Gt+Wtet=ytft,mt=at+Atetft=Ftat,Qt=FtRtFt+VtAt=RtFtQt1,Ct=RtAtQtAt

Übrigens bedeutet die Verteilung von bei den beobachteten Werten bis . Eine einfachere Notation ist aber ich werde mich an Cowpertwaits Notation halten.θt|Dt1θtyt1θt|t1

Der Autor beschreibt die Vorhersage für in Bezug auf die Erwartungen:yt+1|Dt

E[yt+1|Dt]=E[Ft+1θt+1+vt+1|Dt]=Ft+1E[θt+1|Dt]=Ft+1at+1=ft+1

Soweit ich weiß, sind dies die Schritte. Bitte lassen Sie mich wissen, wenn ein Fehler oder eine Ungenauigkeit vorliegt:

  1. Wir beginnen mit , , wir schätzen einen Wert für unsere Schätzungen .m0C0θ0
  2. Wir sagen einen Wert für voraus . Das sollte gleich was . ist bekannt, da .y1|D0f1F1a1a1a1=G1m0
  3. Sobald wir unsere Vorhersage für , berechnen wir den Fehler .y1|D0e1=y1f1
  4. Der Fehler wird verwendet, um die hintere Verteilung zu berechnen , die und erfordert . wird als gewichtete Summe des vorherigen Mittelwerts und des Fehlers angegeben: .e1θ1|D1m1C1m1a1+A1e1
  5. In der folgenden Iteration beginnen wir mit der Vorhersage von wie in Schritt 1. In diesem Fall ist . Da und die Erwartung von , die wir bereits im vorherigen Schritt berechnet haben, können wir mit der Berechnung des Fehlers fortfahren und der Mittelwert der posterioren Verteilung wie zuvor.y2|D1f2=F2a2a2=G2m1m1θ1|D1e2θ2|D2

Ich denke, die Berechnung der posterioren Verteilung ist das, was manche Leute den Aktualisierungsschritt nennen, und die Verwendung der Erwartung von ist der Vorhersageschritt.θt|Dtyt+1|Dt

Der Kürze halber habe ich die Schritte zur Berechnung der Kovarianzmatrizen weggelassen.

Habe ich etwas vergessen? Kennen Sie einen besseren Weg, dies zu erklären? Ich denke, das ist immer noch etwas chaotisch, also gibt es vielleicht einen klareren Ansatz.

Robert Smith
quelle

Antworten:

3

Ich denke, was Sie sagen, ist richtig, und ich denke nicht, dass es chaotisch ist. Eine Art der Formulierung wäre zu sagen, dass der Kalman-Filter ein Fehlerkorrekturalgorithmus ist, der Vorhersagen im Lichte der Diskrepanzen mit aktuellen Beobachtungen modifiziert. Diese Korrektur wird in Schritt 4) unter Verwendung der Verstärkungsmatrix .At

F. Tusell
quelle
Vielen Dank für Ihre Antwort. Vielleicht ist es richtig, aber ich würde gerne eine detailliertere (und natürlichere) Erklärung dazu lesen. Ich habe Beschreibungen in Büchern und Folien gelesen, aber die meisten sind nicht sehr klar und es gibt geringfügige Unterschiede.
Robert Smith