Location of Repository

Approximation of Bayesian inverse problems for PDEs

By Simon L. Cotter, M. Dashti and A. M. Stuart

Abstract

Inverse problems are often ill posed, with solutions that depend sensitively on data. In any numerical approach to the solution of such problems, regularization of some form is needed to counteract the resulting instability. This paper is based on an approach to regularization, employing a Bayesian formulation of the problem, which leads to a notion of well posedness for inverse problems, at the level of probability measures. The stability which results from this well posedness may be used as the basis for quantifying the approximation, in finite dimensional spaces, of inverse problems for functions. This paper contains a theory which utilizes this stability property to estimate the distance between the true and approximate posterior distributions, in the Hellinger metric, in terms of error estimates for approximation of the underlying forward problem. This is potentially useful as it allows for the transfer of estimates from the numerical analysis of forward problems into estimates for the solution of the related inverse problem. It is noteworthy that, when the prior is a Gaussian random field model, controlling differences in the Hellinger metric leads to control on the differences between expected values of polynomially bounded functions and operators, including the mean and covariance operator. The ideas are applied to some non-Gaussian inverse problems where the goal is determination of the initial condition for the Stokes or Navier–Stokes equation from Lagrangian and Eulerian observations, respectively

Topics: QA
Publisher: Society for Industrial and Applied Mathematics
Year: 2010
OAI identifier: oai:wrap.warwick.ac.uk:3313

Suggested articles

Preview

Citations

  1. (2008). A .A p t e ,C .K .R .T .J o n e s ,A .M .S t u a r t ,a n dJ .V o s s , Data assimilation: Mathematical and statistical perspectives,
  2. (2009). A simple proof of uniqueness of the particle trajectories for solutions of the Navier-Stokes equations, doi
  3. (1993). A study of uniqueness for the initialiazation problem for burgers’
  4. (1996). An Introduction to the Mathematical Theory of Inverse Problems, doi
  5. (2007). Approximation errors in nonstationary inverse problems, doi
  6. (2010). Bayesian inverse problems for functions and applications to fluid mechanics, doi
  7. (2007). Convergence rates for sparse chaos approximations of elliptic problems with stochastic coefficients, doi
  8. (2007). Convergence rates for the Bayesian approach to linear inverse problems, doi
  9. (2008). Convergence results for the Bayesian inversion theory, doi
  10. (2003). Data assimilation: Aims and basic concepts, Data Assimilation for the Earth System, doi
  11. (1998). Gaussian Measures, doi
  12. (2007). Introduction to bayesian scientific computing, doi
  13. (2005). Inverse Problem Theory. doi
  14. (2010). Inverse problems: A Bayesian approach, doi
  15. (2006). Karhunen-Loeve approximation of random fields in domains by generalized fast multipole methods, doi
  16. (2006). Kolehmainen,M .S c h w e i g e r ,E .S o m e r s a l o ,a n dT .T a r -vainen, Approximation errors and model reduction with an application to optical diffusion tomography, doi
  17. (1989). Linear inverse problems for generalized random variables, doi
  18. (2001). Monte Carlo strategies in scientific computing, Springer Texts in Statistics, doi
  19. (1995). Navier-Stokes equations and nonlinear functional analysis, doi
  20. (2001). Navier-Stokes equations. AMS doi
  21. (1983). Semigroups of Linear Operators and Applications to Partial Differential Equations, doi
  22. (2006). State estimation approach to nonstationary inverse problems: Discretization error and filtering problem, doi
  23. (2005). Statistical and computational inverse problems, doi
  24. (2007). Statistical inverse problems: Discretization, model reduction and inverse crimes, doi
  25. (1992). Stochastic equations in infinite dimensions, Encyclopedia of Mathematics and its Applications 44, doi
  26. (1989). Stochastic finite element expansion for random media,J doi
  27. (2003). Stochastic Finite Elements: A Spectral Approach, doi
  28. (1997). Unified notation for assimilation: Operational, sequential and variational,
  29. (1970). Well-posed stochastic extensions of ill posed linear problems, doi
  30. (1998). When are quasi-Monte Carlo algorithms efficient for highdimensional integrals? doi

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.