1 / 65

ENGG2012B Lecture 8 Determinant and Cramer’s rule

ENGG2012B Lecture 8 Determinant and Cramer’s rule. Kenneth Shum. Last lecture: Invertible matrix. Given an n n matrix A , if we can find a matrix A ’, such that then A is said to be invertible , or non-singular . We can prove that there is at most one such matrix, if exists.

rory
Download Presentation

ENGG2012B Lecture 8 Determinant and Cramer’s rule

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. ENGG2012BLecture 8Determinant and Cramer’s rule Kenneth Shum ENGG2012B

  2. Last lecture: Invertible matrix • Given an nn matrix A, if we can find a matrix A’, such that then A is said to be invertible, or non-singular. • We can prove that there is at most one such matrix, if exists. • The matrix A’ is called the inverse of A. Multiplied byA Multiplied by A’ Multiplied byIn ENGG2012B

  3. Algorithm for computing inverse • Input: an nn matrix A. • Create an n  2n matrix M = [ A | In ] • The left half is A • The right half is the n by n indentity matrix In • Try to reduce the expanded matrix Msuch that the left half is equal to In. • If succeed, the right half of M is the inverse of A. • If you cannot reduce the left half of M to In, then A is not invertible. ENGG2012B

  4. Elementary matrix • Elementary row operations on a matrix A can be realized by multiplying A by a square matrix, called elementary matrix. • Every matrix A can be reduced to RREF by multiplying from the left a number of elementary matrices. • E1, E2, … are elementary matrices. • R is a matrix in RREF. • We can undo any elementary row operation. The reverse operation is also an elementary row operation. • In matrix term, the inverse of an elementary matrix is also an elementary matrix. • F1, F2,… are inverse of E1, E2…, respectively. ENGG2012B

  5. Vector Notation • We will use two different notations for a point in the 3D space (x,y,z) z z y y x x ENGG2012B

  6. Today’s aim • 2x2 and 3x3 determinants • Cramer’s rule • Extension to nxn determinant • Cramer’s rule • Adjoint formula for matrix inverse ENGG2012B

  7. 2X2 AND 3X3 DETERMINANT ENGG2012B

  8. Area of parallelogram • The sides are vector (a,b) and (c,d). a c b (c,d) d d (a,b) b c a ENGG2012B

  9. 22 determinant as signed area • Definition: Example + – (2,4) Area = 10 Row 1 is the first vector Row 2 is the second vector - bc ad (3,1) ENGG2012B

  10. Notations for 22 determinant Note 1: By viewing a a real number x as a 1x1 matrix, we have det(x) = x . Note 2: For real number x, |x| usually means theabsolute value of x, not the determinant. ENGG2012B

  11. Example Determinant of the 2x2 identity matrix is 1 1 1 ENGG2012B

  12. Example Determinant of a diagonal matrix is the area of a rectangle h w ENGG2012B

  13. Zero column/row  zero determinant Zero area ENGG2012B

  14. Linearly dependent columns/rows  zero determinant If the two columns (or two rows) are constant multiple of each other, the determinant is zero. Zero area ENGG2012B

  15. Geometric meaning of 3x3 determinant • The magnitude of 33 determinant is the volume of a parallelepiped. z y x ENGG2012B

  16. 33 determinant Notation for 33 determinant : Definition: ENGG2012B

  17. Definition of 3x3 determinant – – + + + – ENGG2012B

  18. Determinant of Diagonal matrix • Volume of a rectangular box c b a ENGG2012B

  19. Co-planar  zero determinant • Determinant = 0 Volume = 0  the three vectors lie on the same plane z y A collection of vectorsare said to be co-planarif they lie on the same plane. x ENGG2012B

  20. Application to computing inverse of 2x2 matrix • Given 22 matrix • A formula for A-1: If det A is nonzero, we have • Exchange the two diagonal entries a, d. • Take the negative of the two off-diagonal entries b, c. • Divide by the determinant. ENGG2012B

  21. Solving system of equations by matrix inverse • Want to solve a system of linear equation in matrix form A x= b. • If the inverse of A exists, we can multiply both sides of A x= b by the inverse of A, and get Start with A x= b Multiply by A-1 from the left Matrix product is associative The defining property of inverse I times anything is the same thing ENGG2012B

  22. Example • Solve • If we know the inverse of the 2x2 matrix, we can solve the linear system easily. ENGG2012B

  23. PROPERTIES OF 2X2 AND 3X3 DETERMINANT ENGG2012B

  24. Transpose matrix has the same determinant “The transpose of a matrix” means reflecting the matrixalong the diagonal. Row 1 and row 2 becomecolumn 1 and column 2, andvice versa. We write AT for the transpose of matrix A. – – + + + – Compare with ENGG2012B

  25. Meta rule Any row property of determinant is a column property, and vice versa ENGG2012B

  26. Example of transposing 2x2 matrix (1,4) Area = 10 (2,4) Area = 10 (3,2) (3,1) ENGG2012B

  27. Example of transposing 3x3 matrix • In computing the volume of a parallelepiped, it does not matter whether we write the vector horizontally or vertically in the determinant z Volume of parallelepiped with vertices(0,0,0), (1,2,0), (2,0,1), (–1, 1, 3) equals tothe absolute value of y or x ENGG2012B

  28. Exchanging two rows/columns  sign of determinant is flipped The first kindof elementaryrow operation Example (1,3) (2,1) ENGG2012B

  29. Multiply a row/column by a constant If we multiply a row or a column by a constant c, the value of determinant also increase by a factor of c (Here, c may be zero) (4,4) (0,1) This corresponds to the 2nd kind of elementaryrow operation if c is nonzero (1,1) ENGG2012B

  30. Scale up/down all entries by a constant c • Given an n x n matrix A, if we multiply it by a constant c, the determinant is increased by a factor of cn, i.e., det(cA) = cn det(A). For example, if we multiply the sides of a box by a factor of 3, then the volume increases by a factor of 27. ENGG2012B

  31. Adding a constant multiple of a row/column to another row/column If we add a constant multiple of a row (column) to the other row (column), the determinant does not change. The second row is obtained byadding three times [1 0] to [0 1]. (0,1) (3,1) The 3rd kindof elementaryrow operation (1,0) The rectangle and parallelogramhave the same height and same base. ENGG2012B

  32. Summary ENGG2012B

  33. Additive property If a row (or column) of a determinant is the sum of two rows (or columns), the determinant can be split as the sum of two determinants For example, if the first column is the sum of two column vectors, thenwe have ENGG2012B

  34. Repeated columns/rows  zero determinant If any two columns (or rows) are the identical, then the determinant is zero. Examples: ENGG2012B

  35. Cramer’s rule • If the determinant of a 33 matrix A is non-zero, we can solve the linear system Ax = b by Cramer’s rule. • The solution to is or equivalently A x b Gabriel Cramer (1704-1752) ENGG2012B

  36. Verification of the Cramer’s rule • Claim: The solution to is provided that the determinant of A is nonzero. We verify the formula for x1 in the next slide. Cramer’s rule in wikipedia ENGG2012B

  37. Because x1, x2, x3 satisfy the system of linear equations, we have Then, By substitution By the additive property We can bring outa common factorin a column. det=0 if there arerepeated columns =0 =0 ENGG2012B

  38. Proof continued Hence, There is no division by zero, because det A is non-zero by hypothesis. The verification of the formula for x2 and x3 can be done similarly. QED ENGG2012B

  39. Det is a multi-linear function http://en.wikipedia.org/wiki/Multilinear_map • Combined with the second type of elementary row/column operation, The additive property can be strengthened to the multi-linear property. • If the i-th column of a determinant can be decomposed as then we can write the determinant as the sum of two determinants with all other entries not in the i-th column kept unchanged. i-th column i-th column i-th column ENGG2012B

  40. Useful fact • The determinant of an upper/lower triangular matrix is the product of the diagonal entries. Examples ENGG2012B

  41. GENERALIZATION TO NXN DETERMINANT ENGG2012B

  42. We want a function which … Set of all nxn matrices produces a real number for each n x n matrix, satisfying the analogous properties as in the 2x2 and 3x3 determinants, such as … ? the real number line 4 2 3 -1 1 0 ENGG2012B

  43. Axiom 1 (alternating property) • If any two rows are exchanged, then the determinant is multiplied by -1. ENGG2012B

  44. Axiom 2 (multiplicative property) • A common factor in a row can be brought out of the determinant. ENGG2012B

  45. Axiom 3 (additive property) • If the i-th row is the sum of row vectors ri and si, then the determinant can be decomposed into the sum of two determinants. ENGG2012B

  46. Axiom 4 (identity matrix) • The determinant of the identity matrix is 1. ENGG2012B

  47. A weaker form of Axiom 1 • Axiom 1’: If there are two identical rows, then the determinant is zero. ENGG2012B

  48. A recursive definition of determinant • Expansion on the first row. ENGG2012B

  49. Verification that the recursive definition satisfies the axioms • We can verify that the recursive definition for general nxn matrices satisfies the four axioms. • The verification can be carried out by induction. • Show that the determinant function for 2x2 matrices satisfies the four axioms. • Suppose that the axioms hold for (n-1)x(n-1), (n-2)x(n-2), … Show that they are also true if we define nxn determinant by expansion on the first row. • We illustrate the verification process for n=4. ENGG2012B

  50. Verification of Axiom 4 (identity) True by definition True by induction hypothesis ENGG2012B

More Related