{"title":"An Adaptive Covariance Parameterization Technique for the Ensemble Gaussian Mixture Filter","authors":"Andrey A. Popov, Renato Zanetti","doi":"10.1137/22m1544312","DOIUrl":null,"url":null,"abstract":"SIAM Journal on Scientific Computing, Volume 46, Issue 3, Page A1949-A1971, June 2024. <br/> Abstract. The ensemble Gaussian mixture filter (EnGMF) combines the simplicity and power of Gaussian mixture models with the provable convergence and power of particle filters. The quality of the EnGMF heavily depends on the choice of covariance matrix in each Gaussian mixture. This work extends the EnGMF to an adaptive choice of covariance based on the parameterized estimates of the sample covariance matrix. Through the use of the expectation maximization algorithm, optimal choices of the covariance matrix parameters are computed in an online fashion. Numerical experiments on the Lorenz ’63 equations show that the proposed methodology converges to classical results known in particle filtering. Further numerical results with more advanced choices of covariance parameterization and the medium-size Lorenz ’96 equations show that the proposed approach can perform significantly better than the standard EnGMF and other classical data assimilation algorithms.","PeriodicalId":49526,"journal":{"name":"SIAM Journal on Scientific Computing","volume":"229 1","pages":""},"PeriodicalIF":3.0000,"publicationDate":"2024-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"SIAM Journal on Scientific Computing","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1137/22m1544312","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
SIAM Journal on Scientific Computing, Volume 46, Issue 3, Page A1949-A1971, June 2024. Abstract. The ensemble Gaussian mixture filter (EnGMF) combines the simplicity and power of Gaussian mixture models with the provable convergence and power of particle filters. The quality of the EnGMF heavily depends on the choice of covariance matrix in each Gaussian mixture. This work extends the EnGMF to an adaptive choice of covariance based on the parameterized estimates of the sample covariance matrix. Through the use of the expectation maximization algorithm, optimal choices of the covariance matrix parameters are computed in an online fashion. Numerical experiments on the Lorenz ’63 equations show that the proposed methodology converges to classical results known in particle filtering. Further numerical results with more advanced choices of covariance parameterization and the medium-size Lorenz ’96 equations show that the proposed approach can perform significantly better than the standard EnGMF and other classical data assimilation algorithms.
期刊介绍:
The purpose of SIAM Journal on Scientific Computing (SISC) is to advance computational methods for solving scientific and engineering problems.
SISC papers are classified into three categories:
1. Methods and Algorithms for Scientific Computing: Papers in this category may include theoretical analysis, provided that the relevance to applications in science and engineering is demonstrated. They should contain meaningful computational results and theoretical results or strong heuristics supporting the performance of new algorithms.
2. Computational Methods in Science and Engineering: Papers in this section will typically describe novel methodologies for solving a specific problem in computational science or engineering. They should contain enough information about the application to orient other computational scientists but should omit details of interest mainly to the applications specialist.
3. Software and High-Performance Computing: Papers in this category should concern the novel design and development of computational methods and high-quality software, parallel algorithms, high-performance computing issues, new architectures, data analysis, or visualization. The primary focus should be on computational methods that have potentially large impact for an important class of scientific or engineering problems.