The LU Factorization
[10/100%]Let .
Determine the LU factors of .
Use the LU factors to solve as well as , where
Use the LU factors to determine .
Let and be the matrices
Explain why does not have an LU factorization.
Use partial pivoting and find the permutation matrix as well as the LU factors such that .
Use the information in , , and to solve .
Determine all values of for which fails to have an LU factorization.
If is a nonsingular matrix that possesses an LU factorization, prove that the pivot that emerges after stages of standard Gaussian elimination using only Type III operations is given by
where and are the leading principal submatrices of orders and , respectively. Use this to deduce that all pivots must be nonzero when an LU factorization for exists.
If is a matrix that contains only integer entries and all of its pivots are , explain why must also be an integer matrix.
This fact can be used to construct random integer matrices that possess integer inverses by randomly generating integer matrices and with unit diagonals and then constructing the product .
Consider the tridiagonal matrix .
Assuming that possesses an LU factorization, verify that it is given by
where the 's are generated by the recursion formula
(This holds for tridiagonal matrices of arbitrary size, making the LU factors of such matrices very easy to compute.)
Apply the recursion formula given above to obtain the LU factorization of
is called a band matrix if whenever for some positive integer , called the bandwidth. In other words, the nonzero entries of are constrained to be in a band of diagonal lines above and below the main diagonal. For example, tridiagonal matrices have bandwidth one, and diagonal matrices have bandwidth zero. If is a nonsingular matrix with bandwidth , and if has an LU factorization , then inherits the lower band structure of , and inherits the upper band structure in the sense that has "lower bandwidth" , and has "upper bandwidth" . Illustrate why this is true by using a generic matrix with a bandwidth of .
Construct an example of a nonsingular symmetric matrix that fails to possess an LU (or LDU) factorization.
Construct an example of a nonsingular symmetric matrix that has an LU factorization but is not positive definite.
Recall the LDU factorization: factoring the diagonal pivots out of the upper-triangular factor of an LU factorization produces , where is the diagonal matrix of pivots and is a rightmost upper-triangular matrix with 's on its diagonal (renaming as hereafter, so both and have unit diagonals in the LDU factorization).
Determine the LDU factors for (this is the same matrix used in exercise 3.10.1).
Prove that if a matrix has an LDU factorization, then the LDU factors are uniquely determined.
If is symmetric and possesses an LDU factorization, explain why it must be given by .
Recall: a symmetric matrix possessing an LU factorization in which each pivot is positive is said to be positive definite, and such a matrix can be uniquely factored as , where is an upper-triangular matrix with positive diagonal entries; this is the Cholesky factorization of , and is the Cholesky factor of .
Explain why is positive definite, and then find the Cholesky factor .