inria-00409060, version 1

Large Sample Asymptotics for the Ensemble Kalman Filter

François Le Gland () ^a, 1, Valerie Monbet () ², Vu-Duc Tran () ²

N° RR-7014 (2009)

• a –  INRIA
• 1 :  ASPI (INRIA - IRISA)
• http://www.inria.fr/equipes/aspi
  CNRS : UMR6074 – INRIA – Université de Rennes 1 Campus de Beaulieu 35042 Rennes cedex France
• 2 :  Laboratoire des sciences et techniques de l'information, de la communication et de la connaissance (Lab-STICC)
• http://www.lab-sticc.fr/
  CNRS : UMR3192 – Université de Bretagne Occidentale [UBO] – Université de Bretagne Sud – Institut Mines-Télécom – Télécom Bretagne – PRES Université Européenne de Bretagne [UEB] – Institut Supérieur des Sciences et Technologies de Brest (ISSTB) Technopole Brest Iroise BP 832 29285 BREST CEDEX France

Références bibliographiques

• Type de publication : Rapports
• Domaine : Mathématiques/Probabilités
• Titre : Large Sample Asymptotics for the Ensemble Kalman Filter
• Résumé : The ensemble Kalman filter (EnKF) has been proposed as a Monte Carlo, derivative-free, alternative to the extended Kalman filter, and is now widely used in sequential data assimilation, where state vectors of huge dimension (e.g. resulting from the discretization of pressure and velocity fields over a continent, as considered in meteorology) should be estimated from noisy measurements (e.g. collected at sparse in-situ stations). Even if the state and measurement equations are linear with additive Gaussian white noise, computing and storing the error covariance matrices involved in the Kalman filter is practically impossible, and it has been proposed to represent the filtering distribution with a sample (ensemble) of a few elements and to think of the corresponding empirical covariance matrix as an approximation of the intractable error covariance matrix. Extensions to nonlinear state equations have also been proposed. Surprisingly, very little is known about the asymptotic behaviour of the EnKF, whereas on the other hand, the asymptotic behaviour of many different classes of particle filters is well understood, as the number of particles goes to infinity. Interpreting the ensemble elements as a population of particles with mean-field interactions (and not merely as an instrumental device producing the ensemble mean value as an estimate of the hidden state), we prove the convergence of the EnKF, with the classical rate 1/\sqrt{N}, as the number N of ensemble elements increases to infinity. In the linear case, the limit of the empirical distribution of the ensemble elements is the usual (Gaussian distribution associated with the) Kalman filter, as expected, but in the more general case of a nonlinear state equation with linear observations, this limit differs from the usual Bayesian filter. To get the correct limit in this case, the mechanism that generates the elements in the EnKF should be interpreted as a proposal importance distribution, and appropriate importance weights should be assigned to the ensemble elements.
• Langue du document : Anglais
• Type de rapport : Rapport de recherche
• Nombre de pages : 25
• Date de publication : 2009
• Mots-clés : sequential data assimilation – Kalman filter – ensemble Kalman filter (EnKF) – Bayesian filter – particle filter – mean-field interaction – propagation of chaos
• Référence interne : RR-7014

Liste des fichiers attachés à ce document :

• inria-00409060, version 1
• http://hal.inria.fr/inria-00409060
• oai:hal.inria.fr:inria-00409060
• Contributeur : Francois Le Gland <>
• Soumis le : Mercredi 5 Août 2009, 12:16:55
• Dernière modification le : Jeudi 6 Août 2009, 09:47:30