an Off-Diagonal Block of

Similar documents
Spectral inequalities and equalities involving products of matrices

PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS

arxiv: v3 [math.ra] 22 Aug 2014

b jσ(j), Keywords: Decomposable numerical range, principal character AMS Subject Classification: 15A60

Interlacing Inequalities for Totally Nonnegative Matrices

Linear Operators Preserving the Numerical Range (Radius) on Triangular Matrices

Convexity of the Joint Numerical Range

Maximizing the numerical radii of matrices by permuting their entries

CHI-KWONG LI AND YIU-TUNG POON. Dedicated to Professor John Conway on the occasion of his retirement.

Some inequalities for sum and product of positive semide nite matrices

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

arxiv: v4 [quant-ph] 8 Mar 2013

A Note on Eigenvalues of Perturbed Hermitian Matrices

RANKS OF QUANTUM STATES WITH PRESCRIBED REDUCED STATES

Generalized Interlacing Inequalities

Matrix Inequalities by Means of Block Matrices 1

CANONICAL FORMS, HIGHER RANK NUMERICAL RANGES, TOTALLY ISOTROPIC SUBSPACES, AND MATRIX EQUATIONS

= W(XI - A) = W(Ci(Ai - A)),

Schur complements and matrix inequalities in the Löwner ordering

Extremal Characterizations of the Schur Complement and Resulting Inequalities

Inequalities For Singular Values And Traces Of Quaternion Hermitian Matrices

AMS526: Numerical Analysis I (Numerical Linear Algebra)

arxiv: v2 [quant-ph] 5 Dec 2013

Compound matrices and some classical inequalities

Math 408 Advanced Linear Algebra

Isometries Between Matrix Algebras

Numerical Ranges and Dilations. Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement

Operators with numerical range in a closed halfplane

Some Convexity Features Associated with Unitary Orbits

DAVIS WIELANDT SHELLS OF NORMAL OPERATORS

LINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES

Decomposable Numerical Ranges on Orthonormal Tensors

arxiv: v3 [math.ra] 10 Jun 2016

Singular Value Inequalities for Real and Imaginary Parts of Matrices

arxiv: v1 [math.ra] 8 Apr 2016

Uniqueness of the Solutions of Some Completion Problems

Linear Algebra (Review) Volker Tresp 2017

On Euclidean distance matrices

A Remark on Polynomial Norms and

ACI-matrices all of whose completions have the same rank

Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?

A Note on Product Range of 3-by-3 Normal Matrices

Induced Norms, States, and Numerical Ranges

Matrix operations Linear Algebra with Computer Science Application

3.2 Gaussian Elimination (and triangular matrices)

InequalitiesInvolvingHadamardProductsof HermitianMatrices y

GENERALIZED NUMERICAL RANGES AND QUANTUM ERROR CORRECTION

Fuzhen Zhang s Publication List

Minimum number of non-zero-entries in a 7 7 stable matrix

Eigenvalue and Eigenvector Homework

2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian

ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )

Preservers of eigenvalue inclusion sets

Absolute value equations

Fuzhen Zhang s Publication and Presentation List (updated Dec. 2007)

A new criterion and a special class of k-positive maps

arxiv: v1 [math.ca] 7 Jan 2015

Schur Multiplicative Maps on Matrices

A Maximal Inequality of Non-negative

Algebra Homework, Edition 2 9 September 2010

The minimum rank of matrices and the equivalence class graph

The eigenvalues of an n x n complex matrix A are the roots of the n th degree polynomial det(a-λi) = 0. Further the positive square roots of the

Banach Journal of Mathematical Analysis ISSN: (electronic)

CSL361 Problem set 4: Basic linear algebra

Distribution for the Standard Eigenvalues of Quaternion Matrices

Necessary And Sufficient Conditions For Existence of the LU Factorization of an Arbitrary Matrix.

Clarkson Inequalities With Several Operators

Kantorovich type inequalities for ordered linear spaces

On the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination

Row and Column Distributions of Letter Matrices

Some inequalities for unitarily invariant norms of matrices

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES

Linear Algebra (Review) Volker Tresp 2018

On the second Laplacian eigenvalues of trees of odd order

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

for all A M m and B M n Key words. Complex matrix, linear preserver, spectral norm, Ky Fan k-norm, Schatten p-norm, tensor product.

Matrices. Chapter What is a Matrix? We review the basic matrix operations. An array of numbers a a 1n A = a m1...

SOME INEQUALITIES FOR THE KHATRI-RAO PRODUCT OF MATRICES

Determinantal and Eigenvalue Inequalities for Matrices with Numerical Ranges in a Sector

1 Linear Algebra Problems

On the adjacency matrix of a block graph

Algebra Exam Topics. Updated August 2017

Numerical Linear Algebra Homework Assignment - Week 2

A survey on Linear Preservers of Numerical Ranges and Radii. Chi-Kwong Li

A note on 5 5 Completely positive matrices

Generalized Schur complements of matrices and compound matrices

CHARACTERIZATIONS. is pd/psd. Possible for all pd/psd matrices! Generating a pd/psd matrix: Choose any B Mn, then

New insights into best linear unbiased estimation and the optimality of least-squares

Some Reviews on Ranks of Upper Triangular Block Matrices over a Skew Field

On generalized Schur complement of nonstrictly diagonally dominant matrices and general H- matrices

1 Last time: least-squares problems

Lecture 12 (Tue, Mar 5) Gaussian elimination and LU factorization (II)

Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices

ELA

Positive definite preserving linear transformations on symmetric matrix spaces

Lecture 7: Positive Semidefinite Matrices

CLASSIFICATION OF TREES EACH OF WHOSE ASSOCIATED ACYCLIC MATRICES WITH DISTINCT DIAGONAL ENTRIES HAS DISTINCT EIGENVALUES

Transcription:

J. oflnequal. & Appl., 1999, Vol. 3, pp. 137-142 Reprints available directly from the publisher Photocopying permitted by license only (C) 1999 OPA (Overseas Publishers Association) N.V. Published by license under the Gordon and Breach Science Publishers imprint. Printed in Malaysia. Inequalities on the Singular Values of an Off-Diagonal Block of a Hermitian Matrix* CHI-KWONG LI and ROY MATHIAS Department of Mathematics, College of William & Mary, Williamsburg, VA 23187, USA (Received 26 November 1997," Revised 4 February 1998) A majorization relating the singular values of an off-diagonal block of a Hermitian matrix and its eigenvalues is obtained. This basic majorization inequality implies various new and existing results. Keywords." Hermitian matrix; Eigenvalue; Singular value; Majorization AMS Subject Classification: 15A42; 15A60 1 INTRODUCTION Let A(H) > > A,(H) denote the eigenvalues of an n x n Hermitian matrix H. For an m x n complex matrix X, let ri(x)- V/Ai(X*X) denote the ith singular value for - 1,..., k, where k min{m, n}, and let or(x) (Crl(A),..., r(a)) be the vector of singular values of X. In [2], the following result was obtained as a generalization of a result in [6]. THEOREM Suppose H is an n x n positive definite matrix. Thenfor any n x k matrix X such that X*X I, tr(x*hx (X*HX) 2) <.-(A(H)- A,_j.+ (H)) 2. This research was supported by grants from the NSF. Corresponding author. E-mail: ckli@math.wm.edu. 137

138 C.-K. LI AND R. MATHIAS In fact, this result was conjectured by J. Durbin, and proved in [1] and [3], independently. This result is important in the context of studying the relative performance of the least squares estimator and the best linear unbiased estimator in a linear model [1]. Observe that if U is a unitary matrix of the form IX Y], then X*HY is a k n matrix, and tr(x*hzx (X*HX) 2) rn rj(x*hy) 2, where m min{k, n-k}. In the following, we obtain a majorization result that will allow one to deduce a whole family of inequalities including Theorem 1. In Refs. [1-3], the proof of Theorem was done using partial differentiation to locate the optimal matrix that yields the upper bound of tr(x*h2x (X*HX)2). In our case, we use different approaches to give two proofs for our result Theorem 2 that connect our problem to other subjects. We need some more definitions to state our result. Given two real (row or column) vectors x, y E Rn, we say that x is weakly majorized by y, denoted by x <w y, if the sum of the k largest entries of x is not larger than that of y for each k 1, n. If in addition the sum of the entries of each of the vectors is the same then we say that x is majorized by y. THEOREM 2 Let H be an n n Hermitian matrix. Thenfor any unitary matrix U of the form [X Y], where X is n k matrix, we have cr(x*hy) -<w 1/2(A1 (H) An(H),..., Am(H) An_m+l (n)), where m- min{k, n-k}. Consequently, for any Schur convex increasing functionf: Rm---+I, we have f(a(x*hy)) _< f(1/2(a, (H)- An(H),...,Am(H)- An-m+l(O)) ). Note that if we take f(x)- _,jm=, x2i in Theorem 2, we obtain Theorem 1. In fact, there are many other interesting Schur convex functions (see [5, Chapter 1] for details). For instance,f(x) )P--1 IxJl p with p > and the kth elementary symmetric function Ek(Xl,...,Xm) with _< k <_ m are such examples.

2 PROOFS MAJORIZATION OF AN OFF-DIAGONAL BLOCK 139 We first give a proof of Theorem 2 using the theory of majorization (see [5] for the general background) and a reduction of the problem to the 2 2 case. First proof of Theorem 2 Assume, without loss of generality, that U I and that k < (n-k). Write Hi1 B H B* H2:2 ) where Hl is k k and 922 is (n-k) (n-k). Let B- WEV be a singular value decomposition of B, where V and W are unitary. Then the eigenvalues of the matrix H- 0 H 0 V* are the same as those of H. The 2 in rows and columns and k + is 2 principal submatrix of H lying [i, k + i] ( -lii _ffi(b) ), i-1,...,k. cri(b) hk+i,k+i One easily checks that Theorem 2 is true when n- 2 and k- 1. As a result, if where #i _/i and RRi- 12, then [i, k + i] RT ( #io T]iO ) Ri ffi(b) (i- 1]i)/2. Let R be the n n unitary matrix obtained from In by replacing In[i,k+i] by R for all 1,...,k. Then (R nr)i i and

140 C.-K. LI AND R. MATHIAS (R*R)k+i,k+i- ]i. Since the vector of diagonal entries of R*gR is majorized by the vector of eigenvalues of R*R (e.g., see [5, Chapter 9, B. 1]), for any 1,..., k, we have and Z #i-- Z(g*g)ii <_ Z Ai(H) i=1 i=1 i=1 It follows that Z T]i Z (R nr)k+i,k+i Z )n-i+l(n)" i=1 i=1 i=1 Z o i(b)<_ Z(#i- gli)/2 <_ i=1 i=1 i=1 Z(,i(H)-,n-i+l(H))/2. Next, we give a proof of Theorem 2 using the theory of the C-numerical range (e.g., see [4] and its references for the general background). Second proof of Theorem 2 By the singular value decomposition, one can find unitary matrices V and W of appropriate sizes such that (VX*HYW)jj- crj(x*hy), j- 1,...,m. Thus for any positive integer with _< _< m, if we let Ct -]=1 E+j,j, where {El 1,El2,..., Enn} denotes the standard basis for n then Z crj(x*hy) <_ max {I 7(RX*HYS)jI" R, S unitary} _< max {I Z(Z*HZ)+j,jI" Z unitary} max(ltr(z*hzc,)l" Z unitary} max{itr(hzc,z*)l" Z unitary}, n matrices,

MAJORIZATION OF AN OFF-DIAGONAL BLOCK 141 which can be viewed as the H-numerical radius rh(ct) of Ct (e.g., see [4] for the general background). Moreover, since Ct=( 0, 00) is in the so-called shift block form and (Ct + C)/2 has eigenvalues "i/2,...,1/2,0,...,0,-1/2,...,-1/2, we have (e.g., see [4, (5.1) and (5.2)]) rn(ct) max{[tr(hz(ct + C [)Z*)/2 I" Z unitary} Aj((Ct + Z(/j(H) /n-j+l (H))/2. Remarks Suppose that A1 _> _> An are given real numbers and that k is a positive integer such that < k < n. Let m min{k, n-k}. One can construct 2 2 matrices Hi with eigenvalues,i,/n_m+i, and offdiagonal entries equal to (i-)n_m+ i)/2. Applying a suitable permutation similarity to the matrix will yield a matrix H H1 0"" Hm diag(am+l,..., An-m) where B is k (n-k) such that if < i=j < k, BO 0 otherwise. Clearly,/) has eigenvalues A1,..., An. Thus, we see that our result in Theorem 2 is best possible.

142 C.-K. LI AND R. MATHIAS In the context of statistics one is interested in real symmetric matrices. Since Theorem 2 is true for Hermitian matrices it is afortiori true for real symmetric matrices. It cannot be improved in the case of real symmetric matrices either because the matrix constructed in the example above is a real symmetric matrix. Acknowledgement We thank Professor G. Styan for drawing our attention to [1,2], Professor Z. Jia for sending us a copy of the paper. and References [1] P. Bloomfield and G.S. Watson, The inefficiency of least squares, Biometrika 62 (1975), 121-128. [2] Z. Jia, An extension of Styan s inequality (Chinese), Gongcheng Shuxue Xuebao (J. Engineering Mathematics) 13 (1996), 122-126. [3] M. Knott, On the minimum efficiency of least squares, Biometrika 62 (1975), 129-132. [4] C.K. Li, C-numerical ranges and C-numerical radii, Linear and Multilinear Algebra 37 (1994), 51-82. [5] A.W. Marshall and I. Olkin, Inequalities: The Theory of Majorization and its Applications, Academic Press, 1979. [6] G.P.H. Styan, On some inequalities associated with ordinary least squares and the Kantorovich inequality, Acta Univ. Tampere, Ser. A. 153 (1983), 158-166. Note added in proof X. Zhan has another proof of our main theorem, which is also obtained independently by R. Bhatia, F.C. Silva, P. Assouad and J.A. Dias da Silva.