MN06 - Sistemas de Ecuaciones Lineales - LU

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 19

CM3207 - METODOS NUMERICOS PARA INGENIERIA

SISTEMAS DE ECUACIONES LINEALES – FACTORIZACION LU


Factorización LU
• Tenemos que ; donde y son vectores
Si existen dos matrices y de modo que ,
donde es triangular superior y es triangular
inferior:
• Tomamos un vector de forma que:
y
• Entonces: y

• Como
• Queremos factorizar la siguiente matriz:

• Si hacemos la multiplicación de obtenemos:

• Se obtiene un sistema de ecuaciones indeterminado (6


incógnitas, sólo 4 ecuaciones), por lo tanto se debe
considerar alguna restricción adicional para poder resolver.
• Se pueden asumir dos casos:
(1) Elementos diagonales de como 1 Descomposición
Doolittle
(2) Elementos diagonales de como 1 Desc. Crout
• Por lo general (y para efectos del curso) se toma
Doolittle y de ésta manera:

• Por lo tanto: NOTA: Ojo que el 6/4 es el factor


necesario en EG para eliminar el
6 de la matriz original, y de ahí
se obtiene la matriz U
Ejemplo
• Utilice Factorización LU para determinar la matriz
inversa de:
Solución
• Tenemos que:

• Podríamos decir que:

• Y de manera similar para cada columna de .


Con factorización LU podemos obtener las matrices L
y U para después resolver todos los vectores con sólo
realizar sustitución hacia adelante y hacia atrás.
• Aplicamos EG en orden • Almacenamos el factor
para crear matriz “U” en matriz “L”
• Luego procedemos a solucionar el problema
utilizando la ecuaciones:

;
• Resolvemos con sustitución hacia adelante:

;
Para la 1ra columna de

Resolvemos con sustitución hacia atrás:

Por lo tanto, la primera columna de :


• De similar manera, para las demás columnas de
, se repiten los pasos de sust. hacia adelante y
hacia atrás, cambiando únicamente los vectores
“ ” y “ ” según corresponda, y se obtienen las
siguientes soluciones:
• 2da columna de :

• 3ra columna de :
Ejemplo
• Utilice Factorización LU para determinar la matriz
inversa de:
Solución
Ventajas y desventajas de EG y Factorización LU
• Ambos métodos obtienen la solución exacta
• Si se cambia el vector , EG debe repetirse, Fact LU
permiten hacer la EG en la descomposición LU, y
luego solo se deben repetir los pasos de sustitución
que son más sencillos.
• Esto facilita hallar la matriz inversa y luego se
puede variar el vector de manera sencilla y sólo
multiplicarlo por la inversa para obtener la solución
Sobre la matriz inversa
• Estímulo y respuesta del sistema

• se considera una constante de


proporcionalidad del cambio que sufrirá debido
a un cambio en
Relajación
• Al utilizar métodos iterativos, la misma diferencia entre
una aproximación y otra nos da una indicación de hacia
adonde se está moviendo la respuesta.
• Por ejemplo, si aprox1 = 4 y aprox2 = 5, ¿cuál creen
ustedes que sería el valor de aprox3?

• Como se puede anticipar hacia donde va el resultado,


se puede utilizar el método de relajación (normalmente
SOBRE-relajación) para acelerar la respuesta:
• De esta manera, si :
0

• Pero si :

• En nuestro caso y
:
• Se debe aplicar con precaución pues en algún
momento una solución que venía moviéndose en
una dirección puede cambiar de dirección y aplicar
relajación ocasionaría un mayor error o incluso
problemas de convergencia:

También podría gustarte