Eigenvalue perturbation theory of symplectic,orthogonal, and unitary matrices under generic structured rank one perturbations

We study the perturbation theory of structured matrices under structured rank one perturbations, with emphasis on matrices that are unitary, orthogonal, or symplectic with respect to an indefinite inner product. The rank one perturbations are not necessarily of arbitrary small size (in the sense of norm). In the case of sesquilinear forms, results on selfadjoint matrices can be applied to unitary matrices by using the Cayley transformation, but in the case of real or complex symmetric or skew-symmetric bilinear forms additional considerations are necessary. For complex symplectic matrices, it turns out that generically (with respect to the perturbations) the behavior of the Jordan form of the perturbed matrix follows the pattern established earlier for unstructured matrices and their unstructured perturbations, provided the specific properties of the Jordan form of complex symplectic matrices are accounted for. For instance, the number of Jordan blocks of fixed odd size corresponding to the eigenvalue 1 or ?1 have to be even. For complex orthogonal matrices, it is shown that the behavior of the Jordan structures corresponding to the original eigenvalues that are not moved by perturbations follows again the pattern established earlier for unstructured matrices, taking into account the specifics of Jordan forms of complex orthogonal matrices. The proofs are based on general results developed in the paper concerning Jordan forms of structured matrices (which include in particular the classes of orthogonal and symplectic matrices) under structured rank one perturbations. These results are presented and proved in the framework of real as well as of complex matrices.     

Remarks on Lipschitz properties of matrix groups actions

It is proved that a large class of matrix group actions, including joint similarity and congruence-like actions, as well as actions of the type of matrix equivalence, have local Lipschitz property. Under additional hypotheses, global Lipschitz property is proved. These results are specialized and applied to obtain local Lipschitz property of canonical bases of matrices that are selfadjoint in an indefinite inner product. Real, complex, and quaternionic matrices are considered.     

On the eigenvalues of specially low-rank perturbed matrices

We study the eigenvalues of a matrix A perturbed by a few special low-rank matrices. The perturbation is constructed from certain basis vectors of an invariant subspace of A, such as eigenvectors, Jordan vectors, or Schur vectors. We show that most of the eigenvalues of the low-rank perturbed matrix stayed unchanged from the eigenvalues of A; the perturbation can only change the eigenvalues of A that are related to the invariant subspace. Existing results mostly studied using eigenvectors with full column rank for perturbations, we generalize the results to more general settings. Applications of our results to a few interesting problems including the Google’s second eigenvalue problem are presented.     

On the reduction of the normal Hankel problem to two particular cases

The normal Hankel problem (NHP) is to describe complex matrices that are normal and Hankel at the same time. The available results related to the NHP can be combined into two groups. On the one hand, there are several known classes of normal Hankel matrices. On the other hand, the matrix classes that may contain normal Hankel matrices not belonging to the known classes were shown to admit a parametrization by real 2 × 2 matrices with determinant 1. We solve the NHP for the cases where the characteristic matrix W of the given class has: (a) complex conjugate eigenvalues; (b) distinct real eigenvalues. To obtain a complete solution of the NHP, it remains to analyze two situations: (1) W is the Jordan block of order two for the eigenvalue 1; (2) W is the Jordan block of order two for ?1.     

Stable real invariant semidefinite subspaces and stable factorizations of symmetric rational matrix functions

Stability of certain classes of invariant subspaces of a real matrix which is selfadjoint in a real indefinite inner product, is studied. The study is motivated by, and subsequently applied to, the problem of stability for minimal factorizations of real symmetric rational matrix functions.     

A canonical form for nonderogatory matrices under unitary similarity

A square matrix is nonderogatory if its Jordan blocks have distinct eigenvalues. We give canonical forms for

•

nonderogatory complex matrices up to unitary similarity, and

•

pairs of complex matrices up to similarity, in which one matrix has distinct eigenvalues.

The types of these canonical forms are given by undirected and, respectively, directed graphs with no undirected cycles.     

Eigenvalue perturbation theory of structured real matrices and their sign characteristics under generic structured rank-one perturbations

An eigenvalue perturbation theory under rank-one perturbations is developed for classes of real matrices that are symmetric with respect to a non-degenerate bilinear form, or Hamiltonian with respect to a non-degenerate skew-symmetric form. In contrast to the case of complex matrices, the sign characteristic is a crucial feature of matrices in these classes. The behaviour of the sign characteristic under generic rank-one perturbations is analyzed in each of these two classes of matrices. Partial results are presented, but some questions remain open. Applications include boundedness and robust boundedness for solutions of structured systems of linear differential equations with respect to general perturbations as well as with respect to structured rank perturbations of the coefficients.     

Canonical forms for mixed symmetric-skewsymmetric quaternion matrix pencils

Canonical forms are described for pairs of quaternionic matrices, or equivalently matrix pencils, where one matrix is symmetric and the other matrix is skewsymmetric, under strict equivalence and symmetry respecting congruence. The symmetry is understood in the sense of a fixed involutory antiautomorphism of the skew field of the real quaternions; the involutory antiautomorphism is assumed to be nonstandard, i.e., other than the quaternionic conjugation. Some applications are developed, such as canonical forms for quaternionic matrices under symmetry respecting congruence, and canonical forms for matrices that are skewsymmetric with respect to a nondegenerate symmetric or skewsymmetric quaternion valued inner product.     

Über Matrizen mit Semidefinitem Realteil, Imaginärteil oder Polarem Defekt und Ihre Eigenwerte

We give a new short proof using properties of the field of values to show that

a) a complex matrix with only real eigenvalues is either hermitian or has indefinite imaginary part, and

b) one with only purely imaginary eigenvalues is either skew-hermitian or has indefinite real part, while

c) one whose eigenvalues all have absolute value 1 is either unitary or has indefinite polar defect I—TT^*.

Conversely, every skewsymmetric matrix is the skewsymmetric part of some real matrix that is similar to a real diagonal matrix. The corresponding result for complex matrices is found to be false.     

On two particular cases of solving the normal Hankel problem

The normal Hankel problem is one of characterizing all the complex matrices that are normal and Hankel at the same time. The matrix classes that can contain normal Hankel matrices admit a parameterization by real 2 × 2 matrices with determinant one. Here, the normal Hankel problem is solved in the case where the characteristic matrix of a given class is an order two Jordan block for the eigenvalue 1 or ?1.     

When is a matrix unitary or Hermitian plus low rank?

Hermitian and unitary matrices are two representatives of the class of normal matrices whose full eigenvalue decomposition can be stably computed in quadratic computing complexity once the matrix has been reduced, for instance, to tridiagonal or Hessenberg form. Recently, fast and reliable eigensolvers dealing with low‐rank perturbations of unitary and Hermitian matrices have been proposed. These structured eigenvalue problems appear naturally when computing roots, via confederate linearizations, of polynomials expressed in, for example, the monomial or Chebyshev basis. Often, however, it is not known beforehand whether or not a matrix can be written as the sum of a Hermitian or unitary matrix plus a low‐rank perturbation. In this paper, we give necessary and sufficient conditions characterizing the class of Hermitian or unitary plus low‐rank matrices. The number of singular values deviating from 1 determines the rank of a perturbation to bring a matrix to unitary form. A similar condition holds for Hermitian matrices; the eigenvalues of the skew‐Hermitian part differing from 0 dictate the rank of the perturbation. We prove that these relations are linked via the Cayley transform. Then, based on these conditions, we identify the closest Hermitian or unitary plus rank k matrix to a given matrix A, in Frobenius and spectral norm, and give a formula for their distance from A. Finally, we present a practical iteration to detect the low‐rank perturbation. Numerical tests prove that this straightforward algorithm is effective.     

Eigenvalues of graded matrices and the condition numbers of a multiple eigenvalue

Summary This paper concerns two closely related topics: the behavior of the eigenvalues of graded matrices and the perturbation of a nondefective multiple eigenvalue. We will show that the eigenvalues of a graded matrix tend to share the graded structure of the matrix and give precise conditions insuring that this tendency is realized. These results are then applied to show that the secants of the canonical angles between the left and right invariant of a multiple eigenvalue tend to characterize its behavior when its matrix is slightly perturbed.This work was supported in part by the Air Force Office of Sponsored Research under Contract AFOSR-87-0188     

Eigenvalue perturbation bounds for Hermitian block tridiagonal matrices

We derive new perturbation bounds for eigenvalues of Hermitian matrices with block tridiagonal structure. The main message of this paper is that an eigenvalue is insensitive to blockwise perturbation, if it is well-separated from the spectrum of the diagonal blocks nearby the perturbed blocks. Our bound is particularly effective when the matrix is block-diagonally dominant and graded. Our approach is to obtain eigenvalue bounds via bounding eigenvector components, which is based on the observation that an eigenvalue is insensitive to componentwise perturbation if the corresponding eigenvector components are small. We use the same idea to explain two well-known phenomena, one concerning aggressive early deflation used in the symmetric tridiagonal QR algorithm and the other concerning the extremal eigenvalues of Wilkinson matrices.     

Perturbation of analytic hermitian matrix functions

The behaviour of real eigenvalues of selfadjoint analytic matrix valued functions under small selfadjoint analytic perturbations is studied. Attention is paid mainly to the case when the perturbation is definite (or semidefi-nite). Earlier results of the authors concerning matrix polynomials of first degree are extended to the case of analytic functions.     

Sign patterns that require real, nonreal or pure imaginary eigenvalues

We characterize the n-by-n sign pattern matrices that require all real, all nonreal, and all pure imaginary eigenvalues. Characterization of sign patterns that allow a real eigenvalue and those that allow a nonreal eigenvalue then follow. Some related specialized results and a characterization of sign patterns that allow a positive real eigenvalue are included.     

A Brauer’s theorem and related results

Given a square matrix A, a Brauer’s theorem [Brauer A., Limits for the characteristic roots of a matrix. IV. Applications to stochastic matrices, Duke Math. J., 1952, 19(1), 75–91] shows how to modify one single eigenvalue of A via a rank-one perturbation without changing any of the remaining eigenvalues. Older and newer results can be considered in the framework of the above theorem. In this paper, we present its application to stabilization of control systems, including the case when the system is noncontrollable. Other applications presented are related to the Jordan form of A and Wielandt’s and Hotelling’s deflations. An extension of the aforementioned Brauer’s result, Rado’s theorem, shows how to modify r eigenvalues of A at the same time via a rank-r perturbation without changing any of the remaining eigenvalues. The same results considered by blocks can be put into the block version framework of the above theorem.     

Eigenvalue perturbation theory of classes of structured matrices under generic structured rank one perturbations

udy the perturbation theory of structured matrices under structured rank one perturbations, and then focus on several classes of complex matrices. Generic Jordan structures of perturbed matrices are identified. It is shown that the perturbation behavior of the Jordan structures in the case of singular J-Hamiltonian matrices is substantially different from the corresponding theory for unstructured generic rank one perturbation as it has been studied in [18, 28, 30, 31]. Thus a generic structured perturbation would not be generic if considered as an unstructured perturbation. In other settings of structured matrices, the generic perturbation behavior of the Jordan structures, within the confines imposed by the structure, follows the pattern of that of unstructured perturbations.