Algoritmo de factorización PLU

Objetivos. Estudiar el algoritmo de factorización PLU de una matriz cuadrada invertible.

Requisitos. Factorización LU, matrices de permutación.

1. Definición (factorización PLU). Sea A ∈ Mn (R) una matriz invertible. Una terna
de matrices (P, L, U) se llama factorización PLU de A si P, L, U ∈ Mn (R), PA = LU,
U es una matriz triangular superior con elementos diagonales no nulos, L es una matriz
triangular inferior con elementos diagonales iguales a 1 y P es una matriz de permutación.

2. Ejemplo. Encontrar una factorización PLU de la matriz

 
2 2 3
A =  −4 −4 −3  .
4 8 3

Solución.
 
2 2 3 R2 += 2R1
R3 += −2R1
P1,2,3 A =  −4 −4 −3  −−−−−−→
4 8 3
 
2 2 3
R ↔R
P1,2,3 A =  −2 0 3  −−2−−→
3

2 4 −3
 
2 2 3
P1,3,2 A =  2 4 −3  .
−2 0 3

Factorización PLU en forma explı́cita:

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

Comprobación:
   
2+0+0 2+0+0 3+0+0 2 2 3
LU =  4 + 0 + 0 4+4+0 6−3+0 = 4 8 2  = P1,3,2 A. X
−4 + 0 + 0 −4 + 0 + 0 −6 + 0 + 3 −4 −4 −3

Algoritmo de factorización PLU, página 1 de 6

3. Ejemplo. Usando una factorización PLU resolver el siguiente sistema de ecuaciones
lineales:    
2 2 3 −7
 −4 −4 −3  x =  −1  .
4 8 3 5

Solución. Para aplicar la factorización P1,3,2 A = LU del ejemplo anterior multiplicamos

ambos lados de la ecuación por la matriz de permutación P1,3,2 :
 
−7
LUx = c, donde c = P1,3,2 b =  5  .
−1

Denotamos Ux por y. Primero resolvemos el sistema Ly = c:


 y1 = −7 y1 = −7
Ly = c : 2y1 + y2 = 5 y2 = 5 − 2y1 = 19

−2y1 + y3 = −1 y3 = −1 + 2y1 = −15

Ahora resolvemos el sistema Ux = y:


 2x1 + 2x2 + 3x3 = −7 x1 = (−7 − 2x2 − 3x3 )/2 = 3
Ux = y : 4x2 − 3x3 = 19 x2 = (19 + 3x3 )/4 = 1

3x3 = −15 x3 = −5

Respuesta:  
3
x =  1 .
−5
Comprobación:
      
2 2 3 3 6 + 2 − 15 −7
 −4 −4 −3   1  =  −12 − 4 + 15  =  −1  . X
4 8 3 −5 12 + 8 − 15 5

Algoritmo de factorización PLU, página 2 de 6

4. Ejemplo. Encontrar una factorización PLU de la matriz
 
0 0 4
A= 1 3 2 .
2 8 4

Solución.
 
0 0 4
R ↔R
P1,2,3 A =  1 3 2  −−1−−→
2

2 8 4
 
1 3 2 R2 += 0R1
R += −2R
P2,1,3 A =  0 0 4  3
−− −−−−→
1

2 8 4
 
1 3 2
R ↔R
P2,1,3 A =  0 0 4  −−2−−→
3

2 2 0
 
1 3 2
P2,3,1 A =  2 2 0  .
0 0 4

Factorización PLU en forma explı́cita:

     
0 1 0 0 0 4 1 0 0 1 3 2
 0 0 1  1 3 2  =  2 1 0  0 2 0 .
1 0 0 2 8 4 0 0 1 0 0 4

Comprobación:
   
1+0+0 3+0+0 2+0+0 1 3 2
LU =  2 + 0 + 0 6 + 2 + 0 4 + 0 + 0  =  2 8 4  = P2,3,1 A. X
0+0+0 0+0+0 0+0+4 0 0 4

Algoritmo de factorización PLU, página 3 de 6

5. Ejemplo. Resolver el sistema de ecuaciones lineales:
   
0 0 4 −8
 1 3 2 x =  1 .
2 8 4 6

Solución. El sistema de ecuaciones dado se puede escribir como Ax = b. En el ejemplo

anterior hemos obtenido una factorización PA = LU para esta misma matriz A. Para
aplicar esta factorización multipliquemos ambos lados de la igualdad Ax = b por P por
la izquierda:
   
−8 1
PAx = c, donde c = Pb = P2,3,1  1  =  6  .
6 −8

Para resolver la ecuación LUx = c denotamos Ux por y y resolvemos la ecuación Ly = c:


 y1 = 1 y1 = 1
Ly = c : 2y1 + y2 = 6 y2 = 6 − 2y1 = 4

y3 = −8 y3 = −8

Luego consideramos la ecuación Ux = y:


 x1 + 3x2 + 2x3 = 1 x1 = 1 − 3x2 − 2x3 = −1
Ux = y : 2x2 = 4 x2 = 4/2 = 2

4x3 = −8 x3 = −8/4 = −2

Respuesta:  
−1
x =  2 .
−2
Comprobación:
      
0 0 4 −1 0+0−8 −8
Ax =  1 3 2   2  =  −1 + 6 − 4  =  1  = b. X
2 8 4 −2 −2 + 16 − 8 6

Algoritmo de factorización PLU, página 4 de 6

6. Ejemplo. Construir una factorización PLU de la siguiente matriz A:
 
−3 1 1 0
 −6 2 2 2 
A=  3 3
.
2 −1 
3 7 2 0

Solución.
   
−3 1 1 0 R2 += −2R1 −3 1 1 0
R3 += R1
 −6 2 2 2  R4 += R1  2 0 0 2  R ↔R
P1,2,3,4 A = 
 3
 −−−−−−→ P1,2,3,4 A =   −−2−−→
3

3 2 −1   −1 4 3 −1 
3 7 2 0 8 3 0
   
−3 1 1 0 −3 1 1 0
 −1 4 3 −1  R4 += −2R2 −1 4 3 −1  R ↔R
−−−−→ P1,3,2,4 A =  −−3−−→
 4
P1,3,2,4 A = 
 2
 −− 
0 0 2   2 0 0 2 
8 3 0 −1 2 −3 2
 
−3 1 1 0
 −1 4 3 −1 
P1,3,4,2 A = 
 −1
.
2 −3 2 
2 0 0 2

Comprobación:
  
1 0 0 0 −3 1 1 0
 −1 1 0 0   0 4 3 −1 
LU =  
 −1 2 1 0   0 0 −3

2 
2 0 0 1 0 0 0 2
 
−3 + 0 + 0 + 0 1 + 0 + 0 + 0 1+0+0+0 0+0+0+0
 3 + 0 + 0 + 0 −1 + 4 + 0 + 0 −1 + 3 + 0 + 0 0−1+0+0 
= 3 + 0 + 0 + 0 −1 + 8 + 0 + 0 −1 + 6 − 3 + 0

0−2+2+0 
−6 + 0 + 0 + 0 2 + 0 + 0 + 0 2+0+0+0 0+0+0+2
 
−3 1 1 0
 3 3 2 −1 
= 3 7 2
 = P1,3,4,2 A. X
0 
−6 2 2 2

Algoritmo de factorización PLU, página 5 de 6

7. Ejemplo. Resolver el sistema de ecuaciones lineales a través de una factorización PLU:
   
2 3 −2 5
 2 4 −1  x =  7  .
−2 1 8 1

8. Ejercicio. Resuelva el siguiente sistema de ecuaciones lineales a través de una facto-

rización PLU:    
1 1 2 7
 −2 −2 8  x =  10  .
3 2 4 12

9. Ejercicio. Resuelva el siguiente sistema de ecuaciones lineales a través de una facto-

rización PLU:    
0 0 3 −12
 4 2 1 x =  2 .
4 6 5 6

Algoritmo de factorización PLU, página 6 de 6