Matrices invertibles

Se dice que una matriz cuadrada A es invertible, si existe una matriz B con la propiedad de que

AB = BA = I

Siendo I la matriz identidad. Denominamos a la matriz B la inversa de A y la denotamos por A⁻¹

Ejemplo:

Supongamos A = 25 
13
y B = 3-5 
-12

Entonces:

A·B = 25 · 3-5 = 6 - 5 -10 + 10 = 10 = I
13-123 - 3 -5 + 601
A·B = 3-5 · 25 = 6 - 5 15 - 15 = 10 = I
-1213-2 + 2 -5 + 601

Puesto que AB = BA = I, A y B son invertibles, siendo cada una la inversa de la otra.

Método de Gauss

Sea A = (aₓ) una matriz cuadrada de orden n. Para calcular la matriz inversa de A, que denotaremos como A⁻¹, seguiremos los siguientes pasos:

Paso 1. Construir la matriz n×2·n M = (A ⋮ I) esto es, A está en la mitad izquierda de M y la matriz identidad I en la derecha.

Paso 2. Se deja tal y como está la primera fila de M, y debajo del primer término de la diagonal principal, a₁₁, que llamaremos pivote, ponemos ceros. Luego se opera como se indica en el siguiente ejemplo.

Ejemplo:

Consideremos una matriz 3×3 arbitraria.

A = a₁₁a₁₂a₁₃ 
a₂₁a₂₂a₂₃
a₃₁a₃₂a₃₃

Paso 1

M = (A ⋮ I) = a₁₁a₁₂a₁₃100 ~
a₂₁a₂₂a₂₃010
a₃₁a₃₂a₃₃001

Paso 2

 a₁₁ a₁₂ a₁₃1 0 0 
0 a₁₁·a₂₂ - a₂₁·a₁₂ a₁₁·a₂₃ - a₂₁·a₁₃a₁₁·0 - a₂₁·1 a₁₁·1 - a₂₁·0 a₁₁·0 - a₂₁·0
0 a₁₁·a₃₂ - a₃₁·a₁₂ a₁₁·a₃₃ - a₃₁·a₁₃a₁₁·1 - a₃₁·1 a₁₁·0 - a₃₁·0 a₁₁·1 - a₃₁·0

El siguiente paso es igual que el anterior, pero esta vez se coge como pivote el segundo término de la diagonal principal. Al llegar al último término de la diagonal, se procede igual que antes, pero poniendo los ceros encima del nuevo pivote. Se observa que al coger como pivote el último término de la diagonal, la matriz A se transforma en una matriz triangular.

Una vez realizados todos los pasos, la mitad izquierda de la matriz M se convierte en una matriz diagonal. En este momento hay que proceder a transformar, si es que no lo está, la mitad izquierda en la matriz identidad, dividiendo si fuera necesario las filas de M por un escalar.

Ejemplo: Supongamos que queremos encontrar la inversa de

A = 102 
2-13
418

Primero construimos la matriz M = (A ⋮ I),

M = 102100 
2-13010
418001
~ 1 0 21 0 0 ~
0 -1 - 2,0 3 - 2,20 - 2 1 - 2,0 0
0 1 - 4,0 8 - 4,20 - 4 0 1 - 0
~ 102100 
0-1-1-210
010-401

Luego se toma como pivote a₂₂ = -1,

~ 1 0 21 0 0 ~
0 -1 -1-2 1 0
0 0 0 - (-1)4 - (-2) 0 - 1 -1 - 0
~ 102100 ~
0-1-1-210
0016-1-1

La mitad izquierda de M está en forma triangular, por consiguiente, A es invertible. Si hubiera quedado toda una fila con ceros en la mitad A de M, la operación habría terminado (A no es invertible). A continuación, cogemos como pivote a₃₃, ponemos ceros encima de éste y seguimos operando hasta que nos quede una matriz diagonal.

~ 100-1122 
0-10401
0016-1-1

Ya que la matriz colocada en la mitad izquierda es diagonal, no hay que operar más. Transformamos la matriz diagonal en una matriz identidad; para ello hay que dividir la segunda fila entre -1:

~ 100-1122 
010-401
0016-1-1

La matriz que ha quedado en la mitad derecha de M es precisamente la matriz inversa de A:

A⁻¹ = -1122 
-401
6-1-1

Para comprobar si el resultado es correcto, se procede a multiplicar AA⁻¹, teniendo que dar como resultado la matriz identidad I.

Comprobación:

AA⁻¹ = I

 102 · -1122 =
2-13-401
4186-1-1
= -11 + 0 + 12  2 + 0 - 2  2 + 0 - 2 = 100 = I
-22 + 4 + 18 4 + 0 - 3 4 - 1 - 3010
-44 - 4 + 48 8 + 0 - 8 8 + 1 - 8001

Ejemplo de operaciones con matrices

Ejemplo nº 1

Sean

A = 241 
1-23
50-1
B = 3-1-2 
056
009
C = 201 
0-12
1-25

1) ¿Qué clase de matrices son?

2) Calcular:

-A - B + C.

A + B - C

3·A + C/2

3) Calcular:

(A·B)/C

4) Calcular la inversa de A (A⁻¹) y comprobar el resultado.

Solución

a)

Las tres matrices son cuadradas y de orden tres. A su vez, B es una matriz triangular, ya que todas las entradas debajo de la diagonal principal son ceros, y C es antisimétrica porque los elementos simétricos son opuestos entre sí.

b)

-A - B + C = - 241 - 3-1-2 + 201 =
1-230560-12
50-10091-25
= -2 - 3 + 2 -4 + 1 + 0 -1 + 2 - 1 = -3-30 
-1 - 0 + 02 - 5 - 1-3 - 6 + 2-1-4-7
-5 - 0 + 10 - 0 - 21 - 9 + 5-4-2-3
A + B - C = 241 + 3-1-2 - 201 =
1-230560-12
50-10091-25
= 2 + 3 - 2 4 - 1 - 0 1 - 2 + 1 = 330 
1 + 0 - 0 2 + 5 + 1 3 + 6 - 2147
5 - 0 - 1 0 + 0 + 2 -1 + 9 - 5423
   20-1   
 0-12 
3·A + C/2 = 3· 241 + 1-25 =
1-232
50-1 
= 6123 + 10 =
3-6901
150-3½-15/2
= 6 + 1 12 + 0 3 - ½ = 7125/2 =
3 + 0 -6 - ½ 9 + 13-13/210
15 + ½ 0 - 1 -3 + 5/231/2-1-1/2

c)

Puesto que (A·B)/C = A·B·C⁻¹, calcularemos primero la inversa de C y luego haremos el producto.

C⁻¹ = 20-1100 
0-12010
1-25001
~ 20-1100 ~
0-12010
0-411-102
~ 20-1100 ~
0-12010
00-314-2
~ -600-24-2 ~
030-2-114
00-314-2

Dividimos la primera fila entre -6, la segunda entre 3 y la tercera entre -3 para que en la mitad izquierda quede la matriz identidad,

~ 100-⅓-⅔ 
010-⅔-11/34/3
001-⅓-4/3

Por lo tanto, la matriz inversa de C es:

C⁻¹ = -⅓-⅔ 
-⅔-11/34/3
-⅓-4/3

A continuación, se calcula el producto de las matrices A y B,

A·B = 241 · 3-1-2 = 61829 ,
1-230563-1113
50-100915-5-19

Por último, calculamos (A·B)·C⁻¹

(A·B)·C⁻¹ = 61829 · 1/3-⅔ =
3-1113-⅔-11/34/3
15-5-19-⅓-4/3

Ejemplo del producto de matrices

(A·B)·C⁻¹ = 59/3326/3136/3 
421-5
44/3101/3-43/3

Sacando factor común ⅓, el resultado puede escribirse como:

(A·B)·C⁻¹ = ⅓· 59326136 
1263-15
44101-43

d)

Primero se construye la matriz M = (A ⋮ I) y luego se va desarrollando por Gauss. Así pues:

M = (A ⋮ I) = 241100 ~
1-23010
50-1001
~ 241100 ~
0-85-120
00156-502
~ 241100 ~
0-85-120
00392040-16

Se simplifica un poco para que las operaciones no sean tan costosas, dividiendo la tercera fila entre cuatro. De este modo, se tiene

~ 241100 ~
0-85-120
0039510-4

y se contínua calculando,

~ 78156034-104 ~
0-3120-642820
0039510-4

Se vuelve a simplificar, dividiendo la primera fila entre dos y la segunda entre cuatro,

~ 3978017-52 ~
0-780-1675
0039510-4
~ -3.04200-78-156-546 ~
0-780-1675
0039510-4

Puesto que ya ha quedado una matriz diagonal en la mitad izquierda de M, se procede a transformar esta mitad izquierda en una matriz identidad, dividiendo la primera fila entre -3.042, la segunda entre -78 y la tercera entre 39,

~ 1001/392/397/39 ~
0108/39-7/78-5/78
0015/3910/39-4/39

Así pues, la matriz que ha quedado en la mitad derecha es precisamente la matriz identidad, que sacando factor común 1/78 se puede escribir como:

A⁻¹ = (1/78)· 2414 
16-7-5
1020-8

Para comprobar el resultado, la matriz inversa de A o A⁻¹, tiene que cumplir

A·A⁻¹ = I.

Procedamos a la comprobación:

A·A⁻¹ = 241 · 2414 ·1/78 =
1-2316-7-5
50-11020-8

Ejemplo de cálculo de la inversa de una matriz

A·A⁻¹ = 7800 ·1/78 = 100 = I
0780010
0078001

Editor: Ricardo Santiago Netto (Administrador de Fisicanet).

Éste sitio web usa cookies, si permanece aquí acepta su uso.
Puede leer más sobre el uso de cookies en nuestra política de privacidad.