Fractional Tikhonov regularization for linear discrete ill-posed problems
Tóm tắt
Tikhonov regularization is one of the most popular methods for solving linear systems of equations or linear least-squares problems with a severely ill-conditioned matrix A. This method replaces the given problem by a penalized least-squares problem. The present paper discusses measuring the residual error (discrepancy) in Tikhonov regularization with a seminorm that uses a fractional power of the Moore-Penrose pseudoinverse of AA
T
as weighting matrix. Properties of this regularization method are discussed. Numerical examples illustrate that the proposed scheme for a suitable fractional power may give approximate solutions of higher quality than standard Tikhonov regularization.
Tài liệu tham khảo
Björck, Å.: Numerical Methods for Least Squares Problems. SIAM, Philadelphia (1996)
Calvetti, D., Reichel, L.: Lanczos-based exponential filtering for discrete ill-posed problems. Numer. Algorithms 29, 45–65 (2002)
Calvetti, D., Reichel, L.: Tikhonov regularization with a solution constraint. SIAM J. Sci. Comput. 26, 224–239 (2004)
Golub, G.H., Van Loan, C.F.: Matrix Computations, 3rd edn. Johns Hopkins University Press, Baltimore (1996)
Hansen, P.C.: Rank-Deficient and Discrete Ill-Posed Problems. SIAM, Philadelphia (1998)
Hansen, P.C.: Regularization tools version 4.0 for Matlab 7.3. Numer. Algorithms 46, 189–194 (2007)
Klann, E., Ramlau, R.: Regularization by fractional filter methods and data smoothing. Inverse Probl. 24, 025018 (2008)
Lampe, J., Rojas, M., Sorensen, D., Voss, H.: Accelerating the LSTRS algorithm. Bericht 138, Institute of Numerical Simulation, Hamburg University of Technology, Hamburg, Germany, July 2009
Morozov, V.A.: Methods for Solving Incorrectly Posed Problems. Springer, New York (1984)
Rojas, M., Sorensen, D.C.: A trust-region approach to regularization of large-scale discrete forms of ill-posed problems. SIAM J. Sci. Comput. 23, 1842–1860 (2002)