A wavelet-based method for solving discrete first-kind Fredholm equations with piecewise constant solutions

The inverse problem of nding piecewise constant solutions to discrete Fredholm integral equations of the rst kind arises in many applied elds, e.g. in geophysics. This equation is usually an ill-posed problem when it is considered in a Hilbert space framework, requiring regularization techniques to control arbitrary error ampli cations and to get adequate solutions. In this work, we describe an iterative regularizing method for computing piecewise constant solutions to rst-kind discrete Fredholm integral equations. The algorithm involves two main steps at each iteration: (1) approximating the solution using a new signal reconstruction algorithm from its wavelet maxima which involves a previous step of detecting discontinuities by estimation of its local H older exponents; and (2) obtaining a regularized solution of the original equation using the a priori knowledge and the above approximation. In order to check the behaviour of the proposed technique, we have carried out a statistical study from a high number of simulations obtaining excellent results. Their comparisons with the results coming from using classical Tikhonov regularization by the multiresolution support, total variation (TV) regularization and piecewise polynomial truncated singular value decomposition (PP-TSVD) algorithm, serve to illustrate the advantages of the new method