2,529 research outputs found
Computing the common zeros of two bivariate functions via Bezout resultants
The common zeros of two bivariate functions can be computed by finding the common zeros of their polynomial interpolants expressed in a tensor Chebyshev basis. From here we develop a bivariate rootfinding algorithm based on the hidden variable resultant method and B�ezout matrices with polynomial entries. Using techniques including domain subdivision, B�ezoutian regularization and local refinement we are able to reliably and accurately compute the simple common zeros of two smooth functions with polynomial interpolants of very high degree (� 1000). We analyze the resultant method and its conditioning by noting that the B�ezout matrices are matrix polynomials. Our robust algorithm is implemented in the roots command in Chebfun2, a software package written in object-oriented MATLAB for computing with bivariate functions
Perturbation, extraction and refinement of invariant pairs for matrix polynomials
Generalizing the notion of an eigenvector, invariant subspaces are frequently used in the context of linear eigenvalue problems, leading to conceptually elegant and numerically stable formulations in applications that require the computation of several eigenvalues and/or eigenvectors. Similar benefits can be expected for polynomial eigenvalue problems, for which the concept of an invariant subspace needs to be replaced by the concept of an invariant pair. Little has been known so far about numerical aspects of such invariant pairs. The aim of this paper is to fill this gap. The behavior of invariant pairs under
perturbations of the matrix polynomial is studied and a first-order perturbation expansion is given. From a computational point of view, we investigate how to best extract invariant pairs from a linearization of the matrix polynomial. Moreover, we describe efficient refinement procedures directly based on the polynomial formulation. Numerical experiments
with matrix polynomials from a number of applications demonstrate the effectiveness of our extraction and refinement procedures
Localization theorems for nonlinear eigenvalue problems
Let T : \Omega \rightarrow \bbC^{n \times n} be a matrix-valued function
that is analytic on some simply-connected domain \Omega \subset \bbC. A point
is an eigenvalue if the matrix is singular.
In this paper, we describe new localization results for nonlinear eigenvalue
problems that generalize Gershgorin's theorem, pseudospectral inclusion
theorems, and the Bauer-Fike theorem. We use our results to analyze three
nonlinear eigenvalue problems: an example from delay differential equations, a
problem due to Hadeler, and a quantum resonance computation.Comment: Submitted to SIMAX. 22 pages, 11 figure
- …