350 likes | 1.07k Views
LU Decomposition. Method If an nxn matrix A has an LU-factorization, then the solution of AX = b can be determined by a Forward substitution followed by a Back substitution [ A ] = [ L ][ U ] where [ L ] = lower triangular matrix all of whose diagonal
E N D
LU Decomposition Method If an nxn matrix A has an LU-factorization, then the solution of AX = b can be determined by a Forward substitution followed by a Back substitution [A] = [L][U] where [L] = lower triangular matrix all of whose diagonal entries are different from zero. [U] = upper triangular matrix all of whose diagonal entries are different from zero.
LU Decomposition How can this be used? • Given [A][X] = [b] • Decompose [A] into [L] and [U] • (LU)X = b • L(UX) = b ; Let UX = z • then LZ = b • Use Forward substitution to Solve [L][Z] = [b] for [Z] • Use Back substitution to Solve • [U][X] = [Z] for [X]
[L][Z] = [b] Solve for [Z] using Forward Substitution :
[U][X] = [Z] Solve for [X] using Back Substitution :
Method: [A] Decompose to [L] and [U]Storage of Multipliers Scheme :
Exercise 2.5; Page #136; Qn # 5 Finding the [U] matrix Step 1:
Finding the [U] Matrix Matrix after Step 1: Step 2:
Finding the [L] matrix http://numericalmethods.eng.usf.edu
[L][Z] = [b] Solve for [Z] using Forward Substitution :
[U][X] = [Z] Solve for [X] using Back Substitution :
Limitations to use this Procedure : • Interchanging of any two rows not allowed . • Only elementary row operation permitted is the one that subtract a multiple of one row to another. • In matrix A, if then this procedure fails. • In matrix if , or in if then this procedure fails.