Metodos Matriz Inversa PDF

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

Métodos numéricos para calcular la inversa de una matriz

Estudiante:
Edgar David Hernández Medina
T00054326

Profesor:
Juan Carlos Mantilla Gómez

Universidad Tecnológica de Bolívar


Facultad de Ingeniería
Ingeniería de Sistemas y Computación
Procesamiento numérico
Grupo G
NRC 1642
Cartagena
13 de septiembre de 2020
Matriz inversa
Definición:

En matemáticas, en particular en álgebra lineal, una matriz cuadrada A de orden n se dice


que es invertible, no singular, no degenerada o regular si existe otra matriz cuadrada de
orden n, llamada matriz inversa de 𝐴 y representada como 𝐴−1 , tal que: 𝐴 ∙ 𝐴−1 = 𝐴−1 ∙
𝐴 = 𝐼𝑛 ,donde 𝐼𝑛 es la matriz identidad de orden 𝑛 y el producto utilizado es el producto
de matrices usual.1

Métodos más utilizados para calcularla


1) Método de cofactores:
La matriz inversa puede ser calculada por el método de cofactores cuando necesariamente
el determinante de la matriz sea un valor diferente de 0, cuando esta regla se cumple el
método puede ser aplicado sin ningún problema, el método tiene por definición de matriz
inversa la siguiente fórmula:

𝑎𝑑𝑗(𝐴)𝑡 −1
𝐴 =
𝑑𝑒𝑡(𝐴)
𝑡
Donde 𝑎𝑑𝑗(𝐴) es la matriz de los adjuntos y 𝑑𝑒𝑡(𝐴) es la determinante de la matriz A.

Aplicando método de cofactores para el ejercicio propuesto:

3,027 −2,145 3,336 𝑥1 13,172316


[2,598 3,995 −4,215] ∙ [𝑥2 ] = [−4,953015]
1,005 −2,345 3,874 𝑥3 12,817794

En el ejercicio se requiere resolver un sistema matricial 3x3 donde la matriz de la izquierda es la


matriz de coeficientes (A), la central es la matriz de incógnitas (x) y la de la derecha es la matriz
de términos independientes (b), la solución que se puede dar a este sistema equivale a la
multiplicación de la matriz inversa de la matriz de coeficientes por la matriz de términos
independientes.

De la siguiente manera:

𝑥 = 𝐴−1 ∙ 𝑏

1
https://es.wikipedia.org/wiki/Matriz_invertible
Solución
Primer paso: Se calcula el determinante de la matriz de coeficientes para reconocer si es
posible resolver el sistema por este método.

3,027 −2,145 3,336


𝐴 = [2,598 3,995 −4,215]
1,005 −2,345 3,874
Se procede al calculo del determinante de la matriz por medio de la regla de Sarrus:
det(𝐴) = (3,027)(3,995)(3,874) + (−2,145)(−4,215)(1,005) + (2,598)(−2,345)(3,336)
− (3,336)(3,995)(1,005) − (−2,145)(2,598)(3,874)
− (−4,215)(−2,345)(3,027)
det(𝐴) = 13,88551794
El resultado del determinante resultó un valor diferente a cero así que se puede seguir con la
aplicación del método, como ya está comprobado que existe una matriz inversa se procede con el
calculo de la matriz de cofactores de la matriz de coeficientes. Esta se calcula por medio de un
mecanismo simple, una matriz de cofactores resulta de omitir la fila y la columna de un elemento
de la matriz, esto se repite en cada elemento de la matriz, el signo que tenga el resultado depende
de si la suma de la posición correspondiente, si la suma es par el resultado queda con igual signo,
si es impar el signo debe cambiar por lo que se agrega el signo negativo al resultado, la matriz de
cofactores es del mismo orden de la matriz de coeficientes.

Procedimiento para uno de los elementos:

Elemento 𝐴11 = 3,027


Fila = 1 Columna = 1 Suma = 1+1 = 2 Signo = +
3,027 −2,145 3,336
[2,598 3,995 −4,215]
1,005 −2,345 3,874
Matriz:
3,995 −4,215
[ ]
−2,345 3,874

det(𝐴11 ) = (3,995)(3,874) − (−4,215)(−2,345)


det(𝐴11 ) = 5,592455
Para este elemento se hizo el calculo de su determinante, se hizo la suma del número de fila el
numero de columna y se comprobó que era par por lo que el signo no cambia, el resultado del
determinante corresponde al elemento de la misma posición en la matriz de cofactores.
Para completar la matriz de cofactores se repite el proceso para los demás elementos de la
matriz de coeficientes, lo cual da como resultado la siguiente matriz:
5,592455 −14,300727 −10,107285
𝑐𝑜𝑓(𝐴) = [ 0,48681 8,373918 4,94259 ]
−4,286145 21,425733 17,665575

Para completar el método se debe calcular la matriz adjunta, por lo que se debe hacer la transpuesta
a la matriz de cofactores, una matriz transpuesta consiste en intercambiar las filas y columnas de
una matriz, al realizar este procedimiento la matriz resultante es:

5,592455 0,48681 −4,286145


𝑡
𝑎𝑑𝑗(𝐴) = [−14,300727 8,373918 21,425733 ]
−10,107285 4,94259 17,665575

Teniendo la adjunta ya es posible hacer el cálculo de la matriz inversa, que se hace dividiendo la
adjunta por el determinante de la matriz de coeficientes inicial.
𝑎𝑑𝑗(𝐴)𝑡
𝐴−1 =
det⁡(A)

1 5,592455 0,48681 −4,286145


𝐴−1 = [−14,300727 8,373918 21,425733 ]
13.88551794 −10,107285 4,94259 17,665575
0,402754 0,035058 −0,308677
−1
𝐴 = ⁡ [−1.029902 0,603068 1,543027 ]
−0.727901 0,355952 1,272230
Ahora con la matriz inversa se procede a hallar la solución del sistema de ecuaciones, esto se hace
multiplicando la inversa por la matriz de términos independientes:
𝑥 = 𝐴−1 ∙ 𝑏
𝑥1 0,402754 0,035058 −0,308677 13,172316
[𝑥2 ] = [−1.029902 0,603068 1,543027 ] ∙ [−4,953015]
𝑥3 −0.727901 0,355952 1,272230 12,817794
Esta multiplicación arroja como matriz resultante:
𝒙𝟏 𝟏, 𝟏𝟕𝟓
𝒙
[ 𝟐 ] = [𝟑, 𝟐𝟐𝟓]
𝒙𝟑 𝟒, 𝟗𝟓𝟔
Estos son los valores que le dan solución al sistema de ecuaciones.
2) Matriz inversa por definición, por matriz identidad.

En la definición de una matriz inversa es explicado que existe una matriz que al multiplicarla por
la matriz a la cual se le quiere sacar la matriz inversa el resultado será una matriz conocida como
matriz identidad. Esta definición se puede representar por medio de un ejemplo.

Se tiene la siguiente matriz:


1 2
𝐴=[ ]
3 4
La definición dice que existe una matriz tal que:

1 2 𝑎 𝑏 1 0
[ ]∙[ ]⁡=⁡[ ]
3 4 𝑐 𝑑 0 1

Con esta definición se tiene que lo que se debe hacer es aplicar es un producto matricial entre la
matriz a la cual se le quiere calcular la inversa y otra matriz del mismo orden pero con valores
incógnitos en sus elementos, después de realizar esta multiplicación se debe igualar al número que
resultará por la definición, en el ejemplo si se toma el primer elemento de la matriz identidad se
denota que el valor es 1 por lo que este equivale a la suma de los productos entre la primera
columna de matriz inicial A por la primera fila de matriz de incógnitas B, lo cual resulta en una
ecuación lineal con dos incógnitas. Luego de realizar la multiplicación entre las dos matrices y la
respectiva asociación con el valor correspondiente en la matriz identidad y tener las ecuaciones se
pueden resolver por cualquier método de preferencia para hallar los valores de las incógnitas y
completar la matriz inversa.
Con el ejercicio planteado, la resolución por este método sería:

3,027 −2,145 3,336 𝑎 𝑏 𝑐 1 0 0


[2,598 3,995 −4,215] ∙ [𝑑 𝑒 𝑓 ] = [0 1 0]
1,005 −2,345 3,874 𝑔 ℎ 𝑖 0 0 1

Realizando el producto entre las dos matrices y la asociación las ecuaciones quedan tal que:

3,027𝑎 − 2,145𝑑 + 3,336𝑔 = 1


3,027𝑏 − 2,145𝑒 + 3,336ℎ = 0
3,027𝑐 − 2,145𝑓 + 3,336𝑖 = 0

2,598𝑎 + 3,995𝑑 − 4,215𝑔 = 0


2,598𝑎 + 3,995𝑒 − 4,215ℎ = 1
2,598𝑐 + 3,995𝑓 − 4,215𝑖 = 0

1.005𝑎 − 2.345𝑑 + 3.874𝑔 = 0


1.005𝑏 − 2.345𝑒 + 3.874ℎ = 0
1.005𝑐 − 2.345𝑓 + 3.874𝑖 = 1
En las ecuaciones se puede encontrar una correspondencia entre las incógnitas, así como en la
primera ecuación de grupo de ecuaciones tienen el mismo orden secuencial en sus términos 𝑥𝑎 +
𝑦𝑑 + 𝑧𝑔 = 𝑇𝑖, estas ecuaciones pueden asociarse y tener un sistema de ecuaciones lineales de 3x3
que puede ser resuelto para encontrar cada termino de la primera columna de la matriz inversa.
Representando estas ecuaciones de manera matricial se tiene que:

a d g TI
3,027 -2,145 3,336 1
2,598 3,995 -4,215 0
1,005 -2,345 3,874 0

Este sistema matricial puede ser resuelto por medio de la Regla de Cramer, por lo que se calcula
primero la determinante del sistema.

det(S) = 13,8855179

Luego de esto se intercambia la columna de la variable que se quiere resolver por la columna
correspondiente en la matriz identidad y se halla el determinante:

TI d g
1 -2,145 3,336
0 3,995 -4,215 det(a)= 5,592455
0 -2,345 3,874
a TI g

3,027 1 3,336
2,598 0 -4,215 det(b)= -14,300727
1,005 0 3,874
a d TI

3,027 -2,145 1
2,598 3,995 0 det(c)= -10,107285
1,005 -2,345 0

Ahora teniendo los determinantes de cada matriz se puede hallar el valor de cada una de las
incógnitas, las cuales corresponden esos elementos en la matriz de incógnitas. Estos valores son
calculados dividiendo el determinante de cada variable por el determinante del sistema, dando
como resultados:
a= 0,40275451
d= -1,02990231
g= -0,72790119
Luego se aplica el procedimiento idéntico para las demás ecuaciones y se tienen los siguientes
resultados:

b= 0,03505883
e= 0,60306847
h= 0,35595287

c= -0,30867736
f= 1,54302728
i= 1,27223018

A continuación, se reemplazan estos valores en la matriz de incógnitas planteada inicialmente:

𝑎 𝑏 𝑐 0,40275451 0,03505883 −0,30867736


[𝑑 𝑒 𝑓] = [−1,02990231 0,60306847 1,54302728 ]
𝑔 ℎ 𝑖 −0,72790119 0,35595287 1,27223018

Por lo tanto, ya está calculada la matriz inversa de la matriz de coeficientes, ahora el ultimo paso
es multiplicar esta matriz por la matriz de términos independientes para obtener los valores de las
variables que son solución del sistema de ecuaciones:

𝑥1 0,40275451 0,03505883 −0,30867736 13,172316


[𝑥2 ] = [−1,02990231 0,60306847 1,54302728 ] ∙ [−4,953015]
𝑥3 −0,72790119 0,35595287 1,27223018 12,817794

Al realizar la multiplicación entre estas dos matrices da una matriz resultante:


𝒙𝟏 𝟏, 𝟏𝟕𝟓
𝒙
[ 𝟐 ] = [𝟑, 𝟐𝟐𝟓]
𝒙𝟑 𝟒, 𝟗𝟓𝟔

También podría gustarte