Algoritmo de la factorización LU

1. Objetivo. Estudiar el algoritmo de la factorización LU de una matriz cuadrada in-

vertible.

2. Requisitos:

Matrices elementales y su relación con operaciones elementales.

Matriz inversa, propiedades de multiplicación de matrices.

Criterio de existencia de la factorización LU .

Explicación del método

Definición (factorización LU ). Una factorización LU de una matriz A ∈ Matn (F) es
un par de matrices (L, U ), donde L, U ∈ Matn (F), U es triangular superior, L es triangular
inferior y todos los elementos diagonales de L son iguales a 1.

3. Unicidad de la factorización LU en el caso de matrices invertibles. Nos

restringimos a matrices invertibles. Si una matriz en invertible y admite una factorización
LU , entonces tal factorización es única.

4. Factorización LU en términos de matrices elementales. Dada una matriz A

cuyos menores de esquina todos son no nulos, construyamos las matrices L y U . Vamos a
construirlas paso a paso. Primero ponemos L := I, U := A. En cada paso del algoritmo
será valida la igualdad A = LU . Empezamos a convertir U en una matriz triangular
superior al aplicar operaciones elementales de tipo Ri + = λRj , j < i. Cada vez, cuando
hacemos la operación Ri + = λRj con las filas de U , i.e. multiplicamos U del lado izquiero
por E+ (i, j, λ), tenemos que multiplicar L del lado derecho por E+ (i, j, −λ), i.e. hacer con
L la operación de columnas Cj − = λCi .

5. Ejemplo con razonamientos extensos. Construyamos la factorización LU de la

matriz  
−1 3 2
A =  3 −4 1 .
2 5 −2

Solución. Podemos escribir A en forma A = LU con L = I, U = A:

  
1 0 0 −1 3 2
A =  0 1 0   3 −4 1 .
0 0 1 2 5 −2

página 1 de 6
Ahora vamos a eliminar el elemento U2,1 = 3 usando el elemento U1,1 = −1 como pivote.
Tenemos que hacer con U la operación de filas R2 + = 3R1 . Es lo mismo que multiplicar U
del lado izquierdo por la matriz elemental E+ (2, 1, 3). Para compensar esta multiplicación
y conservar el mismo valor del producto LU , tenemos que multiplicar L del lado derecho
por E− (2, 1, −3):
    
1 0 0 1 0 0 1 0 0 −1 3 2
A =  0 1 0   −3 1 0   3 1 0   3 −4 1 .
0 0 1 0 0 1 0 0 1 2 5 −2

Multipliquemos U por E+ (2, 1, 3) del lado izquierdo, i.e. hagamos con U la operación de
filas R2 + = 3R1 . Multipliquemos L por E− (2, 1, 3) del lado derecho, i.e. hagamos con L
la operación de columnas C1 − = 3C2 :
  
1 0 0 −1 3 2
A =  −3 1 0   0 5 7 .
0 0 1 2 5 −2

Es fácil checar que el producto de las matrices nuevas es igual a la matriz A. Ahora
queremos eliminar el elemento U3,1 = 2 usando como pivote el elemento U1,1 = −1. Para
esto, metemos entre L y U las matrices E+ (3, 1, −2)E+ (3, 1, 2):
    
1 0 0 1 0 0 1 0 0 −1 3 2
A =  −3 1 0   0 1 0   0 1 0   0 5 7 .
0 0 1 −2 0 1 2 0 1 2 5 −2

Hagamos las operaciones elementales correspondientes (R3 + = 2R1 con U , C1 − = 2C2

con L):   
1 0 0 −1 3 2
A = −3
 1 0  0 5 7 .
−2 0 1 0 11 2
Nos falta eliminar U3,2 = 11 usando U2,2 = 5 como pivote. Metemos entre L y U las
matrices elementales E+ (3, 2, 11/5)E+ (3, 2, 11/5):
    
1 0 0 1 0 0 1 0 0 −1 3 2
A =  −3 1 0   0 1 0  0 1 0  0 5 7 .
−2 0 1 0 11/5 1 0 −11/5 1 0 11 2

Apliquemos la operación elemental de filas R3 − = 11

5
a la matriz U y la operación ele-
11
mental de columnas C2 + = 5 a la matriz L:
  
1 0 0 −1 3 2
A =  −3 1 0  0 5 7 .
−2 11/5 1 0 0 −67/5

página 2 de 6
Respuesta:    
1 0 0 −1 3 2
L =  −3 1 0 , U = 0 5 7 .
−2 11/5 1 0 0 −67/5
Comprobación:
   
−1 + 0 + 0 3+0+0 2+0+0 −1 3 2
LU =  3 + 0 + 0 −9 + 5 + 0 −6 + 7 + 0  =  3 −4 1  = A.
2 + 0 + 0 −6 + 11 + 0 −4 + 77/5 − 67/5 2 5 −2

6. Ejemplo sin razonamientos extensos. Construyamos la factorización LU de la

matriz  
1 3 5
A =  4 15 16  .
−3 −7 −10

Solución. Ahora vamos a escribir las matrices L y U juntas y en vez de las matrices
elementales escribimos sólo las operaciones correspondientes que hacemos con U y L:
   
1 0 0 1 3 5 1 0 0 1 3 5
R −= 4R1 R += 3R1
 0 1 0 4 15 16  −−2−−−−→  4 1 0 0 3 −4  −−3−−−−→
C1 += 4C2 C1 −= 3C3
0 0 1 −3 −7 −10 0 0 1 −3 −7 −10
   
1 0 0 1 3 5 R −= R12 1 0 0 1 3 5
 4 1 0 0 3 −4  −−3−−−3−→  4 1 0 0 3 −4  .
C1 += 23 C3
−3 0 1 0 2 5 −3 2/3 1 0 0 23/3

Respuesta:    
1 0 0 1 3 5
L= 4 1 0 , U = 0 3 −4  .
−3 2/3 1 0 0 23/3
Comprobación:
   
1+0+0 3+0+0 5+0+0 1 3 5
LU =  4 + 0 + 0 12 + 3 + 0 20 − 4 + 0  =  4 15 16  = A.
−3 + 0 + 0 −9 + 2 + 0 −15 − 8/3 + 23/3 −3 −1 −10

7. Observación. En el algoritmo de factorización LU es fácil hacer la comprobación en

cada paso del algoritmo: El producto LU siempre debe ser igual a la matriz A.

página 3 de 6
Notación breve para la factorización LU
8. Notación breve para la factorización LU . Primera observación. Cuando hacemos
una operación Rq + = λRp con las filas de U , donde q > p, tenemos que hacer la operación
Cp − = λCq con las columnas de L. Pero en este momento la q-ésima columna de L
coincide con la q-ésima columna de la matriz identidad. Consiguientemente la operación
Cp − = λCq con las columnas de L equivale al poner Lq,p := −λ.
Segunda observación. Después de eliminar el elemento Uq,p , ya no es necesario guardar
su valor nuevo porque sabemos que este valor nuevo es cero. En el mismo lugar podemos
guardar el valor Lq,p = −λ.
Resumen: trabajamos con una sóla matriz B. En su parte superior (incluyendo la
diagonal) construimos paso a paso la matriz U , y en su parte inferior (sin diagonal)
construimos en el mismo tiempo la matriz L.

9. Ejemplo. Construyamos la factorización LU de la matriz

 
−2 4 −1
A =  4 −5 4 .
−6 −3 −14

Solución. Marcamos los elementos de L con otro color.

     
−2 4 −1 R2 += 2R1 −2 4 −1 −2 4 −1
R += −3R1 R += 5R2
 4 −5 4  −−3−−−−−→  −2 3 2  −−3−−−−→  −2 3 2 .
B2,1 :=−2 B3,2 :=−5
−6 −3 −14 B3,1 :=3 3 −15 −11 3 −5 −1

Respuesta:    
1 0 0 −2 4 −1
L =  −2 1 0 , U = 0 3 2 .
3 −5 1 0 0 −1
Comprobación:
   
−2 + 0 + 0 4 + 0 + 0 −1 + 0 + 0 −2 4 −1
LU =  4 + 0 + 0 −8 + 3 + 0 2 + 2 + 0  =  4 −5 4  = A.
−6 + 0 + 0 12 − 15 + 0 −3 − 10 − 1 −6 −3 −14

10. Ejercicios. Para cada una de las siguientes matrices construya la factorización LU
y haga la comprobación:
     
−1 3 −4 1 4 7 1 4 3
 2 −3 5 ,  −2 −8 5  ,  4 11 10  .
3 −5 7 3 2 6 −3 9 5

página 4 de 6
Aplicación de la factorización LU a los sistemas de ecuaciones
lineales
11. Método. Sean A ∈ Matn (F) una matriz invertible, (L, U ) su factorización LU y
b ∈ Fn . Consideremos el sistema de ecuaciones lineales Ax = b, i.e. LU x = b. Denotemos
U x por y. El sistema Ax = b se puede resolver en dos pasos. Primero, calculamos la
solución y de la ecuación Ly = b. Segundo, calculamos la solución x de la ecuación
U x = y.

12. Ejemplo. Usando la factorización LU resolver el sistema de ecuaciones lineales Ax =

b, donde    
2 −3 1 3
A =  −4 9 2 ,b =  4 .
6 −12 −2 −2
Solución. Primero, factoricemos la matriz A:
     
2 −3 1 R2 += 2R1 2 −3 1 2 −3 1
R −= 3R1 R3 += R2
 −4 9 2  −−3−−−−→  −2 3 4  −−−−−→  −2 3 4 .
6 −12 −2 3 −3 −5 3 −1 −1
De allı́    
1 0 0 2 −3 1
L =  −2 1 0 , U = 0 3 4 .
3 −1 1 0 0 −1
Comprobación:
  
1 0 0 2 −3 1
LU =  −2 1 0  0 3 4 
3 −1 1 0 0 −1
   
2 + 0 + 0 −3 + 0 + 0 1+0+0 2 −3 1
=  −4 + 0 + 0 6 + 3 + 0 −2 + 4 + 0  =  −4 9 2  = A. X
6 + 0 + 0 −9 − 3 + 0 3−4−1 6 −12 −2
Resolvamos el sistema de ecuaciones lineales Ly = b:
    
1 0 0 y1 3 y1 = 3;
 −2 1 0   y2  =  4  . =⇒ y2 = 4 + 2y1 = 10;
3 −1 1 y3 −2 y3 = −2 − 3y1 + y2 = −1.
Luego resolvamos el sistema de ecuaciones lineales U x = y (primero despejemos x3 , luego
x2 y x1 ):
3+3x2 −x3

2 −3 1

x1
 
3
 x1 = 2
= 4;
10−4x3
 0 3 4   x2  =  10  . =⇒ x2 = 3
= 2;
0 0 −1 x3 −1 x3 = −1
= 1.
−1

página 5 de 6
Respuesta:  
4
x =  2 .
1
Comprobación:
      
2 −3 1 4 8−6+1 3
Ax =  −4 9 2   2  =  −16 + 18 + 2  =  4  = b. X
6 −12 −2 1 24 − 24 − 2 −2

13. Ejercicios. Usando la factorización LU resuelva los siguientes sistemas de ecuaciones

lineales (haga todas comprobaciones):
         
−3 2 1 x1 −8 3 −1 3 x1 1
 −6 7 0   x2  =  −3  ,  −6 6 −7   x2  =  −3  .
6 −7 −1 x3 5 9 13 3 x3 5

página 6 de 6