<h4>カルマンフィルタ：状態空間モデルの最適フィルタリング</h4><div class='key-point'><div class='key-point-title'>状態空間モデルとカルマンフィルタ</div><p>カルマンフィルタは状態空間モデルにおける最適な状態推定アルゴリズムです。ランダムウォーク・プラス・ノイズモデルは最も基本的な状態空間モデルの一つで、トレンド成分と観測ノイズを分離します。</p></div><h4>モデル設定の確認</h4><p class='step'><strong>Step 1: 状態空間モデルの構造</strong></p><p><strong>状態方程式（システムモデル）：</strong></p><div class='formula'>$\alpha_{t+1} = \alpha_t + \eta_t
lt;/div><p><strong>観測方程式（観測モデル）：</strong></p><div class='formula'>$y_t = \alpha_t + \varepsilon_t
lt;/div><p><strong>ノイズの仮定：</strong></p><ul><li>$\eta_t \sim N(0, Q = 1)$：状態ノイズ（システムノイズ）</li><li>$\varepsilon_t \sim N(0, R = 4)$：観測ノイズ</li><li>$\eta_t$ と $\varepsilon_s$ は独立（すべての $t, s$）</li></ul><p><strong>初期条件：</strong></p><div class='formula'>$\alpha_1|Y_0 \sim N(a_{1|0}, P_1 = 2)
lt;/div><h4>カルマンフィルタのアルゴリズム</h4><p class='step'><strong>Step 2: 一般的なカルマンフィルタの手順</strong></p><p>カルマンフィルタは予測ステップと更新ステップの2段階で構成されます：</p><p><strong>予測ステップ（Prediction Step）：</strong></p><div class='formula'>$\begin{align}a_{t|t-1} &= T a_{t-1|t-1} \\P_{t|t-1} &= T P_{t-1|t-1} T' + Q\end{align}
lt;/div><p><strong>更新ステップ（Update Step）：</strong></p><div class='formula'>$\begin{align}v_t &= y_t - Z a_{t|t-1} \\F_t &= Z P_{t|t-1} Z' + R \\K_t &= P_{t|t-1} Z' F_t^{-1} \\a_{t|t} &= a_{t|t-1} + K_t v_t \\P_{t|t} &= P_{t|t-1} - K_t Z P_{t|t-1}\end{align}
lt;/div><p class='step'><strong>Step 3: 本モデルでのシステム行列の設定</strong></p><p>ランダムウォーク・プラス・ノイズモデルでは：</p><ul><li>$T = 1$（状態遷移行列）</li><li>$Z = 1$（観測行列）</li><li>$Q = 1$（状態ノイズ分散）</li><li>$R = 4$（観測ノイズ分散）</li></ul><h4>$t=1$でのカルマンフィルタ計算</h4><p class='step'><strong>Step 4: 予測ステップの実行</strong></p><p><strong>予測分散の計算：</strong></p><p>初期状態（$t=0$）から$t=1$への予測：</p><div class='formula'>$P_{1|0} = T P_{0|0} T' + Q = 1 \times P_1 \times 1 + Q = 2 + 1 = 3
lt;/div><p>ここで、$P_{0|0} = P_1 = 2$（初期値）を使用しました。</p><div class='key-point'><div class='key-point-title'>予測分散の意味</div><p>予測分散$P_{1|0} = 3$は、観測データ$y_1$を見る前の時点での状態$\alpha_1$の推定不確実性を表します。初期不確実性（2）に状態ノイズ（1）が加わり、不確実性が増大しています。</p></div><p class='step'><strong>Step 5: 更新ステップの実行</strong></p><p><strong>予測誤差分散の計算：</strong></p><div class='formula'>$F_1 = Z P_{1|0} Z' + R = 1 \times 3 \times 1 + 4 = 3 + 4 = 7
lt;/div><p><strong>カルマンゲインの計算：</strong></p><div class='formula'>$K_1 = P_{1|0} Z' F_1^{-1} = 3 \times 1 \times \frac{1}{7} = \frac{3}{7}
lt;/div><p><strong>更新後分散の計算：</strong></p><div class='formula'>$P_{1|1} = P_{1|0} - K_1 Z P_{1|0} = 3 - \frac{3}{7} \times 1 \times 3 = 3 - \frac{9}{7} = \frac{21-9}{7} = \frac{12}{7}
lt;/div><div class='key-point'><div class='key-point-title'>カルマンフィルタの情報統合メカニズム</div><p>更新後分散の減少：$P_{1|1} = 1.2 < P_{1|0} = 3
lt;/p><p>これは観測データ$y_1$が新しい情報をもたらし、状態の推定精度が向上したことを示します。カルマンゲイン$K_1 = \frac{3}{7} \approx 0.43$は、予測値と観測値の重み付けバランスを表しています。</p></div><h4>カルマンフィルタの理論的基礎</h4><p class='step'><strong>Step 7: ベイズ推論としての解釈</strong></p><p>カルマンフィルタはベイズ推論の特殊ケースとして理解できます：</p><p><strong>事前分布：</strong></p><div class='formula'>$\alpha_1|Y_0 \sim N(a_{1|0}, P_{1|0})
lt;/div><p><strong>尤度：</strong></p><div class='formula'>$y_1|\alpha_1 \sim N(\alpha_1, R)
lt;/div><p><strong>事後分布：</strong></p><div class='formula'>$\alpha_1|y_1, Y_0 \sim N(a_{1|1}, P_{1|1})
lt;/div><p>ガウシアン分布の共役性により、事後分布も正規分布となり、平均と分散が解析的に計算できます。</p><h4>拡張カルマンフィルタと粒子フィルタ</h4><p class='step'><strong>Step 8: 非線形・非ガウシアン拡張</strong></p><p><strong>1. 拡張カルマンフィルタ（EKF）</strong></p><p>非線形状態空間モデル：</p><div class='formula'>$\begin{align}\alpha_{t+1} &= f(\alpha_t) + \eta_t \\y_t &= h(\alpha_t) + \varepsilon_t\end{align}
lt;/div><p>線形化による近似：</p><div class='formula'>$\begin{align}F_t &= \frac{\partial f}{\partial \alpha}\bigg|_{\alpha_{t|t}} \\H_t &= \frac{\partial h}{\partial \alpha}\bigg|_{\alpha_{t|t-1}}\end{align}
lt;/div><p><strong>2. 無香料カルマンフィルタ（UKF）</strong></p><p>シグマポイント法による高精度近似：</p><ul><li>決定論的サンプリング</li><li>高次のモーメントマッチング</li><li>EKFよりも精度が高い</li></ul><p><strong>3. 粒子フィルタ</strong></p><p>モンテカルロ法による完全非線形フィルタ：</p><div class='formula'>$p(\alpha_t|y_{1:t}) \approx \sum_{i=1}^{N} w_t^{(i)} \delta(\alpha_t - \alpha_t^{(i)})
lt;/div><table style='width:100%; border-collapse: collapse; margin: 1em 0;'><tr style='background-color: #f5f5f5;'><th style='border: 1px solid #ddd; padding: 8px;'>手法</th><th style='border: 1px solid #ddd; padding: 8px;'>適用範囲</th><th style='border: 1px solid #ddd; padding: 8px;'>利点</th><th style='border: 1px solid #ddd; padding: 8px;'>欠点</th></tr><tr><td style='border: 1px solid #ddd; padding: 8px;'>カルマンフィルタ</td><td style='border: 1px solid #ddd; padding: 8px;'>線形ガウシアン</td><td style='border: 1px solid #ddd; padding: 8px;'>最適、高速</td><td style='border: 1px solid #ddd; padding: 8px;'>制限的</td></tr><tr><td style='border: 1px solid #ddd; padding: 8px;'>拡張カルマンフィルタ</td><td style='border: 1px solid #ddd; padding: 8px;'>弱非線形</td><td style='border: 1px solid #ddd; padding: 8px;'>計算効率良</td><td style='border: 1px solid #ddd; padding: 8px;'>線形化誤差</td></tr><tr><td style='border: 1px solid #ddd; padding: 8px;'>無香料カルマンフィルタ</td><td style='border: 1px solid #ddd; padding: 8px;'>中程度非線形</td><td style='border: 1px solid #ddd; padding: 8px;'>高精度、微分不要</td><td style='border: 1px solid #ddd; padding: 8px;'>計算コスト増</td></tr><tr><td style='border: 1px solid #ddd; padding: 8px;'>粒子フィルタ</td><td style='border: 1px solid #ddd; padding: 8px;'>一般的</td><td style='border: 1px solid #ddd; padding: 8px;'>非常に柔軟</td><td style='border: 1px solid #ddd; padding: 8px;'>計算コスト高</td></tr></table>