Iterative SVD-Based Methods for Ill-Posed Problems
暂无分享,去创建一个
Very large matrices with rapidly decaying singular values commonly arise in the numerical solution of ill-posed problems. The singular value decomposition (SVD) is a basic tool for both the analysis and computation of solutions to such problems. In most applications, it suffices to obtain a partial SVD consisting of only the largest singular values and their corresponding singular vectors. In this paper, two separate approaches—one based on subspace iteration and the other based on the Lanczos method—are considered for the efficient iterative computation of partial SVDs. In the context of ill-posed problems, an analytical and numerical comparison of these two methods is made and the role of the regularization operator in convergence acceleration is explored.
[1] G. Wahba. Practical Approximate Solutions to Linear Operator Equations When the Data are Noisy , 1977 .
[2] Per Christian Hansen,et al. Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..
[3] C. Vogel,et al. Asymptotic theory of filtering for linear operator equations with discrete noisy data , 1987 .
[4] B. Parlett,et al. Lanczos versus subspace iteration for solution of eigenvalue problems , 1983 .