Accurate solutions of polynomial eigenvalue problems

作     者：YOU, YILING RODRIGUEZ, JOSE ISRAEL LIM, LEK-HENG 

作者机构：Department of Statistics University of Chicago ChicagoIL United States Computational and Applied Mathematics Initiative Department of Statistics University of Chicago ChicagoIL United States 

出 版 物：《arXiv》 (arXiv)

年 卷 期：2017年

核心收录：

主　　题：Eigenvalues and eigenfunctions 

摘      要：Quadratic eigenvalue problems (QEP) and more generally polynomial eigenvalue problems (PEP) are among the most common types of nonlinear eigenvalue problems. Both problems, especially the QEP, have extensive applications. A typical approach to solve QEP and PEP is to use a linearization method to reformulate the problem as a higher dimensional linear eigenvalue problem. In this article, we use homotopy continuation to solve these nonlinear eigenvalue problems without passing to higher dimensions. Our main contribution is to show that our method produces substantially more accurate results, and finds all eigenvalues with a certificate of correctness via Smale s α-theory. To explain the superior accuracy, we show that the nonlinear eigenvalue problem we solve is better conditioned than its reformulated linear eigenvalue problem, and our homotopy continuation algorithm is more stable than QZ algorithm | theoretical findings that are borne out by our numerical experiments. Our studies provide yet another illustration of the dictum in numerical analysis that, for reasons of conditioning and stability, it is sometimes better to solve a nonlinear problem directly even when it could be transformed into a linear problem with the same solution mathematically. Copyright © 2017, The Authors. All rights reserved.