Skip to content

Mathematics · Ch 2 — Matrices

Inverse of a Nonsingular Matrix by Elementary Transformation

2.2.1

Inverse of a Nonsingular Matrix by Elementary Transformation

By the definition of inverse, if A−1A^{-1} exists then AA−1=A−1A=IAA^{-1} = A^{-1}A = I. Consider the equation AA−1=IAA^{-1}=I: here AA is the given matrix of order mm, II is the identity matrix of order mm, and the only unknown is A−1A^{-1}. So to find A−1A^{-1}, the strategy is to convert AA into II using elementary transformations -- and whatever happens to AA must happen, in step, to the matrix standing in for A−1A^{-1}.

Why the same row transformation can be applied to both sides. Whenever an elementary row transformation is applied to the product AB=CAB=C of two matrices, it is enough to apply it only to the prefactor AA -- BB stays unchanged -- and to apply the identical transformation to CC as well; the equation remains true. For example, if A=[1234]A=\begin{bmatrix}1&2\\3&4\end{bmatrix} and B=[−1015]B=\begin{bmatrix}-1&0\\1&5\end{bmatrix}, then AB=[110120]=CAB=\begin{bmatrix}1&10\\1&20\end{bmatrix}=C. Transforming CC by R1↔R2R_1\leftrightarrow R_2 gives C∼[120110]C\sim\begin{bmatrix}1&20\\1&10\end{bmatrix}. Applying the same transformation to AA alone (leaving BB unchanged) gives A∼[3412]A\sim\begin{bmatrix}3&4\\1&2\end{bmatrix}, and now [3412][−1015]=[120110]\begin{bmatrix}3&4\\1&2\end{bmatrix}\begin{bmatrix}-1&0\\1&5\end{bmatrix}=\begin{bmatrix}1&20\\1&10\end{bmatrix} -- exactly the already-transformed CC, confirming the shortcut is valid.

Hence the equation AA−1=IAA^{-1}=I can be transformed into IA−1=BIA^{-1}=B by applying the same series of row transformations to both sides of the equation. Symbolically:

A  A−1=I→Row Transformations (both sides)I  A−1=B ⇒ A−1=BA\ \ A^{-1} = I \xrightarrow{\text{Row Transformations (both sides)}} I\ \ A^{-1} = B \ \Rightarrow\ A^{-1}=B

If, instead, one starts from the equally-valid equation A−1A=IA^{-1}A=I, the transformations used must be column transformations, applied to the postfactor AA and to the right-hand II, while the prefactor A−1A^{-1}-slot stays put:

A−1  A=I→Column Transformations (both sides)A−1  I=B ⇒ A−1=BA^{-1}\ \ A = I \xrightarrow{\text{Column Transformations (both sides)}} A^{-1}\ \ I = B \ \Rightarrow\ A^{-1}=B

A row-only derivation and a column-only derivation must never be mixed inside a single computation of A−1A^{-1} -- pick one system and stay in it throughout.

Standard pivoting order for a 3×33\times3 matrix. For A=[a11a12a13a21a22a23a31a32a33]A = \begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix}, reducing to I3=[100010001]I_3=\begin{bmatrix}1&0&0\\0&1&0\\0&0&1\end{bmatrix} by row transformations typically proceeds: (1) reduce a11a_{11} to 11; (2) then reduce a21a_{21} and a31a_{31} to 00; (3) reduce a22a_{22} to 11; (4) then reduce a12a_{12} and a32a_{32} to 00; (5) reduce a33a_{33} to 11; (6) then reduce a13a_{13} and a23a_{23} to 00. A similar (but not identical) working order is used for column transformations. This is a convenient default order, not a rigid law -- any valid sequence of elementary transformations that reaches the identity is an acceptable derivation.

Worked Example (checking invertibility first). Before computing an inverse it is standard to confirm ∣A∣≠0|A|\neq0. For A=[2142]A=\begin{bmatrix}2&1\\4&2\end{bmatrix}: ∣A∣=2(2)−4(1)=0|A|=2(2)-4(1)=0, so AA is singular and not invertible. For B=[cos⁡θsin⁡θ−sin⁡θcos⁡θ]B=\begin{bmatrix}\cos\theta&\sin\theta\\-\sin\theta&\cos\theta\end{bmatrix}: ∣B∣=cos⁡2θ+sin⁡2θ=1≠0|B|=\cos^2\theta+\sin^2\theta=1\neq0, so BB is non-singular and invertible. For C=[132312123]C=\begin{bmatrix}1&3&2\\3&1&2\\1&2&3\end{bmatrix}: expanding along row 1, ∣C∣=1(1⋅3−2⋅2)−3(3⋅3−2⋅1)+2(3⋅2−1⋅1)=1(−1)−3(7)+2(5)=−1−21+10=−12≠0|C|=1(1\cdot3-2\cdot2)-3(3\cdot3-2\cdot1)+2(3\cdot2-1\cdot1)=1(-1)-3(7)+2(5)=-1-21+10=-12\neq0, so CC is non-singular and invertible.

Worked Example -- inverse of a 2×22\times2 matrix by row transformations. For A=[1234]A=\begin{bmatrix}1&2\\3&4\end{bmatrix}: ∣A∣=1(4)−3(2)=−2≠0|A|=1(4)-3(2)=-2\neq0, so A−1A^{-1} exists. Starting from AA−1=IAA^{-1}=I (row transformations only): [1234]A−1=[1001]\begin{bmatrix}1&2\\3&4\end{bmatrix}A^{-1}=\begin{bmatrix}1&0\\0&1\end{bmatrix}. Using R2→R2−3R1R_2\to R_2-3R_1: [120−2]A−1=[10−31]\begin{bmatrix}1&2\\0&-2\end{bmatrix}A^{-1}=\begin{bmatrix}1&0\\-3&1\end{bmatrix}. Using R2→−12R2R_2\to -\tfrac12R_2: [1201]A−1=[1032−12]\begin{bmatrix}1&2\\0&1\end{bmatrix}A^{-1}=\begin{bmatrix}1&0\\\tfrac32&-\tfrac12\end{bmatrix}. Using R1→R1−2R2R_1\to R_1-2R_2: [1001]A−1=[−2132−12]\begin{bmatrix}1&0\\0&1\end{bmatrix}A^{-1}=\begin{bmatrix}-2&1\\\tfrac32&-\tfrac12\end{bmatrix}. Hence A−1=[−2132−12]A^{-1}=\begin{bmatrix}-2&1\\\tfrac32&-\tfrac12\end{bmatrix}. …

Misc 2.2.1aWorked illustration -- why the SAME row transformation must be applied to both factors of a product

Worked out. A concrete numeric product AB=C is transformed by a chosen row operation; applying that operation to A alone (leaving B fixed) and recomputing AB is shown to reproduce exactly the already-transformed C, justifying why the identity-matrix side of the AA^-1=I bookkeeping can be transformed in step with A without breaking the equation. …

Misc 2.2.1bSolved Examples 1-4 -- inverse by elementary row or column transformations, sizes 2x2 and 3x3

Worked out. Four fully worked examples: Example 1 tests invertibility of a 2x2 numeric matrix, a 2x2 trigonometric (cos/sin) matrix, and a 3x3 numeric matrix by evaluating each determinant; Example 2 finds the inverse of a 2x2 matrix by row transformations, ending with a verified closed-form 2x2 inverse; Example 3 finds the inverse of a 3x3 matrix by a sequence of five row transformations; Example 4 finds the inverse of a different 3x3 matrix using column transformations instead, applied to the equation A−1A=IA^{-1}A=I. …