Michiel E. Hochstenbach,Christian Mehl,Bor Plestenjak
Abstract:In Parts I and II of this series of papers, three new methods for the computation of eigenvalues of singular pencils were developed: rank-completing perturbations, rank-projections, and augmentation. It was observed that a straightforward structure-preserving adaption for symmetric pencils was not possible and it was left as an open question how to address this challenge. In this Part III, it is shown how the observed issue can be circumvented by using Hermitian perturbations. This leads to structure-preserving analogues of the three techniques from Parts I and II for Hermitian pencils (including real symmetric pencils) as well as for related structures. It is an important feature of these methods that the sign characteristic of the given pencil is preserved. As an application, it is shown that the resulting methods can be used to solve systems of bivariate polynomials.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
The paper aims to address the challenge of structure - preservation in the singular generalized eigenvalue problem, especially for symmetric and Hermitian pencils. In the first two parts of the research, the author proposes three new methods for calculating the eigenvalues of singular pencils: rank - completing perturbations, rank - projections, and augmentation. However, these methods encounter difficulties when directly applied to symmetric pencils because it is not possible to simply keep the structure unchanged. Specifically, when using symmetric perturbations, random eigenvalues may appear as double real eigenvalues or complex conjugate pairs, which makes the orthogonality constraint unfulfilled.
To solve this problem, this paper proposes a new method, that is, to bypass the above challenges by using Hermitian perturbations. This method can be applied to Hermitian pencils (including real symmetric pencils) and related structures, and importantly, it can preserve the sign characteristic of a given pencil. In addition, the paper shows that these methods can be used to solve binary polynomial systems.
### Formula summary
1. **Singular generalized eigenvalue problem**:
\[
Ax=\lambda Bx
\]
where \(A^{*} = A\) and \(B^{*} = B\) are \(n\times n\) real symmetric or complex Hermitian matrices, and \(\det(A-\lambda B)\equiv0\).
2. **Hermitian perturbation**:
\[
\tilde{A}-\lambda\tilde{B}=(A + \tau UD_{A}U^{*})-\lambda(B+\tau UD_{B}U^{*})
\]
where \(D_{A}\) and \(D_{B}\) are \(k\times k\) Hermitian matrices such that \(D_{A}-\lambda D_{B}\) is regular, \(U\in\mathbb{C}^{n,k}\) is full - rank, and \(\tau\in\mathbb{R}\) is non - zero.
3. **Eigenvalue decomposition**:
\[
P(A-\lambda B)P^{*}=\begin{bmatrix}
R(\lambda)&0\\
0&S(\lambda)
\end{bmatrix}
\]
where \(R(\lambda)\) is the regular part and \(S(\lambda)\) is the singular part.
4. **Perturbed eigenvalue decomposition**:
\[
P(\tilde{A}-\lambda\tilde{B})P^{*}=\begin{bmatrix}
R(\lambda)&0\\
0&R_{\text{new}}(\lambda)
\end{bmatrix}
\]
where \(R_{\text{new}}(\lambda)\) is regular and all of its eigenvalues are different from those of \(R(\lambda)\).
### Conclusion
By introducing Hermitian perturbations, this paper successfully solves the structure - preservation problem in the eigenvalue calculation of symmetric and Hermitian singular pencils. This method not only retains the sign characteristic of the original pencil but also can effectively handle the singular generalized eigenvalue problem, especially showing advantages in applications such as solving binary polynomial systems.