Scribd
Cargar
161 vistas6 páginas

Método de LU

El método LU consiste en descomponer una matriz en el producto de dos matrices triangulares: una matriz triangular inferior (L) y una matriz triangular superior (U). Esto permite resolver sistemas de ecuaciones lineales de manera eficiente al descomponer el problema original en dos etapas: primero resolviendo Ly=b y luego Ux=y. El método se ilustra con un ejemplo de 3 ecuaciones con 3 incógnitas.

Cargado por

Adriano Rodriguez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PPTX, PDF, TXT o lee en línea desde Scribd
161 vistas6 páginas

Método de LU

El método LU consiste en descomponer una matriz en el producto de dos matrices triangulares: una matriz triangular inferior (L) y una matriz triangular superior (U). Esto permite resolver sistemas de ecuaciones lineales de manera eficiente al descomponer el problema original en dos etapas: primero resolviendo Ly=b y luego Ux=y. El método se ilustra con un ejemplo de 3 ecuaciones con 3 incógnitas.

Cargado por

Adriano Rodriguez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PPTX, PDF, TXT o lee en línea desde Scribd

Método de LU

Métodos numéricos
 El método LU consiste en descomponer una matriz regular en el
producto de dos matrices triangulares (una superior y otra inferior)
aplicando el método de Gauss con algunas condiciones. Por esto,
es muy adecuado para aquellas situaciones donde se deben
Descomposición evaluar muchos vectores del lado derecho para un solo valor de
[A] . Aunque hay muchas formas de hacer esto, el análisis se
LU enfocará en mostrar cómo el método de eliminación de Gauss se
implementa como una descomposición LU. Un motivo para
introducir la descomposición LU es que proporciona un medio
eficiente para calcular la matriz inversa.
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:

Funcionamiento 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  Obtener la matriz triangular inferior L y la matriz triangular
superior U.
RESOLVER UN
 Resolver Ly = b (para encontrar y).
SISTEMA DE  El resultado del paso anterior se guarda en una matriz nueva de
ECUACIONES nombre "y".

POR EL  Realizar Ux = y (para encontrar x).


 El resultado del paso anterior se almacena en una matriz nueva
MÉTODO DE llamada "x", la cual brinda los valores correspondientes a las
DESCOMPOSI incógnitas de la ecuación.

CIÓN LU
EJEMPLO 1 DE DESCOMPOSICIÓN LU
Encontrar los valores de x1, x2 y x3 para el siguiente sistema
de ecuaciones:
𝑥 + 4𝑦 − 5𝑧 = −17 Se toman lo coeficientes 1 4 −5
ቐ 3𝑥 − 𝑦 + 𝑧 = 10 A= 3 −1 1
de la variables
2𝑥 + 3𝑦 − 𝑧 = −2 2 3 −1

Crea la matriz triangular superior

1 4 −5 1 4 −5
3 −1 1 𝐹2 = 𝐹2 − 3𝐹1
𝐹3 = 𝐹3 − 2𝐹1 0 −13 16
2 3 −1 0 −5 9

1 4 −5 1 4 −5
0 −13 16 0 −13 16 Se obtiene la matriz U
𝐹3 = 𝐹3 − 5/13𝐹2 U = 37
0 −5 9 0 0
13
1 0 0 𝑎 + 0𝑏 + 0𝑐 = −17
L= 3 1 0 3𝑎 + 𝑏 + 0𝑐 = 10
5 = Se crea el sistema de ecuaciones a partir de la matriz L
5
2 1 2𝑎 + 𝑏 + 𝑐 = −2
13 13
Se obtiene la matriz L
𝑎 + 0𝑏 + 0𝑐 = −17 𝑎 = −17
3𝑎 + 𝑏 + 0𝑐 = 10 𝑏 = 10 − 3𝑎 = 61
5 5 5 111
2𝑎 + 𝑏 + 𝑐 = −2 𝑐 = −2 − 2𝑎 − 𝑏 = −2 − 2(−17) − (61) =
13 13 13 13

Se genera un nuevo sistema de ecuaciones


𝑥 + 4𝑦 − 5𝑧 = −17
0𝑥 − 13𝑦 + 16𝑧 = 61 Los términos independientes son resultado del sistema de
37 111 Ecuaciones de la matriz L y los coeficientes están dados por la matriz U
𝑧=
13 13

Se genera la solución del sistema de ecuaciones

37 111 −13𝑦 + 16𝑧 = 61 𝑥 + 4𝑦 − 5𝑧 = −17


𝑧= 𝑥 + 4(−1) − 5(3) = −17
13 13 −13𝑦 + 16(3) = 61
111 −13𝑦 + 48 = 61 𝑥 − 4 − 15 = −17
𝑧= 𝑥 = −17 + 19
37 13
𝑧=3 𝑦= = −1 𝑥=2
−13

También podría gustarte