Define these matrices in Sage and verify that \(BA =
I\) so that \(B=A^{-1}\text{.}\)
(b)
Find the solution to the matrix equation \(A\xvec =
\threevec{4}{-1}{4}\) using \(A^{-1}\text{.}\)
(c)
Using Sage, multiply \(A\) and \(B\) in the opposite order; that is, what do you find when you evaluate \(AB\text{?}\)
Activity3.1.0.2.
Suppose that \(A\) is an \(n\times n\) invertible matrix with inverse \(A^{-1}\text{.}\) This means that every equation of the form \(A\xvec=\bvec\) has a solution, namely, \(\xvec = A^{-1}\bvec\text{.}\) Which of the following best describes a restatement of this fact?
The columns of \(A\) are linearly independent.
The columns of \(A\) span \(\real^3\text{.}\)
Activity3.1.0.3.
Suppose that \(A\) is an \(n\times n\) invertible matrix with inverse \(A^{-1}\text{.}\) This means that every equation of the form \(A\xvec=\bvec\) has a solution, namely, \(\xvec = A^{-1}\bvec\text{.}\) What can you conclude about the pivot positions of the matrix \(A\text{?}\)
Every column of \(A\) has a pivot position.
Every row of \(A\) has a pivot position.
Every row and every column of \(A\) has a pivot position.
Activity3.1.0.4.
If \(A\) is an invertible \(4\times4\) matrix, what is its reduced row echelon form?
Using the Sage cell below, solve these equations for the columns of \(B\text{.}\)
(b)
What is the matrix \(B\text{?}\) Check that \(AB = I\) and \(BA = I\text{.}\)
(c)
To find the columns of \(B\text{,}\) we solved two equations, \(A\bvec_1=\evec_1\) and \(A\bvec_2=\evec_2\text{.}\) We could do this by augmenting \(A\) two separate times, forming matrices
\begin{equation*}
\begin{aligned}
\left[\begin{array}{r|r} A \amp \evec_1 \end{array}\right] \amp
\\
\left[\begin{array}{r|r} A \amp \evec_2 \end{array}\right] \amp
\\
\end{aligned}
\end{equation*}
and finding their reduced row echelon forms. But instead of solving these two equations separately, we could also solve them together by forming the augmented matrix \(\left[\begin{array}{r|rr} A \amp \evec_1 \amp \evec_2
\end{array}\right]\) and finding the row reduced echelon form. In other words, we augment \(A\) by the matrix \(I\) to form \(\left[\begin{array}{r|r}
A \amp I \end{array}
\right]
\text{.}\)
Form this augmented matrix and find its reduced row echelon form to find \(A^{-1}\text{.}\)
Assuming \(A\) is invertible, we have shown that
\begin{equation*}
\left[\begin{array}{r|r}
A \amp I
\end{array}\right]
\sim
\left[\begin{array}{r|r}
I \amp A^{-1}
\end{array}\right]\text{.}
\end{equation*}
Activity3.1.0.6.
If you have defined a matrix \(A\) in Sage, you can find it's inverse as A.inverse(). Use Sage to find the inverse of the matrix
Suppose that \(n\times n\) matrices \(C\) and \(D\) are both invertible. What do you find when you simplify the product \((D^{-1}C^{-1})(CD)\text{?}\) Explain why this shows the product \(CD\) is invertible and \((CD)^{-1} =
D^{-1}C^{-1}\text{.}\)
We have seen how to use Gaussian elimination to find the inverse of a matrix. Now we want to look at how to use matrix multiplication to perform Gaussian Elimination.
Activity3.1.0.9.
Tweaking the identity matrix slightly allows us to write row operations in terms of matrix multiplication.
(a)
Create a matrix that doubles the third row of \(A\text{:}\)
Express these row operations as matrix multiplication by expressing \(B\) as the product of two matrices and \(A\text{:}\)
\begin{equation*}
B =
\left[\begin{array}{ccc}
\unknown&\unknown&\unknown\\
\unknown&\unknown&\unknown\\
\unknown&\unknown&\unknown
\end{array}\right]
\left[\begin{array}{ccc}
\unknown&\unknown&\unknown\\
\unknown&\unknown&\unknown\\
\unknown&\unknown&\unknown
\end{array}\right]
A
\end{equation*}
Check your work using technology.
We have now seen how to do Gaussian elimination with matrix multiplication. If we are careful to describe our row operation with lower triangular matrices, then we can connect matrix multiplication to finding the inverse. We will see this through the next couple of activities.
When performing Gaussian elimination on \(A\text{,}\) we first apply a row replacement operation in which we multiply the first row by \(-2\) and add to the second row. After this step, we have a new matrix \(A_1\text{.}\)
Find the matrices \(L_2\) and \(L_3\) that perform these row replacement operations so that \(L_3L_2L_1 A = U\text{.}\)
(b)
Explain why the matrix product \(L_3L_2L_1\) is invertible and use this fact to write \(A = LU\text{.}\) What is the matrix \(L\) that you find? Why do you think we denote it by \(L\text{?}\)
Activity3.1.0.13.
Row replacement operations may always be performed by multiplying by a lower triangular matrix. It turns out the other two row operations, scaling and interchange, may also be performed using matrix multiplication. For instance, consider the two matrices