Numerical Computation of the Complex Eigenvalues of a Matrix by solving a Square System of Equations



It is well known that if the largest or smallest eigenvalue of a matrix has been computed by some numerical algorithms and one is interested in computing the corresponding eigenvector, one method that is known to give such good approximations to the eigenvector is inverse iteration with a shift. For complex eigenpairs, instead of using Ruhe’s normalization, we show that the natural two norm normalization for the matrix pencil, yields an underdetermined system of equation and by adding an extra equation, the augmented system becomes square which can be solved by LU factorization at a cheaper rate and quadratic convergence is guaranteed. While the underdetermined system of equations can be solved using QR factorization as shown in an earlier work by the same authors, converting it to a square system of equations has the added advantage that besides using LU factorization, it can be solved by several approaches including iterative methods. We show both theoretically and numerically that both algorithms are equivalent in the absence of roundoff errors.

Full Text: PDF
Download the IISTE publication guideline!

To list your conference here. Please contact the administrator of this platform.

Paper submission email:

ISSN (Paper)2224-3186 ISSN (Online)2225-0921

Please add our address "" into your email contact list.

This journal follows ISO 9001 management standard and licensed under a Creative Commons Attribution 3.0 License.

Copyright ©