4.2.4 Matriz inversa.

DESCOMPOSICIÓN LU Es un método directo para resolver sistemas de ecuaciones de la forma [A] {X}= {B} El principal recurso es que el paso de la eliminación, que toma mucho tiempo, se puede formular de tal manera que solo involucre operaciones con la matriz de coeficientes [A]. Muestra cómo el método de eliminación de Gauss se implementa como una descomposición LU. Una ventaja de este método es que proporciona un medio eficiente para evaluar diversos vectores del lado derecho.

Pasos en la descomposición LU

PROCEDIMIENTO Paso de descomposición LU: [A] se factoriza o descompone en las matrices triangular inferior [L] y superior [U]. Paso de sustitución: [L] y [U] determinan una solución {X} para un lado derecho {B}. De acuerdo a la Ec. [L]{D}={B} genera un vector {D} mediante sustitución hacia delante. El resultado se sustituye en la Ec. [U]{X}={D} mediante sustitución hacia atrás para X.

PROCEDIMIENTO(continuación)

PROCEDIMIENTO(continuación)

PROCEDIMIENTO (continuación)

Ejemplo (Descomposición LU)

EJEMPLO (continuación)

EJEMPLO (continuación)

DESCOMPOSICIÓN CROUT Utiliza una matriz [U] con números 1 sobre la diagonal. Genera [U] y [L] barriendo las columnas y los renglones de la matriz.

FÓRMULAS PARA DESCOMPOSICIÓN CROUT

MATRIZ INVERSA Si una matriz [A] es cuadrada, existe otra matriz [A]-1, conocida como la inversa de [A], en donde: [A][A]-1=[A]-1[A]=[I] En esta sección se estudiará como calcular la inversa por medio del algoritmo de descomposición LU.

PROCEDIMIENTO Se efectúa sustitución hacia delante para calcular la primera columna de la matriz inversa, utilizando un vector unitario (el número 1 en el primer reglón) como el vector del lado derecho.

PROCEDIMIENTO (continuación) De lo anterior se obtiene {D}T, este vector se utilizará como el lado derecho en la siguiente ecuación: Donde se obtiene por sustitución hacia atrás {X}T, estos valores serán la primera columna de la matriz inversa.

PROCEDIMIENTO (continuación) Para calcular la segunda columna se formula la siguiente ecuación: De esta forma se obtiene {D}, y los resultados se usan para calcular {X}T, que es la segunda columna de la inversa.

PROCEDIMIENTO (continuación) Para el cálculo de la tercera columna se realiza de igual manera que las anteriores empleando {B}T=[O,O,1], para obtener {X}T. Finalmente la validez del resultado se comprueba al verificar que [A][A]-1=[I]

No hay comentarios.:
Publicar un comentario