Mixed gradient-Tikhonov methods for solving nonlinear ill-posed problems in Banach spaces. (24th November 2016)

• Record Type: Journal Article
• Authors: Margotti, Fábio
• Is Part Of: Inverse problems. Volume 32:Number 12(2016:Dec.)
• Publication Date: 2016-11-24
• Subjects: Banach spaces -- electrical impedance tomography -- Tikhonov regularization -- gradient methods -- sparsity constraints; Inverse problems (Differential equations) -- Periodicals; 515.357
• DOI: 10.1088/0266-5611/32/12/125012
• Languages: English
• ISSNs: 0266-5611

Abstract

Abstract: Tikhonov regularization is a very useful and widely used method for finding stable solutions of ill-posed problems. A good choice of the penalization functional as well as a careful selection of the topologies of the involved spaces is fundamental to the quality of the reconstructions. These choices can be combined with some a priori information about the solution in order to preserve desired characteristics like sparsity constraints for example. To prove convergence and stability properties of this method, one usually has to assume that a minimizer of the Tikhonov functional is known. In practical situations however, the exact computation of a minimizer is very difficult and even finding an approximation can be a very challenging and expensive task if the involved spaces have poor convexity or smoothness properties. In this paper we propose a method to attenuate this gap between theory and practice, applying a gradient-like method to a Tikhonov functional in order to approximate a minimizer. Using only available information, we explicitly calculate a maximal step-size which ensures a monotonically decreasing error. The resulting algorithm performs only finitely many steps and terminates using the discrepancy principle. In particular the knowledge of a minimizer or even its existence does not need to be assumed. Under standard assumptions, we prove convergence and stability results in relatively general Banach spaces, and subsequently, test its performance numerically, reconstructing conductivities with sparsely located inclusions and different kinds of noise in the 2D electrical impedance tomography.