[Linear Algebra 05] Properties and Applications of Determinants

   Continuing the notes of MIT, this article sorts out the 10 properties of the determinant, of which the first 3 properties serve as the basis for the subsequent properties. Applications are grouped into Cramer's Law, Inverse Matrix and Volume. Note that only square matrices have determinants.


nature

  1. det ⁡ I = 1 \det I =1theI=1
  2. Swap two rows of the determinant, change the sign of the determinant
    Note: exchange once to change number once, exchange iii times, multiply by( − 1 ) i (-1)^i(1)i , so the even-numbered transformation and the odd-numbered transformation are in their own independent transformation categories.
  3. Linearity of "rows":a) Homogeneity: the coefficient of a row of the determinant is proposed as the coefficient of the whole determinant; b) Additivity: when a row of the determinant can be split into two parts, it is equal to the two determinants after splitting of and.Use 2 × 2 2 \times 22×2行列式表示即为:
    ∣ 2 a + 3 e 2 b + 3 f c d ∣ = ∣ 2 a 2 b c d ∣ + ∣ 3 e 3 f c d ∣ = 2 ∣ a b c d ∣ + 3 ∣ e f c d ∣ \begin{vmatrix} 2a + 3e&2b +3f\\ c & d \end{vmatrix}= \begin{vmatrix} 2a &2b \\ c & d \end{vmatrix}+ \begin{vmatrix} 3e&3f\\ c & d \end{vmatrix}= 2\begin{vmatrix} a &b \\ c & d \end{vmatrix}+ 3\begin{vmatrix} e&f\\ c & d \end{vmatrix} 2a _+3 ec2 b+3f _d = 2a _c2 bd + 3 ec3f _d =2 acbd +3 ecfd
    Note: The linearity of the determinant is for the "row", and the linearity of the matrix is ​​for the entire matrix. Therefore ∣ 2 A ∣ = 2 n ∣ A ∣ |2A|=2^n|A|∣2A=2nA

  1. When two rows are equal, the determinant is 0. It can be deduced from property 2.
    A = ∣ abab ∣ swap two lines → − A = ∣ abab ∣ ⇒ A = 0 A = \begin{vmatrix} a &b \\ a & b \end{vmatrix} \underrightarrow{swap two lines} -A =\begin {vmatrix} a &b \\ a & b \end{vmatrix} \Rightarrow A=0A= aabb swap two linesA= aabb A=0
    Note: Further by nature 3, when there are two rows in a multiple relationship, the determinant is also 0.

  2. When one row subtracts k times another row, the determinant does not change.
    Note: The proof similar to property 4 is based on the reasonable extrapolation of property 3. The summary of the determinant is 0 is:When there is a linear combination of other rows in a certain row, the determinant is 0, which is equivalent to requiring the corresponding matrix to be full rank

  3. If one row is 0, the determinant is 0.
    Note: 0 row plus any row is equal to two rows of the original determinant. It can be deduced based on properties 4 and 5.

  4. After the elimination method becomes a triangular matrix, such as A = LUA=LUA=LU transform,Triangular array UUThe product of elements on the main column of U (that is, on the main diagonal) is the value of the determinant.
    Note: The elimination method does not change the value of the determinant.

  5. When the matrix is ​​singular (that is, the rank is not satisfied), the corresponding determinant is 0 .
    Note: It is equivalent to the summary in Property 5.

  6. The determinant of the product of matrices is equal to the product of the determinants , that isdet ⁡ AB = det ⁡ A det ⁡ B \det{AB}=\det{A}\det{B}theAB=theAtheB
    Note: Since all matrices can be transformed into diagonal matrices, if the product of diagonal matrices is considered, obviously the result is valid. Further we can know: det ⁡ A − 1 = 1 / det ⁡ A \det{A^{-1}} = 1/\det{A}theA1=1/theA

  7. The determinant is equal to the transposed determinant. Give a sufficient and concise proof:
    ∣ AT ∣ = ∣ A ∣ ⇒ ∣ UTLT ∣ = ∣ LU ∣ ⇒ ∣ UT ∣ ∣ LT ∣ = ∣ L ∣ ∣ U ∣ (property 9) √ |A^T| = | A| \Rightarrow |U^TL^T| = |LU| \Rightarrow |U^T||L^T| = |L||U| (Property 9) √AT=AUTLT=LUUT∣∣LT=L ∣∣ U (property 9 ) √Note
    : This shows thatThe properties of the rows can be reasonably extrapolated to the columns, for example, when a column is a linear combination of other columns, the determinant will also be 0, and if two columns are exchanged, the sign of the determinant will also change.


  The determinant is equal to the product of the pivots , which implies the importance of these pivots (i.e. eigenvalues).


Expansion of the determinant

  Remember the following formula, namelyThe determinant can be expanded by taking a certain row
det ⁡ A = a i 1 C i 1 + a i 2 C i 2 + ⋯ + a i n C i n \det A = a_{i1}C_{i1}+a_{i2}C_{i2}+\cdots+a_{in}C_{in} theA=ai 1Ci 1+ai2Ci2++ainCin
  Among them, C is called the algebraic cofactor , which is the cofactor after taking the transformed sign . That is,
C ij = ( − 1 ) i + j det ⁡ M ij C_{ij} = (-1)^{i+j}\det M_{ij}Cij=(1)i+jtheMij
  The remainder formula is to delineate aij a_{ij}aijAfter the elements in the row and column, the remaining elements are arranged into a matrix.
∣ + aij − . . . ∣ . . . M ij ∣ \begin{vmatrix} \begin{array}{c:cc} \ +_{a_{ij}}& -_{...} \\ \hdashline |_{...}& M_{ij} \end{array} \end{vmatrix}  +aij......Mij
  Take a 3 × 3 3 \times 33×3. Take the calculation of the determinant as an example. Use the det command in Matlab to find the value of the matrix determinant:
∣ 1 2 3 4 5 6 7 8 9 ∣ = 1 ( − 1 ) 2 ∣ 5 6 8 9 ∣ + 2 ( − 1 ) 3 ∣ 4 6 7 9 ∣ + 3 ( − 1 ) 4 ∣ 4 5 7 8 ∣ = − 3 + 12 − 9 = 0 \begin{vmatrix} 1 & 2 & 3\\ 4& 5 & 6 \\ 7 &8 &9 \\ \end{vmatrix} = 1(-1)^2 \begin{vmatrix} 5 & 6 \\ 8&9 \\ \end{vmatrix} +2(-1)^3\begin{vmatrix} 4 & 6 \\ 7 &9 \\ \end{vmatrix} +3(-1)^4\begin{vmatrix} 4 & 5 \\ 7 & 8 \\ \end{vmatrix} =-3+12-9=0 147258369 =1(1)2 5869 +2(1)3 4769 +3(1)4 4758 =3+129=0

>> A = [1,2,3;4,5,6;7,8,9];
>> det(A) % 行列式

ans =

  -9.5162e-16

>> rank(A) %

ans =

     2

>> rref(A) % 最简行形式

ans =

     1     0    -1
     0     1     2
     0     0     0

  The key to prove the idea is to take a certain line and expand it one by one, it is easy to find that when taking aij a_{ij}aijelements, the elements of the rest of the rows cannot be selected from the rows and columns that have been crossed out, thus becoming a permutation and combination problem , that is, for n × nn \times nn×In terms of the determinant of n , there are a total ofn ! n!n ! sub-items, arranged in a determinant, becomesn × n − 1 × ⋯ 1 n\times n-1 \times \cdots 1n×n1×The level-by-level decomposition of 1 is the level-by-level expansion of the determinant.


application

Inverses (inverse matrix)

  Explicit representation of the inverse matrix:The inverse matrix is ​​the transpose of the algebraic cofactor matrix divided by the determinant of the original matrix。 即
A − 1 = CT / det ⁡ AA^{-1} = C^{T}/\det AA1=CT/theThe derivation of the inverse matrix of A
  starts with the expansion of the determinant, that is to say,ACT = ( det ⁡ A ) I AC^{T} = (\det A) IACT=( itA ) I ratio:
[ a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ] [ C 11 C 21 C 31 C 12 C 22 C 32 C 13 C 23 C 33 ] = [ det ⁡ A 0 0 0 it ⁡ A 0 0 0 it ⁡ A ] \begin{bmatrix} a_{11} & a_{12} & a_{13}\\ a_{21}&a_{22} &a_{23} \\ a_{ 31} & a_{32} &a_{33} \end{bmatrix}\begin{bmatrix}C_{11}&C_{21}&C_{31}\\C_{12}&C_{22}&C_{32}\\ C_{13} & C_{23} &C_{33} \end{bmatrix} = \begin{bmatrix} \the A & 0 & 0\\ 0&\the A &0 \\ 0& 0 &\the A \end{bmatrix } } a11a21a31a12a22a32a13a23a33 C11C12C13C21C22C23C31C32C33 = theA000theA000theA
  Focusing on the elements of the diagonal line, it is not difficult to find that it is the expression of the determinant expansion, and how to explain the remaining zeros? This is a great question. Let's look at an example:
? a 21 C 11 + a 22 C 12 + a 23 C 13 = 0 ? a_{21}C_{11}+a_{22}C_{12}+a_{23}C_{13}=0a21C11+a22C12+a23C13=0
  First quote the explanation in the textbook:

This is the cofactor rule for a new matrix, when the second row of A A A is copied into its first row. The new matrix A ∗ A^* A has two equal rows, so det A ∗ A^* A = 0. Notice that A ∗ A^* A has the same cofactors C 11 C_{11} C11, C 12 C_{12} C12, C 13 C_{13} C13 as A A A-because all rows agree after the first row. Thus the remarkable multiplication is correct.

  In other words, by construction, the algebraic remainders corresponding to the first row of the following two matrices are the same :
A = [ a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ] and A ∗ = [ a 21 a 22 a 23 a 21 a 22 a 23 a 31 a 32 a 33 ] ⇒ det ⁡ A ∗ = a 21 C 11 + a 22 C 12 + a 23 C 13 = 0 A = \begin{bmatrix} a_{11} & a_{12} & a_{13}\\ a_{21}&a_{22} &a_{23} \\ a_{31} & a_{32} &a_{33} \end{bmatrix} \ \ \ and \ \ \ A^* = \begin{bmatrix} a_{21}&a_{22} &a_{23}\\ a_{21}&a_{22} &a_{23} \\ a_{31} & a_{32} &a_{33 } \end{bmatrix} \Rightarrow \det{A^*} =a_{21}C_{11}+a_{22}C_{12}+a_{23}C_{13}=0A= a11a21a31a12a22a32a13a23a33    with A   = a21a21a31a22a22a32a23a23a33 theA=a21C11+a22C12+a23C13=0


Cramer's Rule

  For the equation A x = b Ax=bAx=As far as b is concerned, Cramer's law is a "fancy but not useful" solution formula (not as good as the elimination method). Its main significance lies inexplicitly giving the expression of the solution of the equation, although it will be very cumbersome to do. WhenAAThe determinant of A is not 0 (or when reversible),A x = b Ax = bAx=b字性的解的
x 1 = det ⁡ B 1 det ⁡ A x 2 = det ⁡ B 2 det ⁡ A ⋯ xn = det ⁡ B n det ⁡ A x_1= \frac{\det B_1} {\det A} \ \ x_2= \frac{\det B_2} {\det A} \ \ \cdots \ x_n= \frac{\det B_n} {\det A} \x1=theAtheB1  x2=theAtheB2   xn=theAtheBn 
  Among them B i B_iBilisted as AAA secondiibbfor iThe result after replacement in column b .


Idea 1

  There are two ways of thinking to explain this result. The first is to directly substitute into the expression of the inverse matrix, that is,
x = A − 1 b = CT b det ⁡ A x = A^{-1}b=\frac{C^{T } b}{\det A}x=A1b=theACTb
  take x 1 x_1x1For example, the corresponding CT b C^T bCT bis
b 1 C 11 + b 2 C 21 + b 3 C 31 = 0 b_{1}C_{11}+b_{2}C_{21}+b_{3}C_{31}=0b1C11+b2C21+b3C31=0
  Since the value of the determinant remains unchanged after the transposition, it can be expanded by taking a certain row, and of course it can also be expanded by taking a certain column. And when a certain column is taken to expand, we know that thealgebraic cofactors corresponding to the first column of the following two matrices are the same:
A = [ a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ] And B 1 = A ∗ = [ b 1 a 12 a 13 b 2 a 22 a 23 b 3 a 32 a 33 ] ⇒ det ⁡ A ∗ = b 1 C 11 + b 2 C 21 + b 3 C 31 A = \ begin{bmatrix} a_{11} & a_{12} & a_{13}\\ a_{21}&a_{22} &a_{23} \\ a_{31} & a_{32} &a_{33} \end{ bmatrix} \ \ \ and \ \ \ B_1=A^* = \begin{bmatrix} b_{1} & a_{12} & a_{13}\\ b_{2}&a_{22} &a_{23} \\ b_{3} & a_{32} &a_{33} \end{bmatrix} \Rightarrow \det{A^*} =b_{1}C_{11}+b_{2}C_{21}+b_{3} C_{31}A= a11a21a31a12a22a32a13a23a33    with B   1=A= b1b2b3a12a22a32a13a23a33 theA=b1C11+b2C21+b3C31
  Therefore, for CT b C^TbCCorresponding to xi x_iin T bxiIn terms of the row, it can be used respectively by B i B_iBiSubstitute and get the same result.

Idea 2

  This idea is more ingenious. With the help of the identity matrix, we do the following transformation:
[ a 1 a 2 a 3 ] [ x 1 0 0 x 2 1 0 x 3 0 1 ] = [ ba 2 a 3 ] = B 1 \begin {bmatrix} \\ a_{1} & a_{2} & a_{3} \\ \\ \end{bmatrix} \begin{bmatrix} x_1& 0 &0 \\ x_{2} & 1 & 0\\ x_3 &0 & 1 \\ \end{bmatrix} = \begin{bmatrix} \\ b & a_{2} & a_{3} \\ \\ \end{bmatrix} =B_1 a1a2a3 x1x2x3010001 = ba2a3 =B1
  Taking the determinant on both sides, with the help of properties 7 and 9 of the determinant, there is:
( det ⁡ A ) x 1 = det ⁡ B 1 (\det {A}) x_1 =\det B_1( itA)x1=theB1


Volumes

  We have come into contact with such formulas when learning the cross product of vectors, such as calculating vector a ⃗ = ( 1 , 2 , 3 ) \vec{a}=(1,2,3)a =(1,2,3 ) b ⃗ =( 4 , 5 , 6 ) \vec{b}=(4,5,6)b =(4,5,6 )的叉乘:
a ⃗ × b ⃗ = ∣ 1 2 3 4 5 6 i ⃗ j ⃗ k ⃗ ∣ = ( 12 − 15 ) i ⃗ + ( 12 − 6 ) i ⃗ + ( 5 − 8 ) k ⃗ = ( − 3 , 6 , − 3 ) \vec{a} \times \vec{b} = \begin{vmatrix} 1& 2 &3 \\ 4 & 5 & 6\\ \vec i &\vec j & \vec k \\ \end{vmatrix} =(12-15)\vec{i}+(12-6)\vec{i}+(5-8)\vec{k}=(-3,6,-3 )a ×b = 14i 25j 36k =(1215)i +(126)i +(58)k =(3,6,3)

>> a=[1,2,3];b=[4,5,6];
>> cross(a,b)

ans =

    -3     6    -3

  Further, we find the above result and c ⃗ = ( 0 , 1 , 0 ) \vec{c}=(0,1,0)c =(0,1,0 ) , that is, the triple product, which representsthe volume of the parallelepiped formed by these three vectors:
a ⃗ × b ⃗ ⋅ c ⃗ = ∣ 1 2 3 4 5 6 0 1 0 ∣ = ( 12 − 15 ) + ( 12 − 6 ) 1 + ( 5 − 8 ) 0 = 6 \vec{a} \times \vec{b} \ \cdot \vec{c} = \begin{vmatrix} 1& 2 &3 \\ 4 & 5 & 6\\ 0 &1 & 0 \\ \end{vmatrix} =(12-15)+(12-6)1+(5-8)0=6a ×b  c = 140251360 =(1215)+(126)1+(58)0=6

>> c=[0,1,0];
>> dot(c,cross(a,b))

ans =

     6

  Based on the understanding of the triple product, we can derive the formula for the area of ​​a triangle on the plane, that is, it can be regarded as half the volume of a parallelepiped with a height of 1, so, there is:
1 2 ∣ x 1 y 1 1 x 2 y 2 1 x 3 y 3 1 ∣ = 1 2 [ ( x 2 y 3 − x 3 y 2 ) + ( x 3 y 1 − x 1 y 3 ) + ( x 1 y 2 − x 2 y 1 ) ] = 1 2 [ x 1 ( y 2 − y 3 ) + x 2 ( y 3 − y 1 ) + x 3 ( y 1 − y 2 ) ] \frac{1}{2} \begin{vmatrix} x_1& y_1 &1 \\ x_2 & y_2 & 1\\ x_3 &y_3 & 1 \\ \end{vmatrix} =\frac{1}{2} [(x_2y_3-x_3y_2)+(x_3y_1-x_1y_3)+(x_1y_2-x_2y_1)]=\frac{1} {2}[x_1(y_2-y_3)+x_2(y_3-y_1)+x_3(y_1-y_2)]21 x1x2x3y1y2y3111 =21[(x2y3x3y2)+(x3y1x1y3)+(x1y2x2y1)]=21[x1(y2y3)+x2(y3y1)+x3(y1y2)]
  In particular, when( x 3 , y 3 ) = ( 0 , 0 ) (x_3,y_3)=(0,0)(x3,y3)=(0,0 ) , the area formula can be further simplified as:
1 2 ∣ x 1 y 1 x 2 y 2 ∣ = 1 2 ( x 1 y 2 − x 2 y 1 ) \frac{1}{2} \begin {vmatrix} x_1& y_1 \\ x_2 & y_2 \end{vmatrix} = \frac{1}{2} (x_1y_2-x_2y_1)21 x1x2y1y2 =21(x1y2x2y1)
  It is worth noting that the orthogonal matrixQQQ , sinceQQT = I QQ^T=IQQT=I , so when taking the determinant on both sides, there is( det ⁡ Q ) 2 = 1 (\det Q) ^2=1( itQ)2=1 , which means thatthe rotation matrix Q only changes the angle, but does not change the volume. This knowledge helps to enhance the understanding of operator operations in electromagnetism and phase rotation in signals and systems. In the next lecture on eigenvalues, we will more specifically discover the rich meanings brought by matrix operations.

Guess you like

Origin blog.csdn.net/weixin_47305073/article/details/126207646