Skip to content

Mathematics · Ch 2 — Matrices

Method of Inversion

2.3.1

Method of Inversion

As the name suggests, this method uses the inverse of the coefficient matrix directly. Consider three linear equations

a1x+b1y+c1z=d1,a2x+b2y+c2z=d2,a3x+b3y+c3z=d3.a_1x+b_1y+c_1z=d_1,\qquad a_2x+b_2y+c_2z=d_2,\qquad a_3x+b_3y+c_3z=d_3.

As explained in section 2.3, these can be expressed as

[a1b1c1a2b2c2a3b3c3][xyz]=[d1d2d3]i.e.AX=B,\begin{bmatrix}a_1&b_1&c_1\\a_2&b_2&c_2\\a_3&b_3&c_3\end{bmatrix}\begin{bmatrix}x\\y\\z\end{bmatrix}=\begin{bmatrix}d_1\\d_2\\d_3\end{bmatrix}\quad\text{i.e.}\quad AX=B,

where A,X,BA,X,B are of orders 3×33\times3, 3×13\times1, 3×13\times1 respectively.

The derivation. If the system's solution exists, AA must be non-singular, so A−1A^{-1} exists (found either by the elementary-transformation method or the adjoint method of section 2.2). Pre-multiplying both sides of AX=BAX=B by A−1A^{-1}:

A−1(AX)=A−1(B)A^{-1}(AX)=A^{-1}(B)

i.e.(A−1A)X=A−1B\text{i.e.}\quad(A^{-1}A)X=A^{-1}B

i.e.IX=A−1B\text{i.e.}\quad IX=A^{-1}B

i.e.X=A−1B,\text{i.e.}\quad X=A^{-1}B,

which gives the required solution -- every unknown is read off from the single matrix product A−1BA^{-1}B.

Worked Example -- two equations. Solve 2x+5y=12x+5y=1 and 3x+2y=73x+2y=7 by the method of inversion. Writing AX=BAX=B with A=[2532]A=\begin{bmatrix}2&5\\3&2\end{bmatrix}, X=[xy]X=\begin{bmatrix}x\\y\end{bmatrix}, B=[17]B=\begin{bmatrix}1\\7\end{bmatrix}: ∣A∣=2(2)−3(5)=4−15=−11|A|=2(2)-3(5)=4-15=-11, and adj A=[2−5−32]\text{adj}\,A=\begin{bmatrix}2&-5\\-3&2\end{bmatrix}, so A−1=1−11[2−5−32]A^{-1}=\dfrac{1}{-11}\begin{bmatrix}2&-5\\-3&2\end{bmatrix}. Then

X=A−1B=1−11[2−5−32][17]=1−11[2−35−3+14]=1−11[−3311]=[3−1].X=A^{-1}B=\frac{1}{-11}\begin{bmatrix}2&-5\\-3&2\end{bmatrix}\begin{bmatrix}1\\7\end{bmatrix}=\frac{1}{-11}\begin{bmatrix}2-35\\-3+14\end{bmatrix}=\frac{1}{-11}\begin{bmatrix}-33\\11\end{bmatrix}=\begin{bmatrix}3\\-1\end{bmatrix}.

By equality of matrices, x=3x=3 and y=−1y=-1.

Worked Example -- three equations. Solve x−y+z=4x-y+z=4, 2x+y−3z=02x+y-3z=0, x+y+z=2x+y+z=2 by the method of inversion. The matrix equation is [1−1121−3111][xyz]=[402]\begin{bmatrix}1&-1&1\\2&1&-3\\1&1&1\end{bmatrix}\begin{bmatrix}x\\y\\z\end{bmatrix}=\begin{bmatrix}4\\0\\2\end{bmatrix}, i.e. AX=BAX=B. Here adj A=[422−5051−23]\text{adj}\,A=\begin{bmatrix}4&2&2\\-5&0&5\\1&-2&3\end{bmatrix} and ∣A∣=10|A|=10, so A−1=110[422−5051−23]A^{-1}=\dfrac{1}{10}\begin{bmatrix}4&2&2\\-5&0&5\\1&-2&3\end{bmatrix}. Then

X=A−1B=110[422−5051−23][402]=110[20−1010]=[2−11].X=A^{-1}B=\frac{1}{10}\begin{bmatrix}4&2&2\\-5&0&5\\1&-2&3\end{bmatrix}\begin{bmatrix}4\\0\\2\end{bmatrix}=\frac{1}{10}\begin{bmatrix}20\\-10\\10\end{bmatrix}=\begin{bmatrix}2\\-1\\1\end{bmatrix}.

By equality of matrices, x=2x=2, y=−1y=-1, z=1z=1. …

Misc 2.3.1aSolved Examples 1-2 -- solving a two-variable and a three-variable system by the method of inversion

Worked out. Example 1 sets up two linear equations as AX=B, finds A^-1 by the adjoint formula, and multiplies A^-1 by B to read off both unknowns; Example 2 repeats the same three-step process (set up AX=B, find A^-1 via its adjoint and determinant, compute A^-1 B) for a three-equation, three-unknown system, reading off all three unknowns from the resulting column matrix. …