Stopping criterion for iterative regularization of large-scale ill-posed problems using the Picard parameter
From MaRDI portal
Publication:6288967
arXiv1707.04200MaRDI QIDQ6288967FDOQ6288967
Authors: Eitan Levin, A. Y. Meltzer 
Publication date: 13 July 2017
Abstract: We propose a new stopping criterion for Krylov subspace iterative regularization of large-scale ill-posed inverse problems. Our stopping criterion accurately filters the data using a generalization of the Picard parameter that was originally introduced for direct regularization of small-scale problems. In the one dimension we filter the data in the discrete Fourier transform (DFT) basis using the Picard parameter, which separates noise-dominated Fourier coefficients from the signal-dominated ones. For two-dimensional problems we propose a novel vectorization scheme of the Fourier coefficients of the data based on the Kronecker product structure of the two-dimensional DFT matrix, which effectively reduces the problem to one dimension. At each iteration we compute the distance between the data reconstructed from the iterated solution and the filtered data, terminating the iterations once this distance begins to increase or to level off. The accuracy and robustness of the proposed method is demonstrated by several numerical examples and a MATLAB-based implementation is provided.
Ill-posedness and regularization problems in numerical linear algebra (65F22) Linear equations (linear algebraic aspects) (15A06) Iterative numerical methods for linear systems (65F10)
This page was built for publication: Stopping criterion for iterative regularization of large-scale ill-posed problems using the Picard parameter
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6288967)
