[論文レビュー] On the stability and the uniform propagation of chaos properties of Ensemble Kalman-Bucy filters

Pierre Del Moral, Julian Tugaut|arXiv (Cornell University)|May 30, 2016

Meteorological Phenomena and Simulations参考文献 65被引用数 63

ひとこと要約

本稿では、新規の関数不等式を導入し、Foster-Lyapunov技法、カップリング、スペクトル解析を組み合わせることで、アンサンブルカルマン＝ブチャトフィルタの初めての一様安定性および混沌の伝播に関する結果を確立した。任意の時間スケールにわたる一様L²平均誤差推定を証明し、高次元で悪条件な系におけるフィルタの収束性および安定性の必要十分条件を同定した。

ABSTRACT

The Ensemble Kalman filter is a sophisticated and powerful data assimilation method for filtering high dimensional problems arising in fluid mechanics and geophysical sciences. This Monte Carlo method can be interpreted as a mean-field McKean-Vlasov type particle interpretation of the Kalman-Bucy diffusions. In contrast to more conventional particle filters and nonlinear Markov processes these models are designed in terms of a diffusion process with a diffusion matrix that depends on particle covariance matrices. Besides some recent advances on the stability of nonlinear Langevin type diffusions with drift interactions, the long-time behaviour of models with interacting diffusion matrices and conditional distribution interaction functions has never been discussed in the literature. One of the main contributions of the article is to initiate the study of this new class of models The article presents a series of new functional inequalities to quantify the stability of these nonlinear diffusion processes. In the same vein, despite some recent contributions on the convergence of the Ensemble Kalman filter when the number of sample tends to infinity very little is known on stability and the long-time behaviour of these mean-field interacting type particle filters. The second contribution of this article is to provide uniform propagation of chaos properties as well as Lp-mean error estimates w.r.t. to the time horizon. Our regularity condition is also shown to be sufficient and necessary for the uniform convergence of the Ensemble Kalman filter. The stochastic analysis developed in this article is based on an original combination of functional inequalities and Foster-Lyapunov techniques with coupling, martingale techniques, random matrices and spectral analysis theory.

研究の動機と目的

• 高次元で悪条件な系におけるアンサンブルカルマン＝ブチャトフィルタの長期的安定性および収束性に関する理論的理解の不足に対処すること。
• 任意の時間スケールにわたるサンプル平均および共分散行列に対する一様L²平均誤差推定を確立すること。
• 共分散に依存する拡散行列を有する非線形拡散過程の安定性を分析すること、これは以前に十分に研究されていなかったクラスである。
• 粒子数が無限に近づく際のEnKFの一様収束性のための必要十分な正則性条件を導出すること。
• ジェンセン型不等式およびスペクトル解析を用いて、カルマン＝ブチャトフィルタおよびその関連するリカッチ方程式の安定性を定量化すること。

提案手法

• 共分散に依存する拡散行列を有する非線形McKean-Vlasov型拡散過程の安定性を定量化するための新しい関数不等式のクラスを構築する。
• Foster-Lyapunov技法およびカップリング手法を適用し、Wasserstein距離および相対エントロピーにおける指数的安定性を分析する。
• マルティングール技法、確率的行列理論、およびスペクトル解析を組み合わせて、相互作用粒子系の長期的挙動を研究する。
• カルマン＝ブチャトフィルタの拡散過程の平均場粒子解釈を用い、EnKFを非線形相互作用粒子系としてモデル化する。
• 時間経過に伴う標本統計量（平均および共分散）の真値への収束を分析することで、一様L²平均誤差推定を導出する。
• システムが安定を保つための許容可能なフラクチュエーション行列Qの集合を特定し、不安定化が生じる発散集合を同定する。

実験結果

リサーチクエスチョン

• RQ1 アンサンブルカルマン＝ブチャトフィルタが任意の時間スケールにわたって一様安定性を示す条件は何か？
• RQ2 EnKF粒子系における混沌の伝播性を時間に一様に定量化する方法は何か？
• RQ3 共分散に依存する拡散行列を有する非線形拡散過程の安定性を支配する関数不等式は何か？
• RQ4 システムの正則性と粒子数が無限に近づく際のEnKFの一様収束性との正確な関係は何か？
• RQ5 フラクチュエーション行列Qがどのような条件下でカルマン＝ブチャトフィルタの観測器が指数的発散を示すか？

主な発見

• 本稿では、任意の時間スケールにわたるEnKFの標本平均および標本共分散行列に対する初めての一様L²平均誤差推定を確立した。
• 粒子数が無限に近づく際のEnKFの一様収束性のための必要十分な正則性条件を同定した。
• フラクチュエーション行列Qの発散集合は、Q₁,₂ > 1.5Q₁,₁ − 1.8 で定義される領域と、Q₁,₁ ∈ ]−8.7, −4.7[ で定義される領域の和集合として明示的に特徴付けられた。
• A = [[1,2],[1,3]] である例のシステムでは、Qの安定部分集合は Q₁,₂ > 1.5Q₁,₁ − 1.8 および Q₁,₁ ∈ ]−4.7, ∞[ を満たすことで定義される。
• 解析により、観測器行列の固有値が異符号をとる場合、または複素固有値の実部が正になる場合に不安定化が生じることが明らかになった。
• 許容可能なフラクチュエーション行列Qの集合は、PおよびQを含む行列式および正定値性制約の組み合わせによって定義される。