《カルマンフィルター》

提供: ORWiki
2007年7月19日 (木) 22:42時点におけるOrsjwiki (トーク | 投稿記録)による版 ("《カルマンフィルター》" を保護しました。 [edit=sysop:move=sysop])
ナビゲーションに移動 検索に移動

【かるまんふぃるたー (Kalman filter)】

Sk-0098-b-h-05-1.png
図1:状態空間モデル

スタイル検討


 図1に示すガウス白色雑音を受ける離散時間線形確率システムを考える.



ただし, はn次元状態ベクトル, はp次元観測ベクトル,, および次元ガウス白色雑音ベクトルで, 平均値は0, 共分散行列は



であるとする. ただし, , は非負定値対称行列, はクロネッカーの記号である. は一般に時間の関数であってもよいが,簡単のために添字は省略している. また初期値は平均値,共分散行列のガウス確率ベクトルであり, 雑音とは無相関であるとする.

 カルマンフィルタは, 観測データに基づいて,状態の最小分散推定値(すなわちの条件付き期待値)



を逐次的に計算するアルゴリズムである. にしたがってそれぞれ, 予測, 濾波, 平滑という.以下では, 状態ベクトルの予測推定値, 濾波推定値をと表し, それぞれの推定誤差共分散行列を次のようにおく.



 1960-61年にカルマン(R. E. Kalman)と ビュシー(R. S. Bucy)は式(1), (2)の状態空間モデルに対して, カルマンフィルターと呼ばれる以下のようなアルゴリズムを提案した [1, 2, 3, 4, 5, 6].


(i) フィルタ方程式

(ii) カルマンゲイン


(iii) 推定誤差共分散行列


Sk-0098-b-h-05-2.png
図2:カルマンフィルターのブロック線図

スタイル検討



 図2にカルマンフィルタのブロック線図を示す.カルマンフィルタは観測値を入力とし, 推定値, を逐次的に出力する線形動的システムであり,はイノベーションと呼ばれている.

 がガウス分布でない場合には, 上のアルゴリズムは状態ベクトルおよび信号の線形最小分散推定値を与えるという意味で, 最適なフィルタである.

 式(6), (7)からを消去して,とおくと, 離散時間リッカチ方程式



を得る. また上式のにおける極限を代数リッカチ方程式という. 制御理論の分野には(代数)リッカチ方程式に関する膨大な研究がある [7].

 白色雑音を受ける非線形確率システムに対しても,その線形化モデルにカルマンフィルタを適用することができるので, カルマンフィルタは航空宇宙工学の分野において飛翔体の軌道推定に威力を発揮した. また状態ベクトルだけてなく, モデルに含まれる未知パラメータを同時に推定する拡張カルマンフィルタも提案されており, カルマンフィルタの応用は時系列の推定を始めとして非常に多くの分野に見られる.



参考文献

[1] R. E. Kalman, "A New Approach to Linear Filtering and Prediction Problem," Transactions of American Society of Mechanical Engineers, Journal of Basic Engineering, 82D (1960), 34-45.

[2] R. E. Kalman and R. S. Bucy, "New Results in Linear Filtering and Prediction Theory," Transactions of American Society of Mechanical Engineers, Journal of Basic Engineering, 83D (1961), 95-108.

[3] B. D. O. Anderson and J. B. Moore, Optimal Filtering, Prentice-Hall, 1979.

[4] M. S. Grewal and A. P. Andrews, Kalman Filtering - Theory and Practice, Prentice-Hall, 1993.

[5] 有本 卓, 『カルマンフィルター』, 産業図書, 1977.

[6] 片山 徹, 『新版応用カルマンフィルタ』, 朝倉書店, 2000.

[7] S. Bittanti, A. J. Laub and J. C. Willems (Eds.), The Riccati Equation, Springer, 1991.