Limit points of the iterative scaling procedure

The iterative scaling procedure (ISP) is an algorithm which computes a sequence of matrices, starting from some given matrix. The objective is to find a matrix ’proportional’ to the given matrix, having given row and column sums. In many cases, for example if the initial matrix is strictly positive, the sequence is convergent. It is known that the sequence has at most two limit points. When these are distinct, convergence to these two points can be slow. We give an efficient algorithm which finds the limit points, invoking the ISP only on subproblems for which the procedure is convergent.     

An ordering of dependence for contingency tables

A majorization ordering is defined on matrices with the same row and column sums. This ordering is used as an ordering of dependence for contingency tables. Results are derived for maximal and minimal matrices with respect to the majorization ordering. This theory can be used to maximize and minimize Schur concave functions defined over matrices, when there are row and column sum constraints; in this paper, it is applied to the algorithm of Mehta and Patel (1983) for finding the P-value of Fisher's exact test.     

四元数矩阵方程AXB+CXD=E的广义延拓解

本文在四元数体上讨论矩阵方程AXB+CXD=E的广义行（列）共轭延拓解问题.利用四元数矩阵的复与实分解,以及广义共轭延拓矩阵的结构特点,借助矩阵Kronecker积,把约束四元数矩阵方程转化为实数域上无约束方程,从而得到该方程具有广义行（列）共轭延拓解的充要条件及其通解表达式.最后通过数值算例说明所给算法的可行性.     

Modular Hadamard martrices and related designs

A square matrix with entries ± 1 is called a modular Hadamard matrix if the inner product of each two distinct row vectors is a multiple of some fixed (positive) integer. This paper initiates the study of modular Hadamard matrices and the combinatorial designs associated with them. The related combinatorial designs are the main concern of this paper; some results dealing with the existence and construction of modular Hadamard matrices will be included in a later paper.     

On permanents and the zeros of rook polynomials

The concept of rook polynomial of a “chessboard” may be generalized to the rook polynomial of an arbitrary rectangular matrix. A conjecture that the rook polynomials of “chessboards” have only real zeros is thus carried over to the rook polynomials of nonnegative matrices. This paper proves these conjectures, and establishes interlacing properties for the zeros of the rook polynomials of a positive matrix and the matrix obtained by striking any one row or any one column.     

Fictitious play applied to sequences of games and discounted stochastic games

In this paper, we show that the iterative method of Brown and Robinson, for solving a matrix game, is also applicable to a converging sequence of matrices, where the players choose at staget a row and a column of thet-th matrix in the sequence. As an application of this result, we describe a new solution method for discounted stochastic games with finite state and action spaces.     

关于矩阵的行最简形以及两个正交矩阵之和的几个注记

本文证明了每个ｍ×ｎ矩阵Ａ的行最简形矩阵必唯一．本文还给出了两个正交矩阵之和仍是正交矩阵的两个充分必要条件．     

A PERTURBATION ANALYSIS FOR THE PROJECTION OF A STIFFLY SCALED MATRIX

In this paper we study the perturbation bound of the projection (WA )(WA ), where both the matrices A and W are given with W positive diagonal and severely stiff. When the perturbed matrix A=A δA satisfy several row rank preserving conditions, we derive a new perturbation bound of the projection.     

Recognizing Balanceable Matrices

A 0/±1 matrix is balanced if it does not contain a square submatrix with exactly two nonzero entries per row and per column in which the sum of all entries is 2 modulo 4. A 0/1 matrix is balanceable if its nonzero entries can be signed ±1 so that the resulting matrix is balanced. A signing algorithm due to Camion shows that the problems of recognizing balanced 0/±1 matrices and balanceable 0/1 matrices are equivalent. Conforti, Cornuéjols, Kapoor and Vušković gave an algorithm to test if a 0/±1 matrix is balanced. Truemper has characterized balanceable 0/1 matrices in terms of forbidden submatrices. In this paper we give an algorithm that explicitly finds one of these forbidden submatrices or shows that none exists. Received: October 2004     

Generalized fuzzy matrices

We give a systematic development of fuzzy matrix theory. Many of our results generalize to matrices over the two element Boolean algebra, over the nonnegative real numbers, over the nonnegative integers, and over the semirings, and we present these generalizations. Our first main result is that while spaces of fuzzy vectors do not have a unique basis in general they have a unique standard basis, and the cardinality of any two bases are equal. Thus concepts of row and column basis, row and column rank can be defined for fuzzy matrices. Then we study Green's equivalence classes of fuzzy matrices. New we give criteria for a fuzzy matrix to be regular and prove that the row and column rank of any regular fuzzy matrix are equal. Various inverses are also studied. In the next section, we obtain bounds for the index and period of a fuzzy matrix.     

On the row merge tree for sparse LU factorization with partial pivoting

We consider the problem of structure prediction for sparse LU factorization with partial pivoting. In this context, it is well known that the column elimination tree plays an important role for matrices satisfying an irreducibility condition, called the strong Hall property. Our primary goal in this paper is to address the structure prediction problem for matrices satisfying a weaker assumption, which is the Hall property. For this we consider the row merge matrix, an upper bound that contains the nonzeros in L and U for all possible row permutations that can be performed during the numerical factorization with partial pivoting. We discuss the row merge tree, a structure that represents information obtained from the row merge matrix; that is, information on the dependencies among the columns in Gaussian elimination with partial pivoting and on structural upper bounds of the factors L and U. We present new theoretical results that show that the nonzero structure of the row merge matrix can be described in terms of branches and subtrees of the row merge tree. These results lead to an efficient algorithm for the computation of the row merge tree, that uses as input the structure of A, and has a time complexity almost linear in the number of nonzeros in A. We also investigate experimentally the usage of the row merge tree for structure prediction purposes on a set of matrices that satisfy only the Hall property. We analyze in particular the size of upper bounds of the structure of L and U, the reordering of the matrix based on a postorder traversal and its impact on the factorization runtime. We show experimentally that for some matrices, the row merge tree is a preferred alternative to the column elimination tree. AMS subject classification (2000)  65F50, 65F05, 68R10     

To Solving Multiparameter Problems of Algebra. 5. The ∇V-q Factorization Algorithm and its Applications

The algorithm of ∇V-factorization, suggested earlier for decomposing one- and two-parameter polynomial matrices of full row rank into a product of two matrices (a regular one, whose spectrum coincides with the finite regular spectrum of the original matrix, and a matrix of full row rank, whose singular spectrum coincides with the singular spectrum of the original matrix, whereas the regular spectrum is empty), is extended to the case of q-parameter (q ≥ 1) polynomial matrices. The algorithm of ∇V-q factorization is described, and its justification and properties for matrices with arbitrary number of parameters are presented. Applications of the algorithm to computing irreducible factorizations of q-parameter matrices, to determining a free basis of the null-space of polynomial solutions of a matrix, and to finding matrix divisors corresponding to divisors of its characteristic polynomial are considered. Bibliogrhaphy: 4 titles. __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 309, 2004, pp. 144–153.     

A globally convergent non-interior point algorithm with full Newton step for second-order cone programming

A non-interior point algorithm based on projection for second-order cone programming problems is proposed and analyzed. The main idea of the algorithm is that we cast the complementary equation in the primal-dual optimality conditions as a projection equation. By using this reformulation, we only need to solve a system of linear equations with the same coefficient matrix and compute two simple projections at each iteration, without performing any line search. This algorithm can start from an arbitrary point, and does not require the row vectors of A to be linearly independent. We prove that our algorithm is globally convergent under weak conditions. Preliminary numerical results demonstrate the effectiveness of our algorithm.     

Conley's spectral sequence via the sweeping algorithm

In this article we consider a spectral sequence (Er,dr) associated to a filtered Morse-Conley chain complex (C,Δ), where Δ is a connection matrix. The underlying motivation is to understand connection matrices under continuation. We show how the spectral sequence is completely determined by a family of connection matrices. This family is obtained by a sweeping algorithm for Δ over fields F as well as over Z. This algorithm constructs a sequence of similar matrices Δ⁰=Δ,Δ¹,… , where each matrix is related to the others via a change-of-basis matrix. Each matrix Δr over F (resp., over Z) determines the vector space (resp., Z-module) Er and the differential dr. We also prove the integrality of the final matrix ΔR produced by the sweeping algorithm over Z which is quite surprising, mainly because the intermediate matrices in the process may not have this property. Several other properties of the change-of-basis matrices as well as the intermediate matrices Δr are obtained.     

四元数体上Hermite矩阵的最小化问题

该文建立了四元数矩阵对的标准相关分解(CCD-Q). 借助CCD-Q, GSVD-Q 和有限维内积空间中的投影定理, 该文得到了基于四元数矩阵方程$AXB=C$的Hermite矩阵最小化问题解的表达式.     

Finite iterative solutions to coupled Sylvester-conjugate matrix equations

This paper is concerned with solutions to the so-called coupled Sylveter-conjugate matrix equations, which include the generalized Sylvester matrix equation and coupled Lyapunov matrix equation as special cases. An iterative algorithm is constructed to solve this kind of matrix equations. By using the proposed algorithm, the existence of a solution to a coupled Sylvester-conjugate matrix equation can be determined automatically. When the considered matrix equation is consistent, it is proven by using a real inner product in complex matrix spaces as a tool that a solution can be obtained within finite iteration steps for any initial values in the absence of round-off errors. Another feature of the proposed algorithm is that it can be implemented by using original coefficient matrices, and does not require to transform the coefficient matrices into any canonical forms. The algorithm is also generalized to solve a more general case. Two numerical examples are given to illustrate the effectiveness of the proposed methods.     

An iterative algorithm for solving a class of generalized coupled Sylvester-transpose matrix equations over bisymmetric or skew-anti-symmetric matrices

This paper presents an iterative algorithm to solve a class of generalized coupled Sylvester-transpose matrix equations over bisymmetric or skew-anti-symmetric matrices. When the matrix equations are consistent, the bisymmetric or skew-anti-symmetric solutions can be obtained within finite iteration steps in the absence of round-off errors for any initial bisymmetric or skew-anti-symmetric matrices by the proposed iterative algorithm. In addition, we can obtain the least norm solution by choosing the special initial matrices. Finally, numerical examples are given to demonstrate the iterative algorithm is quite efficient. The merit of our method is that it is easy to implement.     

Deterministic computation of the characteristic polynomial in the time of matrix multiplication

This paper describes an algorithm which computes the characteristic polynomial of a matrix over a field within the same asymptotic complexity, up to constant factors, as the multiplication of two square matrices. Previously, this was only achieved by resorting to genericity assumptions or randomization techniques, while the best known complexity bound with a general deterministic algorithm was obtained by Keller-Gehrig in 1985 and involves logarithmic factors. Our algorithm computes more generally the determinant of a univariate polynomial matrix in reduced form, and relies on new subroutines for transforming shifted reduced matrices into shifted weak Popov matrices, and shifted weak Popov matrices into shifted Popov matrices.     

A sign test for detecting the equivalence of Sylvester Hadamard matrices

In this paper we show that every matrix in the class of Sylvester Hadamard matrices of order 2 ^k under H-equivalence can have full row and column sign spectrum, meaning that tabulating the numbers of sign interchanges along any row (or column) gives all integers 0,1,...,2 ^k  − 1 in some order. The construction and properties of Yates Hadamard matrices are presented and is established their equivalence with the Sylvester Hadamard matrices of the same order. Finally, is proved that every normalized Hadamard matrix has full column or row sign spectrum if and only if is H-equivalent to a Sylvester Hadamard matrix. This provides us with an efficient criterion identifying the equivalence of Sylvester Hadamard matrices.     

Mixed dominating matrices

We characterize the class of matrices for which the set of supports of nonnegative vectors in the null space can be determined by the signs of the entries of the matrix. This characterization is in terms of mixed dominating matrices, which are defined by the nonexistence of square submatrices that have nonzeros of opposite sign in each row. The class of mixed dominating matrices is contained in the class of L-matrices from the theory of sign-solvability, and generalizes the class of S-matrices. We give a polynomial-time algorithm to decide if a matrix is mixed dominating. We derive combinatorial conditions on the face lattice of a Gale transform of a matrix in this class.