Inversa de matrices
ANÁLISIS NUMÉRICO
Descomposición LU
• Un motivo para introducir la descomposición LU es que
proporciona un medio eficiente para calcular la matriz inversa.
• La inversa tiene muchas aplicaciones valiosas en la práctica de la
ingeniería.
• Ésta ofrece también un medio para evaluar la condición de un
sistema.
Descomposición LU
• Como se describió, la eliminación de Gauss sirve para resolver
sistemas de ecuaciones algebraicas lineales,
[A]{X} = {B}
• Aunque la eliminación Gauss representa una forma satisfactoria
para resolver tales sistemas, resulta ineficiente cuando deben
resolverse ecuaciones con los mismos coeficientes [A], pero con
diferentes constantes del lado derecho (las b).
Descomposición LU
• Recuerde que la eliminación de Gauss implica dos pasos:
eliminación hacia adelante y sustitución hacia atrás.
• De ambas, el paso de eliminación hacia adelante es el que
representa la mayor parte del trabajo computacional.
• Esto es particularmente cierto para grandes sistemas de
ecuaciones.
Descomposición LU
• Los métodos de descomposición LU separan el tiempo usado en las
eliminaciones para la matriz [A] de las manipulaciones en el lado
derecho {B}.
• Una vez que [A] se ha “descompuesto”, los múltiples vectores del
lado derecho {B} se pueden evaluar de manera eficiente.
Descomposición LU
• El hecho de que la misma eliminación de Gauss se puede expresar
como una descomposición LU es muy interesante.
• Antes de mostrar cómo se puede realizar esto, demos primero una
demostración matemática de la estrategia de descomposición.
• De manera similar al caso de la eliminación de Gauss, la
descomposición LU requiere de pivoteo para evitar la división
entre cero.
Descomposición LU
• Sin embargo, para simplificar la siguiente descripción,
abordaremos el tema del pivoteo después de que el planteamiento
fundamental se haya elaborado.
• Además, la siguiente explicación se limita a un conjunto de tres
ecuaciones simultáneas.
• Los resultados se pueden extender en forma directa a sistemas n
dimensionales.
Descomposición LU
• La ecuación [A]{X} = {B} se reordena como:
[A] {X} – {B} = 0
• Suponga que la ecuación puede expresarse como un sistema
triangular superior:
Descomposición LU
• Observe que esto es similar a la manipulación que ocurre en el
primer paso de la eliminación de Gauss.
• Es decir, se utiliza la eliminación para reducir el sistema a una
forma triangular superior.
• La ecuación [A] {X} – {B} = 0 también se expresa en notación
matricial y se reordena como:
[U]{X} – {D} = 0
Descomposición LU
• Ahora, suponga que existe una matriz diagonal inferior con números 1
en la diagonal,
que tiene la propiedad de que cuando se premultiplica por la ecuación
[U]{X} – {D} = 0, el resultado es la ecuación [A] {X} – {B} = 0 . Es decir:
Descomposición LU
• [L]{[U]{X} – {D}} = [A]{X} – {B}
• Si esta ecuación se satisface, según las reglas de multiplicación entre
matrices, se obtendrá
[L][U] = [A]
[L]{D} = {B}
Descomposición LU
• Una estrategia de dos pasos para obtener soluciones se basa en las ecuaciones:
[U]{X} – {D} = 0,
[L][U] = [A]
[L]{D} = {B}
1. Paso de descomposición LU. [A] se factoriza o “descompone” en las matrices
triangulares inferior [L] y superior [U].
2. Paso de la sustitución. [L] y [U] se usan para determinar una solución {X} para un
lado derecho {B}. Este paso, a su vez, se divide en dos. Primero, la ecuación[L]{D} =
{B}; se usa para generar un vector intermedio {D} mediante sustitución hacia adelante.
Después, el resultado se sustituye en la ecuación [U]{X} – {D} = 0, la que se resuelve
por sustitución hacia atrás para {X}.
Descomposición LU
• Una estrategia de dos pasos para obtener soluciones se basa en las ecuaciones:
[U]{X} – {D} = 0,
[L][U] = [A]
[L]{D} = {B}
1. Paso de descomposición LU. [A] se factoriza o “descompone” en las matrices
triangulares inferior [L] y superior [U].
2. Paso de la sustitución. [L] y [U] se usan para determinar una solución {X} para un
lado derecho {B}. Este paso, a su vez, se divide en dos. Primero, la ecuación[L]{D} =
{B}; se usa para generar un vector intermedio {D} mediante sustitución hacia adelante.
Después, el resultado se sustituye en la ecuación [U]{X} – {D} = 0, la que se resuelve
por sustitución hacia atrás para {X}.
Descomposición LU
• Aunque a primera vista podría parecer que la eliminación de Gauss no
está relacionada con la eliminación LU, aquélla puede usarse para
descomponer [A] en [L] y [U], lo cual se observa fácilmente para [U],
que es el resultado directo de la eliminación hacia adelante.
• Recuerde que en el paso correspondiente a esta eliminación se
pretende reducir la matriz de coeficientes [A] a la forma:
que es el formato triangular superior deseado
Descomposición LU
• De hecho, esta matriz representa un almacenamiento eficiente de la
descomposición LU de [A].
[A] → [L][U], donde:
Descomposición LU. Ejemplo
• Planteamiento del problema. Obtenga una descomposición LU que
resuelva el siguiente sistema de ecuaciones:
3x1 – 0.1x2 – 0.2x3 = 7.85
0.1x1 + 7x2 – 0.3x3 = –19.3
0.3x1 – 0.2x2 + 10x3 = 71.4
Descomposición LU. Ejemplo
• Solución:
• La matriz de coeficientes [A] es:
Descomposición LU. Ejemplo
• Solución:
• Después de la eliminación hacia adelante, se obtuvo la siguiente
matriz triangular superior:
Descomposición LU. Ejemplo
• Solución:
• Los factores empleados para obtener la matriz triangular superior se
pueden colocar en una matriz triangular inferior.
• Los elementos a21 y a31 se eliminaron al usar los factores:
Descomposición LU. Ejemplo
• Solución:
• Así, la matriz triangular inferior es:
Descomposición LU. Ejemplo
• Solución:
• En consecuencia, la descomposición LU es:
Descomposición LU. Ejemplo
• Solución:
• Este resultado se verifica al realizar la multiplicación de [L][U] que
da:
Descomposición LU. Ejemplo
• Solución:
• Este resultado se verifica al realizar la multiplicación de [L][U] que da:
donde las pequeñas diferencias son debidas a errores de redondeo.
Descomposición LU. Ejemplo
• Solución:
• Como se estableció antes, la intención de la sustitución hacia
adelante es aplicar las operaciones de eliminación al vector {B},
previamente aplicadas a [A].
• El sistema a resolver es:
Descomposición LU. Ejemplo
• Solución:
• La fase de la sustitución hacia adelante se realiza aplicando la
ecuación [L][D] = [A] a nuestro problema
Descomposición LU. Ejemplo
• Solución:
o realizando la multiplicación entre matrices del lado izquierdo e
igualando,
d1 = 7.85
0.0333333d1 + d2 = –19.3
0.1d1 – 0.02713d2 + d3 = 71.4
Descomposición LU. Ejemplo
• Solución:
• Se resuelve la primera ecuación para d1,
d1 = 7.85
• la cual se sustituye en la segunda ecuación y se resuelve para d2
d2 = –19.3 – 0.0333333(7.85) = –19.5617
• Ambas, d1 y d2, se sustituyen en la tercera ecuación para d3
d3 = 71.4 – 0.1(7.85) + 0.02713(–19.5617) = 70.0843
Descomposición LU. Ejemplo
• Solución:
• Así la matriz [D] es:
Descomposición LU. Ejemplo
• Solución:
• Este resultado se sustituye, entonces, en la ecuación, [U]{X} = {D},
para obtener:
Descomposición LU. Ejemplo
• Solución:
• que se resuelve por sustitución hacia atrás para obtener la solución
final:
Descomposición LU. Seudocódigo
Matriz Inversa
• En el estudio de las operaciones con matrices, vimos que si una
matriz [A] es cuadrada, existe otra matriz [A]–1, conocida como la
inversa de [A], para la cual:
[A]= [A] = [I]
• Ahora se enfocará el análisis hacia el modo en que la matriz
inversa se calcula numéricamente.
Matriz Inversa
• La inversa se puede calcular en forma de columna por columna,
generando soluciones con vectores unitarios como las constantes
del lado derecho.
• Por ejemplo, si la constante del lado derecho de la ecuación
tienen un número 1 en la primera posición, y ceros en las otras,
Matriz Inversa
• la solución resultante será la primera columna de la matriz
inversa.
• En forma similar, si se emplea un vector unitario que tiene un
número 1 en el segundo renglón
Matriz Inversa
• la solución resultante será la primera columna de la matriz inversa.
• En forma similar, si se emplea un vector unitario que tiene un
número 1 en el segundo renglón:
• el resultado será la segunda columna de la matriz inversa.
Matriz Inversa
• La mejor forma de realizar un cálculo como éste es con el
algoritmo de descomposición LU.
• Recuerde que una de las ventajas más importantes de la
descomposición LU es que proporciona un medio eficiente para
evaluar diversos vectores del lado derecho.
• Por lo tanto, resulta ideal para evaluar los vectores unitarios
requeridos en el cálculo de la inversa.
Matriz Inversa. Ejemplo
Planteamiento del problema. la descomposición LU para
determinar la matriz inversa de:
Matriz Inversa
Solución:
• Recuerde que la descomposición dio como resultado las siguientes
matrices triangulares inferior y superior:
Matriz Inversa
Solución:
• La primera columna de la matriz inversa puede determinarse al
efectuar el procedimiento de solución por sustitución hacia
adelante, con un vector unitario (con el número 1 en el primer
renglón) como el vector del lado derecho.
• Así, de la ecuación [L]{D} = {B}, el sistema diagonal inferior es:
Matriz Inversa
Solución:
• de donde, por sustitución hacia adelante se obtiene:
= [1 –0.03333 –0.1009].
Este vector se utiliza como el lado derecho de la ecuación [U]{X} = {B},
Matriz Inversa
Solución:
• de donde, por sustitución hacia atrás, se obtiene:
= [0.33249 –0.00518 –0.01008],
• que es la primera columna de la matriz,
Matriz Inversa
Solución:
• Para determinar la segunda columna, la ecuación [L]{D]={B], se formula
como:
Matriz Inversa
Solución:
• De donde se puede obtener {D}, y los resultados se usan con la ecuación
[U]{X}={B} para determinar
= [0.0049440.1429030.00271], que es la segunda columna de la matriz,
Matriz Inversa
Solución:
• Por último, los procedimientos de sustitución hacia adelante y de
sustitución hacia atrás pueden usarse con = ⎣0 0 1⎦, para obtener =
[0.006798 0.004183 0.09988], que es la columna final de la matriz:
• La validez de este resultado se comprueba al verificar que [A] = [I].
Matriz Inversa. Taller
Planteamiento del problema:
a) Resuelva el sistema de ecuaciones siguiente por medio de la
descomposición LU.
8x1 + 4x2 – x3 = 11
–2x1 + 5x2 + x3 = 4
2x1 – x2 + 6x3 = 7
b) Determine la matriz inversa. Compruebe sus resultados por medio de
verificar que [A] = [I].
