An implicit shift bidiagonalization algorithm for ill-posed systems |
| |
Authors: | Åke Björck Eric Grimme Paul Van Dooren |
| |
Institution: | (1) Department of Mathematics, Linköping University, S-581 83 Linköping, Sweden;(2) Coordinated Science Laboratory, University of Illinois at Urbana-Champaign, 1308 W. Main Street, 61801 Urbana, IL |
| |
Abstract: | Iterative methods based on Lanczos bidiagonalization with full reorthogonalization (LBDR) are considered for solving large-scale discrete ill-posed linear least-squares problems of the form min
x
Ax–b2. Methods for regularization in the Krylov subspaces are discussed which use generalized cross validation (GCV) for determining the regularization parameter. These methods have the advantage that no a priori information about the noise level is required. To improve convergence of the Lanczos process we apply a variant of the implicitly restarted Lanczos algorithm by Sorensen using zero shifts. Although this restarted method simply corresponds to using LBDR with a starting vector (AA
T)
p
b, it is shown that carrying out the process implicitly is essential for numerical stability. An LBDR algorithm is presented which incorporates implicit restarts to ensure that the global minimum of the CGV curve corresponds to a minimum on the curve for the truncated SVD solution. Numerical results are given comparing the performance of this algorithm with non-restarted LBDR.This work was partially supported by DARPA under grant 60NANB2D1272 and by the National Science Foundation under grant CCR-9209349. |
| |
Keywords: | 65F20 |
本文献已被 SpringerLink 等数据库收录! |
|