MÉTODOS

NUMÉRICOS
Docente:
Luis Giovanny Guerra, MSc.
lgguerraa@[Link]

Elaboración/Coordinador de Área
2021:
Diego Andrade, MSc.
 La matriz inversa de
una matriz
MATRIZ cuadrada, [A], que
se simboliza por [A-
INVERSA 1
], se define como la
matriz que
multiplicada por la
DE UNA matriz de la cual
[A][A-1] = [A-1][A]
proviene, el
MATRIZ resultado es
matriz unitaria, [I]
la
= [I]

CUADRADA siendo además este

producto el único
caso en el que se
. cumple la Ley
Conmutativa del
producto matricial:
MATRIZ INVERSA DE
UNA MATRIZ
CUADRADA.
• La inversa se obtiene únicamente de
una matriz cuadrada y es una
operación muy útil especialmente
cuando se resuelven sistemas de
ecuaciones lineales [A][X] = [B]
pues obtenida la inversa de la matriz
de los coeficientes, [A] la solución
del sistema se efectúa con el
siguiente procedimiento:
En el sistema [A][X] = [B] , se multiplican los
dos lados de la igualdad por [A-1] y se obtiene:

[A-1][A][X] = [A-1][B] , la igualdad quedará:

[I][X] = [A-1][B] y finalmente:

[X] = [A-1][B]

Los métodos vistos en este capítulo para resolver sistemas de

ecuaciones lineales, se pueden también utilizar para obtener la
matriz inversa de una matriz cuadrada, aplicando el siguiente
concepto:
“Si se resuelve un sistema de ecuaciones
lineales en el que el vector de términos
independientes es una de las columnas de la
matriz unitaria o idéntica, [I], la solución del
sistema es la columna correspondiente (del
mismo número), de la matriz inversa, [A-1] .”
 La demostración de este concepto es muy
Demostración sencilla: supongamos la matriz [A] y su
inversa [A-1]:
 De acuerdo a la definición de matriz
De acuerdo a inversa [A][A-1] = [I]
 Los elementos de la matriz unitaria [I] se
obtienen mediante el producto matricial de
[A][A-1] , es decir, cada elemento Ii,j de la
matriz [I] se obtiene multiplicando la fila i de
la matriz [A] por la columna j de la matriz [A-
MATRIZ 1].

UNITARI
A Apliquemos este concepto en los siguientes
ejercicios de producto matricial
 MATRIZ INVERSA
Si se multiplican las tres filas de la matriz [A] por la primera columna de la
matriz [A-1] se obtendrán los elementos de la primera columna de la matriz
unitaria [I] ; efectuemos esos productos:

• a11 x11 + a12 x21 + a13 x31 = i11 = 1

• a21 x11 + a22 x21 + a23 x31 = i21 = 0
• a31 x11 + a32 x21 + a33 x31 = i31 = 0
Se observa que se ha obtenido un sistema de ecuaciones lineales en el que
el vector de términos independientes es la primera columna de la matriz
unitaria [I] y las incógnitas que se obtienen son la primera columna de
la matriz inversa [A-1].
 MATRIZ INVERSA
De la misma manera, si se multiplican las tres filas de la matriz [A] por la
segunda columna de la matriz [A-1] se obtendrán los elementos de la
segunda columna de la matriz unitaria [I] ; efectuemos esos productos:

• a11 x12 + a12 x22 + a13 x32 = i12 = 0

• a21 x12 + a22 x22 + a23 x32 = i22 = 1
• a31 x12 + a32 x22 + a33 x32 = i32 = 0
Observamos que se ha obtenido un segundo sistema de ecuaciones lineales
en el que el vector de términos independientes es la segunda columna de
la matriz unitaria [I] y las incógnitas que se obtienen son la segunda
columna de la matriz inversa [A-1].
 MATRIZ INVERSA
Finalmente, si se multiplican las tres filas de la matriz [A] por la tercera
columna de la matriz [A-1] se obtendrán los elementos de la tercera
columna de la matriz unitaria [I] ; efectuemos esos productos:
• a11 x13 + a12 x23 + a13 x33 = i13 = 0
• a21 x13 + a22 x23 + a23 x33 = i23 = 0
• a31 x13 + a32 x23 + a33 x33 = i33 = 1
Se ha obtenido un tercer sistema de ecuaciones lineales en el que el vector
de términos independientes es la tercera columna de la matriz unitaria
[I] y las incógnitas que se obtienen son la tercera columna de la matriz
inversa [A-1].
 De esta manera se ha demostrado el concepto
enunciado anteriormente para obtener la matriz inversa

MATRIZ
de una matriz cuadrada a partir de la solución de
sistemas de ecuaciones lineales.

INVERS Al aplicar ese concepto se debe tener en cuenta que

tanto en el método de Eliminación de Gauss, así

A como en el de equivalencia de la matriz [A] en el

producto de las matrices triangulares [L] y [U] , la
primera etapa de esos métodos se realiza una sola vez,
lo que hace que no se tengan resoluciones de sistemas
de ecuaciones totalmente independientes una de otra,
existiendo un ahorro en el procesamiento numérico que
implica obtener la matriz inversa.
 Ejemplos
numéricos.
1. Obtener la
inversa de la
matriz:
 A ) Mediante
Eliminación de
Gauss
 1a Etapa: obtener • La matriz aumentada se forma
la matriz triangular añadiendo a la matriz [A] toda la
superior matriz unitaria, [I]: [A|I]
1er paso de eliminación:
Eliminación de elementos
de la primera columna.
• Fila básica: la primera,
F1
• Operaciones:
• F1 x (-5/4) + F2

• F1 x (-6/4) + F3
Realización de las operaciones:
Resultado después del primer paso
de eliminación:
2o paso de eliminación:
eliminación de elementos
de la segunda columna.
• Fila básica: la segunda,
F2
• Operación:
• F2 x (-2) + F3
Realización de las operaciones:
 Resultado después del segundo paso
de eliminación:

Nótese que se modifica solo el triángulo inferior de la

matriz untaría [I] ; el triángulo superior y la diagonal
principal de la matriz unitaria, no se modifican.
2a Fase: obtención de las
columnas de la inversa:
Se resolverán los Sistemas triangulares
superiores:
 Primera columna de la inversa:

• x31 = -1/4 = -0.25

• x21 = (-1,25 - 3.25*(-0.25)) / (-0.25) =
1.75
• x11 = (1-1.75 - 0.25) / 4 = -1/4 = -0.25

Matriz Segunda columna de la inversa:

• x32 = 2 / 4 = 0.50

inversa:
• x22 = (1- 3.25*0.50) / (-0.25) = 2.5
• x12 = (-2.5 + 0.5) / 4 = -0.5

Tercera columna de la inversa:

• x33 = -1/4 = -0.25

• x23 = -3.25*(-0.25) / (-0.25) = -3.25
• x13 = (3.25 – 0.25) / 4 = 0.75
Matriz inversa:
 b) Mediante
equivalencia
[A] = [L][U]
 1a fila de [U]
• 4, 1, -1

1a Fase: 1a columna de [L] =

• 4/4, 5/4, 6/4 = 1, 1.25, 1.5
obtenció 2a fila de [U]:
n de las • u22 = 1 – 1.25.1 = -0.25
• u23 = 2 – 1.25(-1) = 3.25
matrices 2a columna de [L]:

[L] y [U] • l32 = (1 – 1.5x1) / (-0.25) = 2

3a fila de [U]:
• u33 = 1 – (1.5(-1) + 2x3.25) = -4
planteamiento similar a Gauss
Las matrices [L] y [U] son:
 Primeros sistemas Es decir los términos independientes
[L][Y] = [I] son las tres columnas de la matriz
triangulares inferiores unitaria [I]:
 Solución para [I1] :
y3 = -1.5x1- 2*(-
y1 = 1 y2 = -1.25x1 = -1.25
1.25) = 1

Solucio y1 = 0
Solución para [I2] :
y2 = 1 y3 = -2x1 = -2

nes Solución para [I3] :

y1 = 0 ; y 2 = 0 ; y 3 = 1

Nótese que en este método y con este planteamiento se

modifica solo el triángulo inferior de la matriz unitaria
[I] , en la ejecución de los primeros sistemas de ecuaciones
triangulares inferiores.
Segundos [U][X] = [Y]
sistemas triangulares superiores
 Primera columna de la inversa, solución para
x = -1/4
[Y1]:- 3.25 (- x = (1-1.75 - 0.25) /
x = (-1,25
21 x 11
31
0.25)) / (-0.25) 4 = -1/4
x31 = -0.25
x21= 1.75 x11 = -0.25

Solucio Segunda columna de la inversa, solución

x =2/4
x para [Y0.50)
= (1- 3.25
22 2]: /
x = (-2.5 + 0.5) / 4
x

nes
32
(-0.25) 12
x32= 0.50 x12 = -0.5
x22= 2.5

Tercera columna de la inversa, solución para

[Y3]:
x23 = -3.25x(-0.25) / (-
x33 = -1/4 x13 = (3.25 – 0.25) / 4
0.25)
x33 = -0.25 x13 = 0.75
x23= -3.25
Matriz inversa:
 c) Mediante
equivalencia
[A] = [L][U]
 1a columna de
• [L] = 4, 5, 6

1a Fase: 1a fila de [U] =

• 1, 1/4, -1/4 = 1, 0.25, -0.25

obtención 2a columna de [L]:

de las • l22 = 1 – 5x0.25 = -0.25

• l32 = 1 – 6x0.25 = -0.50

matrices 2a fila de [U]:

[L] y [U] • u23 = (2 – 5(-0.25)/(-0.25) = -13.0

3a columna de [L]:

• l33 = 1 – (6(-0.25) + (-0.50)(-13.0)) = -4.0

 Planteamiento de Crout
Las matrices [L] y [U] son:
 Primeros sistemas Es decir los términos independientes
[L][Y] = [I] son las tres columnas de la matriz
triangulares inferiores unitaria [I]:
 Solución para [I1] :
y2 = 5x0.25 / 0.25 = y3 = (6x0.25 – 0.5x5) /
y1 = 1/4 = 0.25
5 4 = -0.25

Solucio y1 = 0
Solución para [I2] :
y2 = -1/0.25 = - 4 y3 = 0.5x4 / 4 = 0.5

nes Solución para [I3] :

y1 = 0 ; y2 = 0 ; y3 = -1/4 = -0.25

Nótese que en este método y con este planteamiento en la

matriz unitaria [I] no se modifica el triángulo superior, en
la ejecución de los primeros sistemas de ecuaciones
triangulares inferiores.
Segundos [U][X] = [Y]
sistemas triangulares superiores
 Primera columna de la inversa, solución para
x11 = 0.25 -
[Y ]:
x = 5 1+ 13*(- 0.25*1.75 +
x31 = -0.25 21
0.25) = 1.75 0.25*(-0.25)= -
0.25

Solucio Segunda columna de la inversa, solución

para [Y2]:
x = -4 + 13*0.5 x = -0.25*2.5 +

nes
x32 = 0.50 22 12
= 2.5 0.25*0.5= -0.5

Tercera columna de la inversa, solución para

[Y3]:
x13 = -0.25*(-
x23 = 13*(-0.25)
x33 = -0.25 3.25) + 0.25*(-
= -3.25
0.25) = 0.75
Matriz inversa:
2. Obtener la
inversa de la
matriz:
 A ) Mediante
Eliminación de
Gauss
 1a Etapa: obtener • La matriz aumentada se forma
la matriz triangular añadiendo a la matriz [A] toda la
superior matriz unitaria, [I]: [A|I]
1er paso de eliminación:
Eliminación de elementos
de la primera columna.
• Fila básica: la primera,
F1
• Operaciones:
• F1 x(-1/3) + F2
• F1 x(-2/3) + F3
• F1 x(-1/3) + F4
Realización de las operaciones:
Resultado después del primer paso
de eliminación:
2o paso de eliminación:
eliminación de elementos
de la segunda columna.
• Fila básica: la segunda,
F2
• Operación:
• F2 x(1/10) + F3
• F2 x(-1/10) + F4
Realización de las operaciones:
Resultado después del segundo paso
de eliminación:
3er paso de eliminación:
eliminación de elementos
de la tercera columna.
• Fila básica: la tercera, F3
• Operación:
• F3 x(7/27) + F4
Realización de las operaciones:
Resultado después del tercer paso de
eliminación:
2a Fase: obtención de las
columnas de la inversa:
Se resolverán los Sistemas triangulares superiores:
 Primera columna de la inversa:
• x41 = -13/44
• x31 = 10/27(-7/10 + 11/5x(-13/44)) = -1/2
• x21 = 3/10(-1/3 + 1/3x1/2 – 4/3x(-13/44)) = 3/44
• x11 = 1/3(1 – 3/22 – 1/2 + 13/22) = 7/22

Segunda columna de la inversa:

• x42 = -1/22
• x32 = 10/27(1/10 + 11/5x(-1/22)) = 0

Matriz • x22 = 3/10(1 -4/3x(-1/22)) = 7/22

• x12 = 1/3(-7/11 + 1/11) = -2/11

Tercera columna de la inversa:

inversa: • x43 = 7/44

• x33 = 10/27(1+ 11/5x7/44) = 1/2
• x23 = 3/10(-1/3x1/2 – 4/3x7/44) = -5/44
• x13 = 1/3(5/22 + 1/2 – 7/22) = 3/22

Cuarta columna de la inversa:

• x44 = 27/44
• x34 = 10/27x11/5x27/44 = 1/2
• x24 = 3/10(-1/3x1/2 -4/3x27/44) = -13/44
• x14 = 1/3(13/22 + 1/2 – 27/22) = -1/22
Matriz inversa:
 b) Mediante
equivalencia
[A] = [L][U]
 1a fila de [U]
• 3, 2, -1, 2

1a columna de [L] =
• 1, 1/3, 2/3, 1/3

1a Fase: 2a fila de [U]:

• u22 = 4 – 1/3x2 = 4 – 2/3 = 10/3

obtención
• u23 = 0 -1/3(-1) = 1/3
• u24 = 2 – 1/3x2 = 2 – 2/3 = 4/3

2a columna de [L]:

de las • l32 = (1 – 2/3x2) / 10/3 = (1 - 4/3) / 10/3 = -1/10

• l42 = (1 – 1/3x2) / 10/3 = (1 – 2/3) / 10/3 = 1/10

matrices 3a fila de [U]:

• u33 = 2 – (2/3(-1) + (-1/10)1/3) = 2 + 2/3 + 1/30 = 81/30 = 27/10

[L] y [U]
• u34 = -1 – (2/3x2 + (-1/10)4/3) = -1 – 4/3 + 2/15 = -33/15 = -11/5

3a columna de [L]:
• l43 = (-1 – (1/3(-1) + 1/10x1/3) / 27/10 = (-1 + 1/3 - 1/30) / 27/10 =
- 7/27
4a fila de [U]:
• u44 = 3 – (1/3x2 + 1/10x4/3 + (-7/27)(-11/5)) = 3 – 2/3 – 2/15 –
77/135 = 44/27
planteamiento similar a Gauss
Las matrices [L] y [U] son:
Primeros sistemas
[L][Y] = [I] Es decir los términos independientes son las
triangulares tres columnas de la matriz unitaria [I]:
inferiores
 Solución para [I1] : y4 = -1/3 – 1/10x(-
y3 = -2/3 + 1/10x(-
y1 = 1 y2 = -1/3 1/3) + 7/27x(-7/10)
1/3) = -7/10
= -13/27

Solucio y1 = 0 ; y 2 = 1
Solución para [I2] :
y3 = 1/10
y4 = -1/10 + 7/27x1/10 = -
2/27

nes Solución para [I3] :

y1 = 0 ; y2 = 0 ; y3 = 1 ; y4 = 7/27

Solución para [I4] :

y1 = 0 ; y 2 = 0 ; y 3 = 0 ; y 4 = 1
Segundos [U][X] = [Y]
sistemas triangulares superiores
 Primera columna de la inversa, solución para
x = 10/27(-7/10[Y
+ 1]:Z
x = 3/10(-1/3 +
21
x = 1/3(1 – 3/22 – 1/2
x = -13/44 31
1/3 1/2 – 4/3 (-13/44)) 11
41
11/5x(-13/44)) = -1/2 x x + 13/22) = 7/22
= 3/44

Segunda columna de la inversa, solución para

Solucio
x = 10/27(1/10 +[Y2]:x = 3/10(1 -4/3 (- x = 1/3(-7/11 + 1/11)
x = -1/22 32 22 x 12
42
11/5x(-1/22)) = 0 1/22)) = 7/22 = -2/11

nes Tercera columna de la inversa, solución para [Y3]:

x43 = 7/44
x33 = 10/27(1+
11/5x7/44) = 1/2
x23 = 3/10(-1/3x1/2 –
4/3x7/44) = -5/44
x13 = 1/3(5/22 + 1/2 –
7/22) = 3/22

Cuarta columna de la inversa, solución para [Y4]:

x34 = 10/27x11/5x27/44 x24 = 3/10(-1/3x1/2 - x14 = 1/3(13/22 + 1/2 –
x44 = 27/44
= 1/2 4/3x27/44) = -13/44 27/22) = -1/22
Matriz inversa:
 c) Mediante
equivalencia
[A] = [L][U]
 1a columna de [L] =
• 3, 1, 2, 1

1a fila de [U] =
• 1, 2/3, -1/3, 2/3

1a Fase: 2a columna de [L]:

• l22 = 4 – 1x2/3 = 10/3

obtención • l32 = 1 – 2/3x2 = -1/3

• l42 = 1 – 1/3x2 = 1/3

de las
2a fila de [U]:

• u23 = (0 -1/3(-1))/(10/3) = 1/10

• u24 = (2 – 1/3x2)/(10/3) = 2/5

matrices 3a columna de [L]:

• l33 = 2 – (2(-1/3) + (-1/3)(1/10)) = 2 + 7/10 = 27/10

[L] y [U]
• l43 = -1 – (1(-1/3) + 1/10 x1/3) = - 1 + 3/10 = - 7/10

3a fila de [U]:

• u34 = (-1 – (2x2/3 + (-1/3)(2/5))) / (27/10) = (-1 – 4/3 + 2/15) / (27/10) = -

22/27
4a columna de [L]:

• l44 = 3 – (1x2/3 + (1/3)(2/5)+(-7/10)(-22/27)) = 3 - 2/3 - 2/15 – 77/135 =

44/27
 Planteamiento de Crout
Las matrices [L] y [U] son:
Primeros sistemas
[L][Y] = [I] Es decir los términos independientes son las
triangulares tres columnas de la matriz unitaria [I]:
inferiores
 Solución para [I1] : y4 = 27/44(-1/3 +
y3 = 10/27(-2/3 –
y1 = 1/3 y2 = -3/10x1/3 = -1/10 1/3x1/10 – 7/10x7/27)
1/3x1/10) = -7/27
= -13/44

Solucio y1 = 0 ; y2 = 3/10
Solución para [I2] :
y3 = 10/27x1/3x3/10 = 1/27
y4 = 27/44(-1/3x3/10 +
7/10x1/27) = -1/22

nes Solución para [I3] :

y1 = 0 ; y2 = 0 ; y3 = 10/27 y4 = 27/44x7/10x10/27 = 7/44

Solución para [I4] :

y1 = 0 ; y2 = 0 ; y3 = 0 ; y4 = 27/44
Segundos [U][X] = [Y]
sistemas triangulares superiores
 Primera columna de la inversa, soluciónx para [Y1]:-
= 1/3 -2/3 3/44
x31 = -7/27 - x21 = -1/10 + 1/10x1/2 11 x
x41 = -13/44 1/3x1/2 + 2/3x13/44 =
22/27x13/44 = -1/2 + 2/5x13/44 = 3/44
7/22

Segunda columna de la inversa, solución para

Solucio x = -1/22 32
[Y2x]: = 3/10 + 2/5 1/22 =
x = 1/27 - 22/27 1/22 x 22
x = -2/3 7/22 +
x 12 x
42
=0 7/22 2/3x1/22 = -2/11

nes Tercera columna de la inversa, solución para

x43 = 7/44
x33 = 10/27 +
22/27x7/44 = 1/2
[Y ]:
x = 2/3 5/443+
x23 = -1/10x1/2 –
2/5x7/44 = -5/44
13 x
1/3x1/2 – 2/3x7/44 =
3/22

Cuarta columna de la inversa, solución para [Y4]:

x14 = 2/3x13/44 +
x34 = 22/27x27/44 = x24 = -1/10x1/2 –
x44 = 27/44 1/3x1/2 - 2/3x27/44 = -
1/2 2/5x27/44 = -13/44
1/22
Matriz inversa: