On properties and numerical computation of critical points of eigencurves of bivariate matrix pencils

Bor Plestenjak
2024-11-12
Abstract:We investigate critical points of eigencurves of bivariate matrix pencils $A+\lambda B +\mu C$. Points $(\lambda,\mu)$ for which $\det(A+\lambda B+\mu C)=0$ form algebraic curves in $\mathbb C^2$ and we focus on points where $\mu'(\lambda)=0$. Such points are referred to as zero-group-velocity (ZGV) points, following terminology from engineering applications. We provide a general theory for the ZGV points and show that they form a subset (with equality in the generic case) of the 2D points $(\lambda_0,\mu_0)$, where $\lambda_0$ is a multiple eigenvalue of the pencil $(A+\mu_0 C)+\lambda B$, or, equivalently, there exist nonzero $x$ and $y$ such that $(A+\lambda_0 B+\mu_0 C)x=0$, $y^H(A+\lambda_0 B+\mu_0 C)=0$, and $y^HBx=0$. We introduce three numerical methods for computing 2D and ZGV points. The first method calculates all 2D (ZGV) points from the eigenvalues of a related singular two-parameter eigenvalue problem. The second method employs a projected regular two-parameter eigenvalue problem to compute either all eigenvalues or only a subset of eigenvalues close to a given target. The third approach is a locally convergent Gauss--Newton-type method that computes a single 2D point from an inital approximation, the later can be provided for all 2D points via the method of fixed relative distance by Jarlebring, Kvaal, and Michiels. In our numerical examples we use these methods to compute 2D-eigenvalues, solve double eigenvalue problems, determine ZGV points of a parameter-dependent quadratic eigenvalue problem, evaluate the distance to instability of a stable matrix, and find critical points of eigencurves of a two-parameter Sturm-Liouville problem.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the properties of the critical points of the eigencurves of bivariate matrix pencils and their numerical calculation methods. Specifically, the author focuses on bivariate matrix pencils of the form \(A + \lambda B + \mu C\), where \(A, B, C \in \mathbb{C}^{n \times n}\), \(\lambda, \mu \in \mathbb{C}\), and \(x \in \mathbb{C}^n\) is a non - zero vector. The paper mainly studies the Zero - Group - Velocity (ZGV) points on these eigencurves, that is, the points satisfying \(\mu'(\lambda)=0\). ### Main problems 1. **Properties of ZGV points**: - The paper provides a general theory of ZGV points and proves that these points are a special type of 2D points, that is, \(\lambda_0\) is a multiple eigenvalue of the matrix pencil \((A + \mu_0 C)+\lambda B\). - Specific conditions include the existence of non - zero vectors \(x\) and \(y\) satisfying \((A + \lambda_0 B+\mu_0 C)x = 0\), \(y^H(A + \lambda_0 B+\mu_0 C)=0\), and \(y^H B x = 0\). 2. **Numerical calculation methods**: - The paper introduces three numerical methods for calculating 2D and ZGV points: - The first method calculates all 2D (ZGV) points by solving the related singular two - parameter eigenvalue problem. - The second method uses the projected regular two - parameter eigenvalue problem to calculate all or part of the eigenvalues. - The third method is a locally convergent Gauss - Newton - type method, which calculates a single 2D point starting from an initial approximation. ### Application examples - **2D eigenvalues**: Calculate the 2D eigenvalues of bivariate matrix pencils. - **Double eigenvalue problem**: Solve the double eigenvalue problem. - **Parameter - dependent quadratic eigenvalue problem**: Determine the ZGV points of the parameter - dependent quadratic eigenvalue problem. - **Stability distance**: Evaluate the distance from a stable matrix to instability. - **Sturm - Liouville problem**: Find the critical points of the two - parameter Sturm - Liouville problem. ### Mathematical background - **Eigenvalue problem**: Consider the parameter - dependent linear eigenvalue problem of the form \((A + \lambda B+\mu C)x = 0\). - **Biregular matrix pencil**: Define the biregular matrix pencil, that is, for each \((\lambda_0,\mu_0)\in\mathbb{C}^2\), the generalized eigenvalue problems \((A + \lambda_0 B)+\mu C\) and \((A + \mu_0 C)+\lambda B\) are both regular. - **Definition of ZGV points**: If \(\mu'(\lambda_0)=0\), then \((\lambda_0,\mu_0)\) is called a ZGV point. ### Conclusion Through theoretical analysis and numerical methods, the paper systematically studies the ZGV points on the eigencurves of bivariate matrix pencils, provides a complete theoretical framework and calculation methods, and provides important tools and references for research in related fields.