Filtering Dynamical Systems Using Observations of Statistics

Abstract

We consider the problem of filtering dynamical systems, possibly stochastic, using observations of statistics. Thus, the computational task is to estimate a time-evolving density $\rho (v,t)$ given noisy observations of the true density ${\rho }^{†}$; this contrasts with the standard filtering problem based on observations of the state $v$. The task is naturally formulated as an infinite-dimensional filtering problem in the space of densities $\rho$. However, for the purposes of tractability, we seek algorithms in state space; specifically, we introduce a mean-field state-space model, and using interacting particle system approximations to this model, we propose an ensemble method. We refer to the resulting methodology as the ensemble Fokker-Planck filter (EnFPF). Under certain restrictive assumptions, we show that the EnFPF approximates the Kalman-Bucy filter for the Fokker-Planck equation, which is the exact solution to the infinite-dimensional filtering problem. Furthermore, our numerical experiments show that the methodology is useful beyond this restrictive setting. Specifically, the experiments show that the EnFPF is able to correct ensemble statistics, to accelerate convergence to the invariant density for autonomous systems, and to accelerate convergence to time-dependent invariant densities for non-autonomous systems. We discuss possible applications of the EnFPF to climate ensembles and to turbulence modeling.

arXiv comment15 pages, 8 figures

References in corpus (7)

• Climate dynamics and fluid mechanics: Natural variability and related uncertainties
• MAP Estimators and Their Consistency in Bayesian Nonparametric Inverse Problems
• Using machine learning to correct model error in data assimilation and forecast applications
• Accelerating diffusions
• Data Assimilation: The Schrödinger Perspective
• A multi-model ensemble Kalman filter for data assimilation and forecasting
• parasweep: A template-based utility for generating, dispatching, and post-processing of parameter sweeps