Computation of multiple eigenvalues and generalized eigenvectors for matrices dependent on parameters

Abstract

The paper develops Newton's method of finding multiple eigenvalues with one Jordan block and corresponding generalized eigenvectors for matrices dependent on parameters. It computes the nearest value of a parameter vector with a matrix having a multiple eigenvalue of given multiplicity. The method also works in the whole matrix space (in the absence of parameters). The approach is based on the versal deformation theory for matrices. Numerical examples are given. The implementation of the method in MATLAB code is available.

19 pages, 3 figures

References in corpus (1)

• Coupling of eigenvalues of complex matrices at diabolic and exceptional points

Cited by in corpus (11)

• Exceptional points for parameter estimation in open quantum systems: Analysis of the Bloch equations
• Miniversal deformations of matrices under *congruence and reducing transformations
• Robust mode conversion in NV centers using exceptional points
• Linking scalar elastodynamics and non-Hermitian quantum mechanics
• Parameter estimation in atomic spectroscopy using exceptional points
• A high order continuation method to locate exceptional points and to compute Puiseux series with applications to acoustic waveguides
• Experimental observation of exceptional points in coupled pendulums
• Linear level repulsions near exceptional points of non-Hermitian systems
• A Hermitian bypass to the non-Hermitian quantum theory
• Equations of motion governing the dynamics of the exceptional points of parameterically dependent nonhermitian Hamiltonians
• Fast recovery of parametric eigenvalues depending on several parameters and location of high order exceptional points