15A21 Canonical forms, reductions, classification
Refine
Language
- English (12)
Keywords
- Hamiltonian matrix (4)
- canonical form (4)
- symplectic matrix (4)
- matrix polynomial (3)
- perturbation analysis (3)
- Smith form (2)
- indefinite inner product (2)
- matrix pencil (2)
- palindromic eigenvalue problem (2)
- rank one perturbation (2)
- structured linearization (2)
- $H_\infty$ control (1)
- Brunovsky form (1)
- Hamiltonian real Schur form (1)
- Hessenberg matrix (1)
- Jordan form (1)
- Jordan structure (1)
- Kronecker chain (1)
- Kronecker product (1)
- Lagrangian invariant subspace (1)
- Matrix triples (1)
- QR algorithm (1)
- Spark (1)
- Takagi factorization (1)
- alternating matrix polynomial (1)
- complex Hamiltonian Jordan form (1)
- complex bilinear forms (1)
- complex skew-symmetric matrix (1)
- complex symmetric matrix (1)
- compound matrix (1)
- compressed sensing (1)
- conditional stability (1)
- congruence (1)
- eigenvalue cluster (1)
- elementary divisor (1)
- elementary divisors (1)
- even/odd matrix polynomial (1)
- generic perturbation (1)
- ingular value decomposition (1)
- invariant factor (1)
- invariant polynomial (1)
- invariant polynomials (1)
- linear-quadratic control (1)
- linearization (1)
- matrix factorization (1)
- mutual incoherence (1)
- orthogonal matrix (1)
- palindromic matrix polynomial (1)
- palindromic pencil (1)
- restricted isometry property (1)
- singular matrix polynomial (1)
- skew-Hamiltonian matrix (1)
- sparse solution of linear systems (1)
- stability (1)
- staircase form (1)
- structured (1)
- structured SVD (1)
- structured matrices (1)
- trimmed linearization (1)
- unitary congruence (1)
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.
Many applications give rise to matrix polynomials whose coefficients have
a kind of reversal symmetry, a structure we call palindromic.
Several properties of scalar palindromic polynomials are derived,
and together with properties of compound matrices, used to
establish the Smith form of regular and singular T-palindromic matrix polynomials,
over arbitrary fields.
The invariant polynomials are shown to
inherit palindromicity,
and their structure is described in detail.
Jordan structures of palindromic matrix polynomials are characterized,
and necessary conditions for the
existence of structured linearizations established.
In the odd degree case, a constructive procedure for building
palindromic linearizations shows that the necessary conditions are sufficient as well.
The Smith form for *-palindromic polynomials is also analyzed. Finally, results for palindromic matrix polynomials over fields of
characteristic two are presented.
We study 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
is substantially different from the corresponding theory for unstructured generic
rank one perturbations.
Alternating matrix polynomials, that is, polynomials whose coefficients
alternate between symmetric and skew-symmetric matrices,
generalize the notions of even and odd scalar polynomials.
We investigate the Smith forms of alternating matrix polynomials,
showing that each invariant factor is an even or odd scalar polynomial.
Necessary and sufficient conditions
are derived for a given Smith form to be that of an alternating matrix polynomial.
These conditions allow a characterization of the possible Jordan structures
of alternating matrix polynomials,
and also lead to necessary and sufficient conditions
for the existence of structure-preserving strong linearizations.
Most of the results are applicable to singular as well
as regular matrix polynomials.
A generalization of the method of Chu, Liu and Mehrmann
for the computation of the Hamiltonian real Schur form is presented.
The new method avoids some of the difficulties that may arise when
a Hamiltonian matrix has tightly clustered groups of eigenvalues.
A detailed analysis of the method is presented and several numerical examples demonstrate the superior behavior of the method.
Canonical forms for matrix triples $(A,G,\hat G)$, where
$A$ is arbitrary rectangular and $G$, $\hat G$ are either real symmetric
or skew symmetric, or complex Hermitian or skew Hermitian, are derived.
These forms generalize classical product Schur forms as well as
singular value decompositions.
An new proof for the complex case is given, where there is no need to
distinguish whether $G$ and $\hat G$ are Hermitian or skew Hermitian.
This proof is independent from the results in Bolschakov/Reichstein 1995, where
a similar canonical form has been obtained for the complex case,
and it allows generalization to the real case. Here,
the three cases, i.e., that
$G$ and $\hat G$ are both symmetric, both skew symmetric or one each,
are treated separately.
Three properties of matrices: the spark, the mutual incoherence and the restricted isometry property have recently been introduced in the context of compressed sensing. We study these properties for matrices that are Kronecker products and show how these properties relate to those of the factors. For the mutual incoherence we also
discuss results for sums of Kronecker products.
The classical singular value decomposition for a matrix $A\in\Cmn$ is a
canonical form for $A$ that also displays the eigenvalues
of the Hermitian matrices $AA^\ast$ and $A^\ast A$. In this paper, we develop
a corresponding decomposition for $A$ that provides the Jordan canonical forms
for the complex symmetric matrices $AA^T$ and $A^TA$. More generally, we consider
the matrix triple $(A,G_1,G_2)$, where $G_1\in\CC{m}, G_2\in\CC{n}$
are invertible and either complex symmetric and complex skew-symmetric, and we
provide a canonical form under transformations of the form
$(A,G_1,G_2)\mapsto(X^T A Y, X^T G_1X, Y^T G_2Y)$, where $X,Y$ are nonsingular.
We discuss the eigenvalue problem for
general and structured matrix polynomials which may
be singular and may have eigenvalues at infinity.
We derive staircase
condensed forms that allow deflation of the infinite eigenvalue and
singular structure of the matrix polynomial.
The remaining reduced order staircase form leads to
new types of linearizations which determine the finite eigenvalues and
and corresponding eigenvectors. The new linearizations
also simplify the construction of structure preserving linearizations.
In this paper we develop a QR-like algorithm for the palindromic eigenvalue problem $Ax=\lambda A^\adj x$.
We will discuss the two cases that $A^\adj$ denotes the transpose or the conjugate transpose of $A\in\C^{n,n}$.
It is shown that this so-called palindromic QR iteration is equivalent to applying the standard QR algorithm to $A^{-\adj}A$.
Also the concepts of deflation, shifting, and exploiting the invariance of a Hessenberg-type form are adapted.
Moreover, we analyze the problem of reducing a general square matrix to the mentioned Hessenberg-type form
and establish analogies to the Hamiltonian eigenvalue problem.
Finally, we present concrete Hessenberg-type reduction algorithms for special cases.