Solution of Toeplitz Normal Equations by Sine Transform Based Preconditioning
نویسنده
چکیده
The normal equations constructed by a Toeplitz matrix are studied, in order to nd a suitable preconditioner related to the discrete sine transform. New results are given about the structure of the product of two Toeplitz matrices, which allow the CGN method to achieve a superlinear rate of convergence. This preconditioner outperforms the circulant one for the iterative solution of Toeplitz least-squares problems; such strategy can also be applied to nonsymmetric linear systems. A block generalization is discussed.
منابع مشابه
Iterative Solution of Toeplitz Systems by Preconditioning with the Discrete Sine Transform
Solving linear systems or least-squares problems related to Toeplitz matrices is often required in the context of signal and image processing; Conjugate-Gradient-like methods are well-suited for solving such problems. The recent preconditioning technique involving the discrete sine transform is presented: convergence properties are reported and suitable generalizations to block matrices, nonsym...
متن کاملSolving the linearized Navier–Stokes equations using semi-Toeplitz preconditioning
A semi-Toeplitz preconditioner for the linearized Navier– Stokes equation for compressible flow is proposed and tested. The preconditioner is applied to the linear system of equations to be solved in each time step of an implicit method. The equations are solved with flat plate boundary conditions and are linearized around the Blasius solution. The grids are stretched in the normal direction to...
متن کاملSemi-Toeplitz Preconditioning for Linearized Boundary Layer Problems
We have defined and analyzed a semi-Toeplitz preconditioner for timedependent and steady-state convection-diffusion problems. Analytic expressions for the eigenvalues of the preconditioned systems are obtained. An asymptotic analysis shows that the eigenvalue spectrum of the timedependent problem is reduced to two eigenvalues when the number of grid points go to infinity. The numerical experime...
متن کاملPreconditioned Lanczos Methods for the Minimum Eigenvalue of a Symmetric Positive Definite Toeplitz Matrix
In this paper, we apply the preconditioned Lanczos (PL) method to compute the minimum eigenvalue of a symmetric positive definite Toeplitz matrix. The sine transform-based preconditioner is used to speed up the convergence rate of the PL method. The resulting method involves only Toeplitz and sine transform matrix-vector multiplications and hence can be computed efficiently by fast transform al...
متن کاملDisplacement Structure Approach to Discrete-Trigonometric-Transform Based Preconditioners of G.Strang Type and of T.Chan Type
In this paper a displacement structure technique is used to design a class of new preconditioners for the conjugate gradient method applied to the solution of large Toeplitz linear equations. Explicit formulas are suggested for the G.Strang-type and for the T.Chan-type preconditioners belonging to any of 8 classes of matrices diagonalized by the corresponding discrete cosine or sine transforms....
متن کامل