Cuatro ejemplos de factorización P A = LU y dos ejemplos de

factorización de Crout
Ejemplo 1: Una sola permutación de filas al comienzo
Encontrar una factorización de la forma P A = LU para la matriz
 
0 0 2
A = −1 5 −2 ,
3 6 7
utilizando pivoteo parcial.
     
0 0 2 ←
− 3 6 7 1/3 3 6 7
−1 5 −2 −1 5 −2 ←
− + 0 7 1/3
3 6 7 ←
− 0 0 2 0 0 2
Es entonces  
3 6 7
U = 0 7 1/3
0 0 2
La matriz L es  
1 0 0
L = −1/3 1 0
0 0 1
Para matriz P , la matriz de permutación, observamos que cambiamos la fila 1 y 3, es decir
   
1 0 0 ← − 0 0 1
P = 0 1
 0 0 1 0 .
0 0 1 ← − 1 0 0
Se comprueba que
     
0 0 1 0 0 2 1 0 0 3 6 7
0 1 0 −1 5 −2 = −1/3 1 0 0 7 1/3 .
1 0 0 3 6 7 0 0 1 0 0 2

Ejemplo 2: Dos permutaciones de filas

Encontrar una factorización de la forma P A = LU para la matriz
 
0 1 1
A = −1 2 −4 ,
2 −5 1
utilizando pivoteo parcial.

         
0 1 1 ←
− 2 −5 1 1/2 2 −5 1 2 −5 1 2 −5 1
−1 2 −4 −1 2 −4 ←− + 0 −1/2 −7/2 ←
− 0 1 1  1/2 0 1 1 .
2 −5 1 ←
− 0 1 1 0 1 1 ←
− 0 −1/2 −7/2 ←−+ 0 0 −3
Luego  
2 −5 1
U = 0 1 1 .
0 0 −3

1
La matriz L antes del segundo intercambio de filas era
 
1 0 0
−1/2 1 0
0 0 1
Pero al intercambiar las filas 2 y 3 la matriz L cambia a
 
1 0 0
 0 1 0 .
−1/2 0 1
Después de este intercambio queda por agregar la última operación entre filas, y L es finalmente
 
1 0 0
L= 0 1 0 .
−1/2 −1/2 1
Para P , la matriz de permutación, observamos que cambiamos primero la fila 1 y 3, y luego la
fila 2 y 3, es decir      
1 0 0 ← − 0 0 1 0 0 1
P = 0 1 0 0 1 0 ←− 1 0 0 .
0 0 1 ← − 1 0 0 ← − 0 1 0
Se comprueba entonces que
     
0 0 1 0 1 1 1 0 0 2 −5 1
1 0 0 −1 2 −4 =  0 1 0 0 1 1 .
0 1 0 2 −5 1 −1/2 −1/2 1 0 0 −3

Ejemplo 3: Dos permutaciones de filas

Encontrar una factorización de la forma P A = LU para la matriz
 
2 −2 1
A = −8 11 5 ,
4 −13 3
utilizando pivoteo parcial.
     
2 −2 1 ← − −8 11 5 1/4 1/2 −8 11 5
−8 11 5 ← − 2 −2 1 ← −+ 0 3/4 9/4  ←
−.
4 −13 3 4 −13 3 ←−−−− + 0 −15/2 11/2 ←
−
   
−8 11 5 −8 11 5
 0 −15/2 11/2 1/10  0 −15/2 11/2 .
0 3/4 9/4 ←− + 0 0 14/5

Luego  
−8 11 5
U =  0 −15/2 11/2 .
0 0 14/5

La matriz L antes del segundo intercambio de filas era

 
1 0 0
−1/4 1 0
−1/2 0 1

2
Pero al intercambiar las filas 2 y 3 la matriz L cambia a
 
1 0 0
−1/2 1 0 .
−1/4 0 1

Después de este intercambio queda por agregar la última operación entre filas, y L es finalmente
 
1 0 0
L = −1/2 1 0 .
−1/4 −1/10 1

Para P , la matriz de permutación, observamos que cambiamos primero la fila 1 y 2, y luego la

fila 2 y 3, es decir      
1 0 0 ← − 0 1 0 0 1 0
P = 0 1 0 ← − 1 0 0 ←− 0 0 1 .
0 0 1 0 0 1 ← − 1 0 0
Se comprueba entonces que
     
0 1 0 2 −2 1 1 0 0 −8 11 5
0 0 1 −8 11 5 = −1/2 1 0   0 −15/2 11/2 .
1 0 0 4 −13 3 −1/4 −1/10 1 0 0 14/5

Ejemplo 4: Tres permutaciones de filas

Encontrar una factorización de la forma P A = LU para la matriz
 
1 4 0 −4
 5 1 1 −1
A=  3 1 −1 −2 ,


−3 4 6 2

utilizando pivoteo parcial.

     
1 4 0 −4 ← − 5 1 1 −1 −1/5 −3/5 3/5 5 1 1 −1
5
 1 1 −1  ←−  1 4 0 −4 ←
  − + 0

19/5 −1/5 −19/5 ←
 −
3 1 −1 −2  3 1 −1 −2 ←−−−−− + 0 2/5 − /5 − /5 
8 7
−3 4 6 2 −3 4 6 2 ←−−−−−−−−−− + 0 23/5 33/5 7/5 ←
−
     
5 1 1 −1 5 1 1 −1 5 1 1 −1
0 23/5 33/5 7/5  −2/23 −19/23 0 23/5 33/5 7/5  0 23/5 33/5 7/5 
     
0 2/5 − /5 − /5 ←
8 7  −+ 0 0 − /23 −35/23 
50 ←− 0 0 − /23 − /23
 130 114 −5/13
0 19/5 −1/5 −19/5 ←−−−−−− + 0 0 −130/23 −114/23 ←− 0 0 −50/23 −35/23 ←−+
 
5 1 1 −1
0 23/5 33/5 7/5 
0 0 −130/23 −114/23 .
 

0 0 0 5/13

Luego  
5 1 1 −1
0 23/5 33/5 7/5 
U = .
0 0 −130/23 −114/23
0 0 0 5/13

3
La matriz L antes del segundo intercambio de filas era
 
1 0 0 0
 1/5 1 0 0
 .
 3/5 0 1 0
−3/5 0 0 1

Al cambiar las filas 2 y 3, resulta  

1 0 0 0
−3/5 1 0 0
 .
 3/5 0 1 0
1/5 0 0 1
Con las nuevas operaciones elementales, antes del nuevo cambio de filas resulta
 
1 0 0 0
−3/5 1 0 0
 .
 3/5 2/23 1 0
1/5 19/23 0 1

La matriz vuelve a cambiar por la permutación de las filas 3 y 4

 
1 0 0 0
−3/5 1 0 0
 .
 1/5 19/23 1 0
3/5 2/23 0 1

Finalmente, L resulta  
1 0 0 0
−3/5 1 0 0
L=
 1/5 19/23
.
1 0
3/5 2/23 5/23 1
Para la matriz P , observamos
       
1 0 0 0 ← − 0 1 0 0 0 1 0 0 0 1 0 0
 ←−  ←
−
0 1 0 0 1 0 0 0 0 0 0 1 0 0 0 1
      .
0 0 1 0 0 0 1 0 0 0 1 0 ←
− 1 0 0 0
0 0 0 1 0 0 0 1 ←− 1 0 0 0 ←− 0 0 1 0

Se comprueba que
     
0 1 0 0 1 4 0 −4 1 0 0 0 5 1 1 −1
0 0 0 1  5
  1 1 −1 −3/5
= 1 0 0
 0
 23/5 33/5 7/5 
.
1 0 0 0  3 1 −1 −2  1/5 19/23 1 0 0 0 −130/23 −114/23
0 0 1 0 −3 4 6 2 3/5 2/23 5/23 1 0 0 0 5/13

Factorización de Crout a partir de una factorización de Doolittle

Todas las factorizaciones anteriores fueron realizadas según la condición de Doolittle (1’s en la
diagonal principal de L). La factorización según Crout implica 1’s en la diagonal principal de U .

4
Ejemplo 1: La matriz a factorizar es inversible
Sea la matriz  
1 3 5
A = −2 4 6
5 0 3
que se puede expresar según Doolittle como
    
1 3 5 1 0 0 1 3 5
A = −2 4 6 = −2 1 0  0 10 16  .
5 0 3 5 − /2 −1
3 0 0 −2

A partir de esta factorización es fácil construir una factorización de Crout:

A = LDD−1 U
−1
= |{z}
LD D | {z U}
L? U?

donde D es una matriz diagonal que tiene en su diagonal principal los respectos elementos de
la diagonal principal de U , D−1 su inversa y L? y U ? son matrices triangulares que verifican las
condiciones de Crout:
   
1 0 0 1 0 0
D = 0 10 0  , D−1 = 0 1/10 0 .
0 0 −2 0 0 −1/2
   
1 0 0 1 3 5
L? = LD = −2 10 0 , U ? = D−1 U = 0 1 58  .
5 −15 2 0 0 1
Finalmente,     
1 3 5 1 0 0 1 3 5
−2 4 6 = −2 10 0 0 1 58  .
5 0 3 5 −15 2 0 0 1
| {z } | {z } | {z }
A L? U?

Ejemplo 2: La matriz a factorizar no es inversible

Consideramos ahora el caso en que no es posible construir la matriz D (Esto ocurre cuando la
matriz a factorizar no es inversible). Sea la matriz
 
1 2 3
A = 3 2 1
4 4 4

Razonamos sobre el esquema para el producto de matrices

1 * *
U? 0 1 *
L? 0 0 1
* 0 0 1 2 3
* * 0 3 2 1
* * * 4 4 4

5
 ? = 1, porque a
donde los ∗ son entradas a determinar. Observemos que l11 11 = 1 y que de este
? ?
hecho se deduce que u12 = 2 y u13 = 3

1 2 3
U? 0 1 *
L? 0 0 1
1 0 0 1 2 3
* * 0 3 2 1
* * * 4 4 4
? = 3, porque a
Observamos ahora que l21 ?
21 = 3. Además l31 = 4, porque a31 = 4

1 2 3
U? 0 1 *
L? 0 0 1
1 0 0 1 2 3
3 * 0 3 2 1
4 * * 4 4 4
? = −4 porque a
Ahora deducimos que l22 22 = 2, que resulta del producto escalar
   
3 −4 0 · 2 1 0 .

Deducimos el elemento l33 ? : como A es no invertible, L? tampoco lo será. Por ser L? triangular

inferior, el producto de los elementos de su diagonal principal debe ser nulo, de donde resulta
necesariamente l33? = 0.

1 2 3
U? 0 1 *
L? 0 0 1
1 0 0 1 2 3
3 −4 0 3 2 1
4 * 0 4 4 4
Finalmente u?23 = 2 porque a23 = 1 y l32
? = −4 porque a32 = 4. La factorización es
    
1 2 3 1 0 0 1 2 3
3 2 1 = 3 −4 0  0 1 2
4 4 4 4 −4 0 0 0 1
| {z } | {z } | {z }
A L? U?