Matrix inverse problem and its optimal approximation problem for R-skew symmetric matrices

Let R ∈ C^n×n be a nontrivial involution, i.e., R² = I and R ≠ ±I. A matrix A ∈ C^n×n is called R-skew symmetric if RAR = −A. The least-squares solutions of the matrix inverse problem for R-skew symmetric matrices with R∗ = R are firstly derived, then the solvability conditions and the solutions of the matrix inverse problem for R-skew symmetric matrices with R∗ = R are given. The solutions of the corresponding optimal approximation problem with R∗ = R for R-skew symmetric matrices are also derived. At last an algorithm for the optimal approximation problem is given. It can be seen that we extend our previous results [G.X. Huang, F. Yin, Matrix inverse problem and its optimal approximation problem for R-symmetric matrices, Appl. Math. Comput. 189 (2007) 482-489] and the results proposed by Zhou et al. [F.Z. Zhou, L. Zhang, X.Y. Hu, Least-square solutions for inverse problem of centrosymmetric matrices, Comput. Math. Appl. 45 (2003) 1581-1589].     

线性流形上的广义反射矩阵反问题

设 R∈C^m×m 及 S∈C^n×n 是非平凡Hermitian酉矩阵, 即 R^H=R=R^-1≠±I_m ,S^H=S=S^-1≠±I_n.若矩阵 A∈C^m×n 满足 RAS=A, 则称矩阵 A 为广义反射矩阵.该文考虑线性流形上的广义反射矩阵反问题及相应的最佳逼近问题.给出了反问题解的一般表示, 得到了线性流形上矩阵方程AX₂=Z₂, Y₂^H A=W₂^H 具有广义反射矩阵解的充分必要条件, 导出了最佳逼近问题唯一解的显式表示.     

矩阵方程A1Z+ZB1=C1的广义自反最佳逼近解的迭代算法

本文研究了Sylvester复矩阵方程A_1Z+ZB_1=c_1的广义自反最佳逼近解.利用复合最速下降法,提出了一种的迭代算法.不论矩阵方程A_1Z+ZB_1=C_1是否相容,对于任给初始广义自反矩阵Z_0,该算法都可以计算出其广义自反的最佳逼近解.最后,通过两个数值例子,验证了该算法的可行性.     

Matrices A such that AA − AA are nonsingular

In this paper we study the class of square matrices A such that AA^† − A^†A is nonsingular, where A^† stands for the Moore-Penrose inverse of A. Among several characterizations we prove that for a matrix A of order n, the difference AA^† − A^†A is nonsingular if and only if R(A)⊕R(A^∗)=C_n,1, where R(·) denotes the range space. Also we study matrices A such that R(A)^⊥=R(A^∗).     

Oscillation criteria related to integral averaging technique for linear matrix Hamiltonian systems

By employing a generalized Riccati technique and an integral averaging technique, new oscillation criteria are established for the linear matrix Hamiltonian system U′=A(x)U+B(x)V, V′=C(x)U−A∗(x)V under the assumption that all A(x), B(x)=B∗(x)>0, and C(x)=C∗(x) are n×n matrices of real-valued continuous functions on the interval I=[x₀,∞) (−∞<x₀). These criteria extend, improve, and unify a number of existing results and also handle cases not covered by known criteria. Several interesting examples that illustrate the importance of our results are included.     

On the Courant-Fischer theory for Krein spaces

Let J=I_r⊕-I_n-r,0<r<n. An n×n complex matrix A is said to be J-Hermitian if JA=A^∗J. An extension of the classical theory of Courant and Fischer on the Rayleigh ratio of Hermitian matrices is stated for J-Hermitian matrices. Applications to the theory of small oscilations of a mechanical system are presented.     

Minimal ranks of some quaternion matrix expressions with applications

Suppose that p(X, Y) = A − BX − X^(∗)B^(∗) − CYC^(∗) and q(X, Y) = A − BX + X^(∗)B^(∗) − CYC^(∗) are quaternion matrix expressions, where A is persymmetric or perskew-symmetric. We in this paper derive the minimal rank formula of p(X, Y) with respect to pair of matrices X and Y = Y^(∗), and the minimal rank formula of q(X, Y) with respect to pair of matrices X and Y = −Y^(∗). As applications, we establish some necessary and sufficient conditions for the existence of the general (persymmetric or perskew-symmetric) solutions to some well-known linear quaternion matrix equations. The expressions are also given for the corresponding general solutions of the matrix equations when the solvability conditions are satisfied. At the same time, some useful consequences are also developed.     

On (0, 1)-matrices with prescribed row and column sum vectors

Given partitions R and S with the same weight, the Robinson-Schensted-Knuth correspondence establishes a bijection between the class A(R,S) of (0, 1)-matrices with row sum R and column sum S and pairs of Young tableaux of conjugate shapes λ and λ^∗, with S?λ?R^∗. An algorithm for constructing a matrix in A(R,S) whose insertion tableau has a prescribed shape λ, with S?λ?R^∗, is provided. We generalize some recent constructions due to R. Brualdi for the extremal cases λ=S and λ=R^∗.     

Alexandrov’s inequality and conjectures on some Toeplitz matrices

We study determinant inequalities for certain Toeplitz-like matrices over C. For fixed n and N ? 1, let Q be the n × (n + N − 1) zero-one Toeplitz matrix with Q_ij = 1 for 0 ? j − i ? N − 1 and Q_ij = 0 otherwise. We prove that det(QQ^∗) is the minimum of det(RR^∗) over all complex matrices R with the same dimensions as Q satisfying ∣R_ij∣ ? 1 whenever Q_ij = 1 and R_ij = 0 otherwise. Although R has a Toeplitz-like band structure, it is not required to be actually Toeplitz. Our proof involves Alexandrov’s inequality for polarized determinants and its generalizations. This problem is motivated by Littlewood’s conjecture on the minimum 1-norm of N-term exponential sums on the unit circle. We also discuss polarized Bazin-Reiss-Picquet identities, some connections with k-tree enumeration, and analogous conjectured inequalities for the elementary symmetric functions of QQ^∗.     

On single and double Soules matrices

Until now the concept of a Soules basis matrix of sign patternN consisted of an orthogonal matrix R∈R^n,n, generated in a certain way from a positive n-vector, which has the property that for any diagonal matrix Λ = diag(λ₁, … , λ_n), with λ₁ ? ? ? λ_n ? 0, the symmetric matrix A = RΛR^T has nonnegative entries only. In the present paper we introduce the notion of a pair of double Soules basis matrices of sign patternN which is a pair of matrices (P, Q), each in R^n,n, which are not necessarily orthogonal and which are generated in a certain way from two positive vectors, but such that PQ^T = I and such that for any of the aforementioned diagonal matrices Λ, the matrix A = PΛQ^T (also) has nonnegative entries only. We investigate the interesting properties which such matrices A have.As a preamble to the above investigation we show that the iterates, , generated in the course of the QR-algorithm when it is applied to A = RΛR^T, where R is a Soules basis matrix of sign pattern N, are again symmetric matrices generated by the Soules basis matrices R_k of sign pattern N which are themselves modified as the algorithm progresses.Our work here extends earlier works by Soules and Elsner et al.     

矩阵方程AXB+CX~TD=E自反最佳逼近解的迭代算法

本文研究了Sylvester矩阵方程AXB+CXTD=E自反(或反自反)最佳逼近解.利用所提出的共轭方向法的迭代算法,获得了一个结果:不论矩阵方程AXB+CXTD=E是否相容,对于任给初始自反(或反自反)矩阵X1,在有限迭代步内,该算法都能够计算出该矩阵方程的自反(或反自反)最佳逼近解.最后,三个数值例子验证了该算法是有效性的.     

Bifurcation from stability to instability for a free boundary problem modeling tumor growth by Stokes equation

We consider a free boundary problem modeling tumor growth in fluid-like tissue. The model equations include a diffusion equation for the nutrient concentration, and the Stokes equation with a source which represents the proliferation of tumor cells. The proliferation rate μ and the cell-to-cell adhesiveness γ which keeps the tumor intact are two parameters which characterize the “aggressiveness” of the tumor. For any positive radius R there exists a unique radially symmetric stationary solution with radius r=R. For a sequence μ/γ=Mn(R) there exist symmetry-breaking bifurcation branches of solutions with free boundary r=R+εYn,0(θ)+O(ε²) (n even ?2) for small |ε|, where Yn,0 is the spherical harmonic of mode (n,0). Furthermore, the smallest Mn(R), say Mn_∗(R), is such that n_∗=n_∗(R)→∞ as R→∞. In this paper we prove that the radially symmetric stationary solution with R=RS is linearly stable if μ/γ<N_∗(RS,γ) and linearly unstable if μ/γ>N_∗(RS,γ), where N_∗(RS,γ)?Mn_∗(RS), and we prove that strict inequality holds if γ is small or if γ is large. The biological implications of these results are discussed at the end of the paper.     

Unitary part of a contraction

For a contraction A on a Hilbert space H, we define the index j(A) (resp., k(A)) as the smallest nonnegative integer j (resp., k) such that ker(I−Aj∗Aj) (resp., ker(I−Ak*Ak)∩ker(I−AkAk∗)) equals the subspace of H on which the unitary part of A acts. We show that if , then j(A)?n (resp., k(A)?⌈n/2⌉), and the equality holds if and only if A is of class Sn (resp., one of the three conditions is true: (1) A is of class Sn, (2) n is even and A is completely nonunitary with ‖An−2‖=1 and ‖An−1‖<1, and (3) n is even and A=U⊕A^′, where U is unitary on a one-dimensional space and A^′ is of class Sn−1).     

New oscillation criteria for linear matrix Hamiltonian systems

Some new oscillation criteria are established for the matrix linear Hamiltonian system X′=A(t)X+B(t)Y, Y′=C(t)X−A∗(t)Y under the hypothesis: A(t), B(t)=B∗(t)>0, and C(t)=C∗(t) are n×n real continuous matrix functions on the interval [t₀,∞), (−∞<t₀). These results are sharper than some previous results even for self-adjoint second order matrix differential systems.     

Contractive maps on normed linear spaces and their applications to nonlinear matrix equations

In this paper we will give necessary and sufficient conditions under which a map is a contraction on a certain subset of a normed linear space. These conditions are already well known for maps on intervals in R. Using the conditions and Banach’s fixed point theorem we can prove a fixed point theorem for operators on a normed linear space. The fixed point theorem will be applied to the matrix equation X = I_n + A^∗f(X)A, where f is a map on the set of positive definite matrices induced by a real valued map on (0, ∞). This will give conditions on A and f under which the equation has a unique solution in a certain set. We will consider two examples of f in detail. In one example the application of the fixed point theorem is the first step in proving that the equation has a unique positive definite solution under the conditions on A.     

Reynolds operator on functors

Let be an affine R-monoid scheme. We prove that the category of dual functors (over the category of commutative R-algebras) of G-modules is equivalent to the category of dual functors of A^∗-modules. We prove that G is invariant exact if and only if A^∗=R×B^∗ as R-algebras and the first projection A^∗→R is the unit of A. If M is a dual functor of G-modules and w_G?(1,0)∈R×B^∗=A^∗, we prove that M^G=w_G⋅M and M=w_G⋅M⊕(1−w_G)⋅M; hence, the Reynolds operator can be defined on M.     

Products of reflections in the general linear group over a division ring

Let F be a division ring and A?GL_n(F). We determine the smallest integer k such that A admits a factorization A=R₁R₂?R_k?1B, where R₁,…,R_k?1 are reflections and B is such that rank(B?I_n)=1. We find that, apart from two very special exceptional cases, k=rank(A?I_n). In the exceptional cases k is one larger than this rank. The first exceptional case is the matrices A of the form I_m⊕αI_n?m where n?m?2, α≠?1, and α belongs to the center of F. The second exceptional case is the matrices A satisfying (A?I_n)²=0, rank(A?I_n)?2 in the case when char F≠2 only. This result is used to determine, in the case when F is commutative, the length of a matrix A?GL_n(F) with detA=±1 with respect to the set of all reflections in GL_n(F).     

Involutory matrices over finite commutative rings

An n × n matrix A is called involutory iff A²=I_n, where I_n is the n × n identity matrix. This paper is concerned with involutory matrices over an arbitrary finite commutative ring R with identity and with the similarity relation among such matrices. In particular the authors seek a canonical set $\text{C}$ with respect to similarity for the n × n involutory matrices over R—i.e., a set $\text{C}$ of n × n involutory matrices over R with the property that each n × n involutory matrix over R is similar to exactly on matrix in $\text{C}$. Because of the structure of finite commutative rings and because of previous research, they are able to restrict their attention to finite local rings of characteristic a power of 2, and although their main result does not completely specify a canonical set $\text{C}$ for such a ring, it does solve the problem for a special class of rings and shows that a solution to the general case necessarily contains a solution to the classically unsolved problem of simultaneously bringing a sequence A₁,…,A_v of (not necessarily involutory) matrices over a finite field of characteristic 2 to canonical form (using the same similarity transformation on each A_i). (More generally, the authors observe that a theory of similarity fot matrices over an arbitrary local ring, such as the well-known rational canonical theory for matrices over a field, necessarily implies a solution to the simultaneous canonical form problem for matrices over a field.) In a final section they apply their results to find a canonical set for the involutory matrices over the ring of integers modulo 2^m and using this canonical set they are able to obtain a formula for the number of n × n involutory matrices over this ring.     

An alternative approach to unitoidness

In a recent paper [C.R. Johnson, S. Furtado, A generalization of Sylvester’s law of inertia, Linear Algebra Appl. 338 (2001) 287-290], Sylvester’s law of inertia is generalized to any matrix that is ∗-congruent to a diagonal matrix. Such a matrix is called unitoid. In the present paper, an alternative approach to the subject of unitoidness is offered. Specifically, Sylvester’s law of inertia states that a Hermitian n × n matrix of rank r with inertia (p, q, n − r) is ∗-congruent to the direct sum

eⁱ⁰I_p⊕e^iπI_q⊕0I_n-r.