Approximation Accuracy of the Krylov Subspaces for Linear Discrete Ill-Posed Problems
Abstract
For the large-scale linear discrete ill-posed problem or with contaminated by Gaussian white noise, the Lanczos bidiagonalization based Krylov solver LSQR and its mathematically equivalent CGLS, the Conjugate Gradient (CG) method implicitly applied to , are most commonly used, and CGME, the CG method applied to or with , and LSMR, which is equivalent to the minimal residual (MINRES) method applied to , have also been choices. These methods exhibit typical semi-convergence feature, and the iteration number plays the role of the regularization parameter. However, there has been no definitive answer to the long-standing fundamental question: {\em Can LSQR and CGLS find 2-norm filtering best possible regularized solutions}? The same question is for CGME and LSMR too. At iteration , LSQR, CGME and LSMR compute {\em different} iterates from the {\em same} dimensional Krylov subspace. A first and fundamental step towards to answering the above question is to {\em accurately} estimate the accuracy of the underlying dimensional Krylov subspace approximating the dimensional dominant right singular subspace of . Assuming that the singular values of are simple, we present a general theorem for the 2-norm distances between these two subspaces and derive accurate estimates on them for severely, moderately and mildly ill-posed problems. We also establish some relationships between the smallest Ritz values and these distances. Numerical experiments justify the sharpness of our results.
Cite
@article{arxiv.1805.10132,
title = {Approximation Accuracy of the Krylov Subspaces for Linear Discrete Ill-Posed Problems},
author = {Zhongxiao Jia},
journal= {arXiv preprint arXiv:1805.10132},
year = {2020}
}
Comments
33 pages, 11 figures. arXiv admin note: text overlap with arXiv:1701.05708, arXiv:1608.05907