LA FACTORIZACIN LU Y LA MATRIZ INVERSA
MATEMTICAS II GRADO EN INGENIERA DE QUMICA
(16 de marzo de 2012)
La intencin de este complemento (a los apuntes del curso) es comprender por qu el mtodo de
eliminacin gaussiana permite realizar factorizaciones LU . De paso se justifica un mtodo (quizs
conocido) para el clculo de la inversa de una matriz dada.
En todo lo que sigue se entiende que
las matrices son del orden adecuado para que los productos tengan sentido;
si hablamos de la inversa de una matriz es porque sabemos que la inversa existe;
I es la matriz identidad y ser siempre del orden adecuado para que las operaciones tengan
sentido.
1.
Expresin matricial de la eliminacin gaussiana y de la factorizacin LU
Consideremos la matriz
2
1
2
A = 1 2 3 .
3 4 1
Para hacer ceros por debajo de la diagonal aplicamos el mtodo de eliminacin gaussiana (en lo que
sigue Fi indicar la fila i-sima):
1. Restando (en A)
2. Restando (en U1 )
1
2 F1
11
2
3
2
a F2 y 32 F1 a F3 , obtenemos la matriz
2 1
2
2 .
U1 = 0 3
2
0 11
2
2
F2 (es decir,
11
3 F2 )
a F3 , obtenemos la matriz
U2 = 0
0
1
3
2
2 ,
16
3
que ya es triangular superior.
En este ejemplo hemos dado dos pasos que pueden ser expresados mediante el producto de matrices.
1. Restar
1
2 F1
a F2 y 32 F1 a F3 equivale a multiplicar (por la izquierda) por la matriz
1
0 0
1 0 0
1 0 = 12 1 0 .
L1 = 1
2
23 0 1
32 0 1
Ejercicio 1: Comprueba que L1 A = U1 .
1
MATEMTICAS II GRADO EN INGENIERA DE QUMICA
2. Restar
11
3 F2
a F3 equivale a multiplicar (por la izquierda) por la matriz
1
0
0
1
0 .
L2 = 0
0 11
1
3
Ejercicio 2: Comprueba que L2 U1 = U2 .
Es interesante observar que las matrices L1 y L2 se obtienen aplicando, a la matriz identidad, las
operaciones por filas indicadas en cada caso.
De las igualdades L1 A = U1 y L2 U1 = U2 , deducimos que
L2 U1 = L2 (L1 A) = (L2 L1 ) A = U2 .
Pero, teniendo en cuenta que el producto de matrices triangulares inferiores da lugar a otra matriz
triangular inferior, lo que hemos conseguido es la expresin
L A = U,
donde L = L2 L1 es triangular inferior y U = U2 es triangular superior. Finalmente, puesto que la
inversa de una matriz triangular inferior es otra matriz triangular inferior,
A = (L )1 U.
Por tanto, si consideramos L = (L )1 , tenemos una descomposicin LU para A. Adems, si revisamos los clculos, vemos que U es fcil de obtener, pero qu ocurre con L? Veamos que tampoco es
muy difcil de calcular.
Si recordamos que la inversa del producto de dos matrices es el producto de las inversas cambiadas
de orden, tenemos que
L = (L )1 = (L2 L1 )1 = (L1 )1 (L2 )1 .
Para justificar las expresiones de (L1 )1 y (L2 )1 debemos tener en cuenta qu significa que una
matriz sea la inversa de otra. Veamos, si M y N son matrices inversas una de la otra, entonces
(M N ) B = M (N B) = B
para cualquier matriz B. Es decir, tras multiplicar sucesivamente B por N y M tenemos de nuevo
B, o sea, se puede decir que no hemos hecho nada.
Con tal interpretacin de las matrices inversas, parece claro que la inversa de L1 ha de ser una
matriz que haga lo contrario que ella misma, esto es, una matriz que nos permita pasar de U1 a A:
hacemos desaparecer los ceros de la primera columna de U1 para recuperar la primera columna de A.
3
Para ello, necesitamos sumar 1
2 F1 a F2 y 2 F1 a F3 , lo cual equivale a multiplicar (por la izquierda)
por la matriz
1 0 0
1 0 .
(L1 )1 = 1
2
3
0 1
2
Ejercicio 3: Comprueba que (L1 )1 L1 = L1 (L1 )1 = I y que (L1 )1 U1 = A.
Razonando de la misma forma, para la inversa de L2 necesitamos sumar
a multiplicar (por la izquierda) por la matriz
1 0 0
(L2 )1 = 0 1 0 .
0 11
1
3
11
3 F2
a F3 , lo que equivale
Ejercicio 4: Comprueba que (L2 )1 L2 = L2 (L2 )1 = I y que (L2 )1 U2 = U1 .
LA FACTORIZACIN LU Y LA MATRIZ INVERSA
Ahora veamos cul es el resultado de multiplicar las inversas de L1 y L2 (quizs la cuenta ms
sorprendente de este desarrollo):
1 0 0
1
0 0
1 0 0
1 0 0 1 0 = 1
1 0 .
(L1 )1 (L2 )1 = 1
2
2
11
3
3
0 11
1
0 1
1
3
2
2
3
Parece que basta con agrupar en una matriz la primera columna de (L1 )1 y la segunda columna de
(L2 )1 (y aadir la ltima columna de cualquiera de ellas por ser la misma). En realidad, el producto
anterior no tiene nada de sorprendente si se piensa que, al estar (L1 )1 a la izquierda de (L2 )1 , lo
que estamos haciendo es llevar a cabo en (L2 )1 las operaciones por filas indicadas por (L1 )1 (que,
3
recordemos, son sumar 1
2 F1 a F2 y sumar 2 F1 a F3 ).
Tenemos, pues, una manera bastante fcil de construir la factorizacin LU de Doolittle de una
matriz dada. En el caso de la matriz A considerada,
2 1
2
1
0 0
2
1
2
1 0 0 3
2 .
A = 1 2 3 = (L )1 U = (L1 )1 (L2 )1 U2 = 1
2
2
11
3
16
3 4 1
1
0 0
2
3
3
Conviene sealar que el nmero de cuentas para realizar la factorizacin de Doolittle de esta
manera es el mismo que el nmero de cuentas que son necesarias cuando se aplica el algoritmo visto
para implementar en un ordenador. Sin embargo, esta forma de proceder es ms inestable que la del
algoritmo cuando se usan clculos aproximados, esto es, es ms susceptible a los errores de redondeo.
Para acabar esta seccin, realicemos un ejemplo con matrices de orden 4 4. Sea la matriz
2 1 1 1
1 2 1 2
B=
0 1 1 1 .
2 1 1 3
1. Hacemos ceros en la primera columna de B restando 12 F1 a F2 , 20 F1 a F3 y 22 F1 a F4 , quedando
2 1 1 1
0 3 1 3
2
2
2
la matriz U1 =
0 1 1 1 .
0
1
1
2
Observemos que L1 B = U1 si tomamos L1 =
0
1
0
1
0
0
1
2. Hacemos ceros en la segunda columna de U1 restando
3 F2
2
2 1 1 1
0 3
1
32
2
2
a F4 , quedando la matriz U2 =
.
2
0 0
0
3
4
0 0 3 2
1 0
0 1
Observemos que L2 U1 = U2 si tomamos L2 =
2
0 3
0 34
0
0
.
0
1
= 32 F2 a F3 y
0
0
1
0
0
0
1
0
0
0
0
1
2
32
F2 = 43 F2
MATEMTICAS II GRADO EN INGENIERA DE QUMICA
4
3. Por ltimo, hacemos ceros en la tercera columna de U2 restando 23 F3 (= 2F3 ) a F4 , que3
2 1 1 1
0 32 12 32
dando la matriz U3 =
, que ya es triangular superior.
0 0 2 0
3
0 0 0 2
1 0 0 0
0 1 0 0
Observemos que L3 U2 = U3 si tomamos L3 =
0 0 1 0 .
0 0 2 1
Concluimos que la factorizacin LU de Doolittle para la matriz B es
1 0
0 0
2 1 1 1
1 1
1 2 1 2
0 0
2
=LU =
B=
0 1 1 1
0 23
1 0
2 1 1 3
1 34 2 1
la siguiente,
2 1
0 32
0 0
0 0
1
1
3
2 2
2
0
3
0 2
1
Observemos que, para construir L, hemos tomado las tres entradas de la primera columna de L1
que estn debajo de la diagonal principal, las dos entradas de la segunda columna L2 que tambin
estn debajo de la diagonal principal y la nica entrada de la tercera columna de L3 que est, de
nuevo, debajo de la diagonal principal pero, y esto es muy importante, con signo cambiado.
2.
Cmo obtener otras factorizaciones
Supongamos que deseamos obtener la factorizacin LU de Crout. Podemos aprovechar el trabajo
realizado en la seccin anterior? La respuesta es s, como veremos a continuacin.
Recordemos que, en la factorizacin de Crout, es la matriz triangular superior la que tiene unos en
la diagonal principal. Por tanto, si volvemos al caso de la matriz A vista anteriormente, necesitaramos
pasar de la matriz
2 1
2
2
U2 = 0 3
2
16
0 0
3
a otra matriz con unos en la diagonal. Para ello, bastar con multiplicar F1 por 21 , F2 por - 23 y F3
3
por 16
, quedando la matriz
1 12 1
U3 = 0 1 43 .
0 0 1
1
2
Ejercicio 5: Comprueba que U3 = D U2 , siendo D = 0 23
0 0
0 .
3
16
Ejercicio 6: Sea D la matriz del ejercicio 5. Comprueba que D1
2 0
= 0 32
0 0
0 .
16
3
LA FACTORIZACIN LU Y LA MATRIZ INVERSA
Teniendo en cuenta que D D1 = D1 D = I (y que la matriz identidad I acta en el producto
de matrices como el 1 en el producto de nmeros),
2
1
2
A = 1 2 3 = L U = L I U = L D1 D U = LC UC
3 4 1
y, as, LC UC es la factorizacin de Crout al tomar LC = L D1 y UC = D U . Por tanto,
1
0 0
2 0
0
2
0
0
1 0 0 32 0 = 1 32 0 ,
LC = 1
2
3
11
16
0 0 16
3 11
1
2
3
3
2
3
1
1
2 1
2
0
0
1 2 1
2
2 = 0 1 43 .
UC = 0 32 0 0 3
2
16
3
0 0
0 0 1
0 0 16
3
Ejercicio 7: Comprueba que A = LC UC .
Para la matriz B de orden 44 que consideramos en la seccin anterior, puesto que la factorizacin
de Doolittle es
1 0
0 0
2 1 1 1
2 1 1 1
1 1
1
3
3
1 2 1 2
0 0
0 2 2 2
2
=LU =
B=
0 1 1 1
0 23
1 0 0 0 23 0
2 1 1 3
0 0 0 2
1 34 2 1
1
0 0 0
2
0 2 0 0
3
tomando la matriz D =
0 0 3 0 , la descomposicin de Crout ser B = LC UC con
2
0 0 0 21
1 0
0 0
2 0
0 0
2 0 0 0
1 1
3
3
0 0
0 2 0 0 1 2 0 0
LC = L D1 = 2
0 23
1 0 0 0 32 0 0 1 23 0
1 43
1
2
0
0 2
3
UC = D U =
0 0
0 0
0 0 0 2
2 1
0 0
2 1 1 1
3
1
3
0 0
0 2 2 2
3
0 0 32 0
2 0
0 0 0 2
0 21
1 12
0 1
0 0
0 0
43
1
2
13
1
0
2
12
1
0
1
Ejercicio 8: Comprueba que B = LC UC .
A la vista de lo realizado con las matrices A y B para obtener la factorizacin de Crout, es claro
que podramos haber escogido otras matrices
diagonales
y, de esta manera, obtendramos diferentes
1 0 0
factorizaciones. Por ejemplo, con D = 0 2 0 , a partir de la factorizacin de Doolittle de A,
0 0 3
1
0 0
2 1
2
2
1
2
1 0 0 3
2 =
A = 1 2 3 = 1
2
2
3
11
16
3 4 1
1
0 0
2
3
3
MATEMTICAS II GRADO EN INGENIERA DE QUMICA
1
2
3
2
0
1
11
3
0
1 0
0 0 21
1
0 0
1
2
3
2
1
2
11
6
2
0
1 0 0
0 0 2 0 0
1
0 0 3
0
3
1
3
2
2 =
16
3
0
2 1
2
0 0 3 4
1
0 0 16
3
y conseguimos una factorizacin en la que la matriz triangular superior no contiene fracciones.
3.
Un mtodo para calcular la inversa de una matriz
En esta seccin vamos a justificar el funcionamiento de un mtodo (quizs conocido) que est
basado en la eliminacin gaussiana y que permite calcular matrices inversas. La idea de este mtodo
es aplicar eliminacin gaussiana para hacer ceros tanto por debajo como por encima de la diagonal
y, adems, conseguir que todos los elementos de la diagonal sean igual a uno.
Consideremos la matriz A vista en las secciones anteriores para describir dicho mtodo. Adems,
para ser ms prcticos, las operaciones que se vayan realizando sobre A tambin se aplican a la
matriz identidad I. Para simplificar el proceso, se yuxtaponen la matriz A y la matriz identidad de
la siguiente forma,
2
1
2 1 0 0
(A | I) = 1 2 3 0 1 0
3 4 1 0 0 1
Empezamos haciendo ceros en A por debajo de la diagonal. Para ello, primero restamos
y 32 F1 a F3 , operando sobre las filas de A e I conjuntamente, y tenemos
2 1
2 1 0 0
2 12 1 0 .
(M1 | N1 ) = 0 3
2
2 32 0 1
0 11
2
A continuacin, restamos
1
2 F1
a F2
11
3 F2
a F3 para obtener la matriz
2 1
2 1
2 21
(M2 | N2 ) = 0 3
2
16
0 0
10
3
3
0
1
11
3
0 .
1
Ejercicio 9: Comprueba que estas mismas operaciones ya las hicimos antes (aunque slo sobre A y
con los nombres de las matrices cambiados).
1 0 0
1
0
0
1
0 , comprueba que M1 = L1 A,
Ejercicio 10: Tomando L1 = 12 1 0 y L2 = 0
1
0 11
23 0 1
3
V1 = L1 I, M2 = L2 M1 y N2 = L2 N1 .
Ahora vamos a hacer ceros en M2 por encima de la diagonal. Primero en la tercera columna y
despus en la segunda. Por tanto, empezamos restando 83 F3 a F1 y 38 F3 a F2 y obtenemos
11
3
2 1
0 49
8
8
3
3
3
0
(M3 | N3 ) = 0 3
2
4
8
8 .
16
10
11
0 0
3 3
1
3
LA FACTORIZACIN LU Y LA MATRIZ INVERSA
Y restando 23 F2 a F1 queda
9
2 0
0 47
8
3
3
0
(M4 | N4 ) = 0 3
2
4
8
10 11
0 0 16
3
3
3
1 0 83
1 23 0
3
Ejercicio 11: Tomando U1 = 0 1 8 y U2 = 0 1 0
0 0 1
0 0 1
N3 = U1 N2 , M4 = U2 M3 y N4 = U2 N3 .
Para acabar tenemos que hacer unos en la
(M5 | N5 ) = 0
0
18
3
8
, comprueba que M3 = U1 M2 ,
diagonal de M4 . Para ello, hacemos 12 F1 , 23 F2 y
9
1
16
0 0 78
16
1
1 0 12
14 .
4
3
0 1 5 11
8
16
3
16 F3 ,
16
Hemos alcanzado nuestro objetivo pues M5 es la matriz identidad.
1
0
0
2
Ejercicio 12: Tomando D = 0 23 0 , comprueba que M5 = I = D M4 , N5 = D N4 .
3
0 0 16
Ejercicio 13: Comprueba que A N5 = N5 A = I.
Vamos a justificar por qu este mtodo nos proporciona la matriz inversa de A. Por una lado,
I = M5 = D M4 = D U2 M3 = D U2 U1 M2 = D U2 U1 L2 M1 = D U2 U1 L2 L1 A.
Por tanto, D U2 U1 L2 L1 es la inversa de A. Pero, por otro lado,
N5 = D N4 = D U2 N3 = D U2 U1 N2 = D U2 U1 L2 N1 = D U2 U1 L2 L1 I,
es decir, N5 = D U2 U1 L2 L1 es la inversa de A.
Nota 3.1. Con respecto a lo hecho para la matriz A, debemos tener en cuenta una serie de consideraciones
1. En algunos textos, antes de hacer los ceros por encima de la diagonal principal, se convierten los
elementos de la diagonal principal en unos. Como se puede ver en el Ejemplo 3.2, el resultado
final es el mismo y la justificacin de la validez del proceso es anloga a la realizada en el
prrafo precedente.
2. En algunas matrices hay que realizar cambios de filas para poder conseguir elementos distintos
de cero en la diagonal principal. El mtodo de obtencin de la inversa sigue siendo vlido y en
la justificacin del proceso se emplean las llamadas matrices de transposicin. Se ver esta
situacin en el Ejemplo 3.3.
Ejemplo 3.2. Recordemos los dos primeros pasos dados en el clculo de la inversa de A,
2
1
2 1 0 0
(A | I) = 1 2 3 0 1 0
3 4 1 0 0 1
Primero restamos
1
2 F1
a F2 y 32 F1 a F3 ,
(M1 | N1 ) = 0
0
1
3
2
11
2
2
2
2
1
1
2
3
2
0 0
1 0 .
0 1
MATEMTICAS II GRADO EN INGENIERA DE QUMICA
A continuacin, restamos
11
3 F2
a F3 ,
(M2 | N2 ) = 0
0
1
3
2
2
2
16
3
1
1
2
10
3
0
1
11
3
0 .
1
Ahora, en lugar de hacer ceros por encima de la diagonal principal (de M2 ), vamos a hacer unos en
3
esta diagonal. Para ello, consideramos 12 F1 , 23 F2 y 16
F3 ,
0
0
1 12 1 21
M30 | N30 = 0 1 43 13 32 0 .
3
0 0 1 58 11
16
16
Hacemos ceros en la tercera columna tomando F1 F3 y F2 43 F3 ,
11
3
16
1 21 0 89
16
1
M40 | N40 = 0 1 0 21
14 .
4
3
0 0 1 58 11
16
16
Acabamos haciendo ceros en la segunda columna mediante la operacin F1 12 F2 ,
9
1
1 0 0 78
16
16
1
1
M50 | N50 = 0 1 0 12
4
4 .
5
11
3
0 0 1
8 16
16
Por consiguiente, la inversa de A es la matriz N50 que, como no puede ser de otra manera, coincide
con la obtenida anteriormente.
Ejemplo 3.3. Veamos un ejemplo en el que, para calcular la inversa, es necesario realizar un cambio
de filas. Para ello, consideramos la matriz
2 3 2
C = 2 3 3 .
1 2 3
Para empezar escribimos
2 3 2 1 0 0
(C | I) = 2 3 3 0 1 0 .
1 2 3 0 0 1
Hacemos ceros en la primera columna con F2 F1 y F3 12 F1 ,
2 3 2 1 0 0
(M1 | N1 ) = 0 0 5 1 1 0 .
0 12 2 12 0 1
En este momento necesitamos hacer un cambio de filas para poder seguir operando. Concretamente
tenemos que intercambiar la segunda y tercera filas, lo que expresaremos de la forma F2 F3 .
2 3 2 1 0 0
(M2 | N2 ) = 0 12 2 12 0 1 .
0 0 5 1 1 0
LA FACTORIZACIN LU Y LA MATRIZ INVERSA
Hacemos ceros en la tercera columna con F1 + 52 F3 y F2 + 25 F3 ,
2
2 3 0 35
5 0
9
2
(M3 | N3 ) = 0 21 0 10
.
5 1
0 0 5 1 1 0
Hacemos ceros en la segunda columna con F1 6F2 ,
2 0 0 6
9
(M4 | N4 ) = 0 12 0 10
0 0 5 1
2 6
2
1 .
5
1
0
Finalmente, hacemos unos en la diagonal con 12 F1 , 2F2 y 51 F3 ,
1 0 0 3 1 3
4
2 .
(M5 | N5 ) = 0 1 0 95
5
1
0 0 1 5 15 0
3 1 3
4
2 es la inversa de C.
Concluimos que N5 = 95
5
1
51 0
5
Ejercicio 14: Comprueba que C N5 = N5 C = I.
Ejercicio 15: Comprueba la operacin F2 F3 se puede expresar mediante la matriz (de transposicin)
1 0 0
T23 = 0 0 1 .
0 1 0
Comprueba que M2 = T23 M1 y N2 = T23 N1 .
Ejercicio 16: Describe, por medio de matrices, todos los pasos dados para calcular la inversa de C.
Ejercicio 17: Repite el clculo de la inversa de C pero haciendo unos en la diagonal principal antes
de hacer ceros por encima de dicha diagonal.
(Observa que tienes que hacer los clculos de nuevo a partir de (M2 | N2 )).
Ejemplo 3.4. Para acabar, calculemos la inversa
secciones anteriores. Partimos de la expresin
2 1 1
1 2 1
(B | I) =
0 1 1
2 1 1
de la matriz B de orden 4 4 ya utilizada en las
Hacemos ceros (por debajo) en la primera columna
2 1 1
0 32 12
(M1 | N1 ) =
0 1 1
0 2 2
1
2
1
3
1
0
0
0
0
1
0
0
0
0
1
0
restando 12 F1 a
1 1 0
32 21 1
1 0 0
4 1 0
0
0
.
0
1
F2 y F1 a F4 ,
0 0
0 0
.
1 0
0 1
10
MATEMTICAS II GRADO EN INGENIERA DE QUMICA
Hacemos ceros (por debajo) en la segunda columna restando 32 F2 a
2 1 1 1 1
0
0 3
1
3
1
2 2 1
2
2
(M2 | N2 ) =
2
0 0
0 13 23
3
4
0 0 4 2 5
3
Hacemos ceros (por debajo) en la tercera
2
0
(M3 | N3 ) =
0
0
0 1
columna restando 2F3 a
0
1 1 1 1
1
1
3
3
2 2 2 2 1
0 23 0 13 23
0 0 2 1 0
Hacemos ceros (por encima) en la cuarta columna restando 12 F4
2 1 1 0 12
0
0 3 1 0 5 1
2
2
4
(M4 | N4 ) =
0 0 23 0 13 23
0 0 0 2 1 0
Hacemos ceros (por encima) en la tercera columna restando
2 1 0 0 0
0 3 0 0 3
2
2
(M5 | N5 ) =
0 0 23 0 13
0 0 0 2 1
F4 ,
0
0
1
2
0
0
0
1
a F1 y 34 F4 a F2 ,
1 21
3
3
2
4
.
1 0
2
a F1 y 34 F3 a F2 ,
1 12 21
3
3
3
2
4
4
.
23 1 0
0
2 1
3
2 F3
Hacemos ceros (por encima) de la segunda columna restando 23 F2
2 0 0 0 1
0
0 3 0 0 3
3
2
2
2
(M6 | N6 ) =
0 0 23 0 13 23
0 0 0 2 1 0
Finalmente, hacemos unos en la diagonal
1
0
(M7 | N7 ) =
0
0
F3 y 43 F2 a F4 ,
0 0
0 0
.
1 0
a F1 ,
1 0
3
4
1
2
3
4
.
0
1
principal tomando 21 F1 , 23 F2 ,
0 12 0
0 0 0 12
1 0 0 1 1 12 12
0 1 0 12 1 32
0
1
1
0 0 1
2 0
1
2
3
2 F3
y 12 F4 ,
Por tanto, N7 es la inversa de B.
Ejercicio 18: Comprueba que B N7 = N7 B = I.
Ejercicio 19: Describe, por medio de matrices, todos los pasos dados para calcular la inversa de B.
Ejercicio 20: Calcula, usando el mtodo expuesto
0 1
3 0
E=
2 3
1 2
en esta seccin, la inversa de la matriz
2 3
1 2
.
0 1
3 0
