[论文解读] 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技术、耦合方法与谱分析，首次建立了Ensemble Kalman-Bucy滤波器的统一稳定性与混沌传播结果。论文证明了在任意时间范围内样本均值和协方差矩阵的统一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.

研究动机与目标

• 为解决高维、病态系统中Ensemble Kalman-Bucy滤波器长期稳定性与收敛性的理论理解不足问题。
• 在任意时间范围内建立样本均值与协方差矩阵的统一L²-均值误差估计。
• 分析协方差依赖型扩散矩阵的非线性扩散过程的稳定性，此类系统此前未被深入研究。
• 推导当粒子数趋于无穷时EnKF统一收敛性的必要与充分正则性条件。
• 通过Jensen型不等式与谱分析，量化Kalman-Bucy滤波器及其关联Riccati方程的稳定性。

提出的方法

• 提出一类新型函数不等式，以量化具有协方差依赖型扩散矩阵的非线性McKean-Vlasov型扩散过程的稳定性。
• 应用Foster-Lyapunov技术与耦合方法，分析Wasserstein距离与相对熵下的指数稳定性。
• 结合鞅方法、随机矩阵理论与谱分析，研究相互作用粒子系统的长期行为。
• 利用Kalman-Bucy扩散的平均场粒子解释，将EnKF建模为非线性相互作用粒子系统。
• 通过分析样本统计量（均值与协方差）随时间向真实值收敛的过程，推导出统一L²-均值误差估计。
• 表征了使系统保持稳定的可接受波动矩阵Q的集合，识别出发生不稳定的发散区域。

实验结果

研究问题

• RQ1 Ensemble Kalman-Bucy滤波器在何种条件下于任意时间范围内表现出统一稳定性？
• RQ2 如何在时间上统一量化EnKF粒子系统的混沌传播特性？
• RQ3 哪些函数不等式控制具有协方差依赖型扩散矩阵的非线性扩散过程的稳定性？
• RQ4 系统正则性与EnKF在粒子数趋于无穷时的统一收敛性之间存在何种精确关系？
• RQ5 在何种Q波动矩阵条件下，Kalman-Bucy观测器会发生指数发散？

主要发现

• 论文首次建立了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的行列式与正定性约束的组合定义。