Wie lässt sich der stationäre Kalman-Filter-Prädiktor ableiten?

14

In seinem Kapitel über Kalman-Filter heißt es in meinem DSP-Buch, scheinbar aus heiterem Himmel, der stationäre Kalman-Filter für ein System

{x(t+1)=Ax(t)+w(t)y(t)=Cx(t)+v(t)

hat den Prädiktor

x^(t+1|t)=(AAK¯C)x^(t|t1)+AK¯y(t)

und stationäre Zustandsvektorkovarianz und Kalman-Gewinn

ˉ K = ˉ P CT(C ˉ P CT+R)-1

P¯=AP¯ATAP¯CT(CP¯CT+R)1CP¯AT+Q
K¯=P¯CT(CP¯CT+R)1

wobei Q und R die Kovarianzen des Eingangsrauschens w bzw. des Messrauschens v bezeichnen.

Ich kann nicht sehen, wie ich aus dem Minimum-Varianz-Prädiktor zu diesem Ergebnis komme. Könnte es mir jemand erklären oder mich auf eine Ressource verweisen, die den Ausdruck ableitet? Dies ist das zeitvariante Minimum-Varianz-Filter, das ich ableiten kann :

P(t+1|t)=A(P(t|t-1)-P(t|

x^(t+1|t)=(EIN-K(t)C)x^(t|t-1)+K(t)y(t)
K ( t ) = A P ( t | t - 1 ) C T ( C P ( t | t -
P(t+1|t)=EIN(P(t|t-1)-P(t|t-1)CT(CP(t|t-1)CT+R)-1CP(t|t-1))EINT+Q.
K(t)=EINP(t|t-1)CT(CP(t|t-1)CT+R)-1

Ich bin mir nur nicht sicher, wie ich von hier zum stationären Filter weiter oben gehen soll.

Update: Ich kann sehen, dass das Einsetzen von P¯=P(t+1|t)=P(t|t-1) und K(t)=EINK¯ in den zeitvarianten Filter zu führt der stationäre Filter, aber warum mit A multiplizieren EIN? Ist dies nur ein Symptom für eine unglückliche Wahl der Notation, was bedeutet, dass entweder K oder K¯ nicht wirklich den Kalman-Gewinn bedeuten ?

Andreas
quelle
Nein, es ist nicht möglich, den Prädiktor aus den Gleichungen für das System zu "sehen". Ich denke, es wäre besser, wenn Sie ein Lehrbuch über Kalman-Filter lesen, anstatt uns zu bitten, es für Sie abzuleiten (was nur etwas aus einem Lehrbuch herausfordert). Optimales Filtern nach Anderson und Moore könnte ein guter Anfang sein. Es ist in Kapitel 5 abgeleitet, wenn ich mich richtig erinnere.
Lorem Ipsum
@yoda: Danke. Meine Frage war, ob mich jemand auf eine bessere Ressource hinweisen könnte als das Lehrbuch, das mein Kurs empfiehlt. Das ist also eine Antwort.
Andreas
@yoda: Übrigens, für den Fall, dass ich unklar war: Ich bitte nicht um eine Ableitung aus dem Zustandsraumsystem, sondern aus dem Minimum-Varianz-Kalman-Filter. Ich habe die Frage aktualisiert, um zu verdeutlichen, dass ich einen zeitinvarianten Kalman-Filter ableiten kann, nur nicht den stationären.
Andreas
1
Von welchem ​​Text bekommen Sie das Obige? Wenn jemand Zugriff darauf hat, kann dies nützlich sein, damit wir den vollständigen Kontext sehen können.
Jason R

Antworten:

5

Ihre Ableitungen sind korrekt.

K(t)=A ˉ KP¯=P(t|t-1) undK(t)=EINK¯

Ist das deine Verwirrung:

  1. Warum hatten sie in den Kalman Gain und Covariance Matrix Expressions nicht den Term ?t|t-1
  2. Wie kann dies "stationär" sein, wenn Ihre Ableitung zeigt, dass es sich um eine zeitlich variierende Ableitung handelt?

  1. Falsche Wahl der Notation seitens des Buches

Schauen wir uns den Ausdruck: . Die Tatsache, dass eine Funktion von sich selbst ist, zeigt eine rekursive Beziehung. Mit anderen Worten, es verwendet seine vergangenen Werte. Es ist also NICHT für alle Zeitpunkte gleich - es ändert sich bei jeder Iteration.P¯=EINP¯EINT-EINP¯CT(CP¯CT+R)-1CP¯EINT+Q.P¯

  1. Missverständnis des Wortes "stationär".

Als der Autor des Buches "stationär" sagte, meinte er / sie nicht, dass und zu allen Zeitpunkten den gleichen Wert haben. Stattdessen wollte der Autor betonen, dass die Ausdrücke für diese Werte für alle statistischen Realisierungen gleich sind. Stationarität ist ein statistisches Konzept, bei dem die Statistiken des Systems immer gleich sind. Schauen Sie sich die Ausdrücke für und einmal an. Sie hängen nur von \PKP¯K¯

  • Bisherige Werte von sich
  • Übergangsmatrizen und die deterministisch und in Ihrem Fall zeitinvariant sind ( und sind immer gleich)EINCEINC
  • Q. und sind die Rauschkovarianzmatrizen. Diese 2 Matrizen beschreiben die Statistik der Geräusche und sind in allen Realisierungen und Zeitinstanzen gleich.R

Die Kalman-Verstärkung und die Zustands-Kovarianz-Matrix haben für alle Realisierungen dieses Zufallsprozesses den gleichen Wert. ( Randnotiz: Keiner dieser beiden Begriffe hängt von den Messungen ab, . Sie können also vorher berechnet werden. )KPy


Fazit:

Die von Ihnen abgeleiteten "Zeitvarianten" -Gleichungen entsprachen denen im Buch. Abgesehen von den Unterschieden in der Schreibweise gab es bei Ihnen ein leichtes Missverständnis darüber, was sich ändert und was nicht.

ssk08
quelle
1
Ich erinnere mich nicht, welches Problem ich hatte, als ich die Frage stellte, aber jetzt macht es Sinn. Vielen Dank!
Andreas
Ich verstehe das nicht ganz. Wie würden die Gleichungen für ein nicht stationäres Kalman-Filter dann aussehen?
Sandu Ursu