Recursive linearization method for inverse medium scattering problems with complex mixture Gaussian error learning
arXiv:1804.05471 · doi:10.1088/1361-6420/ab08f2
Abstract
This paper is concerned with the modeling errors appeared in the numerical methods of inverse medium scattering problems (IMSP). Optimization based iterative methods are wildly employed to solve IMSP, which are computationally intensive due to a series of Helmholtz equations need to be solved numerically. Hence, rough approximations of Helmholtz equations can significantly speed up the iterative procedure. However, rough approximations will lead to instability and inaccurate estimations. Using the Bayesian inverse methods, we incorporate the modelling errors brought by the rough approximations. Modelling errors are assumed to be some complex Gaussian mixture (CGM) random variables, and in addition, well-posedness of IMSP in the statistical sense has been established by extending the general theory to involve CGM noise. Then, we generalize the real valued expectation-maximization (EM) algorithm used in the machine learning community to our complex valued case to learn parameters in the CGM distribution. Based on these preparations, we generalize the recursive linearization method (RLM) to a new iterative method named as Gaussian mixture recursive linearization method (GMRLM) which takes modelling errors into account. Finally, we provide two numerical examples to illustrate the effectiveness of the proposed method.
arXiv comment26 pages
References in corpus (6)
- MAP Estimators and Their Consistency in Bayesian Nonparametric Inverse Problems
- Well-Posed Bayesian Geometric Inverse Problems Arising in Subsurface Flow
- Iterative Updating of Model Error for Bayesian Inversion
- A Two-stage Method for Inverse Medium Scattering
- MAP Estimators for Piecewise Continuous Inversion
- Infinite-dimensional Bayesian approach for inverse scattering problems of a fractional Helmholtz equation
Cited by in corpus (3)
- Variational Bayes' method for functions with applications to some inverse problems
- Stein variational gradient descent on infinite-dimensional space and applications to statistical inverse problems
- Non-centered parametric variational Bayes' approach for hierarchical inverse problems of partial differential equations