Exercise 1.5.8. This exercise continues the discussion in p. 33-34 regarding factorization of a 3×3 matrix A and the proof that A = LU. An alternative proof begins with the fact that row 3 of U is produced by taking row 3 of A and subtracting times row 1 of U and
times row 2 of U. (Note that for simplicity we assume that no row exchanges are necessary.)
(a) Why do we subtract rows of U rather than rows of A?
(b) The above equation can be rearranged and restated as follows: row 3 of A is equal to the sum of times row 1 of U,
times row 2 of U, and 1 times row 3 of U. Which rule of multiplication (from 1D on p. 25) makes this also equal to row 3 of LU?
Note that applying the results of (a) and (b) for the other rows of U completes the proof.
Answer: (a) By the time we apply elimination to row 3 of A we have already completed the elimination steps for rows 1 and 2 of A. The first step produced row 1 of U, and the second step produced row 2 of U. When we come to the elimination step for row 3 we are therefore subtracting multiples of rows of U, not rows of A.
By similar logic when performing an elimination step on any pivot row we have completed elimination on all the rows above it, and therefore any row subtracted from the pivot row (using the appropriate multiplier) is a row of U.
(b) From 1D part (iii) we have each row of the product matrix AB equal to the product of a row and a matrix; more specifically, row i of AB is equal to row i of A times the matrix B.
So for the product matrix LU we have row 3 of LU equal to row 3 of L times the matrix U.
Now we have
so that row 3 of L consists of the entries ,
, and 1, and since U is upper triangular we have
The product of row 3 of L and the matrix U (which is also row 3 of the product LU) is therefore
which can be expressed as
But this is times row 1 of U plus
times the row 2 of U plus 1 times row 3 of U, and from above we have this equal to row 3 of A. Therefore row 3 of LU is equal to row 3 of A.
By similar arguments we show that row 1 of LU is equal to row 1 of A, and row 2 of LU is equal to row 2 of A. Since each row of LU is equal to the corresponding row of A, we have A = LU.
NOTE: This continues a series of posts containing worked out exercises from the (out of print) book Linear Algebra and Its Applications, Third Edition by Gilbert Strang.
If you find these posts useful I encourage you to also check out the more current Linear Algebra and Its Applications, Fourth Edition, Dr Strang’s introductory textbook Introduction to Linear Algebra, Fourth Edition
and the accompanying free online course, and Dr Strang’s other books
.