Método Por Descomposición Lu

Para encontrar la matriz triangular inferior se busca hacer ceros los valores de arriba de cada pivote, así
como también convertir en 1 cada pivote. Se utiliza el mismo concepto de "factor" explicado
anteriormente y se ubican todos los "factores" debajo de la diagonal según corresponda en cada uno.
Esquemáticamente se busca lo siguiente:

Originalmente se tenía:

Debido a que [A] = [L][U], al encontrar [L] y [U] a partir de [A] no se altera en nada la ecuación y se
tiene lo siguiente:

Por lo tanto, si Ax = b, entonces LUx = b, de manera que Ax = LUx = b.

PASOS PARA RESOLVER UN SISTEMA DE ECUACIONES POR EL MÉTODO DE DESCOMPOSICIÓN LU
1. Obtener la matriz triangular inferior L y la matriz triangular superior U.
2. Resolver Ly = b (para encontrar y).
3. El resultado del paso anterior se guarda en una matriz nueva de nombre "y".
4. Realizar Ux = y (para encontrar x).
5. El resultado del paso anterior se almacena en una matriz nueva llamada "x", la cual brinda
los valores correspondientes a las incógnitas de la ecuación.

EJEMPLO 1 DE DESCOMPOSICIÓN LU
PROBLEMA: Encontrar los valores de x1, x2 y x3 para el siguiente sistema de ecuaciones:

NOTA: Recuérdese que si la matriz es 2x2 se hará 1 iteración; si es 3x3, 2 iteraciones; si es 4x4, 3
iteraciones; y así sucesivamente.
 SOLUCIÓN:

4 -2 -1 9

[A] = 5 1 -1 [B] = 7

1 2 -4 12

ITERACIÓN 1
factor 1 = (a21 / a11) = 5 / 4 = 1.25
factor 2 = (a31 / a11) = 1 / 4 = 0.25
Encontrando [U]
fila 2 = - (factor 1) * (fila 1) + (fila 2)
fila 3 = - (factor 2) * (fila 1) + (fila 3)
a11 = a11
a12 = a12
a13 = a13
a21 = - (1.25) * (4) + (5) = 0
a22 = - (1.25) * (- 2) + (1) = 3.5
a23 = - (1.25) + (- 1) + (- 1) = 0.25
a31 = - (0.25) * (4) + (1) = 0
a32 = - (0.25) * (- 2) + (2) = 2.5
a33 = - (0.25) * (- 1) + (- 1) = - 0.75

4 -2 -1

[U] = 0 3.5 0.25

0 2.5 - 0.75

Encontrando [L]

1 0 0

[L] = 1.25 0 0

0.25 0 0

ITERACIÓN 2
factor 3 = (u32 / u22) = 2.5 / 3.5 = 0.7142857143
Encontrando [U]
fila 3 = - (factor 3) * (fila 2) + (fila 3)
a31 = - (2.5 / 3.5) * (0) + (0) = 0
a32 = - (2.5 / 3.5) * (3.5) + (2.5) = 0
a33 = - (2.5 / 3.5) * (0.25) + (- 0.75) = - 0.9285714286

4 -2 -1

[U] = 0 3.5 0.25

0 0 - 0.9285714286

Encontrando [L]

1 0 0

[L] = 1.25 1 0

0.25 0.7142857143 1

Ahora ya se tiene la matriz [U] y la matriz [L]. El siguiente paso es resolver

Ly = b para encontrar la matriz y. En pocas palabras es como que se pidiera resolver el siguiente
sistema de ecuaciones, encontrando los valores de y1, y2 y y3:

Al resolver el sistema anterior, se obtienen los siguientes valores para y1, y2 y y3:

El último paso es resolver Ux = y para encontrar la matriz x. En otras palabras es como que se pidiera
resolver el siguiente sistema de ecuaciones, encontrando los valores de x1, x2 y x3:

La solución del sistema es:

Este es finalmente el valor de x1, x2 y x3; es decir, la respuesta del ejercicio utilizando la
descomposición LU.
EJEMPLO 2 DE DESCOMPOSICIÓN LU
PROBLEMA: Encontrar los valores de x1, x2 y x3 para el siguiente sistema de ecuaciones:

SOLUCIÓN:

11 -3 -2 18

[A] = 5 -2 -8 [B] = 13

4 -7 2 2

ITERACIÓN 1
factor 1 = (a21 / a11) = 5/11 = 0.4545454545
factor 2 = (a31 / a11) = 4/11 = 0.3636363636

Encontrando [U]
fila 2 = - (factor 1) * (fila 1) + (fila 2)
fila 3 = - (factor 2) * (fila 1) + (fila 3)
a11 = a11
a12 = a12
a13 = a13
a21 = - (0.4545454545) * (11) + (5) = 0
a22 = - (0.4545454545) * (- 3) + (- 2) = - 0.6363636365
a23 = - (0.4545454545) + (- 2) + (- 8) = - 7.0909090919
a31 = - (0.3636363636) * (11) + (4) = 0
a32 = - (0.3636363636) * (- 3) + (- 7) = - 5.909090909
a33 = - (0.3636363636) * (- 2) + (2) = 2.7272727272
 11 -3 -2

[U] = 0 - 0.6363636365 - 7.0909090919

0 - 5.909090909 2.7272727272

Encontrando [L]

1 0 0

[L] = 0.45454545 0 0

0.36363636 0 0

ITERACIÓN 2
factor 3 = (u32/u22) = - 5.909090909 / - 0.6363636365 = 9.285714284

Encontrando [U]
fila 3 = - (factor 3) * (fila 2) + (fila 3)
a31 = - (9.285714284) * (0) + (0) = 0
a32 = - (9.285714284) * (- 0.6363636365) + (- 5.909090909) = 0
a33 = - (9.285714284) * (- 7.0909090919) + (2.7272727272) = 68.57142857

11 -3 -2

-
[U] = 0 - 0.6363636365 7.0909090919

0 0 68.57142857

Encontrando [L]

1 0 0

[L] = 0.4545454545 1 0

0.3636363636 9.285714284 1

Ahora ya se tiene la matriz [U] y la matriz [L]. El siguiente paso es resolver

Ly = b para encontrar la matriz y. En pocas palabras es como que se pidiera resolver el siguiente
sistema de ecuaciones, encontrando los valores de y1, y2 y y3:
Al resolver el sistema anterior, se obtienen los siguientes valores para y1, y2 y y3:

El último paso es resolver Ux = y para encontrar la matriz x. En otras palabras es como que se pidiera
resolver el siguiente sistema de ecuaciones, encontrando los valores de x1, x2 y x3:

La solución del sistema es:

Este es finalmente el valor de x1, x2 y x3; es decir, la respuesta del ejercicio utilizando la
descomposición LU.

EJEMPLO 3 DE DESCOMPOSICIÓN LU

Si en la matriz

1 2 3
A = 2 4 1 , intercambiamos las filas 1ra y 3ra.
obtenemos
3 2 -1

3 2 -1
A’ = 1 2 1
1 2 3
Si de nuevo intercambiamos las mismas filas (1ra. Y 3ra.), obtenemos la matriz A original.

En consecuencia: el intercambio de (las mismas) filas es la operación elemental inversa del

intercambio de filas que es lo que se señala en a) mirando el asunto de una manera más
práctica.

La igualdad algebraica expresada en B se refleja a nivel de filas ya que si una fila se multiplica
por un número c0 y luego se divide por el mismo número (se multiplica por 1/c) esta
termina inalterada. O sea que estas operaciones elementales son inversas, lo cual fue
expresado en b).

Examinemos con un ejemplo lo que sucede cuando a una fila se sustituye por su suma con
múltiplo de otra y luego a la fila resultante se le resta de nuevo el mismo múltiplo de la fila
que se sumó. Por supuesto que obtenemos la fila original.

Si a partir de
1 2 3
A = 2 4 1 , sumamos a la 1ra fila el doble de la segunda.
obtenemos 3 2 -1

Lo cual es equivalente a premultiplicarla por la matriz E 1 + (2) 2 , obtenida de la matriz

idéntica por el mismo cambio elemental, obtenemos

5 10 5
A = 2 4 1
3 2 -1
Si ahora restamos a la 1ra. Fila el doble de la segunda, o sea premultiplicamos por E 1 - (2) 2 ,
retornamos a la matriz original.
Esto que es claro cuando se trata de filas puede probarse algebraicamente ya que
5 10 5 1 2 0 1 2 3

A = 2 4 1 = E 1 + (2) 2 A ‘= 0 1 0 2 4 1
3 2 -1 0 0 1 3 2 -1
Recuerde que la matriz E 1 + (2) 2 se obtiene a partir de la matriz idéntica utilizando la misma
operación elemental (sumandole a la primera fila (1 0 0) el doble de la segunda 2(0 1
0)). Lo supremamente llamativo es que la premultiplicación de una matriz por estas
matrices elementales construidas con operaciones elementales del tipo a) , b) o c), generan
en la matriz el mismo cambio elemental que las definió.

Es claro que E 1 - (2) 2 . E 1 + (2) 2 = I. O sea que E -1 1 + (2) 2 = E 1 - (2) 2.

Revisemos una vez más estas operaciones por filas.

1 2 3 F1 + 2F2 5 10 5 F1 - 2F2 1 2 3
A = 2 4 1 2 4 1 2 4 1 = A
3 2 -1 3 2 -1 3 2 -1

Estos razonamientos nos permiten justificar el siguiente procedimiento para descomponer

una matriz en la forma A = LU. Inversa (sobre I) operación(directa)(sobre A)
F2 + F 1 F2 - F1

1 1 -1 1 0 0 1 1 -1 1 0 0 1 1 -1
A = 1 2 1 = IA = 0 1 0 1 2 1 1 1 0 0 1 2
1 1 1 0 0 1 1 1 1 0 0 1 1 1 1

(operación inversa sobre L)

(-1)
(inversa) L U F3 - F1
(+1) 1 0 0 1 1 -1 Operación directa (sobre A)
1 1 0 0 1 2
F3 + F1 1 0 1 0 0 2
operación inversa
He aquí la descomposición LU de A.

La matriz U es similar a la queda cuando se aplica a A el proceso de Gauss.. Debemos anotar

los siguientes puntos claves con el fin de automatizar la producción de la matriz triangular
inferior L.

a) La matriz L es triangular inferior con 1s en la diagonal

b) En todos los pasos de la producción de L los elementos sobre la diagonal son 0s.
c) Los números que van apareciendo en L, a medida que avanza el proceso,
corresponden al inverso de la operación realizada en las filas de A, en la búsqueda
de la matriz U.

Así: la primera operación sobre las filas de A fue F2 - F1. La operación inversa F2 + F1
genera la matriz
1 0 0
L = 1 1 0 El 1 es el coeficiente de F1 en la operación inversa
0 0 1

donde el elemento enmarcado proviene de la operación inversa realizada en la matriz

idéntica I.

La siguiente operación en búsqueda de la matriz triangular superior U, fue F 3 - F1, la

operación inversa es F3 + F1 . Generándose el nuevo número (en este caso 1 de + F1 , en la
fila 3ra. de la 1ra. Columna).
La nueva matriz L es

1 0 0
L = 1 1 0
1 0 1

Estos pasos pueden consignarse en una sóla matriz lo cual nos permite consignar las
modificaciones producidas a medida que se construye la matriz U y utilizar los espacios de
los 0s en la parte inferior de U para ir consignando los elementos de la matriz triangular
inferior L. La diagonal de L no es necesario calcularla pues sabemos que está formada por
puros 1s. Los elementos de la diagonal que se consignarán serán los de U. Llamaremos a
este método, cálculo de LU con sobreescritura.

Así:
1 1 -1
Descompongamos de nuevo la matriz A = 1 2 1 en su forma LU utilizando el
método de sobreescritura.. 1 1 1
 Operación directa sobre A
Para obtener 0 en la posición señalada de U
F2 - F1 F2 + F1 F3 - F1

1 1 -1 1 1 -1 1 1 -1
A = 1 2 1 Operación inversa en L 1 1 2 1 1 2
1 1 1 1 1 1 F3 + F1 1 0 2
Operación inversa en L

El proceso termina ya que ya apareció un 0 en la 3ra. Fila, 2da. columna y

“milagrosamente” ya apareció la matriz U.

Luego:

1 0 0 1 1 -1
A = 1 1 0 0 1 2
1 0 1 0 0 2
L U