Published on

Notes on Strang Lecture - 4

Lecture link

LU Transformation

A matrix AA as the product of a lower triangular matrix LL and an upper triangular matrix UU:

A=LUA = LU

Basic example

Consider the matrix: A=[2187]A = \begin{bmatrix} 2 & 1 \\ 8 & 7 \end{bmatrix}

We can factorize this as: A=[1041][2103]=LUA = \begin{bmatrix} 1 & 0 \\ 4 & 1 \end{bmatrix} \begin{bmatrix} 2 & 1 \\ 0 & 3 \end{bmatrix} = LU

Where:

  • L=[1041]L = \begin{bmatrix} 1 & 0 \\ 4 & 1 \end{bmatrix} is the lower triangular matrix
  • U=[2103]U = \begin{bmatrix} 2 & 1 \\ 0 & 3 \end{bmatrix} is the upper triangular matrix

Key property: multipliers in L

Important observation: If no row exchange is required during Gaussian elimination, the multipliers used in the elimination process go directly into the lower triangular matrix LL.

In the example above:

  • To eliminate the entry A21=8A_{21} = 8, we use the multiplier m21=82=4m_{21} = \frac{8}{2} = 4
  • This multiplier appears directly in position (2,1)(2,1) of matrix LL
  • The diagonal entries of LL are always 1

General form

For an n×nn \times n matrix AA, the LU decomposition has the form:

  • LL is lower triangular with ones on the diagonal (the inverse of the elimination matrix)
  • UU is upper triangular (the result from Gaussian elimination)