Introducción a Matrices y Determinantes
Introducción a Matrices y Determinantes
1.- GENERALIDADES
Se llama matriz de orden mxn a un conjunto de mxn números (o letras que representan
números) dispuestos en m filas y n columnas y encerrado entre paréntesis o corchetes.
4 2 5
es una matriz de orden 2x3
3 1 2
El conjunto de todas las matrices de orden mxn suele denotarse por Emxn.
La matriz suele denotarse con una letra mayúscula. A cada número que va dentro de la
matriz se le llama elemento. El elemento que está en la fila i y la columna j se denota con
la misma letra que la matriz, pero en minúscula, y dos subíndices que indican la posición
del elemento. Por ejemplo, en una matriz A, el elemento de la fila i y columna j se denota
a ij .
4 2 5
En la matriz el elemento a 13 es 5, el elemento a 23 es 2.
3 1 2
Una matriz suele definirse también a través de sus elementos, de la siguiente manera:
A=( a ij ) mxn . Lo anterior indicaría que tenemos una matriz de orden mxn, llamada A y a
cuyos elementos llamamos a ij .
Los elementos (números) de una matriz pueden ser reales o complejos. En el primer
caso se dice que la matriz es real y en el segundo que es compleja.
2 3 5 2-i 0 2 + i
4 2 1 es una matriz real 4 2 4 − i es una matriz compleja
3 4 8 3 1+ i 8
2 3 2 0 0 1
y son equidimensionales
0 −3 2 0 −2 −5
Se dice que dos matrices A y B son iguales si y sólo si, siendo equidimensionales, los
elementos que ocupan la misma posición son iguales. Es decir:
1
a ij = b ij ∀i = 1, 2,...,m ∧ ∀j = 1, 2, ..., n. (siendo m y n el número de filas y columnas
respectivamente)
Una matriz que tiene todos sus elementos nulos se denomina matriz nula, y se
representa por (0)mxn , siendo mxn el orden de la matriz.
0 0 0
es la matriz nula de 2x3
0 0 0
Se denomina matriz línea a la matriz que consta de una sola línea. Si está formada por
una sola fila se llama matriz fila (su orden es 1xn) y si está formada por una sola columna
se llama matriz columna (su orden es mx1).
2
( −1 2 0 ) es una matriz fila es una matriz columna
−1
−1 3 −2 1 −3 2
Si A= ⇒ -A= es la opuesta de A
2 1 5 −2 −1 −5
2 −1
es una matriz cuadrada de orden 2
4 5
Se llama diagonal principal de una matriz cuadrada a la formada por los elementos a ii
(es decir: a 11, a 22, …, a 33 ).
2
1 2 −1
4 3 −2 los elementos encuadrados forman la diagonal principal.
0 2 2
Una matriz cuadrada se dice triangular superior si todos los términos situados por
debajo de la diagonal principal son ceros (es decir, aij = 0 ∀i > j ). Análogamente, se dice
triangular inferior si todos los términos situados por encima de la diagonal principal son
ceros (es decir, aij = 0 ∀i < j ).
2 4 −1
0 3 4 es una matriz triangular superior de orden 3
0 0 2
a11 0 0 0
a21 a22 0 0
a31 a32 a33 0 es la forma general de una matriz triangular inferior
a ann
n1 an 2 an 3
Se llama matriz diagonal a toda matriz cuadrada cuyos términos situados fuera de la
diagonal principal son nulos (es decir, aij = 0 ∀i ≠ j ). Un caso particular de matriz
diagonal es la matriz escalar en la que todos los elementos de la diagonal principal son
iguales.
a11 0 0 0
2 0 0 0 a22 0 0
0
0 −1 0 0 a33 0 son matrices diagonales
0 0 5
0 ann
0 0
3
a 0 0 0
4 0 0 0 a 0 0
0
0 4 0 0 a 0 son matrices escalares
0 0 4
0 a
0 0
Se denomina matriz unidad o matriz identidad a una matriz escalar en la que los
elementos de la diagonal son todos 1. La matriz unidad de orden n suele denotarse por
In.
1 0 0
I 3 = 0 1 0 matriz unidad de orden 3
0 0 1
1 0 0 0
0 1 0 0
In = 0 0 1 0 matriz unidad de orden n
0 1 n
0 0
A+B= (a ij + b ij )
Es decir, es una matriz cuyos elementos se obtienen sumando uno a uno los elementos
que ocupan la misma posición en A y B. La matriz suma es, por lo tanto, una matriz del
mismo orden que A y B.
1 −2 4 0 3 2 1 1 6
Si A= y B= ⇒ A+B=
2 5 2 −1 4 2 1 9 4
4
a11 + b11 a12 + b12 ... a1n + b1n
a +b a22 + b22 ... a2 n + b2 n
⇒ A+B= 21 21
...
am1 + bm1 am 2 + bm 2 ... amn + bmn
Nota: obsérvese que las propiedades de la suma de matrices son las mismas que las de
la suma de números reales.
1 −2 4 2 −4 8
Si A= y α=2 ⇒ αA=
2 5 2 4 10 4
• α×(A + B) = α ×A + α ×B ∀ α ∧ ∀ A, B
• ( α + β ) ×A = α ×A + β ×A ∀ α, β ∧ ∀A
• ( α×β ) ×A = α ×( β×A ) ∀ α, β ∧ ∀A
• 1×A = A ∀A
5
p
(=
ab)ij ∑a
k =1
=
ik bkj i 1, 2,..., m=
∧ j 1, 2,..., n
-2
2 1 -1 -5 c11 =2 ⋅ (−2) + 1 ⋅ 1 + (−1) ⋅ 2 =−5
. 1 =
2 3 0 2x3 2 -1 2x1 c21 = 2 ⋅ (−2) + 3 ⋅ 1 + 0 ⋅ 2 = -1
3x1
2 1 2 −1
Si A = y B= ⇒ A×B≠ B×A como se comprueba a continuación:
3 2 0 3
2 1 2 −1 4 1
A=
⋅B ⋅ =
3 2 0 3 6 3
2 −1 2 1 1 0
=
B ⋅A ⋅ =
0 3 3 2 9 6
Si se quiere multiplicar ambos miembros de una ecuación matricial X = Y por una matriz
P, es importante que, o bien premultipliquemos ambos miembros por P, o bien
postmultipliquemos ambos miembros por P.
6
El producto de matrices verifica las siguientes propiedades:
• Asociativa
En el producto matricial puede ocurrir que A×B sea nulo siendo A≠(0) y B≠(0) (a diferencia
del producto con escalares donde a×b=0 si y sólo si alguno de los factores es cero).
2 0 0 0 0 0 0 0
Sean A = y B= ; A⋅B =
-3 0 2 4 1 0 0 0
6 2 4 1 3 0
Sean A = , B= y C= ; A ⋅ B = A ⋅ C con A ≠ (0) y B ≠ C
3 1 2 3 5 6
Es fácil comprobar que ∀A∈E mxn , I m ×A=A y A×In=A (es decir, las matrices Im y In juegan el
mismo papel que el número 1 en la multiplicación de números).
A p = (A ⋅ A ⋅ A ⋅ ⋅ A)
p veces
Por lo tanto:
A1 = A
A 2= A ⋅ A
A3 = A ⋅ A ⋅ A
7
Se observa fácilmente que esta operación no puede realizarse con matrices
rectangulares.
Propiedades:
•
ApAq = Ap+q ∀p,q∈+
•
(Ap)q = Apq ∀p,q∈+
1 2 97 -142 97 -142
Para A = , A2A3= A3A2=
-3 2 213 26 213 26
Podría pensarse que las potencias de matrices cuadradas funcionan a todos los efectos
como las potencias de números. Pero esto no es cierto:
(A+B)(A-B) = A2 + BA - AB - B2 ≠ A2 - B2
0 3
0 3 -2
=
A 2×3 ⇒=
t
A 3×2 3 2
3 2 2 -2 2
8
Propiedades:
• (At)t = A ∀A∈Emxn
0 3
0 3 -2 0 3 -2
( )
t
Ejemplo:
= A ⇒
= A 3 2 ⇒
t
=
At
3 2 2 -2 2 3 2 2
0 3 -2 1 -4 1
Ejemplo: A = y B=
3 2 2 -2 -1 0
0 3 1 -2 1 1
=A 3 = t
2 Bt -4 -1 ⇒ =
A + B -1 1
t t
-2 2 1 0 -1 2
1 1
1 -1 -1
⇒ ( A+B )= -1 1 =A t + Bt
t
=
A+B
1 1 2 -1 2
1 2
-8 4
Ejemplo: A = −3 4 y B =
2 1 -7 5
1 −3 2 −8 −7 −22 −4 −23
=At = t
, B ⇒=
Bt ⋅ A t
2 4 1 4 5 14 8 13
-22 14
−22 −4 −23
=A ⋅ B -4 8 ( A ⋅ B)t
⇒= =B ⋅ A
t t
-23 13 14 8 13
Una matriz cuadrada se dice que es simétrica si coincide con su traspuesta, lo cual sólo
se puede cumplir si los elementos dispuestos simétricamente a ambos lados de la
diagonal principal son iguales. Es decir:
1 0 −2
0 6 4 es simétrica
−2 4 3
9
Una matriz cuadrada A se dice que es antisimétrica si coincide con la opuesta de su
traspuesta, lo cual sólo puede ocurrir si los elementos dispuestos simétricamente a
ambos lados de la diagonal principal son opuestos y los elementos de la diagonal
principal son nulos; es decir
0 1 3
−1 0 −5 es antisimétrica
−3 5 0
a11 a12
Orden 2: A = = a11 ⋅ a22 − a12 ⋅ a21
a21 a22
Es más sencillo recordar esta última expresión haciendo uso de la regla de Sarrus, que
“dibuja” esta fórmula mediante un esquema que agrupa los términos que van en cada
sumando. En el dibujo de la regla de Sarrus que se muestra a continuación, los círculos
representan los elementos de la matriz; se han unido con líneas de un mismo color los
elementos que van multiplicados en un mismo sumando (a la izquierda los positivos y a
la derecha los negativos).
10
Términos positivos Términos negativos
11
Cuando un menor de un determinante está formado por las q primeras filas y las q
primeras columnas se denomina menor principal de orden q.
Así, por ejemplo, en un determinante de orden 5 para obtener sus menores principales
de orden 1, 2, 3, 4 y 5 (el propio determinante) vamos seleccionando los elementos (o
filas y columnas) como se indica a continuación:
Para realizar dichas transformaciones deben tenerse en cuenta siempre las propiedades
de los determinantes (no vamos a demostrar ninguna de ellas a pesar de que algunas
demostraciones son muy sencillas):
1.
El determinante de una matriz es igual al de su traspuesta: |A|=|A|t
2. Si una fila (o columna) tiene todos sus elementos nulos, el valor del determinante
es cero.
12
4. Si en un determinante los elementos de una fila o columna son múltiplos de los
de otra, el valor del determinante es cero.
9. Sea el determinante
a11 + b11 a12 + b12 ... a1n + b1n a11 a12 ... a1n b11 b12 ... b1n
a21 a22 ... a2 n a21 a22 ... a2 n a21 a22 ... a2 n
= +
... ... ... ... ... ... ... ... ... ... ... ...
an1 an 2 ... ann an1 an 2 ... ann an1 an 2 ... ann
10. A ⋅ B = A ⋅ B . ∀A,B∈Enxn.
13
a11 a12 ... a1n
a a22 ... a2 n
Así, si D = 21 , desarrollando el determinante por la i-ésima fila se
... ... ... ...
an1 an 2 ... ann
obtiene:
D = a1 j ⋅ A1 j + a2 j ⋅ A2 j + + anj ⋅ Anj
2 1 0 3
4 2 2 1
Ejemplo: Calculemos D =
0 3 2 2
0 0 1 3 ← 2 elementos nulos
1 0 3 2 0 3 2 1 3 2 1 0
D=
−0 ⋅ 2 2 1 + 0 ⋅ 4 2 1 − 1 ⋅ 4 2 1 + 3 ⋅ 4 2 2 =
−66
3 2 2 0 2 2 0 3 2 0 3 2
Hay que tener mucho cuidado en el empleo de las propiedades 7 y 8, ya que alteran el
valor del determinante original, por lo que hay que multiplicar posteriormente el valor
14
1
obtenido por el número adecuado (-1 o , respectivamente) para mantener el valor del
α
determinante inicial.
5.1.- Definición
Se llama rango de una matriz A de orden mxn al orden del mayor menor no nulo que
puede extraerse de ella. También, teniendo en cuenta las propiedades de los
determinantes, se define como el mayor número de filas o columnas que son
linealmente independientes en A.
5.2.- Cálculo
El primer procedimiento que nos puede venir a la mente para calcular el rango de una
matriz es emplear directamente su definición. Así, dada una matriz de orden mxn
comenzaremos extrayendo uno a uno sus menores de mayor orden (supongamos que
éste es n). En cuanto algún menor sea no nulo pararemos el proceso y diremos que el
rango de la matriz será n. Si extraídos todos los menores de orden n no hemos
encontrado ninguno no nulo, sabremos que el rango es menor que n, pero no sabremos
cuánto menor. En este último caso tendríamos que continuar el proceso, extrayendo
uno a uno los menores de orden n-1 hasta que encontremos uno no nulo, momento en
el que podríamos decir que el rango de A es n-1. Si extraídos todos los menores de orden
n-1 no hemos encontrado ninguno no nulo, sabremos que el rango es menor que n-1,
pero no sabremos cuánto menor. Así, tendríamos que pasar a plantear los menores de
orden n-2 … Se continúa con el proceso hasta obtener algún menor no nulo de algún
orden; el rango de la matriz sería este orden.
Puede verse que este proceso es muy laborioso, por lo que si hay algún procedimiento
alternativo es conveniente emplearlo. Uno de estos procedimientos, mucho más
eficiente y comúnmente empleado en la práctica, se muestra en la teoría
correspondiente a sistemas de ecuaciones lineales al explicar cómo se discute la
existencia de solución de este tipo de sistemas cuando aplicamos el método de Gauss.
6.1.- Definición
Una matriz cuadrada A se dice regular si existe otra matriz cuadrada, a la que
llamaremos inversa de A y denotaremos por A-1, tal que A×A-1 =A-1×A=I. En caso de no
existir esta matriz A-1 se dice que A es singular.
0 −1
1 2 -1 1 0
La inversa de A= es A = 1 , ya que A ⋅ A =
-1
1
−1 0 0 1
2 2
15
Cuando una matriz cuadrada tiene inversa, ésta es única.
Pero, además, puede demostrarse que una matriz cuadrada es regular si y sólo si su
determinante es no nulo. De la misma manera, una matriz cuadrada es singular si y sólo
si su determinante es nulo. Recordando la definición de rango de una matriz, lo anterior
es equivalente a decir que una matriz cuadrada de orden n es regular si y sólo si el rango
de A es n, y es singular si el rango de A es menor que n.
1 2
Por ejemplo, para calcular la inversa de A= planteamos que buscamos una
2 1
matriz A-1 que, en caso de existir, ha de ser también cuadrada de orden 2, por lo que será
x y 1 2 x y 1 0
de la forma A −1 = −1
, y ha de cumplir A ⋅ A = I ⇒ ⋅ =
z t 2 1 z t 0 1
Esta expresión representa la igualdad de dos matrices y, puesto que para que dos
matrices sean iguales han de serlo elemento a elemento (según la definición de igualdad
de matrices), de aquí se obtiene el siguiente sistema de ecuaciones lineales:
x + 2z =1
y + 2t =
0
2 x + z = 0
2 y + t =1
16
La discusión de la existencia de solución de este sistema y su resolución se abordarán
por cualquier método conocido. En este ejemplo, el sistema tiene solución y es la
siguiente:
1 2
1 2 2 1 −3 3
x=
− , y= , z= − . Por lo tanto, A −1 =
, t=
3 3 3 3 2 −1
3 3
En caso de que la matriz cuya inversa buscamos no sea regular, llegaremos a un sistema
incompatible (es decir, sin solución). Por ejemplo:
2 8 −1 x y
Si A = ⇒ A =
1 4 z t
2x + 8z = 1
2 8 x y 1 0 2 y + 8t = 0
A ⋅ A −1 =
I ⇒ ⋅ = ⇒ este sistema no tiene solución
1 4 z t 0 1 x + 4z = 0
y + 4t =
1
(no existen valores de x,y,z,t que cumplan todas las ecuaciones simultáneamente). Por
lo tanto, no existe una matriz A-1/ A×A-1=I (es decir, A es singular).
−1 1
6.3.2.- Cálculo de A-1 empleando la fórmula: A= ⋅ Aa
A
Se puede demostrar que si A es regular, su inversa se puede calcular según la expresión
anterior, donde Aa denota a la matriz adjunta de A y es aquella matriz que resulta de
sustituir cada elemento de At por su adjunto.
1 1 −1
Para calcular la inversa de A = 2 3 1 necesitamos su determinante:
1 1 3
A = 1 ⋅ 3 ⋅ 3 + 1 ⋅ 1 ⋅ 1 + 2 ⋅ 1 ⋅ (−1) − 1 ⋅ 3 ⋅ (−1) − 1 ⋅ 1 ⋅ 1 − 2 ⋅ 1 ⋅ 3 = 4 ≠ 0 ⇒ ∃A −1
y su adjunta:
3 1 1 1 1 3
−
1 3 -1 3 -1 1
2 8 −4 4
1 2 1
1 1 1 1 2
At = 1 3 1 −
⇒ Aa = − = −5 4 −3
-1 1 3 1 3 -1 3
-1 1
−1 0 1
2 1
−
1 1 1 2
3 1 3
1 1 1
17
2 −1 1
Entonces, A = -5/4 1 -3 4 .
-1
-1 4 0 1 4
18