TY - JOUR
T1 - Generalization of Strang's preconditioner with applications to Toeplitz least squares problems
AU - Chan, Raymond H.
AU - Ng, Michael K.
AU - Plemmons, Robert J.
PY - 1996/1
Y1 - 1996/1
N2 - In this paper, we propose a method to generalize Strang's circulant preconditioner for arbitrary n-by-n matrices An. The [n/2] th column of our circulant preconditioned Sn is equal to the [n/2 ]th column of the given matrix An. Thus if An is a square Toeplitz matrix, then Sn is just the Strang circulant preconditioner. When Sn is not Hermitian, our circulant preconditioner can be defined as (Sn*Sn) 1/2. This construction is similar to the forward-backward projection method used in constructing preconditioners for tomographic inversion problems in medical imaging. We show that if the matrix An has decaying coefficients away from the main diagonal, then (Sn*Sn) 1/2 is a good preconditioner for An. Comparisons of our preconditioner with other circulant-based preconditioners are carried out for some 1-D Toeplitz least squares problems: min ||b - Ax ||2 . Preliminary numerical results show that our preconditioner performs quite well, in comparison to other circulant preconditioners. Promising test results are also reported for a 2-D deconvolution problem arising in ground-based atmospheric imaging.
AB - In this paper, we propose a method to generalize Strang's circulant preconditioner for arbitrary n-by-n matrices An. The [n/2] th column of our circulant preconditioned Sn is equal to the [n/2 ]th column of the given matrix An. Thus if An is a square Toeplitz matrix, then Sn is just the Strang circulant preconditioner. When Sn is not Hermitian, our circulant preconditioner can be defined as (Sn*Sn) 1/2. This construction is similar to the forward-backward projection method used in constructing preconditioners for tomographic inversion problems in medical imaging. We show that if the matrix An has decaying coefficients away from the main diagonal, then (Sn*Sn) 1/2 is a good preconditioner for An. Comparisons of our preconditioner with other circulant-based preconditioners are carried out for some 1-D Toeplitz least squares problems: min ||b - Ax ||2 . Preliminary numerical results show that our preconditioner performs quite well, in comparison to other circulant preconditioners. Promising test results are also reported for a 2-D deconvolution problem arising in ground-based atmospheric imaging.
KW - Atmospheric imaging
KW - Circulant preconditioned conjugate gradient method
KW - Deconvolution
KW - Image restoration
KW - Medical imaging
KW - Toeplitz least squares problems
UR - http://www.scopus.com/inward/record.url?scp=21844500518&partnerID=8YFLogxK
UR - https://www.scopus.com/record/pubmetrics.uri?eid=2-s2.0-21844500518&origin=recordpage
U2 - 10.1002/(SICI)1099-1506(199601/02)3:1<45::AID-NLA70>3.0.CO;2-T
DO - 10.1002/(SICI)1099-1506(199601/02)3:1<45::AID-NLA70>3.0.CO;2-T
M3 - RGC 21 - Publication in refereed journal
VL - 3
SP - 45
EP - 64
JO - Numerical Linear Algebra with Applications
JF - Numerical Linear Algebra with Applications
SN - 1070-5325
IS - 1
ER -