Abstract
The L-curve is a popular aid for determining a suitable value of the regularization parameter when solving ill-conditioned linear systems of equations with a right-hand side vector, which is contaminated by errors of unknown size. However, for large problems, the computation of the L-curve can be quite expensive, because the determination of a point on the L-curve requires that both the norm of the regularized approximate solution and the norm of the corresponding residual vector be available. Recently, an approximation of the L-curve, referred to as the L-ribbon, was introduced to address this difficulty. The present paper discusses how to organize the computation of the L-ribbon when the matrix of the linear system of equations has many more columns than rows. Numerical examples include an application to computerized tomography.
Similar content being viewed by others
References
M. Bertero and P. Boccacci, Introduction to Inverse Problems in Imaging (Institute of Physics Publ., Bristol, 1998).
D. Calvetti, G.H. Golub and L. Reichel, Estimation of the L-curve via Lanczos bidiagonalization, BIT 39 (1999) 603–619.
D. Calvetti, S. Morigi, L. Reichel and F. Sgallari, Tikhonov regularization and the L-curve for large discrete ill-posed problems, J. Comput. Appl. Math. 123 (2000) 423–446.
D. Calvetti, L. Reichel, F. Sgallari and G. Spaletta, A regularizing Lanczos iteration method for underdetermined linear systems, J. Comput. Appl. Math. 115 (2000) 101–120.
L. Eldén, Algorithms for the regularization of ill-conditioned least squares problems, BIT 17 (1977) 134–145.
G.H. Golub and G. Meurant, Matrices, moments and quadrature, in: Numerical Analysis 1993, eds. D.F. Griffiths and G.A. Watson (Longman, Essex, 1994) pp. 105–156.
G.H. Golub and U. von Matt, Generalized cross-validation for large scale problems, J. Comput. Graph. Stat. 6 (1997) 1–34.
P.C. Hansen, Analysis of discrete ill-posed problems by means of the L-curve, SIAM Review 34 (1992) 561–580.
P.C. Hansen, Regularization tools: a Matlab package for analysis and solution of discrete ill-posed problems, Numer. Algorithms 6 (1994) 1–35.
P.C. Hansen, Rank-Deficient and Discrete Ill-Posed Problems (SIAM, Philadelphia, PA, 1998).
P.C. Hansen and D.P. O'Leary, The use of the L-curve in the regularization of discrete ill-posed problems, SIAM J. Sci. Comput. 14 (1993) 1487–1503.
G.T. Herman, Image Reconstruction from Projections (Academic Press, New York, 1980).
L. Kaufman and A. Neumaier, Regularization of ill-posed problems by envelope guided conjugate gradients, J. Comput. Graph. Stat. 6 (1997) 451–463.
S. Morigi and F. Sgallari, A regularizing L-curve Lanczos method for underdetermined linear systems, Appl. Math. Comput., in press.
C.C. Paige and M.A. Saunders, LSQR: An algorithm for sparse linear equations and sparse least squares, ACM Trans. Math. Software 8 (1982) 43–71.
C.B. Shaw, Jr., Improvements of the resolution of an instrument by numerical solution of an integral equation, J. Math. Anal. Appl. 37 (1972) 83–112.
R.S. Varga, Matrix Iterative Analysis, 2nd ed. (Springer, New York, 2000).
J.H. Wilkinson, The Algebraic Eigenvalue Problem (Oxford Univ. Press, Oxford, 1965).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Calvetti, D., Morigi, S., Reichel, L. et al. An L-ribbon for large underdetermined linear discrete ill-posed problems. Numerical Algorithms 25, 89–107 (2000). https://doi.org/10.1023/A:1016656923184
Issue Date:
DOI: https://doi.org/10.1023/A:1016656923184