Lecture 23: Trace and determinants! (1) (Final lecture)

Similar documents
Lecture 23: Determinants (1)

Lecture 19: Isometries, Positive operators, Polar and singular value decompositions; Unitary matrices and classical groups; Previews (1)

FALL 2011, SOLUTION SET 10 LAST REVISION: NOV 27, 9:45 AM. (T c f)(x) = f(x c).

Linear Algebra II. Ulrike Tillmann. January 4, 2018

Lecture 21: The decomposition theorem into generalized eigenspaces; multiplicity of eigenvalues and upper-triangular matrices (1)

Linear Systems and Matrices

MATH 320: PRACTICE PROBLEMS FOR THE FINAL AND SOLUTIONS

Equality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.

MATRICES ARE SIMILAR TO TRIANGULAR MATRICES

Lecture 1 and 2: Random Spanning Trees

Travis Schedler. Thurs, Oct 27, 2011 (version: Thurs, Oct 27, 1:00 PM)

Here are some additional properties of the determinant function.

Linear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.

MA 1B ANALYTIC - HOMEWORK SET 7 SOLUTIONS

Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1)

Determinants. Beifang Chen

MATH 511 ADVANCED LINEAR ALGEBRA SPRING 2006

Linear algebra I Homework #1 due Thursday, Oct Show that the diagonals of a square are orthogonal to one another.

Determinants An Introduction

Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.

Chapter SSM: Linear Algebra Section Fails to be invertible; since det = 6 6 = Invertible; since det = = 2.

MATH Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product.

Determinants: Uniqueness and more

Lecture 11: Finish Gaussian elimination and applications; intro to eigenvalues and eigenvectors (1)

Linear Algebra: Lecture notes from Kolman and Hill 9th edition.

Matrices and Linear Algebra

Linear Algebra Formulas. Ben Lee

Lecture 19: Polar and singular value decompositions; generalized eigenspaces; the decomposition theorem (1)

1 Last time: least-squares problems

Math 240 Calculus III

0.1 Rational Canonical Forms

5.6. PSEUDOINVERSES 101. A H w.

Math 4153 Exam 3 Review. The syllabus for Exam 3 is Chapter 6 (pages ), Chapter 7 through page 137, and Chapter 8 through page 182 in Axler.

Lecture 2: Eigenvalues and their Uses

Notes on Determinants and Matrix Inverse

Lecture 22: Jordan canonical form of upper-triangular matrices (1)

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

Lecture 18. Ramanujan Graphs continued

Determinants - Uniqueness and Properties

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

The Determinant: a Means to Calculate Volume

Common-Knowledge / Cheat Sheet

Determinants Chapter 3 of Lay

Mathematical Methods wk 2: Linear Operators

MATH 106 LINEAR ALGEBRA LECTURE NOTES

The Determinant. Chapter Definition of the Determinant

Linear Algebra Primer

II. Determinant Functions

Matrix Factorization and Analysis

4. Determinants.

Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

CP3 REVISION LECTURES VECTORS AND MATRICES Lecture 1. Prof. N. Harnew University of Oxford TT 2013

A matrix A is invertible i det(a) 6= 0.

(f + g)(s) = f(s) + g(s) for f, g V, s S (cf)(s) = cf(s) for c F, f V, s S

Math 110 Linear Algebra Midterm 2 Review October 28, 2017

Apprentice Program: Linear Algebra

Foundations of Matrix Analysis

SPRING OF 2008 D. DETERMINANTS

Review for Exam Find all a for which the following linear system has no solutions, one solution, and infinitely many solutions.

The Matrix-Tree Theorem

Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015

Online Exercises for Linear Algebra XM511

Linear Algebra Lecture Notes-II

Review problems for MA 54, Fall 2004.

SYLLABUS. 1 Linear maps and matrices

Properties of the Determinant Function

Chapter 4. Matrices and Matrix Rings

Linear Algebra Primer

Matrix Algebra Determinant, Inverse matrix. Matrices. A. Fabretti. Mathematics 2 A.Y. 2015/2016. A. Fabretti Matrices

REU 2007 Apprentice Class Lecture 8

Differential Topology Final Exam With Solutions

Linear Algebra Primer

HOMEWORK 9 solutions

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

33AH, WINTER 2018: STUDY GUIDE FOR FINAL EXAM

Formula for the inverse matrix. Cramer s rule. Review: 3 3 determinants can be computed expanding by any row or column

MATH 110: LINEAR ALGEBRA PRACTICE MIDTERM #2

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

Eigenvectors. Prop-Defn

Topics in linear algebra

NOTES FOR LINEAR ALGEBRA 133

Quantum Computing Lecture 2. Review of Linear Algebra

CHAPTER 1. Matrices Matrix Algebra

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

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

Chapter 2. Square matrices

MATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants.

Linear Algebra 2 Spectral Notes

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

Linear Algebra. Workbook

MATH 240 Spring, Chapter 1: Linear Equations and Matrices

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

Notes on Linear Algebra

88 CHAPTER 3. SYMMETRIES

1 Positive definiteness and semidefiniteness

Solution of Linear Equations

LINEAR ALGEBRA WITH APPLICATIONS

MH1200 Final 2014/2015

Calculating determinants for larger matrices

det(ka) = k n det A.

Transcription:

Lecture 23: Trace and determinants! (1) (Final lecture) Travis Schedler Thurs, Dec 9, 2010 (version: Monday, Dec 13, 3:52 PM) Goals (2) Recall χ T (x) = (x λ 1 ) (x λ n ) = x n tr(t )x n 1 + +( 1) n det(t ). tr(t ) = sum of eigenvalues. Theorem: equals sum of diagonal entries of M(T ), independent of basis. Deduce: tr(s + T ) = tr(s) + tr(t ). det(t ) = product of eigenvalues. Theorem: equals a unique sum formula on M(T ) using sign of permutations. Real case: Volume of T (R) = det(t ) times volume of R. Theorem: det(st ) = det(s) det(t ). Define χ T (x) := det(xi M(T )) for general F, independent of basis. Formula for A 1 using det(a); Cayley-Hamilton theorem. Trace (3) Definition 1. For a matrix A = (a ij ) over an arbitrary F, define tr(a) := a 11 + + a nn. It is immediate that tr(a + B) = tr(a) + tr(b). Theorem 2 (Corollary 10.10). For all T L(V ), tr(m(t )) does not depend on the basis. For F = C, deduce: tr(m(t )) = λ 1 + + λ n = tr(t ). For general F, define tr(t ) := tr(m(t )). Immediate: tr(s + T ) = tr(s) + tr(t ). Lemma 3 (Proposition 10.9). tr(ab) = tr(ba). 1

The theorem follows immediately from the lemma and the change-of-basis formula: tr(sas 1 ) = tr(s 1 SA) = tr(a). Proof of Lemma. tr(ab) = i,j a ijb ji = i,j b jia ij = tr(ba). Determinant and volume (4) Properties of det(t ) = λ 1 λ n : det(t ) 0 iff T is invertible. Let F = R. Theorem 10.38: vol(t (R)) = det(t ) vol(r). Suppose T is self-adjoint. Then, T dilates in orthogonal directions by λ 1,..., λ n. Taking these as axes, we deduce: for any region R, vol(t (R)) = det(t ) vol(r). For general T, write T = S T T where S is an isometry. By det theorem, det(t ) = det(s) det( T T ). Since S preserves volume, vol(t (R)) = vol( T T (R)). positive ( self-adjoint)! But T T is vol(t (R)) = det( T T ) vol(r) = det(s 1 ) det(t ) vol(r). Finally, det(s 1 ) = ±1 since S 1 is an isometry (its block diagonal matrix has ±1 and rotation matrices as blocks). So vol(t (R)) = det(t ) vol(r). Determinant and alternating forms (5) The determinant has the following essential properties: Multiadditivity: Let A and B be matrices that differ only in a fixed column. Let C be the matrix which is also the same as A and B except in this column, which is the sum of the columns appearing in A and B. Then, det(c) = det(a) + det(b). Interpretation: For R = standard cube (coords 0 and 1), then C(R), A(R), and B(R) are parallelepipeds sharing a common base; the height of C(R) is the sum of the heights of A(R) and B(R). Multihomogeneity: Let B be obtained from A by multiplying a column by λ. Then det(b) = λ det(a). Interpretation: As before, B(R) has the same base as A(R) and λ times the height. Alternation: Let A be a matrix with two identical columns. Then det(a) = 0. 2

A(R) is degenerate, so vol(a(r)) = 0. Alternatively, A is noninvertible det(a) = 0. Theorem: Our sum formula for determinant is the unique function with these properties, up to scaling. Theorem: det(ab) = det(a) det(b) (6) Assume det is the unique function with the above properties up to scaling; scale it so that det(i) = 1. Let A be a fixed matrix. Define F A (B) := det(ab). Claim: F A also has the above properties! Consequence: There exists α such that det(ab) = α det(b) for all B. Plug in B = I. Then det(a) = α det(i) = α. So, det(ab) = det(a) det(b)! Proof of claim: A is linear, so det(a(v 1 v n )) = det(av 1 Av n ) and column operations on v 1,..., v n become the same operations on (Av 1,..., Av n ). Corollary: det(ab) = det(ba). Also, det(sas 1 ) = det(s 1 SA) = det(a). So det(t ) := det(m(t )) makes sense (independent of basis)! So does χ T := det(xi M(T ))! Proof of uniqueness theorem (7) Let dim V = n and suppose that F : V n F satisfies multilinearity (=multiadditivity+multihomogeneity) and alternation. Remark: Alternation implies skew-symmetry: 0 = F (u + v, u + v) F (u, u) F (v, v) = F (u, v) + F (v, u). Let v 1,..., v n be a basis of V. Claim: F is uniquely determined by c := F (v 1,..., v n ). Suppose w j = i a i,jv i. Then F (w 1,..., w n ) = σ:{1,...,n} {1,...,n} F (a σ(1),1v σ(1),..., a σ(n),n v σ(n) ) = σ a σ(1),1 a σ(n),n F (v σ(1),..., v σ(n) ). Unless σ is a bijection, i.e., σ S n := permutations of {1,..., n}, then two of these entries are equal so we get zero. By skew-symmetry, F (v σ(1),..., v σ(n) ) = ( 1) k c where k is the number of swaps of entries needed to reorder (v σ(1),..., v σ(n) ) to (v 1,..., v n ). We get uniqueness! 3

Existence and formula for determinant (8) Set F (v σ(1),..., v σ(n) ) = sign(σ)c. We need sign(σ) = ±1 to be 1 if σ requires an even number of swaps and 1 if it requires an odd number of swaps. If such a sign exists, then sign(σ τ) = sign(σ) sign(τ) automatically. This is the essence of determinant: a multiplicative function on permutation matrices! Existence: define sign(σ) = ( 1) o(σ), where o(σ) = {(i, j) : i < j and σ(i) > σ(j)}. Lemma: If we swap adjacent entries of (σ(1),..., σ(n)), then o(σ) goes up or down by one. So it changes parity (even to odd or vice-versa). Now if we swap entries i and j, this takes 2 j i 1 adacent swaps, so o(σ) still changes parity! We deduce that ( 1) o(σ) = ( 1) k for any sequence of k swaps which equals σ. This proves existence! Matrix formula for determinant (9) Fix a basis v 1,..., v n and let det be the unique multilinear alternating function so that det(v 1,..., v n ) = 1. Given any T L(V ), define det(t ) = det(w 1,..., w n ) where T (v i ) = w i. Write A = M (vj)(t ) = (a i,j ). So w j = i a i,jv i. By our computation, det(t ) = σ S n sign(σ)a σ(1),1 a σ(n),n. Finally, for an upper triangular matrix, det(a) = the product of the diagonal entries. So for F = C, det(t ) = the product of the eigenvalues! Computing the determinant (and χ A (x)) (10) Gaussian elimination: to compute det(a) (or det(xi A)), perform row operations. Multiplying a row by λ multiplies the determinant by λ. If we swap two rows, that multiplies the determinant by 1. If we add a multiple of a row to another row, that does not change the determinant. The determinant of an upper-triangular matrix is the product of the diagonal entries! 4

Inverses of matrices using det (11) Note that det(a) 0 iff A is invertible: if it is invertible, det(a) det(a 1 ) = det(aa 1 ) = 1 so det(a) 0; if it is not invertible, then its columns are linearly dependent so det(a) = 0 as we explained. So can we get a formula for A 1 using det? More precisely, we would like A 1 = det(a) 1 A for some matrix A whose entries are polynomial functions of A. Indeed, let A = (a ij ) where a ij = ( 1)i j times the determinant of the (n 1) (n 1) matrix obtained by striking out the i-th column and the j-th row. Claim: A A = det(a)i = AA. This proves the formula! Proof: (A A) ij = k a ik a kj = det of the matrix obtained from A by replacing the i-th column with the j-th column. This is zero unless i = j, in which case we get det(a). Properties of characteristic polynomial (12) Recall χ T := χ M(T ) = det(xi M(T )). So, χ T (λ) = 0 iff λ is an eigenvalue (i.e., T λi is noninvertible). More generally: Theorem: if f(x) is an irreducible polynomial, then f is a factor of χ T iff f(t ) is noninvertible. See Optional Exercises #2 to get started. Cayley-Hamilton theorem (13) Theorem 4. χ T (T ) = 0. Proof (cf. Wikipedia!) χ A (x)i = det(xi A)I = (xi A) (xi A). Tricky part: Now write (xi A) = B 0 + xb 1 + + x n 1 B n 1 where B 0,..., B n 1 are all matrices with entries in F, not polynomials. We find that det(xi A)I = B 0 A + n 1 i=1 (B ia B i 1 )x i + B n 1 x n. Write χ A (x) = det(xi A) = x n + a n 1 x n 1 + + a 0. All the coefficients of both sides are equal: a 0 I = B 0 A, a 1 I = B 1 A B 0, etc. Now we can plug in: χ A (A) = B 0 A + n 1 i=1 (B ia B i 1 )A i + B n 1 A n. All the terms on the RHS cancel, so χ A (A) = 0. 5

Further properties of characteristic polynomial (14) Theorem: The determinant of a block upper-triangular matrix is the product of the determinants of the diagonal blocks. Proof: Any permutation σ with an entry a σ(i),i above the block diagonal has also an entry a σ(j),j below the block diagonal, which is zero. Consequence: if U V is T -invariant and V = U U, then det(t ) = det(t U ) det(p U,U T U ). Same for χ T. ( ) M(T Proof: In basis for V = U U U ), M(T ) =. 0 M(P U,U T U ) In particular, if U, U are both T -invariant, det(t ) = det(t U ) det(t U ) and tr(t ) = tr(t U ) + tr(t U ). The general decomposition theorem (15) Theorem: V = f V (f), where f ranges over irreducible polynomials such that f(t ) is noninvertible. Each V (f) is a generalized eigentuple, and is T -invariant. Proof: Again, if f(t ) is noninvertible, then V = null f dim V (T ) range f dim V (T ). Apply induction. Restricted to V (f), T has a block (lower)-triangular matrix with diagonal blocks all equal to 0 1 0 a 2, where f(x) = x k + 0 0 0 a 0 1 0 0 a 1....... 0 0 1 a k 1 a k 1 x k 1 + + a 0. Jordan form: We can make it so that all other entries are zero except for a single one in some corners between blocks. 6