Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published February 2018 | Submitted
Journal Article Open

Iterative Updating of Model Error for Bayesian Inversion

Abstract

In computational inverse problems, it is common that a detailed and accurate forward model is approximated by a computationally less challenging substitute. The model reduction may be necessary to meet constraints in computing time when optimization algorithms are used to find a single estimate, or to speed up Markov chain Monte Carlo (MCMC) calculations in the Bayesian framework. The use of an approximate model introduces a discrepancy, or modeling error, that may have a detrimental effect on the solution of the ill-posed inverse problem, or it may severely distort the estimate of the posterior distribution. In the Bayesian paradigm, the modeling error can be considered as a random variable, and by using an estimate of the probability distribution of the unknown, one may estimate the probability distribution of the modeling error and incorporate it into the inversion. We introduce an algorithm which iterates this idea to update the distribution of the model error, leading to a sequence of posterior distributions that are demonstrated empirically to capture the underlying truth with increasing accuracy. Since the algorithm is not based on rejections, it requires only limited full model evaluations. We show analytically that, in the linear Gaussian case, the algorithm converges geometrically fast with respect to the number of iterations when the data is finite dimensional. For more general models, we introduce particle approximations of the iteratively generated sequence of distributions; we also prove that each element of the sequence converges in the large particle limit under a simplifying assumption. We show numerically that, as in the linear case, rapid convergence occurs with respect to the number of iterations. Additionally, we show through computed examples that point estimates obtained from this iterative algorithm are superior to those obtained by neglecting the model error.

Additional Information

© 2018 IOP Publishing Ltd. Published 18 January 2018. The work of D Calvetti is partially supported by NSF grant DMS-1522334. E Somersalo's work is partly supported by the NSF grant DMS-1312424. The research of AM Stuart was partially supported by the EPSRC programme grant EQUIP, by AFOSR Grant FA9550-17-1-0185 and ONR Grant N00014-17-1-2079. M Dunlop was partially supported by the EPSRC MASDOC Graduate Training Program. Both M Dunlop and AM Stuart are supported by DARPA funded program Enabling Quantification of Uncertainty in Physical Systems (EQUiPS), contract W911NF-15-2-0121.

Attached Files

Submitted - 1707.04246.pdf

Files

1707.04246.pdf
Files (5.5 MB)
Name Size Download all
md5:458bb868531baedbe5e799e6a542158f
5.5 MB Preview Download

Additional details

Created:
August 19, 2023
Modified:
March 5, 2024