Determinant Formulas and Cofactors

Determinant formulas and cofactors Now that we know the properties of the determinant, it’s time to learn some (rather messy) formulas for computing it. Formula for the determinant We know that the determinant has the following three properties: 1. det I = 1 2. Exchanging rows reverses the sign of the determinant. 3. The determinant is linear in each row separately. Last class we listed seven consequences of these properties. We can use these ten properties to find a formula for the determinant of a 2 by 2 matrix: � � � � � � � a b � � � � � � � = � a 0 �+� 0 b � � c d � � c d � � c d � � � � � � � � � � a 0 � � a 0 � � 0 b � � 0 b � � � � � � � � � = � + + + c 0 � � 0 d � � c 0 � � 0 d � = 0 + ad + (−cb) + 0 = ad − bc. By applying property 3 to separate the individual entries of each row we could get a formula for any other square matrix. However, for a 3 by 3 matrix we’ll have to add the determinants of twenty seven different matrices! Many of those determinants are zero. The non-zero pieces are: � � a11 � � a21 � � a31 a12 a22 a32 a13 a23 a33 � � � � � � = = � � � � � � � a11 0 0 �� �� a11 0 0 �� �� 0 a12 0 �� � � 0 a22 0 �� + �� 0 0 a23 �� + �� a21 0 0 �� � � 0 0 a33 � � 0 a32 0 � � 0 0 a33 � � � � � � � 0 a12 0 �� �� 0 0 a13 �� �� 0 0 a13 � 0 a23 �� + �� a21 0 0 �� + �� 0 a22 0 + �� 0 � a31 0 0 � � 0 a32 0 � � a31 0 0 a11 a22 a33 − a11 a23 a33 − a12 a21 a33 � � � � � � + a12 a23 a31 + a13 a21 a32 − a13 a22 a31 . Each of the non-zero pieces has one entry from each row in each column, as in a permutation matrix. Since the determinant of a permutation matrix is either 1 or -1, we can again use property 3 to find the determinants of each of these summands and obtain our formula. One way to remember this formula is that the positive terms are products of entries going down and to the right in our original matrix, and the negative terms are products going down and to the left. This rule of thumb doesn’t work for matrices larger than 3 by 3. The number of parts with non-zero determinants was 2 in the 2 by 2 case, 6 in the 3 by 3 case, and will be 24 = 4! in the 4 by 4 case. This is because there are n ways to choose an element from the first row (i.e. a value for α), after which there are only n − 1 ways to choose an element from the second row that avoids a zero determinant. Then there are n − 2 choices from the third row, n − 3 from the fourth, and so on. The big formula for computing the determinant of any square matrix is: det A = ∑ ± a1α a2β a3γ ...anω n! terms where (α, β, γ, ...ω ) is some permutation of (1, 2, 3, ..., n). If we test this on the identity matrix, we find that all the terms are zero except the one corresponding to the trivial permutation α = 1, β = 2, ..., ω = n. This agrees with the first property: det I = 1. It’s possible to check all the other properties as well, but we won’t do that here. Applying the method of elimination and multiplying the diagonal entries of the result (the pivots) is another good way to find the determinant of a matrix. Example In a matrix with many zero entries, many can compute the determinant of: ⎡ 0 0 1 ⎢ 0 1 1 ⎢ ⎣ 1 1 0 1 0 0 terms in the formula are zero. We ⎤ 1 0 ⎥ ⎥ 0 ⎦ 1 by choosing a non-zero entry from each row and column, multiplying those entries, giving the product the appropriate sign, then adding the results. The permutation corresponding to the diagonal running from a14 to a41 is (4, 3, 2, 1). This contributes 1 to the determinant of the matrix; the contribu­ tion is positive because it takes two row exchanges to convert the permutation (4, 3, 2, 1) to the identity (1, 2, 3, 4). Another non-zero term of ∑ ± a1α a2β a3γ a4ω comes from the permutation (3, 2, 1, 4). This contributes −1 to the sum, because one exchange (of the first and third rows) leads to the identity. These are the only two non-zero terms in the sum, so the determinant is 0. We can confirm this by noting that row 1 minus row 2 plus row 3 minus row 4 equals zero. Cofactor formula The cofactor formula rewrites the big formula for the determinant of an n by n matrix in terms of the determinants of smaller matrices. In the 3 × 3 case, the formula looks like: det A = a11 ( a22 a33 − a23 a32 ) + a12 (− a21 a33 + a23 a31 ) + a13 ( a21 a32 − a22 a31 ) � � � � � � � a11 0 0 �� �� 0 a12 0 �� �� 0 0 a13 �� � 0 a23 �� + �� a21 a22 0 �� = �� 0 a22 a23 �� + �� a21 � 0 a32 a33 � � a31 � � 0 a33 a31 a32 0 � This comes from grouping all the multiples of aij in the big formula. Each element is multiplied by the cofactors in the parentheses following it. Note that each cofactor is (plus or minus) the determinant of a two by two matrix. That determinant is made up of products of elements in the rows and columns NOT containing a1j . In general, the cofactor Cij of aij can be found by looking at all the terms in the big formula that contain aij . Cij equals (−1)i+ j times the determinant of the n − 1 by n − 1 square matrix obtained by removing row i and column j. (Cij is positive if i + j is even and negative if i + j is odd.) For n × n matrices, the cofactor formula is: det A = a11 C11 + a12 C12 + · · · + a1n C1n . Applying this to a 2 × 2 matrix gives us: � � � a b � � � � c d � = ad + b(−c). Tridiagonal matrix A tridiagonal matrix is one for which the only non-zero entries lie on or adjacent to the diagonal. For example, the 4 × 4 tridiagonal matrix of 1’s is: ⎡ ⎤ 1 1 0 0 ⎢ 1 1 1 0 ⎥ ⎥ A4 = ⎢ ⎣ 0 1 1 1 ⎦. 0 0 1 1 What is the determinant of an n × n tridiagonal matrix of 1’s? � � � � � 1 1 0 � � 1 1 � � � � = 0, | A3 | = � 1 1 1 � = −1 | A1 | = 1, | A2 | = �� � � 1 1 � � 0 1 1 � � � � � � 1 1 0 � � 1 1 0 � � � � � | A4 | = 1 �� 1 1 1 �� − 1 �� 0 1 1 �� = | A3 | − 1| A2 | = −1 � 0 1 1 � � 0 1 1 � In fact, | An | = | An−1 | − | An−2 |. We get a sequence which repeats every six terms: | A1 | = 1, | A2 | = 0, | A3 | = −1, | A4 | = −1, | A5 | = 0, | A6 | = 1, | A7 | = 1.