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.     

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.     

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 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.     

Some explicit formulas for the polynomial decomposition of the matrix exponential and applications

This paper is devoted to the study of some formulas for polynomial decomposition of the exponential of a square matrix A. More precisely, we suppose that the minimal polynomial M_A(X) of A is known and has degree m. Therefore, e^tA is given in terms of P₀(A),…,P_m−1(A), where the P_j(A) are polynomials in A of degree less than m, and some explicit analytic functions. Examples and applications are given. In particular, the two cases m=5 and m=6 are considered.     

Exact slow-fast decomposition of the singularly perturbed matrix differential Riccati equation

In this paper the Hamiltonian matrix formulation of the Riccati equation is used to derive the reduced-order pure-slow and pure-fast matrix differential Riccati equations of singularly perturbed systems. These pure-slow and pure-fast matrix differential Riccati equations are obtained by decoupling the singularly perturbed matrix differential Riccati equation of dimension n₁+n₂ into the pure-slow regular matrix differential Riccati equation of dimension n₁ and the pure-fast stiff matrix differential Riccati equation of dimension n₂. A formula is derived that produces the solution of the original singularly perturbed matrix differential Riccati equation in terms of solutions of the pure-slow and pure-fast reduced-order matrix differential Riccati equations and solutions of two reduced-order initial value problems. In addition to its theoretical importance, the main result of this paper can also be used to implement optimal filtering and control schemes for singularly perturbed linear time-invariant systems independently in pure-slow and pure-fast time scales.     

A compound matrix algorithm for the computation of the Smith form of a polynomial matrix

In the present paper is presented a numerical method for the exact reduction of a singlevariable polynomial matrix to its Smith form without finding roots and without applying unimodular transformations. Using the notion of compound matrices, the Smith canonical form of a polynomial matrixM(s)^nxn[s] is calculated directly from its definition, requiring only the construction of all thep-compound matricesC _p (M(s)) ofM(s), 1<pn. This technique produces a stable and accurate numerical algorithm working satisfactorily for any polynomial matrix of any degree.     

线性流形上两类矩阵反问题的最小二乘解

给定广义自反矩阵R,S,即R=R=R-1,S=S=S-1,若复矩阵X满足条件RXS=X（或RXS=X）,则称其为（R,S）-对称矩阵（或（R,S）-斜对称矩阵）.分别讨论了线性流形上（R,S）-对称矩阵和（R,S）-斜对称矩阵约束下矩阵方程MZN=E的最小二乘问题,得到了通解表达式.     

The bialternate matrix product revisited

The well known bialternate product of two square matrices is re-examined together with another matrix product defined by means of the permanent function and having similar properties. Old and new results concerning both products are presented in a unified manner. A simple and elegant relation with the Kronecker product of matrices is also given.