Feed aggregator

Generalization of real interval matrices to other fields

WySR RSS Feed - Mon, 2019-08-12 02:19

An interval matrix is a matrix whose entries are intervals in $\R$. This concept, which has been broadly studied, is generalized to other fields. Precisely, a rational interval matrix is defined to be a matrix whose entries are intervals in $\Q$. It is proved that a (real) interval $p \times q$ matrix with the endpoints of all its entries in $\Q$ contains a rank-one matrix if and only if it contains a rational rank-one matrix, and contains a matrix with rank smaller than $\min\{p,q\}$ if and only if it contains a rational matrix with rank smaller than $\min\{p,q\}$; from these results and from the analogous criterions for (real) inerval matrices, a criterion to see when a rational interval matrix contains a rank-one matrix and a criterion to see when it is full-rank, that is, all the matrices it contains are full-rank, are deduced immediately. Moreover, given a field $K$ and a matrix $\al$ whose entries are subsets of $K$, a criterion to find the maximal rank of a matrix contained in $\al$ is described.

Pure PSVD approach to Sylvester-type quaternion matrix equations

WySR RSS Feed - Thu, 2019-07-25 12:36

In this paper, the pure product singular value decomposition (PSVD) for four quaternion matrices is given. The system of coupled Sylvester-type quaternion matrix equations with five unknowns $X_{i}A_{i}-B_{i}X_{i+1}=C_{i}$ is considered by using the PSVD approach, where $A_{i},B_{i},$ and $C_{i}$ are given quaternion matrices of compatible sizes $(i=1,2,3,4)$. Some necessary and sufficient conditions for the existence of a solution to this system are derived. Moreover, the general solution to this system is presented when it is solvable.

Matrix Shanks Transformations

WySR RSS Feed - Mon, 2019-07-01 13:54

Shanks' transformation is a well know sequence transformation for accelerating the convergence of scalar sequences. It has been extended to the case of sequences of vectors and sequences of square matrices satisfying a linear difference equation with scalar coefficients. In this paper, a more general extension to the matrix case where the matrices can be rectangular and satisfy a difference equation with matrix coefficients is proposed and studied. In the particular case of square matrices, the new transformation can be recursively implemented by the matrix $\varepsilon$-algorithm of Wynn. Then, the transformation is related to matrix Pad\'{e}-type and Pad\'{e} approximants. Numerical experiments showing the interest of this transformation end the paper.

Non-sparse Companion Matrices

WySR RSS Feed - Mon, 2019-07-01 13:54

Given a polynomial $p(z)$, a companion matrix can be thought of as a simple template for placing the coefficients of $p(z)$ in a matrix such that the characteristic polynomial is $p(z)$. The Frobenius companion and the more recently-discovered Fiedler companion matrices are examples. Both the Frobenius and Fiedler companion matrices have the maximum possible number of zero entries, and in that sense are sparse. In this paper, companion matrices are explored that are not sparse. Some constructions of non-sparse companion matrices are provided, and properties that all companion matrices must exhibit are given. For example, it is shown that every companion matrix realization is non-derogatory. Bounds on the minimum number of zeros that must appear in a companion matrix, are also given.

Block GLT Sequences: Matrix Functions and Engineering Application

WySR RSS Feed - Tue, 2019-05-28 16:29

The theory of block generalized locally Toeplitz (GLT) sequences is a powerful apparatus for computing the spectral distribution of block-structured matrices arising from the discretization of differential problems, with a special reference to systems of differential equations (DEs) and to the higher-order finite element or discontinuous Galerkin approximation of both scalar and vectorial DEs. In the present paper, the theory of block GLT sequences is extended by proving that $\{f(A_n)\}_n$ is a block GLT sequence as long as $f$ is continuous and $\{A_n\}_n$ is a block GLT sequence formed by Hermitian matrices. It is also provided a relevant application of this result to the computation of the distribution of the numerical eigenvalues obtained from the higher-order isogeometric Galerkin discretization of second-order variable-coefficient differential eigenvalue problems (a topic of interest not only in numerical analysis but also in engineering).

Cone-constrained rational eigenvalue problems

WySR RSS Feed - Tue, 2019-05-28 16:29

This work deals with the eigenvalue analysis of a rational matrix-valued function subject to complementarity constraints induced by a polyhedral cone $K$. The eigenvalue problem under consideration has the general structure \[ \left(\sum_{k=0}^d \lambda^k A_k + \sum_{k =1}^m \frac{p_k(\lambda)}{q_k(\lambda)} \,B_k\right) x = y , \quad K\ni x \perp y\in K^\ast, \] where $K^\ast$ denotes the dual cone of $K$. The unconstrained version of this problem has been discussed in [Y.F. Su and Z.J. Bai. Solving rational eigenvalue problems via linearization. \emph{SIAM J. Matrix Anal. Appl.}, 32:201--216, 2011.] with special emphasis on the implementation of linearization-based methods. The cone-constrained case can be handled by combining Su and Bai's linearization approach and the so-called facial reduction technique. In essence, this technique consists in solving one unconstrained rational eigenvalue problem for each face of the polyhedral cone $K$.

The Determinant and Complex Gershgorin Circles

WySR RSS Feed - Sat, 2019-05-18 05:45

Each connected component of the Gershgorin circles of a matrix contains exactly as many eigenvalues as circles are involved. Thus, the Minkowski (set) product of all circles contains the determinant if all circles are disjoint. In [S.M. Rump. Bounds for the determinant by Gershgorin circles. Linear Algebra and its Applications, 563:215--219, 2019.], it was proved that statement to be true for real matrices whose circles need not to be disjoint. Moreover, it was asked whether the statement remains true for complex matrices. This note answers that in the affirmative. As a by-product, a parameterization of the outer loop of a Cartesian oval without case distinction is derived.

Pairwise Completely Positive Matrices and Conjugate Local Diagonal Unitary Invariant Quantum States

WySR RSS Feed - Sat, 2019-05-18 05:45

A generalization of the set of completely positive matrices called pairwise completely positive (PCP) matrices is introduced. These are pairs of matrices that share a joint decomposition so that one of them is necessarily positive semidefinite while the other one is necessarily entrywise non-negative. Basic properties of these matrix pairs are explored and several testable necessary and sufficient conditions are developed to help determine whether or not a pair is PCP. A connection with quantum entanglement is established by showing that determining whether or not a pair of matrices is pairwise completely positive is equivalent to determining whether or not a certain type of quantum state, called a conjugate local diagonal unitary invariant state, is separable. Many of the most important quantum states in entanglement theory are of this type, including isotropic states, mixed Dicke states (up to partial transposition), and maximally correlated states. As a specific application of these results, a wide family of states that have absolutely positive partial transpose are shown to in fact be separable.

Subscribe to UW Digital aggregator