6 research outputs found
Iteratively regularized Newton-type methods for general data misfit functionals and applications to Poisson data
We study Newton type methods for inverse problems described by nonlinear
operator equations in Banach spaces where the Newton equations
are regularized variationally using a general
data misfit functional and a convex regularization term. This generalizes the
well-known iteratively regularized Gauss-Newton method (IRGNM). We prove
convergence and convergence rates as the noise level tends to 0 both for an a
priori stopping rule and for a Lepski{\u\i}-type a posteriori stopping rule.
Our analysis includes previous order optimal convergence rate results for the
IRGNM as special cases. The main focus of this paper is on inverse problems
with Poisson data where the natural data misfit functional is given by the
Kullback-Leibler divergence. Two examples of such problems are discussed in
detail: an inverse obstacle scattering problem with amplitude data of the
far-field pattern and a phase retrieval problem. The performence of the
proposed method for these problems is illustrated in numerical examples
Regularization of statistical inverse problems and the Bakushinskii veto
In the deterministic context Bakushinskii's theorem excludes the existence of
purely data driven convergent regularization for ill-posed problems. We will
prove in the present work that in the statistical setting we can either
construct a counter example or develop an equivalent formulation depending on
the considered class of probability distributions. Hence, Bakushinskii's
theorem does not generalize to the statistical context, although this has often
been assumed in the past. To arrive at this conclusion, we will deduce from the
classic theory new concepts for a general study of statistical inverse problems
and perform a systematic clarification of the key ideas of statistical
regularization.Comment: 20 page