定常カルマンフィルター予測子の導出方法


14

カルマンフィルターに関する章で、私のDSPの本は、一見すると一見、一見、システムの定常カルマンフィルターであると述べています。

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

予測子があります

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

および定常状態ベクトル共分散とカルマンゲイン

ˉ K = ˉ P CTC ˉ P CT+R-1

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

ここで、とはそれぞれ入力ノイズと測定ノイズの共分散を示します。R w vQRwv

最小分散予測子からこれに到達する方法がわかりません。誰かがそれを説明してもらえますか、または表現を導き出すリソースを教えてくれますか?これは、時変最小分散フィルターであり、これを導出できます。

PT+1|T=APt|t1Pt|

x^(t+1|t)=(AK(t)C)x^(t|t1)+K(t)y(t)
K t = A P t | t 1 C TC P t | t
P(t+1|t)=A(P(t|t1)P(t|t1)CT(CP(t|t1)CT+R)1CP(t|t1))AT+Q
K(t)=AP(t|t1)CT(CP(t|t1)CT+R)1

ここから上記の固定フィルターにどのように移動するかについては確信がありません。

更新:時変フィルターにおよびを代入すると、静止フィルターですが、なぜ乗算するですか?これは、またはいずれかが実際にカルマンゲインを示していないことを意味する、残念な表記法の選択の症状ですか?KT=A ˉ K AK ˉ KP¯=P(t+1|t)=P(t|t1)K(t)=AK¯AKK¯


いいえ、システムの方程式から予測子を「見る」ことはできません。カルマンフィルターに関する教科書を読むように頼むよりも良いと思います(教科書から何かを逆流するだけです)。アンダーソンとムーアによる最適なフィルタリングは、開始するのに適した場所です。正しく覚えていれば、第5章に由来しています。
ロレムイプサム

@yoda:ありがとう。私の質問は、誰かが私のコースで推奨している教科書よりも優れたリソースを教えてくれないかということでした。それが答えです。
アンドレアス

@yoda:ちなみに、私が不明な場合は、状態空間システムからの派生ではなく、最小分散カルマンフィルターからの派生を求めています。質問を更新して、静止フィルターではなく、時不変カルマンフィルターを導出できることを明確にしました。
アンドレアス

1
どのテキストから上記を取得していますか?誰かがそれにアクセスできれば、完全なコンテキストを見ることができるので便利かもしれません。
ジェイソンR

回答:


5

あなたの派生は正しいです。

P¯=P(t|t1)およびK(t)=AK¯

これはあなたの混乱ですか:

  1. なぜカルマンゲインおよび共分散行列式にという用語がないのですか?t|t1
  2. あなたの派生がそれが時変であることを示すとき、これはどのように「静止」することができますか?

  1. 本の一部の表記法の悪い選択

表現で見てみましょう: 。自体が関数であるという事実は、再帰的な関係を示しています。つまり、過去の値を使用します。そのため、すべての時点で同じではありません -反復ごとに変化します。 ˉ PP¯=AP¯ATAP¯CT(CP¯CT+R)1CP¯AT+QP¯

  1. 「静止」という言葉の誤解。

本の著者が「静止した」と言ったとき、彼/彼女はとがすべての瞬間で同じ値を持っていることを意味しませんでした。代わりに、著者は、値に対するそれらの表現がすべての統計的実現に対して同じであることを強調したかった。定常性は、システムの統計が常に同じであることを意味する統計概念です。ともう一度見てください。\にのみ依存しますK ˉ P ˉ KPKP¯K¯

  • 自分自身の以前の値
  • 決定論的であり、あなたの場合は時間不変である遷移行列と( とは常に同じです)C A CACAC
  • RQノイズ共分散行列であるおよびこれらの2つのマトリックスは、ノイズの統計を記述し、すべての実現および時間インスタンスで同じです。R

カルマンゲインおよび状態共分散行列は、このランダムプロセスのすべての実現に対して同じ値を持ちます。(補足:これら2つの項はいずれも測定値依存しないため、事前に計算できます。P yKPy


結論:

あなたが導き出した「時変」方程式は、本の方程式と同等でした。その上、表記上の違いは、何が変化し、何が変化しないかについて、あなたの側にわずかな誤解がありました。


1
私が質問したときに私が抱えていた問題が何であったか覚えていませんが、今では理にかなっています。ありがとう!
アンドレアス

これはよくわかりません。その場合、非定常カルマンフィルターの方程式はどのようになりますか?
サンドゥウルス
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.