Sign Patterns of G-Matrices

Similar documents
Intrinsic products and factorizations of matrices

Sign Patterns that Allow Diagonalizability

MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2

A Simple Proof of Fiedler's Conjecture Concerning Orthogonal Matrices

Matrix Algebra. Matrix Algebra. Chapter 8 - S&B

Foundations of Matrix Analysis

Chapter 3 Transformations

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 13

Algebraic Methods in Combinatorics

MTH 464: Computational Linear Algebra

Review of Linear Algebra

SHORT VITA Frank J. Hall June, B.A. St. Mary's University, San Antonio, Texas, May, 1965

This can be accomplished by left matrix multiplication as follows: I

Chapter 1. Matrix Algebra

3. Vector spaces 3.1 Linear dependence and independence 3.2 Basis and dimension. 5. Extreme points and basic feasible solutions

ELE/MCE 503 Linear Algebra Facts Fall 2018

MAT Linear Algebra Collection of sample exams

Some Results on Generalized Complementary Basic Matrices and Dense Alternating Sign Matrices

Algebraic Methods in Combinatorics

There are two things that are particularly nice about the first basis

Row and Column Distributions of Letter Matrices

Symmetric matrices and dot products

Index. book 2009/5/27 page 121. (Page numbers set in bold type indicate the definition of an entry.)

Central Groupoids, Central Digraphs, and Zero-One Matrices A Satisfying A 2 = J

MAT 2037 LINEAR ALGEBRA I web:

Minimum Ranks and Refined Inertias of Sign Pattern Matrices

Matrix Operations: Determinant

A 2 G 2 A 1 A 1. (3) A double edge pointing from α i to α j if α i, α j are not perpendicular and α i 2 = 2 α j 2

Recall the convention that, for us, all vectors are column vectors.

Math Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88

Lecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Problem Set (T) If A is an m n matrix, B is an n p matrix and D is a p s matrix, then show

Checkered Hadamard Matrices of Order 16

Preliminary Linear Algebra 1. Copyright c 2012 Dan Nettleton (Iowa State University) Statistics / 100

Linear Algebra V = T = ( 4 3 ).

3 (Maths) Linear Algebra

Notes on Linear Algebra and Matrix Theory


1 Linear Algebra Problems

Lecture II: Linear Algebra Revisited

Some bounds for the spectral radius of the Hadamard product of matrices

(v, w) = arccos( < v, w >

Properties of Matrices and Operations on Matrices

Math 240 Calculus III

We first repeat some well known facts about condition numbers for normwise and componentwise perturbations. Consider the matrix

Miscellaneous Results, Solving Equations, and Generalized Inverses. opyright c 2012 Dan Nettleton (Iowa State University) Statistics / 51

orthogonal relations between vectors and subspaces Then we study some applications in vector spaces and linear systems, including Orthonormal Basis,

Row Space and Column Space of a Matrix

Math Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT

Review of Matrices and Block Structures

Lecture 18: The Rank of a Matrix and Consistency of Linear Systems

Lecture Notes in Linear Algebra

Potentially nilpotent tridiagonal sign patterns of order 4

UNIT 3 MATRICES - II

Matrices and systems of linear equations

INVERSE CLOSED RAY-NONSINGULAR CONES OF MATRICES

ORIE 6334 Spectral Graph Theory September 8, Lecture 6. In order to do the first proof, we need to use the following fact.

Notes on Mathematics

Lecture 7: Introduction to linear systems

(v, w) = arccos( < v, w >

Lecture 3: QR-Factorization

INVERSE TRIDIAGONAL Z MATRICES

The Singular Value Decomposition

Introduction to Determinants

Graduate Mathematical Economics Lecture 1

Undergraduate Mathematical Economics Lecture 1

Linear equations in linear algebra

Chap 3. Linear Algebra

Linear algebra and applications to graphs Part 1

ACI-matrices all of whose completions have the same rank

a 11 a 12 a 11 a 12 a 13 a 21 a 22 a 23 . a 31 a 32 a 33 a 12 a 21 a 23 a 31 a = = = = 12

Spectrally arbitrary star sign patterns

The Adjacency Matrix, Standard Laplacian, and Normalized Laplacian, and Some Eigenvalue Interlacing Results

Recitation 8: Graphs and Adjacency Matrices

MAT 1332: CALCULUS FOR LIFE SCIENCES. Contents. 1. Review: Linear Algebra II Vectors and matrices Definition. 1.2.

c c c c c c c c c c a 3x3 matrix C= has a determinant determined by

Matrix Completion Problems for Pairs of Related Classes of Matrices

MATH 315 Linear Algebra Homework #1 Assigned: August 20, 2018

Conceptual Questions for Review

Math 407: Linear Optimization

Sparse spectrally arbitrary patterns

Some New Results on Lyapunov-Type Diagonal Stability

MATRICES. knowledge on matrices Knowledge on matrix operations. Matrix as a tool of solving linear equations with two or three unknowns.

Eigenvalues, Eigenvectors. Eigenvalues and eigenvector will be fundamentally related to the nature of the solutions of state space systems.

QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY

MATH2210 Notebook 2 Spring 2018

CS100: DISCRETE STRUCTURES. Lecture 3 Matrices Ch 3 Pages:

Z-Pencils. November 20, Abstract

Lecture 4 Orthonormal vectors and QR factorization

Chapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors

Linear Algebra: Matrix Eigenvalue Problems

Linear Algebra and Matrix Theory

A matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and

1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )

Advanced Engineering Mathematics Prof. Pratima Panigrahi Department of Mathematics Indian Institute of Technology, Kharagpur

Fall 2016 MATH*1160 Final Exam

Abed Elhashash and Daniel B. Szyld. Report August 2007

EE731 Lecture Notes: Matrix Computations for Signal Processing

Numerical Linear Algebra Homework Assignment - Week 2

DS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.

Transcription:

Department of Mathematics and Statistics Georgia State University Atlanta, GA 30303, USA

G-Matrices Sign Patterns Minimum Rank Open Questions

Outline G-Matrices Sign Patterns Minimum Rank Open Questions

M. Fiedler/F. Hall, G-Matrices, LAA, 2012 All matrices in the paper are real. We call a matrix A a G-matrix if A is nonsingular and there exist nonsingular diagonal matrices D 1 and D 2 such that A T = D 1 AD 2, where A T denotes the transpose of the inverse of A. The matrices D 1 and D 2 are not uniquely determined - we can multiply one by a nonzero number and divide the other by the same number.

A rich class of matrices! Easy to see: All orthogonal matrices are G-matrices. All nonsingular diagonal matrices are G-matrices. If A is a G-matrix, then both A T and A 1 are G-matrices. If A is an n n G-matrix and D is an n n nonsingular diagonal matrix, then both AD and DA are G-matrices. If A is a G-matrix, then both matrices A and A T have the same zero, non-zero structure. Thus the zero, non-zero structures of A and A 1 are symmetric to each other. The direct sum of G-matrices is again a G-matrix.

The product form of the definition leads to the following two facts: Compound matrices of a G-matrix are also G-matrices. The Kronecker product of G-matrices is a G-matrix. Also: A 2 2 matrix is G-matrix if and only if it is nonsingular and has four or two nonzero entries.

Cauchy matrices have form C = [c ij ], where c ij = 1 x i +y j some numbers x i and y j. for We shall restrict to square, say n n, Cauchy matrices - such matrices are defined only if x i + y j 0 for all pairs of indices i, j, and it is well known that C is nonsingular if and only if all the numbers x i are mutually distinct and all the numbers y j are mutually distinct. From M. Fiedler, Notes on Hilbert and Cauchy matrices, LAA, 2010: Every nonsingular Cauchy matrix is a G-matrix.

For generalized Cauchy matrices of order n, additional parameters u 1,..., u n, v 1,..., v n are considered: [ ui v j Ĉ =. x i +y j ] Note that Ĉ = D 1 CD 2 where D 1 = diag(u i ), D 2 = diag(v j ) so that Ĉ is a G-matrix.

Just a couple of results re the diagonal matrices: Theorem. Suppose A is a G-matrix and A T = D 1 AD 2, where D 1 and D 2 are nonsingular diagonal matrices. Then the inertia of D 1 is equal to the inertia of D 2. Proof. We have A T D 1 AD 2 = I and so A T D 1 A = D 1 2. Since A is nonsingular, the result follows from Sylvester s Law of Inertia. Remark. As the example of a diagonal matrix shows, the inertia of both matrices D 1 and D 2 need not be uniquely determined by the given G-matrix.

Theorem. Let A be a symmetric irreducible G-matrix. Then the diagonal matrices D 1, D 2 in A T = D 1 AD 2 can be chosen to be either the same or the negative of each other. If A is positive definite then D 1 and D 2 can be chosen to be the same. If we denote the common matrix as D, then the matrix AD is involutory. Denote the nonnegative square root of the modulus D of D as D 1 2. Then the matrix B = D 1 2 AD 1 2 has the property that all its integral powers are G-matrices. Finally, if S = diag(ε 1, ε 2,..., ε n ) is the matrix of the signs of D, ie. S = D D 1, then the diagonal entries a ii of A and the diagonal entries α ii of A 1 satisfy ε i ( a ii α ii 1) = 0. i

Two papers on Generalized G-Matrices M. Matsuura, A note on generalized G-matrices, LAA, 2012 M. Fiedler/T. Markham, More on G-matrices, LAA, 2013

Outline G-Matrices Sign Patterns Minimum Rank Open Questions

G n denotes the class of all n n sign pattern matrices A that allow a G-matrix, ie, there exists a nonsingular matrix B Q(A) such that B T = D 1 BD 2 for some nonsingular diagonal matrices D 1 and D 2. Theorem The class G n is closed under (i) multiplication (on either side) by a permutation pattern, and (ii) multiplication (on either side) by a diagonal signature pattern. The use of these operations in G n then produces equivalent sign patterns.

Let C n (GC n ) be the class of all sign patterns of the n n nonsingular Cauchy (generalized Cauchy) matrices. It should be clear that C n (GC n ) is closed under operation (i) (operations (i) and (ii)) above. The classes C n and GC n are two particular sub-classes of G n. Since x + y = (x + c) + (y c) and when we permute the rows/columns, we still have a Cauchy matrix C, for a given such nonsingular matrix, we can assume that x 1 > x 2 > > x n > 0, and y 1 > y 2 > > y n. Further, the sign pattern of C is the same as the sign pattern of the matrix [x i + y j ]. So, the following result is clear.

Theorem The class C n is the same as the class of n n sign patterns permutation equivalent to a staircase sign pattern, where the part above (below) the staircase is all + ( ). In this form, whenever there is a minus, then to the right and below there are also minuses. Note that this form includes [ the] all + and all patterns. + However, for example, is excluded as a pattern in + C 2. In fact:

Theorem An n n sign pattern matrix A belongs to C n if and only if it does not contain either one of the patterns [ ] [ ] + + or + + as submatrices. Recall the form of generalized Cauchy matrix: Ĉ = D 1 CD 2. Theorem The class GC n is the same as the class of n n sign patterns that are signature equivalent to a sign pattern in C n. Equivalently, A GC n if and only if there exists a minimum rank one sign pattern matrix M such that A M C n, where denotes the Hadamard product.

Every 2 2 (+, ) sign pattern is a matrix in GC 2. Every 3 3 (+, ) sign pattern is a matrix in GC 3. But, in general GC n does not contain all the (+, ) n n sign patterns, that is, there are forbidden sign patterns. No GC-matrix can have the following pattern (or its equivalent or transpose) as submatrix: + + + S 3 = + + +. + + +

By the way, S 3 = + + + + + + + + + is an L-matrix. Every 3 4 real matrix with the same sign pattern as S 3 has linearly independent rows. R.A. Brualdi, K.L. Chavey, B.L. Shader, Rectangular L-matrices, LAA, 1994.

PO n is the class of n n sign patterns which allow an orthogonal matrix. The class PO n is then a subclass of the class G n. In EHHL, CMJ, 1999, a number of nontrivial sign patterns in PO n (and hence in G n ) are constructed. EG The patterns + + +, + + arising from Householder transformations, are in PO 3. These two patterns are then in the intersection of PO 3 and GC 3.

The sign patterns + + + +, + + + are examples of PO 4 patterns which are not in GC 4, by a forbidden submatrix. The same holds for the corresponding Householder type patterns in PO n, n 4.

Remark Not every sign pattern in G n is in PO n or in GC n. For example, the Kronecker product + + + A = + + + + is in G 9 since the Kronecker product of G-matrices is a G-matrix.

Now, A = + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +. By observing columns 4 and 7 of A, it is clear that A / PO 9. Also, by considering the submatrix at the intersection of rows 1, 3, 5 and columns 2,4,5,6, we know that A is equivalent to a pattern with a forbidden submatrix. So, A / GC 9.

Outline G-Matrices Sign Patterns Minimum Rank Open Questions

Back to C n, the class of all sign patterns of the n n nonsingular Cauchy matrices. Recall: The class C n is the same as the class of n n sign patterns permutation equivalent to a staircase sign pattern, where the part above (below) the staircase is all + ( ). Let C be an n n nonsingular Cauchy matrix. Now, the sign pattern of C is the same as the sign pattern of the reciprocal matrix [x i + y j ]. The reciprocal matrix has rank at most 2, as any 3rd order subdet is 0. But, clearly, excluding the all + and all patterns, any real matrix in the sign pattern class has rank 2. Thus, excluding equivalences of the all + and all patterns, if A C n, then the minimum rank of A is 2.

Recall: The class GC n is the same as the class of n n sign patterns that are signature equivalent to a sign pattern in C n. Thus, excluding equivalences of the all + and all patterns, if A GC n, then the minimum rank of A is 2. Also, if A GC n, then there are signature sign patterns D 1, D 2 and permutation sign patterns P 1, P 2 such that P 1 D 1 AD 2 P 2 is a staircase sign pattern, where the part above (below) the staircase is all + ( ), that is, each row and each column of P 1 D 1 AD 2 P 2 is nonincreasing. This then also holds for the condensed sign pattern A c, where in general condensed means no zero line, no identical rows (or columns), and no two rows (or columns) which are negatives of each other.

This latter result is a special case of the following, from Z. Li et al, Sign patterns with minimum rank 2 and upper bounds on minimum ranks, LAMA, 2013: Theorem An m n sign pattern A has minimum rank 2 iff its condensed sign pattern A c satisfies the conditions (i) A c has at least two rows and two columns, (ii) each row and each column of A c has at most one zero entry, and (iii) there are signature sign patterns D 1, D 2 and permutation sign patterns P 1, P 2 such that each row and each column of P 1 D 1 A c D 2 P 2 is nonincreasing. Note that the conditions in the above theorem imply that A c is equivalent to a certain staircase pattern.

Next, suppose A is an n n (+, ) sign pattern with minimum rank 2. From condition (iii) of the above theorem, there are signature sign patterns D 1, D 2 and permutation sign patterns P 1, P 2 such that each row and each column of P 1 D 1 A c D 2 P 2 is nonincreasing. This then also holds for the sign pattern matrix A, so that A GC n. We thus see that, excluding equivalences of the all + and all patterns, GC n, the class of all sign patterns of the n n nonsingular generalized Cauchy matrices, is the same as the class of all n n (+, ) sign patterns with minimum rank 2!

Recall: The patterns + + +, + + are in the intersection of PO 3 and GC 3. So, the minimum rank of each pattern is 2 (it s also easy to see that each pattern is not SNS).

Recall: The sign pattern + A = + + + is an example of a pattern which is not in GC 4, by a forbidden submatrix. From N. Alon, J.H. Spencer, The Probabilistic Method, 2000: If each row of a (+, ) sign pattern M has at most k sign changes, then mr(m) k + 1. So, mr(a) 3. But, since A / GC 4, mr(a) 2, and clearly, mr(a) 1. Thus, mr(a) = 3, The same holds for the corresponding Householder type pattern in PO n, n 4.

The same can be said for the Householder pattern + +. + Next, what about the 3rd type, up to signature and permutation similarity, irreducible Householder pattern (of course in PO 3 ): 0 + + It s easy to see that this pattern is not SNS and its minimum rank is 2. What about the corresponding Householder type pattern A PO n, n 4??.

From the LAMA paper: The number of polynomial sign changes of a sign pattern M, denoted psc(m), is the largest number of polynomial sign changes of the rows of M (for a sign pattern with zero entries, this entails a more general definition). For every sign pattern M, mr(m) psc(m c ) + 1. Now, for our pattern A, psc(a c ) = 2, so that mr(a) 3. Observe that the 2nd, 3rd, and 4th rows form an L-matrix. So, mr(a) 3. Thus, mr(a) = 3.

A minimum rank 3 reference: Z. Li et al, Sign patterns with minimum rank 3 and point-line configurations

What is the minimum rank of the 9 9 pattern A on pages 20, 21?? A = + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +.

A = A 1 A 2 where mr(a 1 ) = mr(a 2 ) = 2. For real B 1, B 2, rank(b 1 B 2 ) = rank(b 1 )rank(b 2 ). So, mr(a) mr(a 1 )mr(a 2 ) = 4. Rows 3,4,5 of A form an L-matrix. So, mr(a) 3. Thus, mr(a) = 3 or 4. Note that A does not contain an L-matrix with 4 rows.

Outline G-Matrices Sign Patterns Minimum Rank Open Questions

If A G n then there exist a nonsingular matrix B Q(A) and nonsingular diagonal matrices D 1 and D 2 such that B T = D 1 BD 2. Hence, there are signature patterns S 1 and S 2 such that B 1 Q(S 1 A T S 2 ). (1) In particular, we may study the sign patterns A G n for which there exist diagonal matrices D 1, D 2 with positive diagonal entries in the equation B T = D 1 BD 2. Then, (1) becomes B 1 Q(A T ). (2)

In fact, in EHHL, CMJ, 1999, the class T n of all n n sign patterns A for which there exists a nonsingular matrix B Q(A) where (2) holds was studied. There it was asked if the class T n is the same as the subclass PO n. In the same spirit, we may pose the following questions: (i) Suppose A is an n n sign pattern matrix and there exists B Q(A) such that (2) holds. Do there then exist diagonal matrices D 1, D 2 with positive diagonal entries such that B T = D 1 BD 2? (ii) More generally, suppose that (1) holds for some A and B. Is B then necessarily a G-matrix, or, does there exist some ˆB Q(A) that is a G-matrix?

Let A be an n n sign pattern. Suppose there exist a nonsingular matrix B Q(A) and nonsingular diagonal matrices D 1, D 2 with + diagonal entries such that B T = D 1 BD 2. Let F = D 1, G = D 2. Then (FBG) 1 = (FBG) T. So, FBG is an orthogonal matrix in Q(A). Conversely, if C is an orthogonal matrix in Q(A), then C T = C = I n CI n. Thus, A allows a G-matrix with associated diagonal matrices with + diagonal entries iff A PO n. For such a sign pattern A, we have that A T n. Does the reverse hold???

Fiedler/Hall, Some inheritance properties for complementary basic matrices, LAA, 2010. Fiedler/Hall, G-matrices, LAA, 2012. Fiedler/Hall, A note on permanents and generalized complementary basic matrices, LAA, 2012. Fiedler/Hall, Some graph theoretic properties of generalized complementary basic matrices, LAA, 2013. Fiedler/Hall/Marsli, Gersgorin discs revisited, AMM, 2013. Fiedler/Hall, Combinatorial aspects of generalized complementary basic matrices, CEJM, 2013. Fiedler/Hall, Factorizable matrices, Special Matrices journal, 2013. Fiedler/Hall/Stroev, Dense alternating sign matrices and extensions, LAA, 2014. Fiedler/Hall, Max algebraic complementary matrices, LAA, 2014. Fiedler/Hall/Stroev, Permanents, determinants, and generalized complementary basic matrices, OaM, to appear.

Fiedler/Hall/Stroev, Discriminant matrices, in preparation. Fiedler/Hall, More on generalized complementary basic matrices, in preparation. Maybe hamburger matrices...