Metodos Numericos

3.4 METODO DE DESCOMPOSICION LU.
La factorizaci´
on LU de una matriz es una factorización que resume el proceso de eliminaci´ on gaussiana
umero total de operaciones de punto flotante cuando
aplicado a la matriz y que es conveniente en terminos del n´
se desea calcular la inversa de una matriz o cuando se resolverá una serie de sistemas de ecuaciones con una
misma matriz de coeficientes. En la lectura, primeramente consideraremos la factorización LU sin intercambio
basada en matrices elementales y que es conocida como de Doolittle y posteriormente veremos el algoritmo
que da la factorización PA = LU.
on LU
Factorizaci´
Suponga que la matriz A es una matriz m × n se puede escribir como el producto de dos matrices:
A = LU
donde L es una matriz triangular inferior m × m y U es una matriz escalonada m × n. Entonces para resolver
el sistema:
A x = b,
escribimos
A x = (L U) x = L (U x) .
Una posible estrategia de solución consiste en tomar y = U x y resolver para y:
L y = b.
Como la matriz L es triangular superior este sistema puede resolverse mediante sustitución hacia abajo, lo
cual se hace fácilmente en m2 FLOPS. Una vez con los valores encontrados de y, las incógnitas al sistema
inicial se resuelve despejando x de
U x = y.
C. Chapra Raymond, P. Canale.

Métodos numéricos para ingenieros.
Quinta edición; McGRAW-HILL.
Nuevamente, como U es escalonada, este sistema puede resolverse en caso de tener soución mediante susti-
tución hacia atrás, lo cual es sencillo. Estas observaciones nos dan la pauta para ver la conveniencia de una
factorización como la anterior, es decir factorizar A como el producto de una matriz L triangular superior, por
otra U la cual es escalonada. Esta factorización se llama usualmente Descomposición LU.
Uso de la factorizaci´
on LU
Ejemplo:
on LU de A:
Use la factorizaci´
    
4 −2 1 1 0 0 4 −2 1
A =  20 −7 12  =  5 1 0   0 3 7  = LU
−8 13 17 −2 3 1 0 0 −2
para despejar x del sistema:  

11
A x =  70  = b
17
Solución
Sea y = (y1 , y2 , y3 ) un nuevo vector de incógnitas. Primero resolveremos el sistema triangular inferior L y = b:
   
1 0 0 11
 5 1 0  y =  70 
−2 3 1 17
Este sistema escrito en su forma de ecuaciones queda:
y1 = 11
5 y1 + y2 = 70
−2 y1 + 3 y2 + y3 = 17
Por eliminación directa de la:
primera ecuación:
y1 = 11,
segunda ecuación:
y2 = 70 − 5 y1 = 70 − 5 (11) = 15,
y de la tercera:
y3 = 17 + 2y1 − 3 y2 = 17 + 2 (11) − 3 (15) = −6.
Ahora el sistema U x = y:    
4 −2 1 11
 0 3 7  x =  15 
0 0 −2 −6
El cual escrito en su forma de ecuaciones queda:
4 x1 − 2 x2 + x3 = 11
3 x2 + 7 x3 = 15
− 2 x3 = −6
El cual al ser resuelto por sustituci´

on hacia atrás queda:
de la última ecuación:
x3 = 3,
segunda ecuación:
x2 = 5 − 7/3 x3 = 5 − 7/3 (3) = −2,
y de la primera:
x1 = 11/4 + 1/2x2 − 1/4 x3 = 11/4 + 1/2 (−2) − 1/4 (−3) = 1
C. Chapra Raymond, P. Canale.

Métodos numéricos para ingenieros.
Quinta edición; McGRAW-HILL.

Metodos Numericos

Enviado por

Direitos autorais:

Formatos disponíveis

Metodos Numericos

Enviado por

Dados do documento

Descrição original:

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Metodos Numericos

Enviado por

Direitos autorais:

Formatos disponíveis

3.4 METODO DE DESCOMPOSICION LU.

C. Chapra Raymond, P. Canale.

para despejar x del sistema:  

Este sistema escrito en su forma de ecuaciones queda:

Por eliminación directa de la:

El cual al ser resuelto por sustituci´

C. Chapra Raymond, P. Canale.

Você também pode gostar