Solving a structured quadratic eigenvalue problem by a structure-preserving doubling algorithm

dc.contributor.authorGuo, Chun-Hua
dc.contributor.authorLin, Wen-Wei
dc.date.accessioned2014-04-27T23:16:01Z
dc.date.available2014-04-27T23:16:01Z
dc.date.issued2010
dc.description.abstractIn studying the vibration of fast trains, we encounter a palindromic quadratic eigenvalue problem (QEP) $(\lambda^2 A^T + \lambda Q + A)z = 0$, where $A, Q \in \mathbb{C}^{n \times n}$ and $Q^T = Q$. Moreover, the matrix $Q$ is block tridiagonal and block Toeplitz, and the matrix $A$ has only one nonzero block in the upper-right corner. So most of the eigenvalues of the QEP are zero or infinity. In a linearization approach, one typically starts with deflating these known eigenvalues, for the sake of efficiency. However, this initial deflation process involves the inverses of two potentially ill-conditioned matrices. As a result, large error might be introduced into the data for the reduced problem. In this paper we propose using the solvent approach directly on the original QEP, without any deflation process. We apply a structure-preserving doubling algorithm to compute the stabilizing solution of the matrix equation $X+A^TX^{-1}A=Q$, whose existence is guaranteed by a result on the Wiener--Hopf factorization of rational matrix functions associated with semi-infinite block Toeplitz matrices and a generalization of Bendixson's theorem to bounded linear operators on Hilbert spaces. The doubling algorithm is shown to be well defined and quadratically convergent. The complexity of the doubling algorithm is drastically reduced by using the Sherman--Morrison--Woodbury formula and the special structures of the problem. Once the stabilizing solution is obtained, all nonzero finite eigenvalues of the QEP can be found efficiently and with the automatic reciprocal relationship, while the known eigenvalues at zero or infinity remain intact.en_US
dc.description.authorstatusFacultyen_US
dc.description.peerreviewyesen_US
dc.description.sponsorshipNSERC, NSC (Taiwan), NCTS (Taiwan)en_US
dc.identifier.citationSIAM J. Matrix Anal. Appl.en_US
dc.identifier.urihttps://hdl.handle.net/10294/5264
dc.language.isoenen_US
dc.publisherSIAMen_US
dc.titleSolving a structured quadratic eigenvalue problem by a structure-preserving doubling algorithmen_US
dc.typeArticleen_US

Files

Original bundle

Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
GL10.pdf
Size:
309.36 KB
Format:
Adobe Portable Document Format

License bundle

Now showing 1 - 1 of 1
No Thumbnail Available
Name:
license.txt
Size:
2.24 KB
Format:
Item-specific license agreed upon to submission
Description:

Collections