AN INVERSE EIGENVALUE PROBLEM FOR JACOBI MATRICES

In this paper, we discuss an inverse eigenvalue problem for constructing a 2n × 2n Jacobi matrix T such that its 2n eigenvalues are given distinct real values and its leading principal submatrix of order n is a given Jacobi matrix. A new sufficient and necessary condition for the solvability of the above problem is given in this paper. Furthermore, we present a new algorithm and give some numerical results.     

Relationship Between the Stiffly Weighted Pseudoinverse and Multi-Level Constrained Rseudoinverse

It is known that for a given matrix A of rank r, and a set D of positive diagonal matrices, supw∈D‖（W^1/2A) W^1/2‖ = (miniσ (A^(i))^-1, in which (A^(i) is a submatrix of A formed with r = (rank(A)) rows of A, such that (A^(i) has full row rank r. In many practical applications this value is too large to be used. In this paper we consider the case that both A and W(∈D) are fixed with W severely stiff. We show that in this case the weighted pseudoinverse (W^1/2‖A) W^1/2‖ is close to a multilevel constrained weighted pseudoinverse therefore ‖(W^1/2A) W^1/‖2 is uniformly bounded.We also prove that in this case the solution set the stiffly weighted least squares problem is close to that of corresponding multi-level constrained least squares problem.     

用Levenberg-Marquardt类的投影收缩方法解运输问题

For solving linear variational inequalities (LVI), the projection and contraction method of Levenberg-Marquardt type needs less iterations than an elementary projection and contraction method. However, the method of Levenberg-Marquardt type has to calculate the inverse of a matrix and hence it is unsuitable for large problems. In this paper, using the special structure of the constraint matrix, we present a PC method of Levenberg-Marquardt type for LVI arising from transportation problem without calculating any inverse matrices.Several computational experiments are presentded to indicate that the methods is good for solving the transportation problem.     

The Inverse Problem of Centrosymmetric Matrices with a Submatrix Constraint

By using Moore-Penrose generalized inverse and the general singular value decomposition of matrices, this paper establishes the necessary and sufficient conditions for the existence of and the expressions for the centrosymmetric solutions with a submatrix constraint of matrix inverse problem AX = B. In addition, in the solution set of corresponding problem, the expression of the optimal approximation solution to a given matrix is derived.     

Matrix Completions and Chordal Graphs

In a matrix-completion problem the aim is to specifiy the missing entries of a matrix in order to produce a matrix with particular properties. In this paper we survey results concerning matrix-completion problems where we look for completions of various types for partial matrices supported on a given pattern. We see that the existence of completions of the required type often depends on the chordal properties of graphs associated with the pattern.     

三对角逆M-矩阵

In this paper we study a class of inverse M-matrices:tridiagonal inverse M-matrices,Graph theory is used to discuss the structure and properties of tridiagonal inverse M-matrices,A sufficient and necessary condtion for a nonnegative tridiagonal matrix to be an inverse M-matrix is given.Finally,it is proved that the set of the inverses of M-matrices with unipathic is closed under Hadamard product.     

ON HERMITIAN POSITIVE DEFINITE SOLUTION OF NONLINEAR MATRIX EQUATION X＋A^＊X^-2A=Q

Based on the fixed-point theory, we study the existence and the uniqueness of the maximal Hermitian positive definite solution of the nonlinear matrix equation X＋A^＊X^-2A=Q, where Q is a square Hermitian positive definite matrix and A＊ is the conjugate transpose of the matrix A. We also demonstrate some essential properties and analyze the sensitivity of this solution. In addition, we derive computable error bounds about the approximations to the maximal Hermitian positive definite solution of the nonlinear matrix equation X＋A^＊X^-2A=Q. At last, we further generalize these results to the nonlinear matrix equation X＋A^＊X^-nA=Q, where n≥2 is a given positive integer.     

The Core-EP,Weighted Core-EP Inverse of Matrices and Constrained Systems of Linear Equations

We study the constrained systemof linear equations Ax=b,x∈R(A^k)for A∈C^n×nand b∈Cn,k=Ind(A).When the system is consistent,it is well known that it has a unique A^Db.If the system is inconsistent,then we seek for the least squares solution of the problem and consider min_x∈R(A^k)||b?Ax||2,where||·||2 is the 2-norm.For the inconsistent system with a matrix A of index one,it was proved recently that the solution is A^■b using the core inverse A^■of A.For matrices of an arbitrary index and an arbitrary b,we show that the solution of the constrained system can be expressed as A^■b where A^■is the core-EP inverse of A.We establish two Cramer’s rules for the inconsistent constrained least squares solution and develop several explicit expressions for the core-EP inverse of matrices of an arbitrary index.Using these expressions,two Cramer’s rules and one Gaussian elimination method for computing the core-EP inverse of matrices of an arbitrary index are proposed in this paper.We also consider the W-weighted core-EP inverse of a rectangular matrix and apply the weighted core-EP inverse to a more general constrained system of linear equations.     

反中心对称矩阵的广义特征值反问题

Given matrix X and diagonal matrix A , the anti-centrosymmetric solutions (A, B) and its optimal approximation of inverse generalized eigenvalue problem AX = BXA have been considered. The general form of such solutions is given and the expression of the optimal approximation solution to a given matrix is derived. The algorithm and one numerical example for solving optimal approximation solution are included.     

Norm estimates of w-circulant operator matrices and isomorphic operators for w-circulant algebra

An n × n ω-circulant matrix which has a specific structure is a type of important matrix. Several norm equalities and inequalities are proved for ω-circulant operator matrices with ω = e~(iθ)(0≤θ 2π) in this paper. We give the special cases for norm equalities and inequalities, such as the usual operator norm and the Schatten p-norms. Pinching type inequality is also proposed for weakly unitarily invariant norms. Meanwhile,we present that the set of ω-circulant matrices with complex entries has an idempotent basis. Based on this basis, we introduce an automorphism on the ω-circulant algebra and then show different operators on linear vector space that are isomorphic to the ω-circulant algebra. The function properties, other idempotent bases and a linear involution are discussed for ω-circulant algebra. These results are closely related to the special structure of ω-circulant matrices.     

The solvability conditions for inverse eigenproblem of symmetric and anti‐persymmetric matrices and its approximation

The problem of generating a matrix A with specified eigen‐pair, where A is a symmetric and anti‐persymmetric matrix, is presented. An existence theorem is given and proved. A general expression of such a matrix is provided. We denote the set of such matrices by ??????_Eⁿ. The optimal approximation problem associated with ??????_Eⁿ is discussed, that is: to find the nearest matrix to a given matrix A* by A∈??????_Eⁿ. The existence and uniqueness of the optimal approximation problem is proved and the expression is provided for this nearest matrix. Copyright © 2002 John Wiley & Sons, Ltd.     

The representation and perturbation of theW-weighted Drazin inverse

LetA andE bem x n matrices andW an n xm matrix, and letA _d,W denote the W-weighted Drazin inverse ofA. In this paper, a new representation of the W-weighted Drazin inverse ofA is given. Some new properties for the W-weighted Drazin inverseA _d,W and B_d,W are investigated, whereB =A+E. In addition, the Banach-type perturbation theorem for the W-weighted Drazin inverse ofA andB are established, and the perturbation bounds for ∥B_d,W∥ and ∥B_{d, W}, -A_d,W∥/∥A_d,W∥ are also presented. WhenA andB are square matrices andW is identity matrix, some known results in the literature related to the Drazin inverse and the group inverse are directly reduced by the results in this paper as special cases.     

Shorted operators of partitioned matrices and application

In this paper, our main objective is to study the effect of appending/deleting a column/row on the shorted operators. It turns out that for matrices A and B for which the shorted operator S(A|B) exists, S(A₁|B₁) of the matrix A₁=[A:a] with respect to the matrix B₁=[B:b], when it exists, is obtained by appending a suitable column to S(A|B). Moreover, if S(A₁|B₁) exists, then S(A|B) exists and is obtained from S(A₁|B₁) by dropping its last column. In the process, we study the effect of appending/deleting a column/row on the space pre-order and the parallel sum of parallel summable matrices. Finally, we specialize to the case of and matrices and study the effect of bordering (by an additional column and a row) on the shorted operator. We conclude the paper with an application to Linear Models with singular dispersion structure.     

Computing the nearest stable matrix pairs

In this paper, we study the nearest stable matrix pair problem: given a square matrix pair (E,A), minimize the Frobenius norm of (Δ_E,Δ_A) such that (E+Δ_E,A+Δ_A) is a stable matrix pair. We propose a reformulation of the problem with a simpler feasible set by introducing dissipative Hamiltonian matrix pairs: A matrix pair (E,A) is dissipative Hamiltonian if A=(J−R)Q with skew‐symmetric J, positive semidefinite R, and an invertible Q such that Q^TE is positive semidefinite. This reformulation has a convex feasible domain onto which it is easy to project. This allows us to employ a fast gradient method to obtain a nearby stable approximation of a given matrix pair.     

The error bound of the perturbation of the Drazin inverse

Let A and E be n×n matrices and B = A + E. Denote the Drazin inverse of A by A^D. In this paper we give an upper bound for the relative error ∥B^D ? A^D∥/∥A^D∥₂ and a lower bound for ∥B^D∥₂ under certain circumstances. The continuity properties and the derivative of the Drazin inverse are also considered.     

The Lyapunov order for real matrices

The real Lyapunov order in the set of real n×n matrices is a relation defined as follows: A?B if, for every real symmetric matrix S, SB+B^tS is positive semidefinite whenever SA+A^tS is positive semidefinite. We describe the main properties of the Lyapunov order in terms of linear systems theory, Nevenlinna-Pick interpolation and convexity.     

Amalgams of free inverse semigroups

We study inverse semigroup amalgams of the formS * _U T whereS andT are free inverse semigroups andU is an arbitrary finitely generated inverse subsemigroup ofS andT. We make use of recent work of Bennett to show that the word problem is decidable for any such amalgam. This is in contrast to the general situation for semigroup amalgams, where recent work of Birget, Margolis and Meakin shows that the word problem for a semigroup amalgamS * _U T is in general undecidable, even ifS andT have decidable word problem,U is a free semigroup, and the membership problem forU inS andT is decidable. We also obtain a number of results concerning the structure of such amalgams. We obtain conditions for theD-classes of such an amalgam to be finite and we show that the amalgam is combinatorial in such a case. For example every one-relator amalgam of this type has finiteD-classes and is combinatorial. We also obtain information concerning when such an amalgam isE-unitary: for example every one relator amalgam of the formInv<A ∪B :u =v > whereA andB are disjoint andu (resp.v) is a cyclically reduced word overA ∪A ⁻¹ (resp.B ∪B ⁻¹) isE-unitary. Research of all authors supported by a grant from the Italian CNR. The first and third authors’ research was partially supported by MURST. The second author’s research was also partially supported by NSF and the Center for Communication and Information Science of the University of Nebraska at Lincoln.     

The inverse problem of bisymmetric matrices with a submatrix constraint

An n×n real matrix A is called a bisymmetric matrix if A=A^T and A=S_nAS_n, where S_n is an n×n reverse unit matrix. This paper is mainly concerned with solving the following two problems: Problem I Given n×m real matrices X and B, and an r×r real symmetric matrix A₀, find an n×n bisymmetric matrix A such that where A([1: r]) is a r×r leading principal submatrix of the matrix A. Problem II Given an n×n real matrix A^*, find an n×n matrix  in S_E such that where ∥·∥ is Frobenius norm, and S_E is the solution set of Problem I. The necessary and sufficient conditions for the existence of and the expressions for the general solutions of Problem I are given. The explicit solution, a numerical algorithm and a numerical example to Problem II are provided. Copyright © 2003 John Wiley & Sons, Ltd.     

Exponential sums and rectangular partitions

Let F denote a finite field with q=p^f elements, and let σ(A) equal the trace of the square matrix A. This paper evaluates exponential sums of the form S(E,X₁,…,X_n)e{?σ(CX₁?X_nE)}, where S(E,X₁,…,X_n) denotes a summation over all matrices E,X₁,…,X_n of appropriate sizes over F, and C is a fixed matrix. This evaluation is then applied to the problem of counting ranked solutions to matrix equations of the form U₁?U_αA+DV₁?V_β=B where A,B,D are fixed matrices over F.     

On the representability of totally unimodular matrices on bidirected graphs

We present a polynomial time algorithm to construct a bidirected graph for any totally unimodular matrix B by finding node-edge incidence matrices Q and S such that QB=S. Seymour’s famous decomposition theorem for regular matroids states that any totally unimodular (TU) matrix can be constructed through a series of composition operations called k-sums starting from network matrices and their transposes and two compact representation matrices B₁,B₂ of a certain ten element matroid. Given that B₁,B₂ are binet matrices we examine the k-sums of network and binet matrices. It is shown that thek-sum of a network and a binet matrix is a binet matrix, but binet matrices are not closed under this operation for k=2,3. A new class of matrices is introduced, the so-called tour matrices, which generalise network, binet and totally unimodular matrices. For any such matrix there exists a bidirected graph such that the columns represent a collection of closed tours in the graph. It is shown that tour matrices are closed under k-sums, as well as under pivoting and other elementary operations on their rows and columns. Given the constructive proofs of the above results regarding the k-sum operation and existing recognition algorithms for network and binet matrices, an algorithm is presented which constructs a bidirected graph for any TU matrix.