Multilevel Ensemble Kalman-Bucy Filters
Abstract
In this article we consider the linear filtering problem in continuous-time. We develop and apply multilevel Monte Carlo (MLMC) strategies for ensemble Kalman-Bucy filters (EnKBFs). These filters can be viewed as approximations of conditional McKean-Vlasov-type diffusion processes. They are also interpreted as the continuous-time analogue of the \textit{ensemble Kalman filter}, which has proven to be successful due to its applicability and computational cost. We prove that an ideal version of our multilevel EnKBF can achieve a mean square error (MSE) of $\mathcal{O}(\epsilon^2), \ \epsilon>0$ with a cost of order $\mathcal{O}(\epsilon^{-2}\log(\epsilon)^2)$. In order to prove this result we provide a Monte Carlo convergence and approximation bounds associated to time-discretized EnKBFs. This implies a reduction in cost compared to the (single level) EnKBF which requires a cost of $\mathcal{O}(\epsilon^{-3})$ to achieve an MSE of $\mathcal{O}(\epsilon^2)$. We test our theory on a linear problem, which we motivate through high-dimensional examples of order $\sim \mathcal{O}(10^4)$ and $\mathcal{O}(10^5)$.
- Publication:
-
arXiv e-prints
- Pub Date:
- November 2020
- DOI:
- 10.48550/arXiv.2011.04342
- arXiv:
- arXiv:2011.04342
- Bibcode:
- 2020arXiv201104342C
- Keywords:
-
- Mathematics - Numerical Analysis;
- Statistics - Computation