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

Size: px
Start display at page:

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

Transcription

1 Determinantal and Eigenvalue Inequalities for Matrices with Numerical Ranges in a Sector Raymond Nung-Sing Sze The Hong Kong Polytechnic University July 29, 2014 Based on a joint work with Chi-Kwong Li 2014 Workshop on Numerical Ranges and Numerical Radii Tsinghua Sanya International Mathematics Forum (TSIMF) Sanya, China Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 1 / 25

2 Introduction - Fischer Inequality Fischer inequality A11 A 12 Suppose A = M A 21 A n, where A 11 M m, is positive semi-definite. 22 Then det(a) det(a 11) det(a 22). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 2 / 25

3 Introduction - Fischer Inequality Fischer inequality A11 A 12 Suppose A = M A 21 A n, where A 11 M m, is positive semi-definite. 22 Then det(a) det(a 11) det(a 22). Outline of Proof. Suppose A 22 is invertible. Then A11 A 12 Im A A = = 12A 1 22 A11 A 12A 1 22 A21 0. A 21 A 22 0 I n m A 21 A 22 Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 2 / 25

4 Introduction - Fischer Inequality Fischer inequality A11 A 12 Suppose A = M A 21 A n, where A 11 M m, is positive semi-definite. 22 Then det(a) det(a 11) det(a 22). Outline of Proof. Suppose A 22 is invertible. Then A11 A 12 Im A A = = 12A 1 22 A11 A 12A 1 22 A21 0. A 21 A 22 0 I n m A 21 A 22 Because A 11 A 12A 1 22 A21 is positive semi-definite, det(a) = det(a 11 A 12A 1 22 A 21) det(a 22) det(a 11) det(a 22). By continuity, the general result holds. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 2 / 25

5 Introduction - Fischer Inequality Fischer inequality A11 A 12 Suppose A = M A 21 A n, where A 11 M m, is positive semi-definite. 22 Then det(a) det(a 11) det(a 22). Outline of Proof. Suppose A 22 is invertible. Then A11 A 12 Im A A = = 12A 1 22 A11 A 12A 1 22 A21 0. A 21 A 22 0 I n m A 21 A 22 Because A 11 A 12A 1 22 A21 is positive semi-definite, det(a) = det(a 11 A 12A 1 22 A 21) det(a 22) det(a 11) det(a 22). By continuity, the general result holds. Remark: This proof reveals that estimating the determinant and eigenvalues of the Schur complement of A w.r.t. A 22, i.e., A 11 A 12A 1 22 A 21. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 2 / 25

6 3-7 June, ILAS Conference In 2013 ILAS Conference, Providence, RI, USA, 3-7 June, 2013, Minghua Lin of University of Victoria gave a talk with title Fischer type determinantal inequalities for accretive-dissipative matrices In his[ talk, he considered ] the following problem. Given an n n matrix A11 A 12 A = M A 21 A n, where A 11 M m with m n/2. Find optimal 22 (smallest) γ > 0 such that det(a) γ det(a 11) det(a 22), which is connected to the study of growth factor in Gaussian elimination. By Fischer inequality, γ = 1 if A is positive semi-definite. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 3 / 25

7 3-7 June, ILAS Conference Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 4 / 25

8 Some early results before 2013 A matrix A is called accretive-dissipative if A + A 0 and i(a A) 0. Ikramov showed that if A is accretive-dissipative, then det(a) γ det(a 11) det(a 22) with γ = 3 m. [JMS 121: (2004)] Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 5 / 25

9 Some early results before 2013 A matrix A is called accretive-dissipative if A + A 0 and i(a A) 0. Ikramov showed that if A is accretive-dissipative, then det(a) γ det(a 11) det(a 22) with γ = 3 m. Lin improved the bound to γ = { 2 3m/2 if m n/3, 2 n/2 if n/3 < m n/2. [JMS 121: (2004)] [LAA 438: (2013)] (Submitted 15 October 2012) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 5 / 25

10 Some early results before 2013 A matrix A is called accretive-dissipative if A + A 0 and i(a A) 0. Ikramov showed that if A is accretive-dissipative, then det(a) γ det(a 11) det(a 22) with γ = 3 m. Lin improved the bound to γ = { 2 3m/2 if m n/3, 2 n/2 if n/3 < m n/2. [JMS 121: (2004)] Lin further proposed the following. Conjecture 1 Suppose A is accretive-dissipative. Then [LAA 438: (2013)] (Submitted 15 October 2012) det(a) 2 m det(a 11) det(a 22). Im I m Note that the above equality holds if A = I I m I n 2m. m Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 5 / 25

11 Numerical range The numerical range of a matrix A M n is defined by W (A) = {x Ax : x C n, x x = 1}. W (A) is always nonempty and convex. W (A) = {λ} if and only if A = λi. W (A) [0, ) if and only if A 0. W (A) R if and only if A is Hermitian. W (A) = conv σ(a) if A is normal. W (A) = {µ C : t [0,2π) eit µ + e it µ λ 1(e it A + e it A )}. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 6 / 25

12 Numerical range in a Sector For any α [0, π/2), let S α = {z C : Im z (Re z) tan α}. A subset of C is a sector of half angle α if it is of the form {e iϕ z : z S α} for some ϕ [0, 2π). 2α For any matrix A, one can always find α and θ in [0, 2π) such that W (A) e iθ S α W (e iθ A) S α. W (A) (WLOG, we may assume θ = 0.) If A is positive semi-definite, then W (A) S α with α = 0. 2α If A is accretive-dissipative, then W (e iπ/4 A) S α with α = π 4. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 7 / 25

13 Results and conjectures of Drury Drury proved that if W (A) is a subset of a sector of half angle α such that 0 α < π/(2m), then det(a) sec 2 (mα) det(a 11) det(a 22). [LAA 439: (2013)] (Submitted 6 June 2013) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 8 / 25

14 Results and conjectures of Drury Drury proved that if W (A) is a subset of a sector of half angle α such that 0 α < π/(2m), then det(a) sec 2 (mα) det(a 11) det(a 22). Drury further proposed the following. Conjecture 2 [LAA 439: (2013)] (Submitted 6 June 2013) If W (A) is a subset of a sector of half angle α such that 0 α < π/2, then det(a) sec 2m (α) det(a 11) det(a 22). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 8 / 25

15 Results and conjectures of Drury Drury proved that if W (A) is a subset of a sector of half angle α such that 0 α < π/(2m), then det(a) sec 2 (mα) det(a 11) det(a 22). Drury further proposed the following. Conjecture 2 [LAA 439: (2013)] (Submitted 6 June 2013) If W (A) is a subset of a sector of half angle α such that 0 α < π/2, then det(a) sec 2m (α) det(a 11) det(a 22). Moreover, if A 11 is 1 1 nonzero and A 22 is (n 1) (n 1) invertible, then lies in the set { R = det(a) det(a 11) det(a = det(a) 22) A 11 det(a 22) re i2φ C : 0 r } 2(cos(2φ) cos(2α)), α φ α. sin 2 (2α) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 8 / 25

16 Results and conjectures of Drury y R with α = π 12 x Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 9 / 25

17 Results and conjectures of Drury y y R with α = π 12 R with α = π 6 x x Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 9 / 25

18 Results and conjectures of Drury y y R with α = π 12 R with α = π 6 x x y R with α = π 4 x Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 9 / 25

19 Results and conjectures of Drury y y R with α = π 12 R with α = π 6 x x y y R with α = π 4 R with α = π 3 x x Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 9 / 25

20 8-14 Maryland, USA After the ILAS meeting, I visited Professor Li at Maryland for a week (8-14 June) with Zejun Huang. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 10 / 25

21 8-14 Maryland, USA After the ILAS meeting, I visited Professor Li at Maryland for a week (8-14 June) with Zejun Huang. We discussed these two conjectures. After trying several different approaches, finally we were able to obtain some results. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 10 / 25

22 8-14 Maryland, USA After the ILAS meeting, I visited Professor Li at Maryland for a week (8-14 June) with Zejun Huang. We discussed these two conjectures. After trying several different approaches, finally we were able to obtain some results. At the same time, we were also exploring how to make Zongzi for celebrating Dragon Boat Festival (12th June). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 10 / 25

23 Our approach - A generalized eigenvalue problem Let A = A11 A 12 M A 21 A n, where A 11 M m with m n/2. 22 Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 11 / 25

24 Our approach - A generalized eigenvalue problem Let A = A11 A 12 M A 21 A n, where A 11 M m with m n/2. 22 We consider the following equality for determinant det(a) = det(a 22) det(a 11 A 12A 1 22 A 21) = det(a 22) det(a 11) det(i m A 1 11 A 12A 1 22 A 21). and affirm (and improve) the conjectures of Drury and Lin. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 11 / 25

25 Our approach - A generalized eigenvalue problem Let A = A11 A 12 M A 21 A n, where A 11 M m with m n/2. 22 We consider the following equality for determinant det(a) = det(a 22) det(a 11 A 12A 1 22 A 21) = det(a 22) det(a 11) det(i m A 1 11 A 12A 1 22 A 21). and affirm (and improve) the conjectures of Drury and Lin. Determine the optimal containment region of the generalized eigenvalue λ satisfying A A12 λ x = x for some nonzero x C n. 0 A 22 A 21 0 Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 11 / 25

26 Our approach - A generalized eigenvalue problem Let A = A11 A 12 M A 21 A n, where A 11 M m with m n/2. 22 We consider the following equality for determinant det(a) = det(a 22) det(a 11 A 12A 1 22 A 21) = det(a 22) det(a 11) det(i m A 1 11 A 12A 1 22 A 21). and affirm (and improve) the conjectures of Drury and Lin. Determine the optimal containment region of the generalized eigenvalue λ satisfying A A12 λ x = x for some nonzero x C n. 0 A 22 A 21 0 We can then obtain the optimal eigenvalue containment region of the matrix I m A 1 11 A 12A 1 22 A 21 in case A 11 and A 22 are invertible. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 11 / 25

27 Connection to the conjectures of Drury and Lin Let B 1 = invertible. A11 0 and B 0 A 2 = 22 0 A12 and assume that B A is Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 12 / 25

28 Connection to the conjectures of Drury and Lin A11 0 Let B 1 = and B 0 A 2 = 22 invertible. Then λ satisfies λb 1x = B 2x if and only if λx = B 1 1 B2x, 0 A12 and assume that B A is for some nonzero x C n Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 12 / 25

29 Connection to the conjectures of Drury and Lin A11 0 Let B 1 = and B 0 A 2 = 22 invertible. Then λ satisfies λb 1x = B 2x if and only if λx = B A12 and assume that B A is for some nonzero x C n B2x, i.e., λ is an eigenvalue of the matrix [ ] O A 1 11 A12, B = B 1 1 B 2 = A 1 22 A21 O Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 12 / 25

30 Connection to the conjectures of Drury and Lin A11 0 Let B 1 = and B 0 A 2 = 22 invertible. Then λ satisfies λb 1x = B 2x if and only if λx = B A12 and assume that B A is for some nonzero x C n B2x, i.e., λ is an eigenvalue of the matrix [ ] O A 1 11 A12, B = B 1 1 B 2 = A 1 22 A21 O whose eigenvalues have the form ±λ 1,..., ±λ m, 0,..., 0. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 12 / 25

31 Connection to the conjectures of Drury and Lin A11 0 Let B 1 = and B 0 A 2 = 22 invertible. Then λ satisfies λb 1x = B 2x if and only if λx = B A12 and assume that B A is for some nonzero x C n B2x, i.e., λ is an eigenvalue of the matrix [ ] O A 1 11 A12, B = B 1 1 B 2 = A 1 22 A21 O whose eigenvalues have the form ±λ 1,..., ±λ m, 0,..., 0. The matrix B 2 = [ A 1 11 A12A 1 O 22 A21 O A 1 22 A21A 1 11 A12 has eigenvalues: λ 2 1,..., λ 2 m, and λ 2 1,..., λ 2 m, 0,..., 0. ] Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 12 / 25

32 Connection to the conjectures of Drury and Lin A11 0 Let B 1 = and B 0 A 2 = 22 invertible. Then λ satisfies λb 1x = B 2x if and only if λx = B A12 and assume that B A is for some nonzero x C n B2x, i.e., λ is an eigenvalue of the matrix [ ] O A 1 11 A12, B = B 1 1 B 2 = A 1 22 A21 O whose eigenvalues have the form ±λ 1,..., ±λ m, 0,..., 0. The matrix B 2 = [ A 1 11 A12A 1 O 22 A21 O A 1 22 A21A 1 11 A12 has eigenvalues: λ 2 1,..., λ 2 m, and λ 2 1,..., λ 2 m, 0,..., 0. Thus, we need only estimate the eigenvalues of B to get information about those of B 2, and det(a) = det(a 11) det(a 22) det(i m A 1 11 A 12A 1 22 A 21). ] Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 12 / 25

33 Connection to the conjectures of Drury and Lin A11 0 Let B 1 = and B 0 A 2 = 22 invertible. Then λ satisfies λb 1x = B 2x if and only if λx = B A12 and assume that B A is for some nonzero x C n B2x, i.e., λ is an eigenvalue of the matrix [ ] O A 1 11 A12, B = B 1 1 B 2 = A 1 22 A21 O whose eigenvalues have the form ±λ 1,..., ±λ m, 0,..., 0. The matrix B 2 = [ A 1 11 A12A 1 O 22 A21 O A 1 22 A21A 1 11 A12 has eigenvalues: λ 2 1,..., λ 2 m, and λ 2 1,..., λ 2 m, 0,..., 0. Thus, we need only estimate the eigenvalues of B to get information about those of B 2, and det(a) = det(a 11) det(a 22) det(i m A 1 11 A 12A 1 22 A 21). The rest is hammering out the details. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 12 / 25 ]

34 Main result Consider the generalized eigenvalue problem A A12 λ x = x for some nonzero x C 0 A 22 A 21 0 n. (1) Li and Sze [JMAA 410: (2014)] (Submitted 5 July 2013) A11 A 12 Let A = M A 21 A n with A 11 M m be such that m n/2, and 22 W (A) be a subset of a sector of half angle α [0, π/2). (a) Suppose A 11 A 22 is singular, and x C n is a nonzero vector in the kernel. Then Ax = 0 and (1) holds for every λ C with this nonzero vector x. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 13 / 25

35 Main result Consider the generalized eigenvalue problem A A12 λ x = x for some nonzero x C 0 A 22 A 21 0 n. (1) Li and Sze [JMAA 410: (2014)] (Submitted 5 July 2013) A11 A 12 Let A = M A 21 A n with A 11 M m be such that m n/2, and 22 W (A) be a subset of a sector of half angle α [0, π/2). (a) Suppose A 11 A 22 is singular, and x C n is a nonzero vector in the kernel. Then Ax = 0 and (1) holds for every λ C with this nonzero vector x. (b) Suppose A 11 and A 22 are invertible. If λ C satisfies (1), then { } 1 λ 2 R = re i2φ 2(cos(2φ) cos(2α)) : 0 r, α φ α. sin 2 (2α) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 13 / 25

36 Optimality of the containment region Let λ C be such that λ 2 = 1 re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). y re i2φ x

37 Optimality of the containment region Let λ C be such that λ 2 = 1 re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). y re i2φ x

38 Optimality of the containment region Let λ C be such that λ 2 = 1 re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). Define a = cot θ sin φ and b = tan θ cos φ and y e iφ a + ib A = a + ib e iφ. re i2φ x

39 Optimality of the containment region Let λ C be such that λ 2 = 1 re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). Define a = cot θ sin φ and b = tan θ cos φ and y e iφ a + ib re i2φ A = a + ib e iφ. r 1e iθ Notice that A is normal and has eigenvalues e iφ + (a + ib) = r 1e iθ, r 1 0, e iφ (a + ib) = r 2e iθ, r 2 0. r 2e iθ x

40 Optimality of the containment region Let λ C be such that λ 2 = 1 re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). Define a = cot θ sin φ and b = tan θ cos φ and y e iφ a + ib re i2φ A = a + ib e iφ. Notice that A is normal and has eigenvalues Therefore, e iφ + (a + ib) = r 1e iθ, r 1 0, e iφ (a + ib) = r 2e iθ, r 2 0. r 2e iθ r 1e iθ W (A) x W (A) = {αr 1e iθ + (1 α)r 2e iθ : α [0, 1]} S θ. Finally, λ = ±(a + ib)e iθ satisfy e iφ 0 0 a + ib λ 0 e iφ x = x. a + ib 0

41 Optimality of the containment region Let λ C be such that λ 2 = 1 re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). Define a = cot θ sin φ and b = tan θ cos φ and y e iφ a + ib re i2φ A = a + ib e iφ. Notice that A is normal and has eigenvalues Therefore, e iφ + (a + ib) = r 1e iθ, r 1 0, e iφ (a + ib) = r 2e iθ, r θ r 2e iθ r 1e iθ W (A) x W (A) = {αr 1e iθ + (1 α)r 2e iθ : α [0, 1]} S θ. Finally, λ = ±(a + ib)e iθ satisfy e iφ 0 0 a + ib λ 0 e iφ x = x. a + ib 0 Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 14 / 25

42 Optimality of the containment region Example Let λ C be such that 1 λ 2 = re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 15 / 25

43 Optimality of the containment region Example Let λ C be such that 1 λ 2 = re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 15 / 25

44 Optimality of the containment region Example Let λ C be such that 1 λ 2 = re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). Let ( ) e iφ a + ib A = I m a + ib e iφ (e iφ I n 2m) with a = cot θ sin φ and b = tan θ cos φ. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 15 / 25

45 Optimality of the containment region Example Let λ C be such that 1 λ 2 = re i2φ with r 2(cos 2φ cos 2α)/ sin 2 (2α). Suppose θ ( φ, α] satisfies r = 2(cos(2φ) cos(2θ))/ sin 2 (2θ). Let ( ) e iφ a + ib A = I m a + ib e iφ (e iφ I n 2m) with a = cot θ sin φ and b = tan θ cos φ. Then W (A) lying in a sector of half angle θ α, and A has eigenvalues: e iφ + (a + ib) = r 1e iθ, r 1 0, with multiplicity m, e iφ (a + ib) = r 2e iθ, r 2 0, with multiplicity m, and e iφ with multiplicity n 2m, all in S α. Moreover, λ = ±(a + ib)e iφ satisfy the generalized eigenvalue problem (1), and 1 λ 2 = re 2iφ. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 15 / 25

46 Optimality of the containment region All eigenvalues of the matrix C = I m A 1 11 A12A 1 22 A21 lie in the set { } re i2φ 2(cos(2φ) cos(2α)) : 0 r, φ α. sin 2 (2α) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 16 / 25

47 Optimality of the containment region All eigenvalues of the matrix C = I m A 1 11 A12A 1 22 A21 lie in the set { } re i2φ 2(cos(2φ) cos(2α)) : 0 r, φ α. sin 2 (2α) Thus, the spectral radius of m m matrix C is bounded by 2(cos(2φ) cos(2α)) 2(1 cos(2α)) max = = sec 2 (α), φ α sin 2 (2α) sin 2 (2α) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 16 / 25

48 Optimality of the containment region All eigenvalues of the matrix C = I m A 1 11 A12A 1 22 A21 lie in the set { } re i2φ 2(cos(2φ) cos(2α)) : 0 r, φ α. sin 2 (2α) Thus, the spectral radius of m m matrix C is bounded by 2(cos(2φ) cos(2α)) 2(1 cos(2α)) max = = sec 2 (α), φ α sin 2 (2α) sin 2 (2α) and hence det(c) sec 2m (α). So det(a) = det(a 11) det(a 22)det(C) sec 2m (α) det(a 11) det(a 22). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 16 / 25

49 Optimality of the containment region All eigenvalues of the matrix C = I m A 1 11 A12A 1 22 A21 lie in the set { } re i2φ 2(cos(2φ) cos(2α)) : 0 r, φ α. sin 2 (2α) Thus, the spectral radius of m m matrix C is bounded by 2(cos(2φ) cos(2α)) 2(1 cos(2α)) max = = sec 2 (α), φ α sin 2 (2α) sin 2 (2α) and hence det(c) sec 2m (α). So det(a) = det(a 11) det(a 22)det(C) sec 2m (α) det(a 11) det(a 22). By continuity, one can remove the invertibility assumption on A 11 A 22. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 16 / 25

50 Consequences - confirmation of Conjectures 1 and 2 Corollary 1 [Conjecture 2] A11 A 12 Let A = M A 21 A n with A 11 M m such that m n/2. 22 Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 17 / 25

51 Consequences - confirmation of Conjectures 1 and 2 Corollary 1 [Conjecture 2] A11 A 12 Let A = M A 21 A n with A 11 M m such that m n/2. 22 Suppose W (A) is a subset of a sector of half angle α [0, π/2). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 17 / 25

52 Consequences - confirmation of Conjectures 1 and 2 Corollary 1 [Conjecture 2] A11 A 12 Let A = M A 21 A n with A 11 M m such that m n/2. 22 Suppose W (A) is a subset of a sector of half angle α [0, π/2). Then det(a) sec 2m (α) det(a 11) det(a 22). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 17 / 25

53 Consequences - confirmation of Conjectures 1 and 2 Corollary 1 [Conjecture 2] A11 A 12 Let A = M A 21 A n with A 11 M m such that m n/2. 22 Suppose W (A) is a subset of a sector of half angle α [0, π/2). Then det(a) sec 2m (α) det(a 11) det(a 22). If A 11 and A 22 are invertible, then the eigenvalues of the matrix lies in the region { R = C = I m A 1 11 A 12A 1 22 A 21 M m re i2φ : 0 r } 2(cos(2φ) cos(2α)), α φ α. sin 2 (2α) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 17 / 25

54 Consequences - confirmation of Conjectures 1 and 2 Corollary 1 [Conjecture 2] A11 A 12 Let A = M A 21 A n with A 11 M m such that m n/2. 22 Suppose W (A) is a subset of a sector of half angle α [0, π/2). Then det(a) sec 2m (α) det(a 11) det(a 22). If A 11 and A 22 are invertible, then the eigenvalues of the matrix lies in the region { R = C = I m A 1 11 A 12A 1 22 A 21 M m re i2φ : 0 r } 2(cos(2φ) cos(2α)), α φ α. sin 2 (2α) In particular, if A 11 M 1 and A 22 M n 1 are invertible, then det(a)/[det(a 11) det(a 22)] = det(c) R. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 17 / 25

55 Consequences - confirmation of Conjectures 1 and 2 Corollary 2 [Conjecture 1] Suppose A satisfies the hypothesis of Corollary 1 with α = π/4. Then det(a) 2 m det(a 11) det(a 22). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 18 / 25

56 Consequences - confirmation of Conjectures 1 and 2 Corollary 2 [Conjecture 1] Suppose A satisfies the hypothesis of Corollary 1 with α = π/4. Then det(a) 2 m det(a 11) det(a 22). If A 11 and A 22 are invertible, then the eigenvalues of the matrix C = I m A 1 11 A 12A 1 22 A 21 M m lies in the set {z C : z 1 1}. We thank Dr. Zejun Huang for some discussion. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 18 / 25

57 17-23 Tokyo, Japan When we were preparing the final version of our paper, we learned that Drury had independently affirmed Conjecture 2 for α [0, π/4]. [LAMA, accepted] (Submitted 18 June, 2013) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 19 / 25

58 17-23 Tokyo, Japan When we were preparing the final version of our paper, we learned that Drury had independently affirmed Conjecture 2 for α [0, π/4]. [LAMA, accepted] (Submitted 18 June, 2013) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 19 / 25

59 17-23 Tokyo, Japan When we were preparing the final version of our paper, we learned that Drury had independently affirmed Conjecture 2 for α [0, π/4]. [LAMA, accepted] (Submitted 18 June, 2013) Finally, we finished our manuscript and uploaded to ArXiv on 22th June, Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 19 / 25

60 17-23 Tokyo, Japan When we were preparing the final version of our paper, we learned that Drury had independently affirmed Conjecture 2 for α [0, π/4]. [LAMA, accepted] (Submitted 18 June, 2013) Finally, we finished our manuscript and uploaded to ArXiv on 22th June, Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 19 / 25

61 Singular value inequalities Let σ 1(X) σ n(x) 0 be the singular values of an n n matrix X. Drury and Lin [OAM, accepted] (Submitted 10 July 2013) A11 A 12 Let A = M A 21 A n with A 11 M m. Suppose W (A) is a subset of a 22 sector of half angle α [0, π/2). Then σ j ( A11 A 12A 1 22 A 21 ) sec 2 (α) σ j(a 11) j = 1,..., m. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 20 / 25

62 Singular value inequalities Let σ 1(X) σ n(x) 0 be the singular values of an n n matrix X. Drury and Lin [OAM, accepted] (Submitted 10 July 2013) A11 A 12 Let A = M A 21 A n with A 11 M m. Suppose W (A) is a subset of a 22 sector of half angle α [0, π/2). Then σ j ( A11 A 12A 1 22 A 21 ) sec 2 (α) σ j(a 11) j = 1,..., m. Notice that det ( A11 A 12A 1 22 A 21 ) = m j=1 σ j ( A11 A 12A 1 22 A 21 ) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 20 / 25

63 Singular value inequalities Let σ 1(X) σ n(x) 0 be the singular values of an n n matrix X. Drury and Lin [OAM, accepted] (Submitted 10 July 2013) A11 A 12 Let A = M A 21 A n with A 11 M m. Suppose W (A) is a subset of a 22 sector of half angle α [0, π/2). Then σ j ( A11 A 12A 1 22 A 21 ) sec 2 (α) σ j(a 11) j = 1,..., m. Notice that det ( A11 A 12A 1 22 A 21 ) = m j=1 σ j ( A11 A 12A 1 22 A 21 ) m sec 2 (α) σ j(a 11) = sec 2m (α) det(a 11). j=1 Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 20 / 25

64 Singular value inequalities Let σ 1(X) σ n(x) 0 be the singular values of an n n matrix X. Drury and Lin [OAM, accepted] (Submitted 10 July 2013) A11 A 12 Let A = M A 21 A n with A 11 M m. Suppose W (A) is a subset of a 22 sector of half angle α [0, π/2). Then σ j ( A11 A 12A 1 22 A 21 ) sec 2 (α) σ j(a 11) j = 1,..., m. Notice that det ( A11 A 12A 1 22 A 21 ) = Then m j=1 σ j ( A11 A 12A 1 22 A 21 ) m sec 2 (α) σ j(a 11) = sec 2m (α) det(a 11). j=1 det A = det ( A11 A 12A 1 22 A 21 ) det(a22) sec 2m (α) det(a 11) det(a 22), which also confirmed Conjecture 2. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 20 / 25

65 More inequalities A norm is unitarily invariant on M n if A = UAV for any A M n and any unitaries U and V M n. Corollary - Norm Inequality A11 A 12 Let A = M A 21 A n with A 11 M m. Suppose W (A) is a subset of a 22 sector of half angle α [0, π/2). Then A/A 11 sec 2 (α) A 22 for any unitarily invariant norm. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 21 / 25

66 More inequalities A norm is unitarily invariant on M n if A = UAV for any A M n and any unitaries U and V M n. Corollary - Norm Inequality A11 A 12 Let A = M A 21 A n with A 11 M m. Suppose W (A) is a subset of a 22 sector of half angle α [0, π/2). Then A/A 11 sec 2 (α) A 22 for any unitarily invariant norm. One of the major steps: If W (A) S α, then A = Xdiag ( e iθ 1,..., e iθn) X with θ j α, j = 1,..., n for some invertible X M n. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 21 / 25

67 More inequalities Fan and Hoffman Inequality For any A M n, λ j(re A) σ j(a) j = 1,..., n. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 22 / 25

68 More inequalities Fan and Hoffman Inequality For any A M n, λ j(re A) σ j(a) j = 1,..., n. Reversed Inequality Let A M n be such that W (A) S α. Then λ j(a) sec 2 (α)(re A) j = 1,..., n. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 22 / 25

69 More inequalities Fan and Hoffman Inequality For any A M n, λ j(re A) σ j(a) j = 1,..., n. Reversed Inequality Let A M n be such that W (A) S α. Then λ j(a) sec 2 (α)(re A) j = 1,..., n. F. Zhang [LAMA, accepted] (submitted 28 October 2013) Let A M n be such that W (A) S α. Then λ j(a) w sec(α)(re A). Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 22 / 25

70 More norm inequalities Lin and Zhou [JMAA 407: (2013)] (submitted 27 February 2013) Let A M n be accretive-dissipative. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } 4 A 11 A 22. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 23 / 25

71 More norm inequalities Lin and Zhou [JMAA 407: (2013)] (submitted 27 February 2013) Let A M n be accretive-dissipative. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } 4 A 11 A 22. Conjecture Let A M n be accretive-dissipative. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } 2 A 11 A 22. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 23 / 25

72 More norm inequalities Lin and Zhou [JMAA 407: (2013)] (submitted 27 February 2013) Let A M n be accretive-dissipative. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } 4 A 11 A 22. Conjecture Let A M n be accretive-dissipative. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } 2 A 11 A 22. The conjecture is confirmed by Y. Zhang. [JMAA: 412: (2014)] (Submitted 26 August 2013) Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 23 / 25

73 More norm inequalities Lin and Zhou [JMAA 407: (2013)] (submitted 27 February 2013) Let A M n be accretive-dissipative. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } 4 A 11 A 22. Conjecture Let A M n be accretive-dissipative. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } 2 A 11 A 22. The conjecture is confirmed by Y. Zhang. [JMAA: 412: (2014)] Let Then A = [ ] max{ A 12 2, A 21 2 } = 2 = 2 A 11 A 22. y (Submitted 26 August 2013) W (A) x Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 23 / 25

74 More norm inequalities F. Zhang [LAMA, accepted] (submitted 28 October 2013) Let A M n be such that W (A) S α. Then for any unitarily invariant norm, max{ A 12 2, A 21 2 } sec 2 (α) A 11 A 22. Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 24 / 25

75 Raymond Nung-Sing Sze WONRA2014, Sanya, China P. 25 / 25

arxiv: v1 [math.ra] 8 Apr 2016

arxiv: v1 [math.ra] 8 Apr 2016 ON A DETERMINANTAL INEQUALITY ARISING FROM DIFFUSION TENSOR IMAGING MINGHUA LIN arxiv:1604.04141v1 [math.ra] 8 Apr 2016 Abstract. In comparing geodesics induced by different metrics, Audenaert formulated

More information

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation

More information

University, Fort Lauderdale, FL, USA. Published online: 30 Jun To link to this article:

University, Fort Lauderdale, FL, USA. Published online: 30 Jun To link to this article: This article was downloaded by: [172.4.212.45] On: 12 August 2014, At: 11:32 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer

More information

Linear and Multilinear Algebra. Linear maps preserving rank of tensor products of matrices

Linear and Multilinear Algebra. Linear maps preserving rank of tensor products of matrices Linear maps preserving rank of tensor products of matrices Journal: Manuscript ID: GLMA-0-0 Manuscript Type: Original Article Date Submitted by the Author: -Aug-0 Complete List of Authors: Zheng, Baodong;

More information

Eigenvalues of the sum of matrices. from unitary similarity orbits

Eigenvalues of the sum of matrices. from unitary similarity orbits Eigenvalues of the sum of matrices from unitary similarity orbits Department of Mathematics The College of William and Mary Based on some joint work with: Yiu-Tung Poon (Iowa State University), Nung-Sing

More information

Math Homework 8 (selected problems)

Math Homework 8 (selected problems) Math 102 - Homework 8 (selected problems) David Lipshutz Problem 1. (Strang, 5.5: #14) In the list below, which classes of matrices contain A and which contain B? 1 1 1 1 A 0 0 1 0 0 0 0 1 and B 1 1 1

More information

Elementary linear algebra

Elementary linear algebra Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The

More information

HIGHER RANK NUMERICAL RANGES OF NORMAL MATRICES

HIGHER RANK NUMERICAL RANGES OF NORMAL MATRICES HIGHER RANK NUMERICAL RANGES OF NORMAL MATRICES HWA-LONG GAU, CHI-KWONG LI, YIU-TUNG POON, AND NUNG-SING SZE Abstract. The higher rank numerical range is closely connected to the construction of quantum

More information

LINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES

LINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES LINEAR MAPS PRESERVING NUMERICAL RADIUS OF TENSOR PRODUCTS OF MATRICES AJDA FOŠNER, ZEJUN HUANG, CHI-KWONG LI, AND NUNG-SING SZE Abstract. Let m, n 2 be positive integers. Denote by M m the set of m m

More information

Math 108b: Notes on the Spectral Theorem

Math 108b: Notes on the Spectral Theorem Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator

More information

Workshop on Matrices and Operators. July 13, 2007 T5, Meng Wah Complex The University of Hong Kong

Workshop on Matrices and Operators. July 13, 2007 T5, Meng Wah Complex The University of Hong Kong Workshop on Matrices and Operators July 13, 2007 T5, Meng Wah Complex The University of Hong Kong Morning Session (Chaired by Yiu-Tung Poon) 09:00 09:45 Bit-Shun Tam, Tamkang University Exponents of K-primitive

More information

HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION)

HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) PROFESSOR STEVEN MILLER: BROWN UNIVERSITY: SPRING 2007 1. CHAPTER 1: MATRICES AND GAUSSIAN ELIMINATION Page 9, # 3: Describe

More information

Math 408 Advanced Linear Algebra

Math 408 Advanced Linear Algebra Math 408 Advanced Linear Algebra Chi-Kwong Li Chapter 4 Hermitian and symmetric matrices Basic properties Theorem Let A M n. The following are equivalent. Remark (a) A is Hermitian, i.e., A = A. (b) x

More information

Abstract. In this article, several matrix norm inequalities are proved by making use of the Hiroshima 2003 result on majorization relations.

Abstract. In this article, several matrix norm inequalities are proved by making use of the Hiroshima 2003 result on majorization relations. HIROSHIMA S THEOREM AND MATRIX NORM INEQUALITIES MINGHUA LIN AND HENRY WOLKOWICZ Abstract. In this article, several matrix norm inequalities are proved by making use of the Hiroshima 2003 result on majorization

More information

Singular Value Inequalities for Real and Imaginary Parts of Matrices

Singular Value Inequalities for Real and Imaginary Parts of Matrices Filomat 3:1 16, 63 69 DOI 1.98/FIL16163C Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Singular Value Inequalities for Real Imaginary

More information

An Even Order Symmetric B Tensor is Positive Definite

An Even Order Symmetric B Tensor is Positive Definite An Even Order Symmetric B Tensor is Positive Definite Liqun Qi, Yisheng Song arxiv:1404.0452v4 [math.sp] 14 May 2014 October 17, 2018 Abstract It is easily checkable if a given tensor is a B tensor, or

More information

Properties of Linear Transformations from R n to R m

Properties of Linear Transformations from R n to R m Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation

More information

Recent study in preserver problems. Chi-Kwong Li Department of Mathematics College of William and Mary Williamsburg, VA 23185

Recent study in preserver problems. Chi-Kwong Li Department of Mathematics College of William and Mary Williamsburg, VA 23185 Recent study in preserver problems Chi-Kwong Li Department of Mathematics College of William and Mary Williamsburg, VA 23185 1 1. Preserver Problems on Matrices/Operators Let M be a matrix space or matrix

More information

Linear Algebra: Matrix Eigenvalue Problems

Linear Algebra: Matrix Eigenvalue Problems CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given

More information

PRODUCT OF OPERATORS AND NUMERICAL RANGE

PRODUCT OF OPERATORS AND NUMERICAL RANGE PRODUCT OF OPERATORS AND NUMERICAL RANGE MAO-TING CHIEN 1, HWA-LONG GAU 2, CHI-KWONG LI 3, MING-CHENG TSAI 4, KUO-ZHONG WANG 5 Abstract. We show that a bounded linear operator A B(H) is a multiple of a

More information

MINIMAL NORMAL AND COMMUTING COMPLETIONS

MINIMAL NORMAL AND COMMUTING COMPLETIONS INTERNATIONAL JOURNAL OF INFORMATION AND SYSTEMS SCIENCES Volume 4, Number 1, Pages 5 59 c 8 Institute for Scientific Computing and Information MINIMAL NORMAL AND COMMUTING COMPLETIONS DAVID P KIMSEY AND

More information

Numerical Ranges of the Powers of an Operator

Numerical Ranges of the Powers of an Operator Numerical Ranges of the Powers of an Operator Man-Duen Choi 1 Department of Mathematics, University of Toronto, Toronto, Ontario, Canada M5S 2E4 E-mail: choi@math.toronto.edu Chi-Kwong Li 2 Department

More information

Analysis Preliminary Exam Workshop: Hilbert Spaces

Analysis Preliminary Exam Workshop: Hilbert Spaces Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H

More information

Math 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination

Math 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column

More information

AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)

AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) AMS526: Numerical Analysis (Numerical Linear Algebra for Computational and Data Sciences) Lecture 14: Eigenvalue Problems; Eigenvalue Revealing Factorizations Xiangmin Jiao Stony Brook University Xiangmin

More information

Lecture 5. Ch. 5, Norms for vectors and matrices. Norms for vectors and matrices Why?

Lecture 5. Ch. 5, Norms for vectors and matrices. Norms for vectors and matrices Why? KTH ROYAL INSTITUTE OF TECHNOLOGY Norms for vectors and matrices Why? Lecture 5 Ch. 5, Norms for vectors and matrices Emil Björnson/Magnus Jansson/Mats Bengtsson April 27, 2016 Problem: Measure size of

More information

Review of Some Concepts from Linear Algebra: Part 2

Review of Some Concepts from Linear Algebra: Part 2 Review of Some Concepts from Linear Algebra: Part 2 Department of Mathematics Boise State University January 16, 2019 Math 566 Linear Algebra Review: Part 2 January 16, 2019 1 / 22 Vector spaces A set

More information

Matrix Inequalities by Means of Block Matrices 1

Matrix Inequalities by Means of Block Matrices 1 Mathematical Inequalities & Applications, Vol. 4, No. 4, 200, pp. 48-490. Matrix Inequalities by Means of Block Matrices Fuzhen Zhang 2 Department of Math, Science and Technology Nova Southeastern University,

More information

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

Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.

More information

Numerical radius, entropy and preservers on quantum states

Numerical radius, entropy and preservers on quantum states 2th WONRA, TSIMF, Sanya, July 204 Numerical radius, entropy and preservers on quantum states Kan He (å Ä) Joint with Pro. Jinchuan Hou, Chi-Kwong Li and Zejun Huang Taiyuan University of Technology, Taiyuan,

More information

AMS526: Numerical Analysis I (Numerical Linear Algebra)

AMS526: Numerical Analysis I (Numerical Linear Algebra) AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 16: Eigenvalue Problems; Similarity Transformations Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 18 Eigenvalue

More information

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 )

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 ) Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O

More information

Compression, Matrix Range and Completely Positive Map

Compression, Matrix Range and Completely Positive Map Compression, Matrix Range and Completely Positive Map Iowa State University Iowa-Nebraska Functional Analysis Seminar November 5, 2016 Definitions and notations H, K : Hilbert space. If dim H = n

More information

Convexity of the Joint Numerical Range

Convexity of the Joint Numerical Range Convexity of the Joint Numerical Range Chi-Kwong Li and Yiu-Tung Poon October 26, 2004 Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement. Abstract Let A = (A 1,..., A m ) be an

More information

LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM

LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator

More information

Part IA. Vectors and Matrices. Year

Part IA. Vectors and Matrices. Year Part IA Vectors and Matrices Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2018 Paper 1, Section I 1C Vectors and Matrices For z, w C define the principal value of z w. State de Moivre s

More information

Therefore, A and B have the same characteristic polynomial and hence, the same eigenvalues.

Therefore, A and B have the same characteristic polynomial and hence, the same eigenvalues. Similar Matrices and Diagonalization Page 1 Theorem If A and B are n n matrices, which are similar, then they have the same characteristic equation and hence the same eigenvalues. Proof Let A and B be

More information

Yimin Wei a,b,,1, Xiezhang Li c,2, Fanbin Bu d, Fuzhen Zhang e. Abstract

Yimin Wei a,b,,1, Xiezhang Li c,2, Fanbin Bu d, Fuzhen Zhang e. Abstract Linear Algebra and its Applications 49 (006) 765 77 wwwelseviercom/locate/laa Relative perturbation bounds for the eigenvalues of diagonalizable and singular matrices Application of perturbation theory

More information

determinantal conjecture, Marcus-de Oliveira, determinants, normal matrices,

determinantal conjecture, Marcus-de Oliveira, determinants, normal matrices, 1 3 4 5 6 7 8 BOUNDARY MATRICES AND THE MARCUS-DE OLIVEIRA DETERMINANTAL CONJECTURE AMEET SHARMA Abstract. We present notes on the Marcus-de Oliveira conjecture. The conjecture concerns the region in the

More information

MATRICES ARE SIMILAR TO TRIANGULAR MATRICES

MATRICES ARE SIMILAR TO TRIANGULAR MATRICES MATRICES ARE SIMILAR TO TRIANGULAR MATRICES 1 Complex matrices Recall that the complex numbers are given by a + ib where a and b are real and i is the imaginary unity, ie, i 2 = 1 In what we describe below,

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors Sec. 6.1 Eigenvalues and Eigenvectors Linear transformations L : V V that go from a vector space to itself are often called linear operators. Many linear operators can be understood geometrically by identifying

More information

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

CANONICAL FORMS, HIGHER RANK NUMERICAL RANGES, TOTALLY ISOTROPIC SUBSPACES, AND MATRIX EQUATIONS CANONICAL FORMS, HIGHER RANK NUMERICAL RANGES, TOTALLY ISOTROPIC SUBSPACES, AND MATRIX EQUATIONS CHI-KWONG LI AND NUNG-SING SZE Abstract. Results on matrix canonical forms are used to give a complete description

More information

Here is an example of a block diagonal matrix with Jordan Blocks on the diagonal: J

Here is an example of a block diagonal matrix with Jordan Blocks on the diagonal: J Class Notes 4: THE SPECTRAL RADIUS, NORM CONVERGENCE AND SOR. Math 639d Due Date: Feb. 7 (updated: February 5, 2018) In the first part of this week s reading, we will prove Theorem 2 of the previous class.

More information

THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR

THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR WEN LI AND MICHAEL K. NG Abstract. In this paper, we study the perturbation bound for the spectral radius of an m th - order n-dimensional

More information

REPRESENTATION THEORY WEEK 7

REPRESENTATION THEORY WEEK 7 REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable

More information

Symmetric and anti symmetric matrices

Symmetric and anti symmetric matrices Symmetric and anti symmetric matrices In linear algebra, a symmetric matrix is a square matrix that is equal to its transpose. Formally, matrix A is symmetric if. A = A Because equal matrices have equal

More information

arxiv: v3 [math.ra] 22 Aug 2014

arxiv: v3 [math.ra] 22 Aug 2014 arxiv:1407.0331v3 [math.ra] 22 Aug 2014 Positivity of Partitioned Hermitian Matrices with Unitarily Invariant Norms Abstract Chi-Kwong Li a, Fuzhen Zhang b a Department of Mathematics, College of William

More information

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

b jσ(j), Keywords: Decomposable numerical range, principal character AMS Subject Classification: 15A60 On the Hu-Hurley-Tam Conjecture Concerning The Generalized Numerical Range Che-Man Cheng Faculty of Science and Technology, University of Macau, Macau. E-mail: fstcmc@umac.mo and Chi-Kwong Li Department

More information

Symmetric Norm Inequalities And Positive Semi-Definite Block-Matrices

Symmetric Norm Inequalities And Positive Semi-Definite Block-Matrices Symmetric Norm Inequalities And Positive Semi-Definite lock-matrices Antoine Mhanna To cite this version: Antoine Mhanna Symmetric Norm Inequalities And Positive Semi-Definite lock-matrices 15

More information

Functional Analysis Review

Functional Analysis Review Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all

More information

Matrix Theory. A.Holst, V.Ufnarovski

Matrix Theory. A.Holst, V.Ufnarovski Matrix Theory AHolst, VUfnarovski 55 HINTS AND ANSWERS 9 55 Hints and answers There are two different approaches In the first one write A as a block of rows and note that in B = E ij A all rows different

More information

Operators with numerical range in a closed halfplane

Operators with numerical range in a closed halfplane Operators with numerical range in a closed halfplane Wai-Shun Cheung 1 Department of Mathematics, University of Hong Kong, Hong Kong, P. R. China. wshun@graduate.hku.hk Chi-Kwong Li 2 Department of Mathematics,

More information

Linear Algebra - Part II

Linear Algebra - Part II Linear Algebra - Part II Projection, Eigendecomposition, SVD (Adapted from Sargur Srihari s slides) Brief Review from Part 1 Symmetric Matrix: A = A T Orthogonal Matrix: A T A = AA T = I and A 1 = A T

More information

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.

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. LINEAR MAPS PRESERVING KY FAN NORMS AND SCHATTEN NORMS OF TENSOR PRODUCTS OF MATRICES AJDA FOŠNER, ZEJUN HUANG, CHI-KWONG LI, AND NUNG-SING SZE Abstract. For a positive integer n, let M n be the set of

More information

MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian.

MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. Spanning set Let S be a subset of a vector space V. Definition. The span of the set S is the smallest subspace W V that contains S. If

More information

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 2

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 2 EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 2 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 5, 2012 Andre Tkacenko

More information

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

Math Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88 Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors Chapter 1 Eigenvalues and Eigenvectors Among problems in numerical linear algebra, the determination of the eigenvalues and eigenvectors of matrices is second in importance only to the solution of linear

More information

Lecture 5. Complex Numbers and Euler s Formula

Lecture 5. Complex Numbers and Euler s Formula Lecture 5. Complex Numbers and Euler s Formula University of British Columbia, Vancouver Yue-Xian Li March 017 1 Main purpose: To introduce some basic knowledge of complex numbers to students so that they

More information

Math 4A Notes. Written by Victoria Kala Last updated June 11, 2017

Math 4A Notes. Written by Victoria Kala Last updated June 11, 2017 Math 4A Notes Written by Victoria Kala vtkala@math.ucsb.edu Last updated June 11, 2017 Systems of Linear Equations A linear equation is an equation that can be written in the form a 1 x 1 + a 2 x 2 +...

More information

OHSx XM511 Linear Algebra: Solutions to Online True/False Exercises

OHSx XM511 Linear Algebra: Solutions to Online True/False Exercises This document gives the solutions to all of the online exercises for OHSx XM511. The section ( ) numbers refer to the textbook. TYPE I are True/False. Answers are in square brackets [. Lecture 02 ( 1.1)

More information

1. General Vector Spaces

1. General Vector Spaces 1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule

More information

Extremal numbers of positive entries of imprimitive nonnegative matrices

Extremal numbers of positive entries of imprimitive nonnegative matrices Extremal numbers of positive entries of imprimitive nonnegative matrices Xingzhi Zhan East China Normal University Shanghai, China e-mail: zhan@math.ecnu.edu.cn A square nonnegative matrix A is said to

More information

Diagonalizing Matrices

Diagonalizing Matrices Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B

More information

Math 489AB Exercises for Chapter 2 Fall Section 2.3

Math 489AB Exercises for Chapter 2 Fall Section 2.3 Math 489AB Exercises for Chapter 2 Fall 2008 Section 2.3 2.3.3. Let A M n (R). Then the eigenvalues of A are the roots of the characteristic polynomial p A (t). Since A is real, p A (t) is a polynomial

More information

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent. Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u

More information

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear

More information

The eigenvalue map and spectral sets

The eigenvalue map and spectral sets The eigenvalue map and spectral sets p. 1/30 The eigenvalue map and spectral sets M. Seetharama Gowda Department of Mathematics and Statistics University of Maryland, Baltimore County Baltimore, Maryland

More information

NOTES ON THE NUMERICAL RANGE

NOTES ON THE NUMERICAL RANGE NOTES ON THE NUMERICAL RANGE JOEL H. SHAPIRO Abstract. This is an introduction to the notion of numerical range for bounded linear operators on Hilbert space. The main results are: determination of the

More information

University of Colorado at Denver Mathematics Department Applied Linear Algebra Preliminary Exam With Solutions 16 January 2009, 10:00 am 2:00 pm

University of Colorado at Denver Mathematics Department Applied Linear Algebra Preliminary Exam With Solutions 16 January 2009, 10:00 am 2:00 pm University of Colorado at Denver Mathematics Department Applied Linear Algebra Preliminary Exam With Solutions 16 January 2009, 10:00 am 2:00 pm Name: The proctor will let you read the following conditions

More information

Review of similarity transformation and Singular Value Decomposition

Review of similarity transformation and Singular Value Decomposition Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm

More information

Generalized Interlacing Inequalities

Generalized Interlacing Inequalities Generalized Interlacing Inequalities Chi-Kwong Li, Yiu-Tung Poon, and Nung-Sing Sze In memory of Professor Ky Fan Abstract We discuss some applications of generalized interlacing inequalities of Ky Fan

More information

R, 1 i 1,i 2,...,i m n.

R, 1 i 1,i 2,...,i m n. SIAM J. MATRIX ANAL. APPL. Vol. 31 No. 3 pp. 1090 1099 c 2009 Society for Industrial and Applied Mathematics FINDING THE LARGEST EIGENVALUE OF A NONNEGATIVE TENSOR MICHAEL NG LIQUN QI AND GUANGLU ZHOU

More information

4. Determinants.

4. Determinants. 4. Determinants 4.1. Determinants; Cofactor Expansion Determinants of 2 2 and 3 3 Matrices 2 2 determinant 4.1. Determinants; Cofactor Expansion Determinants of 2 2 and 3 3 Matrices 3 3 determinant 4.1.

More information

1 Linear Algebra Problems

1 Linear Algebra Problems Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and

More information

Ranks and determinants of the sum of matrices from unitary orbits

Ranks and determinants of the sum of matrices from unitary orbits Ranks and determinants of the sum of matrices from unitary orbits Chi-Kwong Li, Yiu-Tung Poon and Nung-Sing Sze Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his 70th birthday. Abstract The

More information

j=1 x j p, if 1 p <, x i ξ : x i < ξ} 0 as p.

j=1 x j p, if 1 p <, x i ξ : x i < ξ} 0 as p. LINEAR ALGEBRA Fall 203 The final exam Almost all of the problems solved Exercise Let (V, ) be a normed vector space. Prove x y x y for all x, y V. Everybody knows how to do this! Exercise 2 If V is a

More information

Some inequalities for unitarily invariant norms of matrices

Some inequalities for unitarily invariant norms of matrices Wang et al Journal of Inequalities and Applications 011, 011:10 http://wwwjournalofinequalitiesandapplicationscom/content/011/1/10 RESEARCH Open Access Some inequalities for unitarily invariant norms of

More information

Introduction to Matrix Algebra

Introduction to Matrix Algebra Introduction to Matrix Algebra August 18, 2010 1 Vectors 1.1 Notations A p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the line. When p

More information

Eigenvalue and Eigenvector Problems

Eigenvalue and Eigenvector Problems Eigenvalue and Eigenvector Problems An attempt to introduce eigenproblems Radu Trîmbiţaş Babeş-Bolyai University April 8, 2009 Radu Trîmbiţaş ( Babeş-Bolyai University) Eigenvalue and Eigenvector Problems

More information

Decomposable Numerical Ranges on Orthonormal Tensors

Decomposable Numerical Ranges on Orthonormal Tensors Decomposable Numerical Ranges on Orthonormal Tensors Chi-Kwong Li Department of Mathematics, College of William and Mary, P.O.Box 8795, Williamsburg, Virginia 23187-8795, USA. E-mail: ckli@math.wm.edu

More information

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 4

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 4 EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 4 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 12, 2012 Andre Tkacenko

More information

Linear Algebra Highlights

Linear Algebra Highlights Linear Algebra Highlights Chapter 1 A linear equation in n variables is of the form a 1 x 1 + a 2 x 2 + + a n x n. We can have m equations in n variables, a system of linear equations, which we want to

More information

Isometries Between Matrix Algebras

Isometries Between Matrix Algebras Abstract Isometries Between Matrix Algebras Wai-Shun Cheung, Chi-Kwong Li 1 and Yiu-Tung Poon As an attempt to understand linear isometries between C -algebras without the surjectivity assumption, we study

More information

On a decomposition lemma for positive semi-definite block-matrices

On a decomposition lemma for positive semi-definite block-matrices On a decomposition lemma for positive semi-definite bloc-matrices arxiv:10.0473v1 [math.fa] Feb 01 Jean-Christophe ourin, Eun-Young Lee, Minghua Lin January, 01 Abstract This short note, in part of expository

More information

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

c c c c c c c c c c a 3x3 matrix C= has a determinant determined by Linear Algebra Determinants and Eigenvalues Introduction: Many important geometric and algebraic properties of square matrices are associated with a single real number revealed by what s known as the determinant.

More information

Characterization of half-radial matrices

Characterization of half-radial matrices Characterization of half-radial matrices Iveta Hnětynková, Petr Tichý Faculty of Mathematics and Physics, Charles University, Sokolovská 83, Prague 8, Czech Republic Abstract Numerical radius r(a) is the

More information

arxiv: v1 [math.ra] 11 Aug 2014

arxiv: v1 [math.ra] 11 Aug 2014 Double B-tensors and quasi-double B-tensors Chaoqian Li, Yaotang Li arxiv:1408.2299v1 [math.ra] 11 Aug 2014 a School of Mathematics and Statistics, Yunnan University, Kunming, Yunnan, P. R. China 650091

More information

1. Select the unique answer (choice) for each problem. Write only the answer.

1. Select the unique answer (choice) for each problem. Write only the answer. MATH 5 Practice Problem Set Spring 7. Select the unique answer (choice) for each problem. Write only the answer. () Determine all the values of a for which the system has infinitely many solutions: x +

More information

An inverse numerical range problem for determinantal representations

An inverse numerical range problem for determinantal representations An inverse numerical range problem for determinantal representations Mao-Ting Chien Soochow University, Taiwan Based on joint work with Hiroshi Nakazato WONRA June 13-18, 2018, Munich Outline 1. Introduction

More information

MATH 211 review notes

MATH 211 review notes MATH 211 review notes Notes written by Mark Przepiora 1 Determinants 1. Let A = 3 1 1 2 1 2 5 2 1 (a) Compute det(a). Solution: To find the determinant, we expand along the first row. det(a) = 2 1 1 2

More information

MATH 320: PRACTICE PROBLEMS FOR THE FINAL AND SOLUTIONS

MATH 320: PRACTICE PROBLEMS FOR THE FINAL AND SOLUTIONS MATH 320: PRACTICE PROBLEMS FOR THE FINAL AND SOLUTIONS There will be eight problems on the final. The following are sample problems. Problem 1. Let F be the vector space of all real valued functions on

More information

Computational math: Assignment 1

Computational math: Assignment 1 Computational math: Assignment 1 Thanks Ting Gao for her Latex file 11 Let B be a 4 4 matrix to which we apply the following operations: 1double column 1, halve row 3, 3add row 3 to row 1, 4interchange

More information

MATHEMATICS: SPECIALIST UNITS 3C AND 3D FORMULA SHEET 2015

MATHEMATICS: SPECIALIST UNITS 3C AND 3D FORMULA SHEET 2015 MATHEMATICS: SPECIALIST UNITS 3C AND 3D FORMULA SHEET 05 Copyright School Curriculum and Standards Authority, 05 This document apart from any third party copyright material contained in it may be freely

More information

Spectral inequalities and equalities involving products of matrices

Spectral inequalities and equalities involving products of matrices Spectral inequalities and equalities involving products of matrices Chi-Kwong Li 1 Department of Mathematics, College of William & Mary, Williamsburg, Virginia 23187 (ckli@math.wm.edu) Yiu-Tung Poon Department

More information

q n. Q T Q = I. Projections Least Squares best fit solution to Ax = b. Gram-Schmidt process for getting an orthonormal basis from any basis.

q n. Q T Q = I. Projections Least Squares best fit solution to Ax = b. Gram-Schmidt process for getting an orthonormal basis from any basis. Exam Review Material covered by the exam [ Orthogonal matrices Q = q 1... ] q n. Q T Q = I. Projections Least Squares best fit solution to Ax = b. Gram-Schmidt process for getting an orthonormal basis

More information

Lecture 7: Positive Semidefinite Matrices

Lecture 7: Positive Semidefinite Matrices Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.

More information

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

Numerical Ranges and Dilations. Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement Numerical Ranges and Dilations Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement Man-Duen Choi Department of Mathematics, University of Toronto, Toronto, Ontario, Canada M5S 3G3.

More information

Linear Algebra Review. Fei-Fei Li

Linear Algebra Review. Fei-Fei Li Linear Algebra Review Fei-Fei Li 1 / 51 Vectors Vectors and matrices are just collections of ordered numbers that represent something: movements in space, scaling factors, pixel brightnesses, etc. A vector

More information

Norm Inequalities of Positive Semi-Definite Matrices

Norm Inequalities of Positive Semi-Definite Matrices Norm Inequalities of Positive Semi-Definite Matrices Antoine Mhanna To cite this version: Antoine Mhanna Norm Inequalities of Positive Semi-Definite Matrices 15 HAL Id: hal-11844 https://halinriafr/hal-11844v1

More information