Circulant preconditioners for discrete ill-posed Toeplitz systems

被引:0
|
作者
L. Dykes
S. Noschese
L. Reichel
机构
[1] Kent State University,Department of Mathematical Sciences
[2] University School,Dipartimento di Matematica “Guido Castelnuovo”
[3] SAPIENZA Università di Roma,undefined
来源
Numerical Algorithms | 2017年 / 75卷
关键词
Ill-posed problem; Deconvolution; FFT; Image deblurring; 65F10; 65F15; 65F30;
D O I
暂无
中图分类号
学科分类号
摘要
Circulant preconditioners are commonly used to accelerate the rate of convergence of iterative methods when solving linear systems of equations with a Toeplitz matrix. Block extensions that can be applied when the system has a block Toeplitz matrix with Toeplitz blocks also have been developed. This paper is concerned with preconditioning of linear systems of equations with a symmetric block Toeplitz matrix with symmetric Toeplitz blocks that stem from the discretization of a linear ill-posed problem. The right-hand side of the linear systems represents available data and is assumed to be contaminated by error. These kinds of linear systems arise, e.g., in image deblurring problems. It is important that the preconditioner does not affect the invariant subspace associated with the smallest eigenvalues of the block Toeplitz matrix to avoid severe propagation of the error in the right-hand side. A perturbation result indicates how the dimension of the subspace associated with the smallest eigenvalues should be chosen and allows the determination of a suitable preconditioner when an estimate of the error in the right-hand side is available. This estimate also is used to decide how many iterations to carry out by a minimum residual iterative method. Applications to image restoration are presented.
引用
收藏
页码:477 / 490
页数:13
相关论文
共 50 条