Aplicaciones Metodos
Aplicaciones Metodos
Aplicaciones Metodos
FACTORIZACION LU
Si tenemos el sistema de ecuaciones Ax=b, podemos cambiarlo por LUx=b. Pero Ux será
también un vector columna, por lo que nos queda por un lado que Ly=b, y por otro que
Ux=y. Ambos son sistemas triangulares. Resolvemos el primero para obtener y, y a
continuación el segundo para obtener la solución buscada x
En ocasiones podemos tener varios sistemas lineales en los que la matriz del sistema A es
siempre la misma, y en los que el vector de términos independiente bk k ∈ [1,M]es
diferente para cada sistema.
Los métodos exactos para resolver sistemas numéricamente realizan del orden de N3
operaciones para resolver un sistema, donde N es el orden de la matriz A. Por tanto, para
resolver M sistemas necesitaríamos MN3 operaciones en total. El método de Doolittle
también necesita N3 operaciones para obtener los factores. Pero la resolución de sistemas
lineales solo necesita del orden de N2 operaciones.
Por tanto, si usamos la factorización de LU para obtener los factores de la matriz A, una vez
hallados, podremos resolver cada uno de los sistemas triangulares siguientes:
1. Lyk=bk
2. Uxk=yk
El cálculo de la matriz inversa es similar al caso anterior de resolver varios sistemas lineales
con la misma matriz del sistema A. En este caso, los vectores de términos independientes son
cada una de las columnas de la matriz inversa ik. Cada una de las soluciones xk obtenidas de
esta manera es la columna correspondiente de la matriz inversa de A.
MÉTODO DE JACOBI
METODO DE GAUSS-SEIDEL
AX=b
En esta ecuación:
Bibliografía
PARA EL PREZI
Resolución de un sistema lineal
La aplicación más directa de la factorización LU es la resolución de sistemas lineales. Al
transformar la matriz del sistema en un producto de matrices triangulares, el sistema de
ecuaciones puede transformarse en dos sistemas triangulares, que son muy sencillos de
resolver en un ordenador usando el algoritmo de sustitución hacia atrás.
Los métodos exactos para resolver sistemas numéricamente realizan del orden de N3
operaciones para resolver un sistema, donde N es el orden de la matriz A. Por tanto, para
resolver M sistemas necesitaríamos MN3 operaciones en total. El método de Doolittle
también necesita N3 operaciones para obtener los factores. Pero la resolución de sistemas
lineales solo necesita del orden de N2 operaciones.
En este caso, los vectores de términos independientes son cada una de las columnas de la
matriz inversa ik. Cada una de las soluciones xk obtenidas de esta manera es la columna
correspondiente de la matriz inversa de A.
MÉTODO DE JACOBI
METODO DE GAUSS-SEIDEL
AX=b