On the Lanczos and Golub–Kahan reduction methods applied to discrete ill‐posed problems. Issue 1 (12th October 2015)

• Record Type: Journal Article
• Authors: Gazzola, Silvia; Onunwor, Enyinda; Reichel, Lothar; Rodriguez, Giuseppe
• Is Part Of: Numerical linear algebra with applications. Volume 23:Issue 1(2016:Jan.)
• Page Start: 187
• Page End: 204
• Publication Date: 2015-10-12
• Subjects: discrete ill‐posed problems -- Lanczos decomposition -- Golub–Kahan bidiagonalization -- LSQR -- TSVD; Algebras, Linear -- Periodicals; 512.5
• DOI: 10.1002/nla.2020
• Languages: English
• ISSNs: 1070-5325

Abstract

Summary: The symmetric Lanczos method is commonly applied to reduce large‐scale symmetric linear discrete ill‐posed problems to small ones with a symmetric tridiagonal matrix. We investigate how quickly the nonnegative subdiagonal entries of this matrix decay to zero. Their fast decay to zero suggests that there is little benefit in expressing the solution of the discrete ill‐posed problems in terms of the eigenvectors of the matrix compared with using a basis of Lanczos vectors, which are cheaper to compute. Similarly, we show that the solution subspace determined by the LSQR method when applied to the solution of linear discrete ill‐posed problems with a nonsymmetric matrix often can be used instead of the solution subspace determined by the singular value decomposition without significant, if any, reduction of the quality of the computed solution. Copyright © 2015 John Wiley & Sons, Ltd.