A perturbation result for the asymptotic behavior of matrix powers

We prove a perturbation result for the asymptotic behavior of the sequence (A ⁿ c) _nN , whereAG|(d), the space of invertibled×d matrices, andc ^d .     

The Rayleigh-Ritz method, refinement and Arnoldi process for periodic matrix pairs

We extend the Rayleigh-Ritz method to the eigen-problem of periodic matrix pairs. Assuming that the deviations of the desired periodic eigenvectors from the corresponding periodic subspaces tend to zero, we show that there exist periodic Ritz values that converge to the desired periodic eigenvalues unconditionally, yet the periodic Ritz vectors may fail to converge. To overcome this potential problem, we minimize residuals formed with periodic Ritz values to produce the refined periodic Ritz vectors, which converge under the same assumption. These results generalize the corresponding well-known ones for Rayleigh-Ritz approximations and their refinement for non-periodic eigen-problems. In addition, we consider a periodic Arnoldi process which is particularly efficient when coupled with the Rayleigh-Ritz method with refinement. The numerical results illustrate that the refinement procedure produces excellent approximations to the original periodic eigenvectors.     

The perturbed laplacian matrix of a graph

For a graph G, we define its perturbed Laplacian matrix as D?A(G) where A(G) is the adjacency matrix of G and D is an arbitrary diagonal matrix. Both the Laplacian matrix and the negative of the adjacency matrix are special instances of the perturbed Laplacian. Several well-known results, contained in the classical work of Fiedler and in more recent contributions of other authors are shown to be true, with suitable modifications, for the perturbed Laplacian. An appropriate generalization of the monotonicity property of a Fiedler vector for a tree is obtained. Some of the results are applied to interval graphs.     

On the matrix algebra of elliptic biquaternions

In this study, we introduce the concept of elliptic biquaternion matrices. Firstly, we obtain elliptic matrix representations of elliptic biquaternion matrices and establish a universal similarity factorization equality for elliptic biquaternion matrices. Afterwards, with the aid of these representations and this equality, we obtain various results on some basic topics such as generalized inverses, eigenvalues and eigenvectors, determinants, and similarity of elliptic biquaternion matrices. These valuable results may be useful for developing a perfect theory on matrix analysis over elliptic biquaternion algebra in the future.     

Optimal multilevel matrix algebra operators

We study the optimal Frobenius operator in a general matrix vector space and in particular in the multilevel trigonometric matrix vector spaces, by emphasizing both the algebraic and geometric properties. These general results are used to extend the Korovkin matrix theory for the approximation of block Toeplitz matrices via trigonometric vector spaces. The abstract theory is then applied to the analysis of the approximation properties of several sine and cosine based vector spaces. Few numerical experiments are performed to give evidence of the theoretical results.     

Commuting pairs and triples of matrices and related varieties

In this note, we show that the set of all commuting d-tuples of commuting n×n matrices that are contained in an n-dimensional commutative algebra is a closed set, and therefore, Gerstenhaber's theorem on commuting pairs of matrices is a consequence of the irreduciblity of the variety of commuting pairs. We show that the variety of commuting triples of 4×4 matrices is irreducible. We also study the variety of n-dimensional commutative subalgebras of M_n(F), and show that it is irreducible of dimension n²−n for n4, but reducible, of dimension greater than n²−n for n7.     

Computing eigenspaces with specified eigenvalues of a regular matrix pair (A,B) and condition estimation: theory,algorithms and software

Theory, algorithms and LAPACK-style software for computing a pair of deflating subspaces with specified eigenvalues of a regular matrix pair (A, B) and error bounds for computed quantities (eigenvalues and eigenspaces) are presented. Thereordering of specified eigenvalues is performed with a direct orthogonal transformation method with guaranteed numerical stability. Each swap of two adjacent diagonal blocks in the real generalized Schur form, where at least one of them corresponds to a complex conjugate pair of eigenvalues, involves solving a generalized Sylvester equation and the construction of two orthogonal transformation matrices from certain eigenspaces associated with the diagonal blocks. The swapping of two 1×1 blocks is performed using orthogonal (unitary) Givens rotations. Theerror bounds are based on estimates of condition numbers for eigenvalues and eigenspaces. The software computes reciprocal values of a condition number for an individual eigenvalue (or a cluster of eigenvalues), a condition number for an eigenvector (or eigenspace), and spectral projectors onto a selected cluster. By computing reciprocal values we avoid overflow. Changes in eigenvectors and eigenspaces are measured by their change in angle. The condition numbers yield bothasymptotic andglobal error bounds. The asymptotic bounds are only accurate for small perturbations (E, F) of (A, B), while the global bounds work for all (E, F.) up to a certain bound, whose size is determined by the conditioning of the problem. It is also shown how these upper bounds can be estimated. Fortran 77software that implements our algorithms for reordering eigenvalues, computing (left and right) deflating subspaces with specified eigenvalues and condition number estimation are presented. Computational experiments that illustrate the accuracy, efficiency and reliability of our software are also described.     

A topological invariant for pairs of maps

In this paper we develop the notion of contact orders for pairs of continuous self-maps (f, g) from ℝⁿ, showing that the set Con(f, g) of all possible contact orders between f and g is a topological invariant (we remark that Con(f, id) = Per(f)). As an interesting application of this concept, we give sufficient conditions for the graphs of two continuous self-maps from ℝ intersect each other. We also determine the ordering of the sets Con(f, 0) and Con(f, h), for h ∈ Hom(ℝ) such that f ∘ h = h ∘ f. For this latter set we obtain a generalization of Sharkovsky’s theorem.     

Function of a square matrix with repeated eigenvalues

An analytical function f(A) of an arbitrary n×n constant matrix A is determined and expressed by the “fundamental formula”, the linear combination of constituent matrices. The constituent matrices Z_kh, which depend on A but not on the function f(s), are computed from the given matrix A, that may have repeated eigenvalues. The associated companion matrix C and Jordan matrix J are then expressed when all the eigenvalues with multiplicities are known. Several other related matrices, such as Vandermonde matrix V, modal matrix W, Krylov matrix K and their inverses, are also derived and depicted as in a 2-D or 3-D mapping diagram. The constituent matrices Z_kh of A are thus obtained by these matrices through similarity matrix transformations. Alternatively, efficient and direct approaches for Z_kh can be found by the linear combination of matrices, that may be further simplified by writing them in “super column matrix” forms. Finally, a typical example is provided to show the merit of several approaches for the constituent matrices of a given matrix A.     

线性互补问题中几类特殊矩阵与半正定矩阵之间的关系

线性互补问题中特殊矩阵M 的性质是线性互补问题中研究的重要部分之一,本文深入研究了Cf0矩阵与半正定矩阵、子正定矩阵与半正定矩阵之间的关系,并且得到了特殊矩阵是半正定矩阵的一些充分条件。     

A CLASS OF SINGULARLY PERTURBED INITIAL BOUNDARY PROBLEM FOR REACTION DIFFUSION EQUATION

The singularly perturbed initial boundary value problem for a class of reaction diffusion equation is considered. Under appropriate conditions, the existence-uniqueness and the asymptotic behavior of the solution are showed by using the fixed-point theorem.     

一类半线性椭圆型方程奇摄动Robin边值问题

The singularly perturbed Robin boundary value problems for the semilinear elliptic equation are considered. Under suitable conditions and by using the fixed point theorem the existence ,uniqueness and asymptotic behavior of solution for the boundary value problems are studied.     

Perturbation, extraction and refinement of invariant pairs for matrix polynomials

Generalizing the notion of an eigenvector, invariant subspaces are frequently used in the context of linear eigenvalue problems, leading to conceptually elegant and numerically stable formulations in applications that require the computation of several eigenvalues and/or eigenvectors. Similar benefits can be expected for polynomial eigenvalue problems, for which the concept of an invariant subspace needs to be replaced by the concept of an invariant pair. Little has been known so far about numerical aspects of such invariant pairs. The aim of this paper is to fill this gap. The behavior of invariant pairs under perturbations of the matrix polynomial is studied and a first-order perturbation expansion is given. From a computational point of view, we investigate how to best extract invariant pairs from a linearization of the matrix polynomial. Moreover, we describe efficient refinement procedures directly based on the polynomial formulation. Numerical experiments with matrix polynomials from a number of applications demonstrate the effectiveness of our extraction and refinement procedures.     

The efficient evaluation of the hypergeometric function of a matrix argument

We present new algorithms that efficiently approximate the hypergeometric function of a matrix argument through its expansion as a series of Jack functions. Our algorithms exploit the combinatorial properties of the Jack function, and have complexity that is only linear in the size of the matrix.

     

On a class of matrices generated by certain generalized permutation matrices and applications

ABSTRACT

In this paper, we study a particular class of matrices generated by generalized permutation matrices corresponding to a subgroup of some permutation group. As applications, we first present a technique from which we can get closed formulas for the roots of many families of polynomial equations with degree between 5 and 10, inclusive. Then, we describe a tool that shows how to find solutions to Fermat's last theorem and Beal's conjecture over the square integer matrices of any dimension. Finally, simple generalizations of some of the concepts in number theory to integer square matrices are presented.     

On perturbations of matrix pencils with real spectra. II

A well-known result on spectral variation of a Hermitian matrix due to Mirsky is the following: Let and be two Hermitian matrices, and let and be their eigenvalues arranged in ascending order. Then for any unitarily invariant norm . In this paper, we generalize this to the perturbation theory for diagonalizable matrix pencils with real spectra. The much studied case of definite pencils is included in this.

     

The contragredient equivalence: Application to solve some matrix systems

In this paper we apply the contragredient equivalence to solve two matrix systems. Firstly, we characterize and build all possible solutions of the matrix system P = XY, Q = YX, giving a recursive formula for the number of contragrediently nonequivalent solutions. And, secondly, we find the solution of the matrix system AX = YC, BY = XD.