Cramer's Rule, Volume, and Linear Transformations

Cramer’s Rule

  • Theorem 7 (Cramer’s Rule): Let AA be an invertible n×nn \times n matrix. For any b\mathbf{b} in Rn\mathbb{R}^n, the unique solution x\mathbf{x} of Ax=bA\mathbf{x} = \mathbf{b} has entries given by:

    • xi=det(Ai(b))det(A)x_i = \frac{\det(A_i(\mathbf{b}))}{\det(A)}, for i=1,2,,ni = 1, 2, \dots, n
    • Where Ai(b)A_i(\mathbf{b}) is the matrix formed by replacing column ii of AA with the vector b\mathbf{b}.
  • Proof Logic for Cramer's Rule:

    • Denote the columns of AA as a1,,an\mathbf{a}_1, \dots, \mathbf{a}_n and the columns of the n×nn \times n identity matrix II as e1,,en\mathbf{e}_1, \dots, \mathbf{e}_n.
    • If Ax=bA\mathbf{x} = \mathbf{b}, then by the definition of matrix multiplication:
      • AIi(x)=A[e1xen]=[Ae1AxAen]=[a1ban]=Ai(b)A \cdot I_i(\mathbf{x}) = A[\mathbf{e}_1 \dots \mathbf{x} \dots \mathbf{e}_n] = [A\mathbf{e}_1 \dots A\mathbf{x} \dots A\mathbf{e}_n] = [\mathbf{a}_1 \dots \mathbf{b} \dots \mathbf{a}_n] = A_i(\mathbf{b})
    • By the multiplicative property of determinants:
      • det(A)det(Ii(x))=det(Ai(b))\det(A) \cdot \det(I_i(\mathbf{x})) = \det(A_i(\mathbf{b}))
    • The determinant of Ii(x)I_i(\mathbf{x}) is simply xix_i (by cofactor expansion along row ii).
    • Hence, det(A)xi=det(Ai(b))\det(A) \cdot x_i = \det(A_i(\mathbf{b})).
    • Since AA is invertible, det(A)0\det(A) \neq 0, allowing us to solve for xix_i.
  • Example 1: Solving a System via Cramer’s Rule:

    • Consider the system:
      • 3x12x2=63x_1 - 2x_2 = 6
      • 5x1+4x2=8-5x_1 + 4x_2 = 8
    • Matrix form Ax=bA\mathbf{x} = \mathbf{b}:
      • A=(3254)A = \begin{pmatrix} 3 & -2 \\ -5 & 4 \end{pmatrix}, b=(68)\mathbf{b} = \begin{pmatrix} 6 \\ 8 \end{pmatrix}
    • Calculations:
      • det(A)=(3×4)(2×5)=1210=2\det(A) = (3 \times 4) - (-2 \times -5) = 12 - 10 = 2
      • A1(b)=(6284)    det(A1(b))=(6×4)(2×8)=24+16=40A_1(\mathbf{b}) = \begin{pmatrix} 6 & -2 \\ 8 & 4 \end{pmatrix} \implies \det(A_1(\mathbf{b})) = (6 \times 4) - (-2 \times 8) = 24 + 16 = 40
      • A2(b)=(3658)    det(A2(b))=(3×8)(6×5)=24+30=54A_2(\mathbf{b}) = \begin{pmatrix} 3 & 6 \\ -5 & 8 \end{pmatrix} \implies \det(A_2(\mathbf{b})) = (3 \times 8) - (6 \times -5) = 24 + 30 = 54
    • Solutions:
      • x1=402=20x_1 = \frac{40}{2} = 20
      • x2=542=27x_2 = \frac{54}{2} = 27

A Formula for the Inverse Matrix

  • The Adjugate Matrix: The adjugate matrix (also called the classical adjoint) is the transpose of the matrix of cofactors formed from matrix AA.

  • Theorem 8: Let AA be an invertible n×nn \times n matrix. Then:

    • A1=1det(A)adj(A)A^{-1} = \frac{1}{\det(A)} \text{adj}(A)
  • Example 3: Finding the Inverse of a 3x3 Matrix:

    • Matrix A=(213111142)A = \begin{pmatrix} 2 & 1 & 3 \\ 1 & -1 & 1 \\ 1 & 4 & -2 \end{pmatrix}
    • The nine cofactors (CijC_{ij}) are calculated as follows:
      • C11=2C_{11} = -2, C12=3C_{12} = 3, C13=5C_{13} = 5
      • C21=14C_{21} = 14, C22=7C_{22} = -7, C23=7C_{23} = -7
      • C31=4C_{31} = 4, C32=1C_{32} = 1, C33=3C_{33} = -3
    • The adjugate matrix is the transpose of the cofactor matrix:
      • adj(A)=(2144371573)\text{adj}(A) = \begin{pmatrix} -2 & 14 & 4 \\ 3 & -7 & 1 \\ 5 & -7 & -3 \end{pmatrix}
    • To find det(A)\det(A), we can check the product (adj(A))A(\text{adj}(A))A:
      • (adj(A))A=(140001400014)=14I(\text{adj}(A))A = \begin{pmatrix} 14 & 0 & 0 \\ 0 & 14 & 0 \\ 0 & 0 & 14 \end{pmatrix} = 14I
    • Since (adj(A))A=det(A)I(\text{adj}(A))A = \det(A)I, then det(A)=14\det(A) = 14.
    • The inverse is:
      • A1=114(2144371573)A^{-1} = \frac{1}{14} \begin{pmatrix} -2 & 14 & 4 \\ 3 & -7 & 1 \\ 5 & -7 & -3 \end{pmatrix}

Determinants as Area or Volume

  • Theorem 9:

    • If AA is a 2×22 \times 2 matrix, the area of the parallelogram determined by the columns of AA is det(A)|\det(A)|.
    • If AA is a 3×33 \times 3 matrix, the volume of the parallelepiped determined by the columns of AA is det(A)|\det(A)|.\n
  • Geometric Observations for the Proof:

    • The theorem is clearly true for any diagonal matrix. For example, in the 2×22 \times 2 case with a diagonal matrix having entries aa and dd, the area is ad|ad|.
    • Column operations: It is shown that column replacement does not change the area or volume.
    • Specifically, if a1,a2\mathbf{a}_1, \mathbf{a}_2 are nonzero vectors, for any scalar cc, the area of the parallelogram determined by a1\mathbf{a}_1 and a2\mathbf{a}_2 equals the area of the parallelogram determined by a1\mathbf{a}_1 and a2+ca1\mathbf{a}_2 + c\mathbf{a}_1.
    • Proof reasoning: This is because both parallelograms share the same base (from 00 to a1\mathbf{a}_1) and the same altitude (since a2+ca1\mathbf{a}_2 + c\mathbf{a}_1 lies on a line parallel to the base).
    • Any matrix AA can be transformed into a diagonal matrix using column operations that do not change its determinant (up to a sign) or the volume of its associated shape.
  • Definition of Volume for Parallelepipeds (n=3n=3):

    • The volume is the area of the base in the plane Span{a1,a2}\text{Span}\{\mathbf{a}_1, \mathbf{a}_2\} multiplied by the altitude of a3\mathbf{a}_3 above that plane.
    • A column replacement operation where a3\mathbf{a}_3 is changed to a3+ca1\mathbf{a}_3 + c\mathbf{a}_1 does not change the altitude because the new vector still has the same perpendicular distance to the base plane.
  • Example 4: Calculating Area for a Parallelogram not at the Origin:

    • Vertices of the parallelogram: (1,1),(3,6),(7,2),(9,7)(1, 1), (3, 6), (7, 2), (9, 7).
    • To use the determinant method, translate the parallelogram to the origin by subtracting one vertex from all others.
    • Subtract $(1, 1)$ from each vertex:
      • (1,1)(1,1)=(0,0)(1, 1) - (1, 1) = (0, 0)
      • (3,6)(1,1)=(2,5)(3, 6) - (1, 1) = (2, 5)
      • (7,2)(1,1)=(6,1)(7, 2) - (1, 1) = (6, 1)
      • (9,7)(1,1)=(8,6)(9, 7) - (1, 1) = (8, 6)
    • The parallelogram is now determined by the vectors v1=(25)\mathbf{v}_1 = \begin{pmatrix} 2 \\ 5 \end{pmatrix} and v2=(61)\mathbf{v}_2 = \begin{pmatrix} 6 \\ 1 \end{pmatrix}.
    • The matrix is A=(2651)A = \begin{pmatrix} 2 & 6 \\ 5 & 1 \end{pmatrix}.
    • det(A)=(2×1)(6×5)=230=28\det(A) = (2 \times 1) - (6 \times 5) = 2 - 30 = -28.
    • The area is det(A)=28=28|\det(A)| = |-28| = 28.

Linear Transformations and Geometric Change

  • Theorem 10: Let T:RnRnT: \mathbb{R}^n \rightarrow \mathbb{R}^n be the linear transformation determined by a matrix AA.

    • If n=2n=2 and SS is a parallelogram in R2\mathbb{R}^2, then:
      • area(T(S))=det(A)area(S)\text{area}(T(S)) = |\det(A)| \cdot \text{area}(S)
    • If n=3n=3 and SS is a parallelepiped in R3\mathbb{R}^3, then:
      • volume(T(S))=det(A)volume(S)\text{volume}(T(S)) = |\det(A)| \cdot \text{volume}(S)
  • Proof Summary:

    • A parallelogram SS starting at the origin determined by b1\mathbf{b}_1 and b2\mathbf{b}_2 has the form S={x1b1+x2b2:0x1,x21}S = \{x_1\mathbf{b}_1 + x_2\mathbf{b}_2 : 0 \leq x_1, x_2 \leq 1\}.
    • The image T(S)T(S) consists of points T(x1b1+x2b2)=x1T(b1)+x2T(b2)T(x_1\mathbf{b}_1 + x_2\mathbf{b}_2) = x_1T(\mathbf{b}_1) + x_2T(\mathbf{b}_2).
    • The image T(S)T(S) is a parallelogram determined by the columns of the matrix ABAB, where B=[b1bn]B = [\mathbf{b}_1 \dots \mathbf{b}_n].
    • The area of the image is det(AB)|\det(AB)|, which equals det(A)det(B)|\det(A) \cdot \det(B)|.
    • Since det(B)|\det(B)| is the area of SS, then area(T(S))=det(A)area(S)\text{area}(T(S)) = |\det(A)| \cdot \text{area}(S).
    • This logic extends to arbitrary parallelograms p+Sp + S because translation does not change area.
  • Expansion to Other Regions:

    • Theorem 10 generalizes to any region with a finite area or volume, as these shapes can be approximated by unions of small squares or cubes.
  • Example 5: Area of an Ellipse:

    • Standard equation for an ellipse EE: x12a2+x22b2=1\frac{x_1^2}{a^2} + \frac{x_2^2}{b^2} = 1, where a,b>0a, b > 0.
    • The ellipse EE is the image of the unit disk DD (u12+u221u_1^2 + u_2^2 \leq 1) under the transformation T(u)=AuT(\mathbf{u}) = A\mathbf{u} where:
      • A=(a00b)A = \begin{pmatrix} a & 0 \\ 0 & b \end{pmatrix}
    • Area calculation:
      • area(E)=det(A)area(D)\text{area}(E) = |\det(A)| \cdot \text{area}(D)
      • det(A)=ab\det(A) = ab
      • area(D)=π(1)2=π\text{area}(D) = \pi(1)^2 = \pi
      • area(E)=abπ\text{area}(E) = ab\pi