Capítulo 3. Matrices y Determinantes

Apuntes de Álgebra Lineal
Capítulo 3
Matrices y Determinantes
3.1.
Operaciones con matrices
3.1.1.
Suma, resta y multiplicación por escalares
Las matrices de un tamaño fijo m × n se pueden sumar entre sí y esta operación de sumar se
puede definir operando elemento a elemento o bien operando columna a columna. Por ejemplo,
consideremos las siguientes matrices de tamaño 2 × 3:
A=
1
4
2
5
3
6
,
B=
−1
0
.
1 −2
7
−4
Operar elemento a elemento es poner:
A+B =
1+7 2−1 3+0
4−4 5+1 6−2
=
8
0
3
,
4
1
6
mientras que operar columna a columna es poner A = [a1 a2 a3 ], B = [b1 b2 b3 ] con:
a1 =
1
,
4
a2 =
2
,
5
a3 =
3
,
6
b1 =
7
,
−4
b2 =
−1
,
1
b3 =
0
,
−2
de forma que
a1 + b1 =
1
7
8
+
=
,
4
−4
0
a2 + b2 =
2
−1
1
+
=
,
5
1
6
a3 + b3 =
3
0
3
+
=
,
6
−2
4
y entonces
8
0
1
6
3
.
4
La definición de suma de matrices elemento a elemento es fácil de entender y de aplicar,
pero tediosa para la demostración de propiedades. La segunda es una definición un poco más
sutil y parece más complicada, pero bien mirada es incluso más sencilla que la primera y tiene la
gran ventaja de que permite hacer demostraciones más sencillas de las propiedades de matrices.
En cualquier caso, usaremos la definición primera sólo para demostrar el número más reducido
posible de propiedades, de las cuales se puedan deducir todas las demás sin tener que recurrir
constantemente a la definición.
1
Versión de 8 de noviembre de 2016, 10:36 h.
A + B = [ a1 + b1 | a2 + b2 | a3 + b3 ] =
3.1. Operaciones con matrices
3. Matrices y Determinantes
Igual que ocurre con la suma de matrices, la multiplicación de números por matrices se puede definir elemento a elemento o por columnas. Por ejemplo, con la matriz B anterior, podemos
poner:
3×7
3 × (−1) 3 × 0
21 −3
0
3B =
=
3 × (−4) 3 × 1
3 × (−2)
−12
3 −6
o podemos calcular:
21
−3
0
3b1 =
, 3b2 =
, 3b3 =
,
−12
3
−6
y
3B = [3b1 3b2 3b3 ] =
21
−12
−3
0
.
3 −6
Usando cualquiera de esas definiciones de suma y multiplicación por números no tiene
ninguna dificultad el demostrar que estas operaciones de las matrices reales de un tamaño fijo
m × n cumplen los axiomas de espacio vectorial real donde la matriz cero m × n, denotada 0, es
aquella todos cuyos elementos son cero y la matriz opuesta, − A, de una matriz A es aquella
cuyos elementos son los opuestos de los de A:
(a) Propiedades de la suma de matrices (“Grupo conmutativo”):
1. Propiedad asociativa de la suma: A + ( B + C ) = ( A + B) + C.
2. Matriz cero: 0 + A = A + 0 = A.
3. Opuesta de una matriz: − A + A = A + (− A) = 0.
(abreviación: en lugar de A + (− B) se puede escribir para mayor brevedad A − B).
4. Propiedad conmutativa de la suma: B + A = A + B.
(b) Propiedades distributivas de la multiplicación por números:
1. Propiedad distributiva para la suma de matrices: x ( A + B) = xA + xB.
2. Propiedad distributiva para la suma de números: ( x + y) A = xA + yA.
(c) “Acción de los escalares”:
1. Propiedad asociativa del producto de números por matrices: x (yA) = ( xy) A.
2. Ley de identidad: 1 A = A.
(El número 1 es neutro para la multiplicación de números por matrices.)
Además, se cumplen las siguiente propiedades, las cuales ya sabemos que se deducen de las
anteriores pero que en este caso es más fácil demostrar directamente a partir de la definición
del producto de un número por una matriz:
(a) Toda matriz multiplicada por cero da la matriz cero: 0A = 0.
(b) Todo múltiplo de la matriz cero es la matriz cero: c 0 = 0.
(c) Multiplicar una matriz por −1 da la matriz opuesta: (−1) A = − A.
3.1.2.
Producto de matrices
Motivación
Sabemos multiplicar una matriz por un vector y sabemos que el resultado es otro vector:
b = Ax. Supongamos ahora que tenemos una segunda matriz, B, que se puede multiplicar
por b y que la usamos para calcular el vector Bb = B( Ax). Problema: ¿Qué matriz hay que
multiplicar por x para obtener directamente el vector B( Ax)?
Para contestar a esto es necesario recordar que Ax es una combinación lineal de las columnas
de A: Ax = x1 a1 + · · · + xn an y que la operación de producto de una matriz por un vector
2
3. Matrices y Determinantes
3.1. Operaciones con matrices
respeta las combinaciones lineales, por lo que B( Ax) = B( x1 a1 + · · · + xn an ) es igual a una
combinación lineal de los vectores Ba1 , . . . Ban obtenidos al multiplicar B por cada columna de
A. Estos vectores son matrices columna con el mismo número de elementos que cada columna
de B y forman una nueva matriz que llamaremos C. La combinación lineal de estas columnas
usando los elementos de x como coeficientes es a la vez el producto de C por x y también el
vector B( Ax), por tanto C es la matriz que multiplicada por x nos da B( Ax).
Definición del producto de matrices
El producto de una matriz A por una matriz B es la matriz, denotada AB, cuyas columnas
son el resultado de multiplicar A por cada columna de B. Es decir: cada columna del producto AB
es el producto de A por la correspondiente columna de B:
Definición del
col j ( AB) = A · col j ( B).
producto de
(3.1) matrices.
Por tanto cada columna de AB es una combinación lineal de las columnas de A
La propiedad característica del producto así definido es que para cualquier vector x que propiedad
característica
tenga tantos elementos como columnas tiene B, se cumple
del producto
( AB)x = B( Ax).
La regla “fila por columna”: Suponiendo que las columnas de B son b1 , . . . , b p y para calcular
AB calculamos cada uno de los productos Ab j por la regla “fila×columna” llegamos a la conclusión de que cada fila del producto AB es igual a la correspondiente fila de A multiplicada
por B:
filai ( AB) = filai ( A) · B
(3.2)
1 Ejercicio de tarea. Usa las ecuaciones (3.1) y (3.2) como corresponda para calcular la fila 3 y la
columna 2 del producto AB y lo mismo para el producto BA, donde:




1 1 1
2 0 0
A = 1 2 3 ,
B = 0 3 0 .
1 4 5
0 0 5
Propiedades del producto de matrices
A continuación se presenta la lista de las propiedades de la multiplicación de matrices. Aquí
λ es un escalar cualquiera e Ik denota la matriz identidad k × k.
(a) A( BC ) = ( AB)C (propiedad asociativa).
(b) A( B + C ) = AB + BC (propiedad distributiva por la izquierda).
(c) ( B + C ) A = BA + CA (propiedad distributiva por la derecha).
En las propiedades (d) y (e) se supone que A tiene m filas y n columnas.
(d) Im A = A = AIn .
(e) (λIm ) A = λA = A(λIn ).
Toda matriz de la forma λIm (o sea, que sea un múltiplo escalar de una matriz identidad)
se llama una matriz numérica o matriz escalar.
De ésta junto con la propiedad asociativa se puede deducir la siguiente, de la cual ésta es
un caso particular:
(f) (λA) B = λ( AB) = A(λB).
2 Ejercicio de tarea. Demuestra que la propiedad (f) se deduce de las propiedades (a) y (e).
3
3.1. Operaciones con matrices
3. Matrices y Determinantes
Sobre la no-conmutatividad del producto de matrices y otras advertencias.
(a) No conmutatividad: En general AB 6= BA. Por ejemplo:
1
3 4
3 4 =
mientras que
2
6 8
3
4
1
= 11.
2
Pero la razón de esto no es que los resultados sean matrices de distinto tamaño, pues en
general AB 6= BA incluso para matrices cuadradas:
0 1
4 0
∗ 5
4 0
0 1
∗ 4
=
mientras que
=
.
2 3
0 5
∗ ∗
0 5
2 3
∗ ∗
6= 0 tales que A2 = 0. Por ejemplo,
1
0 1
0 0
=
.
0
0 0
0 0
(b) Matrices nilpotentes: Existen matrices A
2 0 1
0
=
0 0
0
(c) Divisores de cero: Un producto de matrices igual a cero no implica que algún factor sea
cero, es decir, si AB = 0, en general no se puede concluir que A = 0 o B = 0. Por ejemplo:
0
1 0
= 0.
1
(d) No se cumple la Ley de cancelación: De AB = AC, en general no se puede deducir B = C.
Por ejemplo:
0 1
0 2
0 1
0 3
=
.
0 0
0 0
0 0
0 0
3.1.3.
Traspuesta de una matriz y sus propiedades.


a1
 
Para cualquier vector (o matriz de una sola columna) a =  ... , la matriz traspuesta es la
am
matriz fila que tiene los mismos elementos que a:
aT = a1
...
am .
Para una matriz general de varias columnas A = [a1 . . . an ], la traspuesta es la matriz cuyas
filas son las traspuestas de las columnas de A:
 T
a1
 .. 
T
A =  . .
aTn
De esta definición se deduce que la traspuesta de una matriz m × n es una matriz n × m.
Por ejemplo:
1
4
2
5
3
6
T

1
= 2
3

4
5.
6
La operación de calcular la matriz traspuesta tiene las siguientes propiedades cuya demostración es un ejercicio que debe hacer el estudiante:
T
(a) AT = A
(b) ( A + B)T = AT + BT
(c) (λA)T = λ AT
(d) ( AB)T = BT AT (¡cambian de orden!)
4
3. Matrices y Determinantes
3.2. Matrices elementales
Enlaces a todos los ejercicios de tarea de esta sección
Usa los siguientes enlaces para visualizar cada uno de los ejercicios de tarea que aparecen
en esta sección:
Enlaces: Ejercicio 1, Ejercicio 2.
3.2.
Matrices elementales
Ejecución de operaciones elementales mediante multiplicación de matrices
Consideremos las siguientes matrices:

1
E1 =  0
3

0
0,
1
0
1
0

0
E2 =  1
0
1
0
0

0
0,
1

1
E3 =  0
0

0
0.
5
0
1
0
Obsérvese que todas ellas son el resultado de realizar una operación elemental de filas sobre la
matriz identidad 3 × 3. Ahora calculemos los productos E1 A, E2 A, E3 A, donde A es una matriz
cualquiera de tres filas. Por ejemplo:

a
A=c
e

b
d
f

con lo cual:

b
d,
f + 3b
a
c
E1 A = 
e + 3a

c
E2 A =  a
e

d
b,
f

a
E3 A =  c
5e

b
d
5f
Vemos que cada uno de los productos nos da como resultado el mismo que si hubiésemos hecho
sobre A la operación elemental con la que se consiguió la correspondiente matriz E. Esto es un
hecho general y se cumple el siguiente teorema:
Teorema: Sea A una matriz de m filas. Para realizar una operación elemental sobre A, basta realizar
dicha operación elemental sobre las filas de la matriz identidad Im y multiplicar el resultado por A.
Definición de matrices elementales
Definición: Se llama matriz elemental a toda matriz obtenida realizando una operación elemental sobre las filas de una matriz identidad.
Ejemplos: (1) Sea E la matriz obtenida al realizar una operación elemental de intercambio de
filas sobre la matriz identidad 2 × 2:
I2 =
1
0
0
1
F ↔F
2
1
−−−
→
E=
0
1
1
0
E es la matriz elemental correspondiente a la operación elemental
“intercambiar las filas 1 y 2”.
.
Sea ahora A una matriz cualquiera con 2 filas. Entonces:
A=
a
d
b
e
c
f
F1 ↔ F2
−−−→
d
a
e
b
f
c
y también:
EA =
0
1
1
0
a
d
b
e
c
f
=
d
a
e
b
f
c
El producto EA da como resultado el mismo que realizar sobre A la operación elemental de
“intercambiar las filas 1 y 2 de A”.
5
.
3.3. Matriz inversa
3. Matrices y Determinantes
(2) Sea E la matriz obtenida al realizar una operación elemental de reemplazo de filas sobre
la matriz identidad 3 × 3:




1 0 0
1 0 0
E es la matriz elemental correspondiente
F2 +λF3
a la operación elemental “sumar
I3 =  0 1 0  −−−−→ E =  0 1 λ  .
a la fila 2 la 3 multiplicada por λ”.
0 0 1
0 0 1
Sea ahora A una matriz cualquiera con 3 filas. Entonces:
!
!
A=
a
c
e
b
d
f
F2 +λF3
−−−−→
a
c + λe
e
b
d + λf
f
y también:
EA =
1
0
0
0
1
0
0
λ
1
!
a
c
e
b
d
f
!
=
a
c + λe
e
b
d + λf
f
!
El producto EA da como resultado el mismo que realizar sobre A la operación elemental de
“suma a la fila 2 la fila 3 multiplicada por λ”.
3.3.
Matrices inversibles y cálculo de la matriz inversa
inversa por la
izquierda
Sea A una matriz de m filas y n columnas. Se llama matriz inversa por la izquierda de A a
toda matriz B tal que BA = In . En tal caso B será necesariamente n × m. Además, el sistema
homogéneo Ax = 0 es determinado (el vector 0 es la única solución porque si hubiera otra, x,
entonces x = BAx = B0 = 0) y por por ser determinado no tiene variables libres. Por tanto, si
A tiene una inversa por la izquierda, A tiene un pivote en cada columna y no puede tener menos filas que
columnas: m ≥ n.
inversa por la
derecha
Se llama matriz inversa por la derecha de A a toda matriz C tal que AC = Im . En tal caso C
será necesariamente n × m. Además, todo sistema Ax = b será compatible (ya que al menos el
vector Cb será una solución). Por tanto, si A tiene una inversa por la derecha, entonces A tiene un
pivote en cada fila y no puede tener más filas que columnas: m ≤ n.
matriz
inversible
Una matriz A se llama inversiblesi tiene una inversa por la derecha y una inversa por la
izquierda. En tal caso, se puede demostrar que la inversa por la izquierda es la misma matriz
que la inversa por la derecha1 (llamándose entonces simplemente la matriz inversa de A). Por lo
dicho antes, la matriz A es necesariamente una matriz cuadrada (al igual que su inversa). Por
tanto una matriz es inversible si y sólo si existe una matriz B tal que
BA = I
y
AB = I.
Toda matriz inversible es una matriz cuadrada. La matriz inversa de A se denota por A−1 .
matriz singular
Se llama matriz singular a toda matriz cuadrada que no sea inversible.
3 Ejercicio de tarea. Suponiendo que M es una inversa de A por la izquierda y que N es una
inversa de A por la derecha, demuestra:
(a) M y N tienen el mismo número de filas y columnas que la traspuesta de A.
(b) M = N.
(c) A es una matriz cuadrada.
4 Ejercicio de tarea. Demostrar que si A es una matriz cuadrada entonces basta que tenga una
inversa por un lado para que sea inversible.
1 Si
AC = Im y BA = In entonces C = ( BA)C = B( AC ) = B.
6
3. Matrices y Determinantes
3.3. Matriz inversa
Propiedades de la inversa
−1
(a)
A −1
(b)
( AB)−1
(c)
T −1
A
= A.
= B −1 A −1 .
T
= A −1 .
De la segunda de estas propiedades se deduce que el producto de varias matrices inversibles
de tamaño n × n es de nuevo una matriz inversible y su inversa es el producto de las inversas
en el orden opuesto:
−1
1
−1
A1 · · · A k
= A−
k · · · A1 .
Inversa de una matriz
2 × 2. La inversa (suponiendo que existe) de una matriz cuadrada
a b
de orden 2, A = c d , se obtiene multiplicando el inverso de ad − bc por la matriz obtenida al
intercambiar las posiciones de los elementos de la diagonal y cambiar de signo a los otros dos:
a
c
b
d
−1
=
1
ad − bc
d
−c
La comprobación de esto es un simple cálculo:
a b
d −b
ad − bc
=
c d
−c
a
0
−b
.
a
0
ad − bc
.
En el caso c = 0 es útil expresar la fórmula de la inversa de esta forma:
a
0
b
d
−1
=
a −1
0
− a−1 bd−1
.
d −1
Uso de la inversa en la resolución de sistemas
Para aquellos sistemas de ecuaciones lineales Ax = b en los que la matriz de coeficientes, A,
sea una matriz inversible existe automáticamente solución única y esa solución está dada por
x = A−1 b.
Evidentemente esto sólo se aplica a sistemas con el mismo número de ecuaciones que incógnitas.
Los casos en los que esta fórmula para la solución tiene mayor utilidad son los de sistemas de
dos ecuaciones con dos incógnitas, en los cuales es muy sencillo ver si la matriz de coeficientes
tiene inversa y apenas es nada costoso el calcular la inversa.
Ejemplo:
El sistema
3x1 + 4x2 = 3
5x1 + 6x2 = 7
tiene matriz de coeficientes inversible y su inversa es la matriz
3
5
4
6
−1
1
=
18 − 20
6
−5
−4
3
por tanto:
7
x=A
−1
1
b=
2
−6
4
5 −3
3
5
=
.
7
−3
3.3. Matriz inversa
3. Matrices y Determinantes
Caracterizaciones de las matrices inversibles
Teorema de la matriz inversible Una matriz cuadrada real A de orden n es una matriz inversible si y sólo si se cumple cualquiera de las siguientes condiciones:
(a) A tiene una inversa por la derecha.
(b) A tiene una inversa por la izquierda.
(c) El sistema homogéneo representado por la ecuación Ax = 0 es determinado (no tiene más
solución que la trivial).
(d) A tiene n columnas pivote.
(e) La forma escalonada reducida de A es la matriz identidad de orden n.
(f) Para todo b ∈ Rn el sistema representado por la ecuación Ax = b es compatible.
(g) Las columnas de A generan Rn .
(h) AT es una matriz inversible.
Lo mismo es cierto para las matrices complejas cambiando R por C en los enunciados en los
que aparece.
Las inversas de las matrices elementales
Toda matriz elemental tiene inversa y su inversa es otra matriz elemental: La correspondiente
a la operación elemental inversa. En consecuencia, es muy sencillo escribir la matriz inversa de
una matriz elemental. Por ejemplo:




1 0 0
1
0
0
0 1
0
1
1 −λ .
Si E =
entonces E−1 =
; si E =  0 1 λ  entonces E−1 = 0
1 0
1 0
0 0 1
0
0
1
Algoritmo para averiguar si una matriz es inversible y calcular su inversa
Teorema. Una matriz cuadrada A es inversible si y sólo si es equivalente por filas a la correspondiente
matriz identidad. En ese caso, cualquier sucesión de operaciones elementales de filas que reducen A a la
identidad I también reducen I a la inversa de A.
Demostración: Cada paso de la reducción corresponde a la multiplicación por la izquierda por
una matriz elemental:
A ∼ E1 A ∼ E2 E1 A ∼ · · · ∼ E p ( E p−1 · · · E1 ) A = I
por tanto ( E p · · · E1 ) A = I, lo que significa que A es inversible y su inversa es:
A−1 = E p · · · E1 .
Corolario. Si una matriz cuadrada A es inversible entonces su forma escalonada reducida de ( A| I ) es
( I | A−1 ). Recíprocamente, si la forma escalonada reducida de ( A| I ) es ( I | M) entonces A es inversible y
su inversa es A−1 = M
De esto se deduce que la forma escalonada reducida de la matriz por bloques formada por
una fila de dos bloques en la que el primero es la matriz A y el segundo la matriz identidad, es
la que tiene como primer bloque la identidad y segundo la inversa de A:
( A | I ) ∼ · · · ∼ ( I | A −1 )
En otras palabras: Las mismas operaciones elementales que transforman una matriz inversible
en la matriz identidad, transforman la matriz identidad en la inversa de la matriz.
8
3. Matrices y Determinantes
3.4. Factorización LU
Cómo calcular solamente una columna (o fila) particular de la matriz inversa.
la columna j de A−1 basta resolver el sistema
Para calcular
Ax = e j
donde e j es la columna j de la matriz identidad.
Para calcular la fila i de A−1 basta calcular la columna i de AT
dicho antes, basta resolver el sistema
AT x = e i .
−1
, para lo cual, según lo
Enlaces a todos los ejercicios de tarea de esta sección
Usa los siguientes enlaces para visualizar cada uno de los ejercicios de tarea que aparecen
en esta sección:
Enlaces: Ejercicio 3, Ejercicio 4.
3.4.
Un caso especial de la forma escalonada de una matriz: La
Factorización LU
Motivación
Supongamos que necesitamos resolver varios sistemas de ecuaciones que tienen la misma
matriz de coeficientes A y sólo se diferencian en los términos independientes. Si se conocen dos
matrices L, U tales que L es una matriz unitriangular inferior, es decir, triangular inferior con
unos en la diagonal, U es una matriz escalonada y son tales que LU = A, entonces el sistema
Ax = b o L(Ux) = b se puede resolver en dos pasos muy sencillos: Primero resolvemos el
sistema
Ly = b
y después el sistema
Ux = y.
Una descomposición de una matriz, A = LU, como producto de dos matrices como las descritas
se conoce como una factorización LU (o factorización “Lower-Upper") de A (leído: « factorización
“L”-“U” »).
Existencia y unicidad
No toda matriz admite una factorización LU, pero, como se verá más abajo, una matriz A
la admitirá siempre que sea posible obtener una forma escalonada de A sin realizar operaciones
elementales de intercambio de filas (ni, evidentemente, operaciones de reemplazo que se puedan
combinar para conseguir un intercambio de filas, lo cual excluye las operaciones de reemplazo
“regresivas” en que se suma un múltiplo de una fila a otra fila anterior ya que combinando tales
operaciones con operaciones de reemplazo “progresivas” es posible conseguir un intercambio
de filas). En otras palabras: La matriz A admite una factorización LU si se puede poner A
en forma escalonada utilizando únicamente operaciones de reemplazo progresivas (sumar un
múltplo de una fila a otra posterior a ella).
Recíprocamente, si A admite una factorización LU, entonces es evidente que L, por ser una
matriz triangular inferior, se puede poner en forma escalonada utilizando únicamente operaciones de reemplazo progresivas (¡obteniéndose una matriz identidad!). Esas operaciones transformarán A en U y por tanto A se puede poner en forma escalonada utilizando únicamente
operaciones de reemplazo progresivas. En consecuencia:
9
Una matriz se
llama
unitriangular si
es una matriz
triangular y
todos los
elementos de
su diagonal
son iguales a 1.
3.4. Factorización LU
3. Matrices y Determinantes
La factorización LU de una matriz A existe si y sólo si A se puede poner en forma escalonada utilizando
únicamente operaciones elementales de reemplazo progresivas (sumar un múltplo de una fila a otra
posterior a ella).
En general, la factorización LU de una matriz no es única, como muestra el siguiente ejemplo.
Sean






1 0 0
1 0 0
1 2
L1 = 2 1 0 ,
L2 = 2 1 0 ,
U = 0 0 .
3 0 1
3 1 1
0 0
Evidentemente L1 y U constituyen una factorización LU de su producto L1 U e igualmente L2 y
U constituyen una factorización LU de su producto L2 U. Por otra parte, es fácil comprobar que
L1 U = L2 U, por consiguiente tenemos aquí dos factorizaciones LU de la misma matriz.
Algoritmo de la factorización LU
(a) Se reduce A a una forma escalonada, U, mediante una sucesión de operaciones de reemplazo progresivas si es posible. Si no es posible, la matriz A no tiene factorización LU.
(b) Se colocan los elementos bajo la diagonal de unos de L de forma que la misma sucesión de
operaciones elementales reduzcan L a la matriz identidad.
Entonces
( E p · · · E1 ) A = U
( E p · · · E1 ) L = I.
y
De aquí se deduce
A = ( E p · · · E1 )−1 U = LU.
Ejemplo
Calcular la factorización LU de la matriz

2
4
 −4 −5
A=
 2 −5
−6
0

−1
5 −2
3 −8
1

−4
1
8
7 −3
1
La primera columna de L es la primera columna pivote de A dividida por el primer elemento:

1
 −2
L=
 1
−3
0
1
0
0
1

0
0

0
1
Claramente, las operaciones elementales que produzcan ceros en la primera columna de A
también producirán ceros en la primera columna de L. Ahora intentamos hallar una forma
escalonada de A utilizando solamente operaciones de reemplazo. Después de crear ceros en la
primera columna, el problema ha sido reducido a uno con una fila y una columna menos:

2
 −4
A=
 2
−6
4
−5
−5
0
 

−1
5 −2
2
4 −1
5 −2

3 −8
1
3
1
2 −3
 ∼ 0

−4
1
8 0 −9 −3 −4 10
7 −3
1
0 12
4 12 −5
10
3. Matrices y Determinantes

3
A 0 =  −9
12
3.5. Factorización de rango máximo


1
2 −3
−3 −4 10
4 12 −5
de donde
1
L 0 =  −3
4

1
 −2
L=
 1
−3

0
1
0
0
1
y
0
1
−3
4
0
0
1

0
0

0
1
En el siguiente paso:

3
A 0 =  −9
12
1
−3
4
 

−3
3 1 2 −3
10 ∼ 0 0 2
1
−5
0 0 4
7
2
−4
12

A00 =
y, puesto que
A00
2
4
1
7
=
2
4

de donde
1
7
2
 −4
A=
2
−6
3.5.
∼
4
−5
−5
0
2
0
−1
3
−4
7
L00 =
1
2
0
1
1
 −2
L=
 1
−3
y
0
1
−3
4
0
0
1
2

0
0

0
1
1
, tenemos:
5
5
−8
1
−3
−2
2
1
0
∼ 0
8
1
0


4
3
0
0
−1
1
0
0
5
2
2
0
−2
−3
=U
1
5

Rango de una matriz y factorización de rango máximo
Se llama rango de una matriz A, y se denota ran( A), al número de pivotes de A (es decir, el
número de filas no nulas o de elementos principales en cualquier forma escalonada de A):
Rango de una
matriz.
ran( A) = número de pivotes de A.
Como los distintos pivotes están en distintas filas, el rango de una matriz es siempre menor
o igual al número de filas de la matriz. Análogamente, como los distintos pivotes están en
distintas columnas, el rango de una matriz es siempre menor o igual al número de columnas
de la matriz. En consecuencia, si A es una matriz m × n de rango r, su rango es siempre menor
o igual que el mínimo de m y n: r ≤ mı́n(m, n). En consecuencia, decimos que A tiene rango
máximo si su rango verifica:
Matriz de
rango máximo.
r = mı́n(m, n)
En la matriz (ampliada) de un sistema de ecuaciones lineales compatible no hay ningún
pivote en la columna de los términos independientes por tanto su rango es igual al de la matriz
de coeficientes. Por otra parte, el sistema será determinado si y sólo si ese rango es igual al
número de incógnitas.
Se llama factorización de rango máximo de una matriz A a toda descomposición de A como
producto de dos matrices, A = MN, de tal forma que M y N sean matrices de rango máximo y
ambas con el mismo rango que A. Veremos a continuación que el cálculo de la forma escalonada
reducida nos ofrece sin más cálculos una factorización de rango máximo de cualquier matriz no
nula.
Evidentemente la factorización de rango máximo de una matriz no es única, pues si A = MN
es una factorización de rango máximo de una matriz no nula de rango r, para cualquier matriz
inversible r × r, B, A = ( MB)( B−1 N ) es también una factorización de rango máximo de A.
11
3.5. Factorización de rango máximo
3. Matrices y Determinantes
La factorización de rango máximo proporcionada por la forma escalonada
reducida
Para una matriz cualquiera A, de m filas, n columnas y rango r < mı́n(m, n), sea M la matriz
formada por sus columnas pivote y sea N la matriz formada por las filas no nulas de su forma
escalonada reducida. El número de columnas de M es igual al rango de A y lo mismo ocurre
con el número de filas de N ya que la forma escalonada reducida tiene exactamente una fila no
nula por cada pivote de A. En consecuencia se puede multiplicar M por N y el resultado es una
matriz m × n.
Ejemplo.
Sea A la matriz


−2 −6
6
4
6
 0
0
1
2
0
.
A=
 0
0
0
3
9
1
3 −3 −2 −3
Calculamos su forma escalonada reducida:

1
0
A −−−−−−−→ 
0
F4 − F1
0
1
3 F3
− 12 F1
3
0
0
0
−3
1
0
0
−2
2
1
0
F2 − 2F3
1
−3
0  F1 + 2F3  0
−−−−−−−→ 
0
3  F + 3F
2
1
0
0


3
0
0
0
0
1
0
0
0
0
1
0
−15
−6 
.
3
0

Las columnas pivote de A son la primera, tercera y cuarta, por tanto las matrices M y N son




−2
6
4
1 3 0 0 −15

 0
1
2
−6 .
, y N = 0 0 1 0
M=
 0
0
3
0 0 0 1
3
1 −3 −2
El producto de M y N es:


−2
6
4 
1 3 0 0
 0
1
2
 0 0 1 0
MN = 
 0
0
3
0 0 0 1
1 −3 −2



−2 −6
6
4
6
−15
 0
0
1
2
0
.
−6 = 
 0
0
0
3
9
3
1
3 −3 −2 −3
¿Por qué es el resultado de este producto igual a la matriz A de partida? La explicación es
muy sencilla. Por un lado, al multiplicar M por una cualquiera de las columnas pivote de N (que
son columnas de una matriz identidad) da como resultado la correspondiente columna pivote de
A. Esto explica por qué se recuperan las columnas pivote y colocadas en su sitio. Pero también
es claro que cada columna no pivote de N consiste en los coeficientes que necesitamos en una
combinación lineal de las columnas pivote de A para que el resultado sea la correspondiente
columna no pivote de A. En consecuencia, el cálculo de la forma escalonada reducida de una
matriz nos proporciona una factorización de rango máximo de esa matriz:
Toda matriz de rango r > 0 es igual al producto de la matriz formada por sus columnas pivote por la
matriz formada por las filas no nulas de su forma escalonada reducida.
Ejemplo.
La matriz

1
A = 0
1
2
2
2
12
3
4
3

4
6
4
3. Matrices y Determinantes
3.6. Matrices por bloques
tiene claramente rango 2 y sus dos primeras columnas son linealmente independientes, por lo
tanto se tiene

 

1 2 3 4
1 2 1 0 x1 y1
A = 0 2 4 6 = 0 2
0 1 x2 y2
1 2 3 4
1 2
x1
y1
donde las columnas
e
son respectivamente las soluciones (¡únicas!) de los sistemas
x2
y2

1
0
1

 
3
2 x
1
2
= 4 ,
x2
2
3

Resolviéndolos:
3.6.
1
0
1
2
2
2
3
4
3


 
2 4
y
1
2
= 6 .
y2
2
4
1
0
1
 
4
1
6 = 0
4
1

2 1
2
0
2
0
1
−1 −2
.
2
3
Matrices por bloques
En la introducción de matrices como tablas de números que representaban sistemas de ecuaciones y en el método de las operaciones elementales para su resolución predominaba la idea
de matriz como conjunto de filas, una para cada ecuación. Por otra parte, la definición del producto de matrices y de la matriz traspuesta se ha basado en considerar una matriz como un
conjunto ordenado de columnas y no simplemente como un rectángulo de números. En otras
palabras, hemos estado considerando una matriz como una fila de matrices columna. Esto se
podría llamar “una descomposición de una matriz en bloques-columna”. En esta sección vamos
a ver que la descomposición de las matrices en otros tipos de bloques puede ser enormemente
útil. En general, podemos partir tanto las filas como las columnas de una matriz m × n, A, para
obtener una descomposición de A en submatrices A I J , donde I = 1, 2, . . . , p y J = 1, 2, . . . , q:

A=




A11
A21
..
.
A p1
A12
A22
..
.
A p2
...
A1q
A2q
..
.
A pq
n1
n2
...
nq
...
...





m1
m2
..
.
mp
donde la submatriz A I J tiene m I filas y n J columnas, es decir tamaño m I × n J .
Esta descomposición se llama descomposición por bloques de la matriz A. Los números m1 ,
m2 , . . . , m p y n1 , n2 , . . . , nq que indican los tamaños de los bloques verifican, evidentemente, que
m1 + m2 + · · · + m p = m
y
n1 + n2 + . . . nq = n.
La submatriz A I J es el bloque I J de esta descomposición por bloques de A.
Operaciones con matrices por bloques
Lo interesante de operar con matrices por bloques es que podemos actuar como si cada
bloque fuese un número. Entonces las operaciones con matrices por bloques se reducen a las
operaciones con matrices que ya conocemos, excepto que hay que asegurarse de que cuando
vayamos a sumar o multiplicar dos bloques sus tamaños sean compatibles para la operación
que queremos realizar.
13
3.6. Matrices por bloques
3. Matrices y Determinantes
Multiplicación por escalares El caso de multiplicación de una matriz por bloques por un
escalar no ofrece dificultad ya que todo bloque, independientemente de su tamaño, se puede
multiplicar por cualquier número.
Suma de matrices por bloques Para poder sumar dos matrices por bloques, es necesario no
sólo que sean matrices del mismo tamaño sino también que estén divididas en bloques de la
misma forma de tal manera que bloques correspondientes sean del mismo tamaño y se puedan
sumar.
Producto de matrices por bloques Para hallar el producto AB de dos matrices por bloques
se puede usar la regla usual de “fila por columna” siempre que el número de bloques en cada
“fila de bloques” de A sea igual al número de bloques en cada “columna de bloques” de B y
además que los bloques en esa fila de bloques de A sean compatibles para multiplicación por
los bloques de la columna de bloques de B.
Una de las consecuencias de la multiplicación de matrices por bloques es la siguiente alternativa a la regla “fila por columna” para la multiplicación de matrices:
Regla columna por fila: Sean a1 , . . . , an las columnas de A y sean b1 , . . . , bn las filas de B.
Entonces el producto matricial a1 b1 es una matriz con tantas filas como A y tantas columnas
como B. Lo mismo ocurre con los demás productos ai bi y el producto de matrices AB es igual
a la suma:
AB = a1 b1 + · · · + an bn
Ejemplo.
Sea
A=
A1
A3
A2
A4
una matriz cuadrada de orden 2n descompuesta en cuatro bloques n × n, A1 , A2 , A3 , A4 . Se
trata de encontrar una matriz
P1 P2
P=
P3 P4
con P1 , P2 , P3 , P4 bloques n × n, verificando que el producto PA intercambie la primera y la
segunda fila de bloques de la matriz A, es decir, tal que
PA =
A3
A1
A4
.
A2
Dado que el intercambio de filas es una operación elemental, basta realizar esta operación
sobre la matriz identidad por bloques:
I 0
0 I
, P=
0 I
I 0
Efectuando el producto PA por bloques obtenemos
P1 A1 + P2 A3 = 0A1 + I A3 = A3
P1 A2 + P2 A4 = 0A2 + I A4 = A4
P3 A1 + P4 A3 = I A1 + 0A3 = A1
P3 A2 + P4 A4 = I A2 + 0A4 = A2
y se verifica lo pedido.
14
3. Matrices y Determinantes
3.6. Matrices por bloques
5 Ejercicio de tarea. En las condiciones del ejemplo anterior hállense:
(a) Una matriz P tal que PA sea igual al resultado de multiplicar la primera fila de bloques
de A por la izquierda por una matriz inversible X de orden n.
(b) Una matriz P tal que PA sea igual al resultado de sumar a la segunda fila de bloques de
A la primera fila de bloques multiplicada por la izquierda por una matriz cuadrada X de
orden n.
Inversa de una matriz por bloques
En general no es sencillo calcular la inversa de una matriz por bloques realizando solamente
operaciones por bloques, pero hay un caso especial en el que sí es posible: Es el caso de una
matriz partida en 2 × 2 bloques, que sea triangular por bloques y tal que los bloques en la
diagonal sean cuadrados.
Sean A y C matrices cuadradas (no necesariamente del mismo tamaño) y sea B una matriz
con el mismo número de filas que A y el mismo número de columnas que C, de forma que se
puede formar la matriz por bloques A0 CB , que resulta ser una matriz cuadrada. Entonces
A
0
B
C
−1
=
A −1
0
− A−1 BC −1
.
C −1
(3.3)
Es sencillo demostrar la fórmula (3.3) por medio de la multiplicación directa, es decir, calculando el producto
−1
A
− A−1 BC −1
A B
.
0 C
0
C −1
6 Ejercicio de tarea. Usando las propiedades de la traspuesta deduce la siguiente fórmula a partir
de (3.3):
−1 A 0
A −1
0
=
.
(3.4)
B C
−C −1 BA−1 C −1
7 Ejercicio de tarea. Usa las fórmulas (3.3) y (3.4) según convenga para calcular la siguiente matriz inversa:
 −1

1 2 −2 0 0
3 7
3 0 0


0 0
1 0 0

 .
1 2 −1 2 1
3 2
5 3 2

7

 −3
0
 13
−27
Solución: 

−2
1
0
−4
8
20
−9
1
45
−91
0
0
0
2
−3

0

0
0

−1
2
Enlaces a todos los ejercicios de tarea de esta sección
Usa los siguientes enlaces para visualizar cada uno de los ejercicios de tarea que aparecen
en esta sección:
Enlaces: Ejercicio 5, Ejercicio 6, Ejercicio 7.
15
3.7. Determinantes
3. Matrices y Determinantes
3.7.
Determinantes
3.7.1.
Definición de determinante
Para calcular el determinante de una matriz cuadrada de orden n tenemos que saber elegir
n elementos de la matriz de forma que tomemos solo un elemento de cada fila y de cada columna como
se ilustra en la siguiente figura:







2
7
7
8
6
1
6
9
0
3
3
8
5
6
9
4
9
3
4
7
5
0
1
2
4




.


(3.5)
Además, hacer una elección de elementos de esta forma no es suficiente. Hay, en general, muchas formas de elegir n elementos de forma que se cumpla dicha condición y es necesario
realizar todas las elecciones posibles.
Ejercicio (a) ¿Cuál es el número total de dichos productos en una matriz de tamaño 3 × 3?
(b) ¿Cuál es la fórmula general que nos da el número de dichos productos para una matriz de
orden n?
Solución: (a) 6. (b) n!
Para cada elección de n elementos elegidos como se ha dicho, tenemos que hallar el producto
todos ellos y volver a elegir otros n elementos tomando sólo uno de cada fila y uno de cada
columna y hallar su producto y así sucesivamente.
Definición de
determinante.
Definición. El determinante de una matriz cuadrada es la suma equilibrada de todos esos posibles
productos que se pueden formar eligiendo un elemento de cada fila de forma que estén todos
en distintas columnas.
Significado de
“suma
equilibrada”.
¿Qué es una “suma equilibrada”? La expresión “suma equilibrada” en la definición de determinante significa que cada uno de esos productos va multiplicado por +1 o por −1 según cómo
se hayan elegido los factores. Concretamente, cada uno de esos productos va multiplicado por
(−1) p siendo p el número de intercambios de filas y columnas que sean necesarios para colocar
todos los factores de ese producto en la diagonal.
Por ejemplo, para el producto de los elementos elegidos en (3.5), el número de intercambios
de filas o columnas es p = 3 ya que se necesitan tres intercambios:

1
2
3
4
5
6

9

0
7
8
9
0
7
5
3
8
6
3
6
9


6
7
1

9

0
3
4
7
5
3
4

 F1 ↔ F2
 −−−−−−→

2 
2
1
8
6
4
7
4
3
6
9
7

8
9
0

5




2 
1
4

6
7
8
9
0
1
F3 ↔ F5 
−−−−−−→ 
3
0
2
3
4
5
6
9
7
8
6
4

 F4 ↔ F5
4 
 −−−−−−→
2 
9
7
5
3
1
16


6
7
8
9
0
1

3

9
2
3
4
5
6
9
7
7
5
3


4 

1 
0
8
6
4
2
3. Matrices y Determinantes
3.7. Determinantes
por lo tanto el sumando correspondiente a esta elección de elementos es:
(−1)3 × 6 × 2 × 9 × 3 × 2.
Si todos los factores se han elegido sobre la diagonal, p = 0 y el producto va multiplicado
por (−1)0 = 1. El ejemplo típico de esto es el de cualquier matriz identidad. En ella hay un único
producto de elementos elegidos de distintas filas y columnas que no da cero: es el producto de
los elementos de la diagonal. Por tanto tenemos:
Matriz identidad.
orden n,
El determinante de una matriz identidad es 1. Si In es la matriz identidad de
det In = 1 ,
det(− In ) = (−1)n .
Matriz diagonal. Igual que en una matriz identidad, en una matriz diagonal el producto de
los elementos de la diagonal es el único producto de elementos elegidos de distintas filas y
columnas que no da cero. Por lo tanto:

d1

0
det 
.
 ..
0

··· 0
.. 
..
.
.
 = d1 · · · d n .

..
..
. 0
.
· · · 0 dn
0
..
.
Determinante de una matriz 2 × 2. El ejemplo general no trivial más sencillo de cálculo de
determinantes es el del determinante de una matriz 2 × 2, A = ac db . En este caso sólo hay dos
productas posibles: ad y bc. El primero llevará signo positivo y el segundo negativo, por tanto:
a
det
c
3.7.2.
b
d
= ad − bc.
Consecuencias inmediatas de la definición
Determinante de la matriz traspuesta. Dado que la definición de determinante es simétrica respecto
a las filas y columnas, el valor del determinante de una matriz es igual al de su traspuesta.
det( AT ) = det A.
Consecuencia: Toda propiedad de los determinantes enunciada en términos de filas y columnas se cumple también al sustituir cada ocurrencia de la palabra “fila” por “columna” y cada
ocurrencia de la palabra “columna” por “fila”.
Matriz con una fila o columna de ceros. Si todos los elementos de una fila o de una columna son
cero el determinante es cero. (Puesto que en cada uno de los productos hay un factor igual a cero.)
Por ejemplo:

a
det d
g
b
e
h
17

0
0 = 0.
0
3.7. Determinantes
3. Matrices y Determinantes
Efecto de un intercambio de filas o columnas. Si se intecambian las posiciones de dos filas o de dos
columnas de una matriz, se cambia el signo de su determinante (ya que se cambia el signo de cada
sumando en la “suma equilibrada”).
Consecuencia 1: Si Pjk es la matriz elemental que intercambia las filas j y k,
(a) det( Pjk A) = − det A.
(b) det Pjk = −1.
(c) det( Pjk A) = det Pjk · det A.
Consecuencia 2: Si una matriz cuadrada tiene dos filas o dos columnas iguales, su determinante es cero
(ya que intercambiando esas dos filas o columnas el determinante no cambia y a la vez cambia
de signo). Por ejemplo:


a b a
det d e d  = 0.
g h g
Trasladar una fila arriba o abajo o una columna a derecha o izquierda. Trasladar una fila
hacia arriba atravesando k filas es colocarla justo encima del bloque formado por las k filas
inmediatamente anteriores sin alterar el orden de esas k filas. Esto es equivalente a realizar
k intercambios de filas y por tanto el determinante queda multiplicado por (−1)k . Lo mismo
ocurre al trasladar una fila atravesando k filas hacia abajo o cuando esto se hace con columnas
en lugar de filas.
Si se traslada una fila a derecha o izquierda atravesando k filas o se tralada una columna arriba o abajo
atravesando k columnas el determinante queda multiplicado por (−1)k .
8 Ejercicio de tarea. Usa la propiedad anterior para demostrar que si en una matriz se traslada
todo un bloque de h filas a derecha o izquierda atravesando un bloque de k filas el determinante
queda multiplicado por (−1)hk . Usa este resultado para calcular
0 In
det
.
Im 0
Efecto de un reescalado de una fila o columna. Si se multiplican todos los elementos de una fila o
de una columna de una matriz por un número, se multiplica su determinante por ese número.
Consecuencia 1: Si Eλ es una matriz elemental de reescalado por el escalar λ,
(a) det( Eλ A) = λ det A.
(b) det Eλ = λ.
(c) det( Eλ A) = det Eλ · det A.
Consecuencia 2: Si una matriz cuadrada tiene dos filas o dos columnas que son una múltiplo de la otra,
su determinante es cero (ya que es un múltplo del determinante de una matriz con dos filas o
columnas iguales).
Por ejemplo:

a
det d
g
b
e
h


3a
a
3d  = 3 det d
3g
g
b
e
h

a
d  = 0.
g
Múltiplo escalar de una matriz. Si una matriz cuadrada de orden n la multiplicamos por un número
p, todas las columnas quedan multiplicadas por p, luego el determinante de la matriz queda multiplicado
por p n veces, es decir, el determinante queda multiplicado por pn .
Esta propiedad nos permite contestar a la siguiente pregunta:
18
3. Matrices y Determinantes
Ejercicio:
5
 0
Sea A = 
 −5
0
−7
3
−8
5
2
0
0
0

2
−4
.
3
−6
Sabiendo que det A = 20
calcular
det(3A).
Solución: det(3A) = 34 × 20 = 1620

3.7. Determinantes
Matriz triangular. Si todos los elementos encima o debajo de la diagonal son cero, todos los productos
en los que intervenga un elemento fuera de la diagonal son cero y el determinante es igual al producto de
los elementos de la diagonal.

3.7.3.
a11
0

 a21
det 
 .
 ..
a22
..
.
am1
am2


···
0
a11

.. 
..

.
. 
 = a11 · · · amm = det  0
 ..

..
 .
.
0 
0
· · · amm

a21
···
a22
..
.
···

· · · am2 

..  .
..
.
. 
0 amm
am1
Cálculo del determinante por desarrollo de cofactores de una fila o
columna
Desarrollo por la primera columna. Si queremos calcular un determinante, necesitamos formar todos los posibles productos que se pueden formar eligiendo un elemento de cada fila y de
cada columna. Esto se puede hacer ordenadamente de la siguiente forma: Primero formamos
todos los productos en los que aparece el primer elemento de la primera columna, a11 . La suma
de todos estos productos (cada uno con su signo) es igual a a11 multiplicado por el determinante
de la matriz A11 que se obtiene al eliminar en la original la primera fila y la primera columna, es
decir (ver más abajo), el determinante del menor del elemento (1, 1). Después formamos todos
los productos en los que interviene el segundo elemento de la primera columna, a21 . La suma
equilibrada de estos productos es igual a − a21 multiplicado por el determinante de la matriz
A21 que se obtiene al eliminar en la original la segunda fila y la primera columna (el menor del
elemento (2, 1)). Continuando de esta manera, vemos que el determinante de la matriz se puede
expresar como una suma de productos de los elementos ai1 de la primera columna, cada uno
de ellos multiplicado por el determinante de una matriz de un orden menor que la dada:
det A = a11 det A11 − a21 det A21 + · · · ± an1 det An1
(3.6)
Para entender perfectamente esta forma de calcular un determinante necesitamos introducir
algunos conceptos:
Menor de un elemento. Dada una matriz A, se llama menor del elemento que ocupa la posición
(i, j) (es decir, fila i, columna j) y se denota Aij a la matriz obtenida al eliminar toda la fila i y toda la
columna j de la matriz dada.
Cofactor de un elemento. Dada una matriz cuadrada A, se llama cofactor del elemento que ocupa
la posición (i, j) al determinante det Aij del menor de ese elemento multiplicado por +1 o −1 dependiendo
de si i + j es par o impar.
Así, el cofactor del elemento (i, j) de la matriz A se calcula por la fórmula:
Cij = (−1)i+ j det Aij .
Usando cofactores, la fórmula (3.6) se puede escribir:
det A = a11 C11 + · · · + an1 Cn1
19
3.7. Determinantes
3. Matrices y Determinantes
Los cofactores de los elementos de una fila o columna nos permiten calcular el determinante de
una matriz cuadrada mediante la fórmula de expansión del determinante por los cofactores de una
fila o columna. La expansión del determinante de una matriz n × n, A = ( aij ) por la columna j
es:
det A = a1j C1j + · · · + anj Cnj
La expansión del determinante de la misma matriz por la fila i es:
det A = ai1 Ci1 + · · · + ain Cin
3.7.4.
Consecuencias no tan inmediatas de la definición
Efecto de descomponer una fila o columna como suma de dos. Si en una matriz cuadrada se
descompone una fila o columna como suma de dos, su determinante se descompone en suma de dos (esto
se demuestra desarrollando el determinante por los elementos de dicha fila o columna).
Por ejemplo:

a
det d
g
b
e
h

c+p
f + q = (c + p)C13 + ( f + q)C23 + (k + r )C33
k+r
= (cC13 +

a
= det d
g
f C23 + kC33 ) + ( pC13 + qC23 + rC33 )



b c
a b p
e f  + det d e q  .
h k
g h r
Efecto de sumar o restar a una fila otra fila o a una columna otra columna. Si en una matriz
cuadrada se le suma o resta a una fila otra fila o a una columna otra columna, su determinante no cambia.
Por ejemplo:

a
det d
g
b
e
h


c+b
a
f + e  = det d
k+h
g
b
e
h


c
a
f  + det d
k
g
b
e
h


b
a
e  = det d
h
g
b
e
h

c
f.
k
Efecto de una operación de reemplazo de una fila o columna. Si a una fila o columna se le suma
otra multiplicada por un escalar el determinante de la matriz no cambia.
Por ejemplo:

a
det d
g
b
e
h


c + 3a
a
f + 3d = det d
k + 3g
g
b
e
h

c
f.
k
Consecuencia: Si E es una matriz elemental de reemplazo y A es una matriz cuadrada del mismo
tamaño que E,
(a) det( EA) = det A.
(b) det E = 1.
(c) det( EA) = det E · det A.
Con lo anterior se completan las consecuencias de los efectos que tiene sobre el determinante
de una matriz el realizar operaciones elementales de filas. En el siguiente ejercicio se obtiene
una importante propiedad de la que se obtendrán consecuencias en el ejercicio 11.
20
3. Matrices y Determinantes
9 Ejercicio de tarea. (a) Usando el hecho de que toda matriz unitriangular inferior (triangular
inferior con unos en la diagonal) multiplicada por la izquierda por cualquier matriz efectúa
sobre ésta una sucesión de operaciones de reemplazo de filas que no cambian el determinante,
justifica que si M es una matriz unitriangular inferior entonces
det( MA) = det A.
(3.7)
(b) Usa las propiedades de la traspuesta para deducir, a partir de (3.7), que si N es una matriz
unitriangular superior
det( AN ) = det A.
(3.8)
Pista: (a) M es el producto de matrices elementales de reemplazo de filas.
matriz
unitriangular
3.7. Determinantes
3.7.5.
Cálculo de un determinante por reducción a forma escalonada
De las propiedades enunciadas en la sección anterior se deduce que si una matriz A se transforma, mediante operaciones elementales de filas, en una matriz escalonada U y solamente se
han usado operaciones de reemplazo y de intercambio (o sea, sin usar operaciones de reescalado, lo cual, por otra parte, siempre es posible), entonces el determinante de la matriz escalonada
U es igual al determinante de A multiplicado por ±1 dependiendo de si el número de operaciones de intercambio ha sido par o impar. En otras palabras, si U es una forma escalonada de
A obtenida sin operaciones de reescalado y con r operaciones de intercambio, entonces
det A = (−1)r det U.
Al aplicar esta técnica de cálculo de un determinante no es necesario limitarse a operaciones
elementales de filas. Se pueden realizar operaciones elementales de filas y de columnas mezcladas según convenga, teniendo cuidado de contabilizar en r el número total de intercambios de
filas y columnas.
Puesto que toda matriz cuadrada escalonada es triangular, el determinante de U es igual
al producto de todos los elementos de su diagonal y en consecuencia, el determinante de A es
igual a (−1)r multiplicado por todos los elementos de la diagonal de U. Si A no es inversible
(es decir, es singular) entonces la matriz escalonada U tiene algún elemento de la diagonal igual
a cero y en consecuencia su determinante es cero. Recíprocamente si U tiene determinante cero,
algún elemento de su diagonal es cero y su número de pivotes es menor que el número de
columnas. Esto implica que A es necesariamente singular, Por tanto tenemos:
Proposición.
cero:
Una matriz es singular (no tiene inversa) si y sólo si su determinante es igual a
det A = 0
3.7.6.
si y sólo si
A es singular.
Teorema fundamental del cálculo de determinantes
Teorema. (a) De los efectos que tiene sobre el determinante de una matriz el realizar operaciones elementales de filas (véase más arriba) se deduce que para toda matriz elemental E y toda
matriz A del mismo tamaño se verifica
det( EA) = det E · det A.
(3.9)
(b) Del apartado (a) se deduce que si A y B son dos matrices cuadradas del mismo tamaño,
det( AB) = det A · det B.
21
(3.10)
3.7. Determinantes
3. Matrices y Determinantes
Demostración: Basta demostrar (b) ya que el apartado (a) se ha justificado ampliamente más
arriba. En el caso de que A no es inversible entonces AB tampoco lo es y, por la proposición de
la sección anterior, ambos miembros de la ecuación son iguales a cero.
Si A es inversible entonces es igual a un producto de matrices elementales: A = E1 · · · Ek y
en consecuencia:
det( AB) = det( E1 ) det( E2 · · · Ek B) = · · · = det E1 · · · det Ek det B = det( E1 · · · Ek ) det B = det A det B.
Corolario 1.
Si A es una matriz inversible, entonces
det( A−1 ) =
1
.
det A
(3.11)
Corolario 2. Si A y B son dos matrices semejantes (esto es, existe una matriz inversible, M, tal
que A = MBM−1 ) entonces,
det A = det B.
3.7.7.
Ejemplos de determinantes especiales
(a) Si en una matriz cuadrada de orden n todos los elementos de una columna son iguales a p, su
determinante es igual a p multiplicado por un determinante de orden n − 1. Lo mismo ocurre si
todos los elementos de una fila son iguales. Se resta la primera fila de cada una de las demás,
con lo que la columna en cuestión se convierte en ( p, 0, . . . , 0) y ahora se desarrolla el
determinante por esa columna.

a
det  d
g
b
e
h


p
a
p  = det  d − a
p
g−a
b
e−b
h−b

p
d−a
0  = (−1)1+3 p · det
g−a
0
e−b
h−b
.
(b) Si en una matriz cuadrada de orden n la suma de los elementos de una fila es igual a la de los de
otra fila cualquiera (todas las filas tienen la misma suma) entonces el determinante es igual a esa
suma multiplicada por un determinante de orden n − 1.
Para verlo, realizamos sobre la matriz original las siguientes operaciones de reemplazo:
Sumamos a la primera columna todas las demás columnas, con lo cual la primera columna
tiene todos los elementos iguales y estamos en la situación del ejemplo anterior.

a
det  b
c
3.7.8.
b
c
b


c
a+b+c
a  = det  a + b + c
a
a+b+c
b
c
b

c
c−b
a  = ( a + b + c) det
b−b
a
a−c
a−c
.
Determinantes de matrices por bloques
Sean A y C matrices cuadradas (no necesariamente del mismo tamaño) y sea B una matriz
con el mismo número de filas que A y el mismo número de columnas que C, de forma que se
puede formar la matriz por bloques A0 CB . Entonces
A
det
0
B
C
= det A · det C.
22
(3.12)
3. Matrices y Determinantes
3.7. Determinantes
Esto es una sencilla consecuencia del teorema fundamental (ver ecuación (3.10)) combinado
con la siguiente identidad:
A B
A 0
I A −1 B
=
,
0 C
0 I
0
C
pero también puede demostrarse de forma independiente del teorema fundamental usando
solamente la definición de determinante. La demostración se basa en la observación
de que al
formar los productos de elementos que forman el determinante de A0 CB , si en un producto
determinado apareciese un elemento de B, también habría uno del bloque nulo simétrico y ese
producto sería cero. Cada uno de los productos no nulos es, pues, un producto de dos factores:
uno es uno de los productos que forman det A y el otro es uno de los productos que forman
det C.
Usando las propiedades de la traspuesta es muy sencillo deducir de (3.12) su versión “traspuesta”:
A 0
det
= det A · det C.
(3.13)
B C
La fórmula (3.12) (y análogamente la (3.13)) sigue siendo cierta en un caso más general.
Supongamos que A es una matriz n × n partida en k × k bloques (no necesariamente del mismo
tamaño) y tal que los bloques de la diagonal, Aii para i = 1, . . . k, son cuadrados (siendo el
tamaño de Aii , pi × pi ) y los bloques debajo de la diagonal son todos matrices nulas,


A11 A12 . . . A1k
 0 A22 . . . A2k 


A=
.. 
..
..
..

.
.
.
.
0 ...
0 Akk
Entonces el determinante de A es el producto de los determinantes de las matrices de la diagonal:
det( A) = det( A11 ) · · · det( Akk ).
(3.14)
Esta fórmula puede deducirse de (3.12) por el método de inducción.
10 Ejercicio de tarea. Si A es una matriz inversible y M =
tra
I
0
A
A B
=
C D
−CA−1 I
0
A B
C D
es una matriz cuadrada, demues
B
D − CA−1 B
y usa el ejercicio 9-ecuación (3.7) para deducir de ello
A B
det
= det A · det( D − CA−1 B).
C D
Las operaciones de matrices por bloques son una poderosa herramienta para demostrar
algunas fórmulas que aparentemente no tienen nada que ver con las matrices por bloques.
Como ejemplo, el siguiente ejercicio demuestra la fórmula (3.10).
11 Ejercicio de tarea. El objetivo de este ejercicio es demostrar el teorema fundamental del cálculo
de determinantes,
det( AB) = det A · det B,
usando solamente las propiedades de las operaciones de matrices por bloques. Para ello hacemos los siguientes pasos donde A, B, C y D son matrices cuadradas del mismo tamaño e I es la
matriz identidad del mismo tamaño que A, B, etc.:
23
3.8. Regla de Cramer y fórmula de la matriz inversa
3. Matrices y Determinantes
(a) Demuestra que
−C − D
.
A
B
I
I
0
I
I
0
−I
I
A
C
B
D
= det
−I
I
I
0
B
D
A
C
Pista: Calcula el producto
det
(b) Usa el resultado anterior combinado con (3.13) (supuesto demostrado independientemente
del teorema fundamental) para deducir
A B
det
= det(−C ) det B.
C 0
(c) Usa la multiplicación de matrices por bloques para comprobar la identidad
A 0
I B
A AB
=
−I B
0 I
−I
0
y deduce de ella, del apartado anterior y del ejercicio 9-ecuación (3.8), que
A 0
det( AB) = det
.
−I B
(d) Combina (3.13) con el resultado anterior para deducir la fórmula buscada.
Enlaces a todos los ejercicios de tarea de esta sección
Usa los siguientes enlaces para visualizar cada uno de los ejercicios de tarea que aparecen
en esta sección:
Enlaces: Ejercicio 8, Ejercicio 9, Ejercicio 10, Ejercicio 11.
3.8.
Regla de Cramer y fórmula de la matriz inversa
La regla de Cramer es un método que nos permite escribir una fórmula para cada incógnita
de un sistema de ecuaciones lineales cuya matriz de coeficientes es una matriz inversible. Dichas
fórmulas nos permiten calcular la solución de dicho sistema por medio de determinantes. Para
explicar la fórmula que corresponde a cada incógnita es necesario comprender una operación
especial que se puede realizar sobre una matriz de m filas y que consiste en sustituir una
columna determinada por un vector dado de Rm .
La operación de sustituir una columna de una matriz por un vector dado
Supongamos que A es una matriz de m filas y n columnas, y que b es un vector de Rm .
Denotamos A( j←b) la matriz obtenida a partir de A al sustituir la columna j de A por el vector
b. Por ejemplo, si


 
1 0
a
A = 2 4
y
b = b
5 9
c
entonces

A (1← b )
a
= b
c

0
4,
9

A (2← b )
1
= 2
5

a
b.
c
En relación con esta operación de matrices el cálculo de determinantes tiene las siguientes
propiedades:
24
3. Matrices y Determinantes
3.8. Regla de Cramer y fórmula de la matriz inversa
(a) Si x es el vector de componentes ( x1 , . . . , xn ), entonces
det I( j←x) = x j .
Por ejemplo:

Si

x1
x =  x2  ,
x3

det I(3←x)
entonces
1 0
= det  0 1
0 0

x1
1
x2  = x3 det
0
x3
0
1
= x3 .
(b) Para cualquier vector x ∈ Rn , el producto de A por la matriz I( j←x) es:
A I( j←x) = A( j← Ax) .
Por ejemplo:

A I( j←x)
1
= 2
5

0 1
4
0
9
x1
x2


1
1
x1
= A
A
= 2 Ax = A(2← Ax)
0
x2
5
Combinando esta propiedad con la anterior deducimos:
det A( j← Ax) = x j det A.
(3.15)
(c) Si A es una matriz cuadrada de orden n y e j es la columna j de la matriz identidad In ,
entonces el cofactor del elemento ( j, i ) de A es:
Cji = det A(i←e j ) .
Por ejemplo,

1
A = 2
5
Si

3
7
6
0
4
9

entonces
A (3← e2 )
1
= 2
5
0
4
9

0
1
0
y por tanto

det A(3←e2 )
1
= det  2
5
0
4
9

0
1
2+3

1 = (−1)
det
5
0
0
9
= − det A2 3 = C2 3 .
La Regla de Cramer
En términos de la operación estudiada más arriba, la regla de Cramer se enuncia así:
Regla de Cramer. Si A es una matriz n × n inversible y b ∈ Rn , entonces el vector x =
( x1 , . . . , xn ) cuyos elementos están definidos por:
det A( j←b)
xj =
det A
es la solución única del sistema compatible determinado Ax = b.
Demostración: Si A es una matriz inversible y Ax = b entonces A( j←b) = A( j← Ax) . Tomando
determinantes y usando (3.15): det A( j←b) = x j det A.
25
3.9. Ejercicios adicionales
3. Matrices y Determinantes
Una fórmula para la matriz inversa
Si A es una matriz inversible n × n, la columna j de A−1 es un vector x j que satisface el
sistema de ecuaciones Ax j = e j (donde, igual que antes, usamos la notación e1 , . . . , en para las
columnas de la matriz identidad de orden n). En consecuencia, por la regla de Cramer, x j está
dado por




det A(1←e j )
Cj1
1 
.. 
.

= 1 
xj =
 ..  .
.

det A
det A

Cjn
det A(n←e j )
Ahora bien, la columna de cofactores que aparece en esta fórmula tiene por elementos los de la
fila j de la matriz de cofactores de A. En conclusión,
La matriz inversa de A es igual al inverso del determinante de A multiplicado por la traspuesta de
la matriz de cofactores de A:
T

C11 · · · C1n
1  . .
..  .
..
A −1 =
 ..
.
det A
Cn1 · · · Cnn
3.9.
Ejercicios adicionales
Operaciones con matrices
I 1. Dado que los vectores en Rn pueden ser considerados comomatrices n × 1, las propiedades
1 −3
5
de las traspuestas también se aplican a vectores. Sean A =
y x =
Calcula
−2
4
3
( Ax)T , x T A T , xx T , y x T x. ¿Está definido el producto A T x T ?
I 2. Sean A una matriz 4 × 4 y x un vector en R4 . ¿Cuál es la forma más rápida de calcular A2 x:
Haciendo A( Ax) o haciendo ( A · A)x?. Cuenta las multiplicaciones que hay que hacer en cada
caso.
En los ejercicios 3 y 4, calcula cada suma o producto si la matriz está definida. Si alguna
expresión no está definida, explica por qué.
Sean
1 2
C=
,
−2 1
3 5
D=
,
−1 4
−1
A=
,
2
7 −5
1
B=
,
1 −4 −3
I 3. −2A,
2
4
0
−5
B − 2A,
AC,
CD.
I 4. A + 2B,
E=
3C − E,
−5
.
3
CB,
EB.
En el resto de esta serie de ejercicios y en las series que siguen, debe suponerse que cada
expresión de matrices está definida. Esto es, los tamaños de las matrices (y de los vectores)
involucrados “se corresponden” de manera apropiada.
26
3. Matrices y Determinantes
I 5. Calcula 3I2 − A y (3I2 ) A con
A=
4
5
−1
.
−2
3.9. Ejercicios adicionales
I 6. Calcula A − 5I3 y (5I3 ) A, con


9 −1
3
7 −6 .
A =  −8
−4
1
8
En los ejercicios 7 y 8, calcula el producto AB en dos formas: (a) mediante la definición, donde Ab1 y Ab2 se calculan por separado, y (b) mediante la regla fila-por-columna para calcular
AB.




−1
2
4 −2
3 −2
1
3
4, B =
0, B =
I 7. A =  5
I 8. A = −3
−2
1
2 −1
2 −3
3
5
I 9. Si una matriz A es de orden 5 × 3 y el producto AB es de orden 5 × 7, ¿cuál es el orden de B?
I 10. ¿Cuántas filas tiene B si BC es una matriz de orden 3 × 4?.
I 11. Sean A
2
−3
I 12. Sean A =
5
4
,yB=
1
3
2
−4
−5
. ¿Qué valor(es) de k, si hay, hacen que AB = BA?.
k
−3
8 4
5 −2
,B=
,yC=
6
5 5
3
1
Comprueba que AB = AC y que sin embargo B 6= C.




1 1 1
2 0 0
I 13. Sean A = 1 2 3 y D = 0 3 0. Calcula AD y DA. Explica cómo cambian las filas
1 4 5
0 0 5
o columnas de A cuando se multiplica por D a la derecha o a la izquierda. Halla una matriz B
de orden 3 × 3, que no sea la matriz identidad o la matriz cero, tal que AB = BA.
3 −6
I 14. Sea A =
. Construye una matriz B de orden 2 × 2 tal que AB sea igual a la matriz
−1
2
cero. Las columnas de B no deben ser iguales entre sí y deben ser distintas de cero.
I 15. Sean r1 , . . . , r p vectores en Rn , y sea Q una matriz de orden m × n. Escribe la matriz
[ Qr1 . . . Qr p ] como un producto de dos matrices sin usar una matriz identidad.
Los ejercicios 16 y 17 tratan de matrices arbitrarias A, B y C para las cuales las sumas y
productos indicados están definidos. Indica para cada una de las siguientes afirmaciones si es
verdadera o falsa. Justifica tus respuestas.
I 16.
(a) Si A y B son matrices de orden 2 × 2 con columnas a1 , a2 y b1 , b2 , respectivamente,
entonces AB = [a1 b1 a2 b2 ].
(b) Toda columna de AB es una combinación lineal de las columnas de B usando como coeficientes los elementos de la columna correspondiente de A.
(c) AB + AC = A( B + C )
27
3.9. Ejercicios adicionales
3. Matrices y Determinantes
(d) A T + B T = ( A + B) T
(e) La traspuesta de un producto de matrices es igual al producto de sus traspuestas en el
mismo orden.
I 17.
(a) Si A y B son matrices 3 × 3 y B = [b1 b2 b3 ], entonces AB = [ Ab1 + Ab2 + Ab3 ].
(b) La segunda fila de AB es la segunda fila de A multiplicada a la derecha por B.
(c) ( AB)C = ( AC ) B
(d) ( AB) T = A T B T
(e) La traspuesta de una suma de matrices es igual a la suma de sus traspuestas.
I 18. Si A =
1
−2
−2
5
y AB =
−1
6
2
−9
−1
, halla la primera y la segunda columna de B.
3
I 19. Supongamos que las dos primeras columnas de B son iguales. ¿Qué puede decirse acerca de
las columnas de AB (suponiendo que este producto está definido)?. ¿Por qué?
I 20. Supongamos que la tercera columna de B es la suma de las primeras dos columnas. ¿Qué
puede decirse acerca de la tercera columna de AB? ¿Por qué?
I 21. Supongamos que la segunda columna de B es toda cero. ¿Qué puede decirse acerca de la
segunda columna de AB?
I 22. Supongamos que la última columna de AB es completamente cero, pero B por sí sola no
tiene ninguna columna de ceros. ¿Qué puede decirse acerca de las columnas de A?
I 23. Demuestra que si las columnas de B son linealmente dependientes, también lo son las columnas de AB.
I 24. Supongamos que CA = In (la matriz identidad n × n). Demuestra que la ecuación Ax = 0
tiene únicamente la solución trivial. Explica por qué A no puede tener más columnas que filas.
I 25. Supongamos que AD = Im , (la matriz identidad m × m). Demuestra que para todo b en Rm ,
la ecuación Ax = b tiene una solución. [Sugerencia: Piensa en la ecuación ADb = b.] Explica
por qué A no puede tener más filas que columnas.
I 26. Supongamos que A es una matriz de orden m × n y que existen matrices n × m, C y D, tales
que CA = In y AD = Im . Demuestra que m = n y C = D. [Sugerencia: Piensa en el producto
CAD).]
I 27. Supongamos que A es una matriz de orden 3 × n cuyas columnas generan R3 . Explica cómo
construir una matriz D de orden n × 3 tal que AD = I3 .
En los ejercicios 28 y 29, considera los vectores en Rn como matrices n × 1. Para u y v en
el producto de matrices u T v es una matriz 1 × 1, llamada producto escalar, o producto interno,
de u y v. Por lo general, se escribe como un único número real sin paréntesis o corchetes. El
producto de matrices uv T es una matriz de orden n × n, llamada producto exterior de u y v.
Rn ,
28
3. Matrices y Determinantes
3.9. Ejercicios adicionales


 
−2
a
I 28. Sean u =  3 y v = b . Calcula u T v, v T u, u v T , v u T .
−4
c
I 29. Si u y v están en Rn , ¿qué relación hay entre u T v y v T u? ¿Y entre u v T y v u T ?
I 30. Demuestra que Im A = A cuando A es una matriz de orden m × n. Puedes utilizar el hecho
de que Im x = x para todo x en Rm .
I 31. Demuestra que AIn = A cuando A es una matriz de orden m × n. [Sugerencia: Usa la
definición (de columnas) del producto de matrices AIn .]
I 32. Halla una fórmula para ( ABx)T , donde x es un vector y A y B son matrices con los tamaños
apropiados.
I 33. Calcula Sk para k = 2, . . . , 6 donde

0
0

S=
0
0
0
1
0
0
0
0
0
1
0
0
0
0
0
1
0
0

0
0

0
.
1
0
I 34. (Para hacer con Mathematica en una práctica de laboratorio.) Describe con palabras qué
pasa al calcular A5 , A10 , A20 y A30 para


1/6 0,5 1/3
A =  0,5 1/4 1/4
1/3 1/4 5/12
Matrices inversas


1 2
I 1. Sea A =  1 3 . Construye una matriz C de 2 × 3 (mediante ensayo y error) usando sólo
1 5
los números 1, −1 y 0 como elementos, de tal forma que CA = I2 . Calcula AC y observa que
AC 6= I3 .
1 1 1 0
. Construye una matriz D de 4 × 2 usando sólo los números 1 y 0
0 1 1 1
como elementos, de tal forma que AD = I2 . ¿Es posible que CA = I4 para alguna matriz C de
orden 4 × 2?. ¿Por qué sí o por qué no?.
I 2. Sea A =
En los ejercicios 3 a 6 halla las inversas de las matrices dadas.
8 6
3 2
8
5
I 3.
I 4.
I 5.
5 4
7 4
−7 −5
I 6.
I 7. Usa la inversa de la matriz del ejercicio 3 para resolver el sistema
8x1 + 6x2 = 2
5x1 + 4x2 = −1
I 8. Usa la inversa de la matriz del ejercicio 5 para resolver el sistema
8x1 + 5x2 = −9
−7x2 − 5x2 = 11
I 9. Sean A =
1
5
2
−1
1
2
3
, b1 =
, b2 =
, b3 =
, b4 =
.
12
3
−5
6
5
29
3
7
−4
−8
3.9. Ejercicios adicionales
3. Matrices y Determinantes
(a) Halla A−1 y utilízala para resolver las cuatro ecuaciones
Ax = b1 , Ax = b2 , Ax = b3 , Ax = b4
(b) Las cuatro ecuaciones del apartado (a) pueden resolverse con el mismo conjunto de operaciones de fila, puesto que la matriz de coeficientes es la misma en cada caso. Resuelve las cuatro ecuaciones del apartado (a) reduciendo por filas la matriz aumentada
[ A b1 b2 b3 b4 ] para hallar su forma escalonada reducida.
I 10. Usa el álgebra de matrices para mostrar que si A es inversible y D satisface AD = I, entonces
D = A −1 .
En los ejercicios 11 y 12, indica para cada afirmación si es verdadera o falsa. Justifica tus
respuestas.
I 11.
(a) Para que una matriz B sea inversa de A, ambas ecuaciones AB = I y BA = I deben ser
ciertas.
(b) Si A y B son matrices n × n inversibles, entonces A−1 B−1 es la inversa de AB.
a b
(c) Si A =
, y ab − cd = 0, entonces A es inversible.
c d
(d) Si A es una matriz inversible n × n, entonces la ecuación A x = b es consistente para toda
b en Rn .
(e) Toda matriz elemental es inversible.
I 12.
(a) Un producto de matrices n × n inversibles es inversible, y la inversa del producto es el
producto de sus inversas en el mismo orden.
(b) Si A es inversible, entonces la inversa de A−1 es la propia A.
a b
(c) Si A =
, y ad = bc, entonces A no es inversible.
c d
(d) Si A se puede reducir por filas a la matriz identidad, entonces A es inversible.
(e) Si A es una matriz inversible n × n, entonces las operaciones elementales de filas que
reducen A a la identidad In también reducen A−1 a In .
I 13. Sea A una matriz inversible de n × n y sea B una matriz n × p. Demuestra que la ecuación
AX = B tiene una única solución A−1 B.
I 14. Sea A una matriz inversible n × n, y sea B una matriz n × p. Explica por qué A−1 B puede
calcularse mediante reducción por filas:
Si
[ A B] ∼ · · · ∼ [ I X ] entonces X = A−1 B
Observación: Si A tiene orden más grande que 2 × 2, entonces la reducción por filas de [ A B] es mucho más rápida
(requiere menos operaciones) que calcular A−1 y luego el producto A−1 B.
30
3. Matrices y Determinantes
3.9. Ejercicios adicionales
I 15. Supongamos que AB = AC, donde B y C son matrices n × p y A es inversible. Demuestra
que B = C. ¿Es esto cierto en general si A no es inversible?
I 16. Supongamos ( B − C ) D = 0, donde B y C son matrices m × n y D es inversible. Demuestra
que B = C.
I 17. Supongamos que A, B y C son matrices inversibles n × n. Demuestra que ABC también es
inversible construyendo una matriz D tal que ( ABC ) D = I y D ( ABC ) = I.
I 18. Supongamos que A y B son matrices n × n, y que AB es inversible. Demuestra que A y B
son inversibles.
Sugerencia: Pon C = AB deduce qué se debe multiplicar a derecha/izquierda de C −1 para obtener las inversas requeridas
y luego demuestra que efectivamente son las inversas requeridas...
I 19. Despeja A en la ecuación AB = BC suponiendo que A, B y C son cuadradas y que B es
inversible.
I 20. Supongamos que P es inversible y A = PBP−1 . Halla B en términos de A.
I 21. Si A, B y C son matrices inversibles n × n, la ecuación C −1 ( A + X ) B−1 = In ¿tiene alguna
solución para X? Si es así, hállala.
I 22. Supongamos que A, B y X son matrices n × n con A, X, y A − AX inversibles, y supongamos
que
( A − AX )−1 = X −1 B
(3.16)
(a) Explica por qué B es inversible.
(b) Resuelve la ecuación (3.16) en X. Si es necesario usar la inversa de una matriz, explique
por qué dicha matriz es inversible.
I 23. Explica por qué las columnas de una matriz A de n × n son linealmente independientes
cuando A es inversible.
I 24. Explica por qué las columnas de una matriz A de orden n × n generan Rn cuando A es
inversible.
I 25. Supongamos que A es n × n y que la ecuación A x = 0 tiene solamente la solución trivial.
Explica por qué A tiene n columnas pivote y es equivalente por filas a In . Observación: Esto implica
que A debe ser inversible.
I 26. Supongamos que para una matriz cuadrada A de orden n la ecuación Ax = b tiene una
solución para todo b en Rn . Explica por qué A debe ser inversible.
Sugerencia: Piensa si A es equivalente por filas a In
Los ejercicios 27 y 28 demuestran el teorema que dice:
a b
Si A =
entonces A es inversible sólo si ad − bc 6= 0 y en ese caso su inversa es
c d
1
d −b
−1
A =
.
a
ad − bc −c
I 27. Demuestra que si ad − bc = 0, entonces la ecuación Ax = 0 tiene más de una solución. ¿Por
qué implica esto que A no es inversible?
Sugerencia: Primero, considera el caso a = b = 0. Después, si a y b no son ambos cero, considera el vector x =
31
−b
a
.
3.9. Ejercicios adicionales
3. Matrices y Determinantes
I 28. Demuestra que si ad − bc 6= 0, la fórmula para A−1 es correcta.
Halla las inversas de las matrices dadas en los ejercicios 31 a 34, caso de que existan. Usa
el algoritmo explicado en clase consistente en hallar la forma escalonada reducida de la matriz
ampliada [ A I ].
I 29.
1
4

2
7
1
I 31.  −3
2
I 30.
0
1
−3

−2
4
4
5
4
10
7

1
I 32.  4
−2

−2
1
−7
3
6 −4
I 33. Usa el algoritmo explicado en clase para hallar las inversas de




1 0 0 0
1 0 0


1 1 0 y 1 1 0 0.
1 1 1 0
1 1 1
1 1 1 1
Para cada una de las inversas halladas verifica que efectivamente es la inversa comprobando
que cumple AB = I.


−2 −7 −9
5
6 . Halla la tercera columna de A−1 sin calcular las otras columnas.
I 34. Sea A =  2
1
3
4

1
I 35. Si existe, halla la inversa de la matriz A =  −1
5
3.9.1.

−2 −1
5
6
−4
5
Operaciones elementales realizadas mediante matrices elementales
Los siguientes dos ejercicios demuestran casos especiales de las propiedades de las matrices
elementales. Aquí A es una matriz 3 × 3 e I = I3 .
I 1. Usa la ecuación filai ( AB) = filai ( A) · B para demostrar que para i = 1, 2, 3,
filai ( A) = filai ( I ) · A.
I 2. Demuestra que si las filas 1 y 2 de A se intercambian, entonces el resultado es igual a EA,
donde E es la matriz elemental formada al intercambiar las filas 1 y 2 de I.
I 3. Demuestra que si la fila 3 de A se multiplica por 5, entonces el resultado es igual a EA, donde
E es la matriz elemental formada al multiplicar la fila 3 de I por 5.
I 4. Demuestra que si la fila 3 de A es reemplazada por fila3 ( A) − 4 fila1 ( A), el resultado es igual
a EA, donde E es la matriz elemental formada a partir de I al reemplazar la fila 3 de I por
fila3 ( I ) − 4 fila1 ( I ).
32
3. Matrices y Determinantes
3.9. Ejercicios adicionales
Determinantes
En los siguientes enunciados se utilizan las barras verticales para representar un determinante, es decir, se usa la notación | A| = det A. En el caso del determinante de una matriz dada
explícitamente por sus elementos las barras verticales sustituyen a los paréntesis que normalmente delimitan los elementos de la matriz. Cada una de las ecuaciones que aparecen en los
ejercicios 1 a 4 ilustra una propiedad de los determinantes. Enuncia la propiedad que corresponda.
1 −3
2 −6
1 −3
0
2 4 6 5 −2 5 −2 5 −2 = 2 3
5 −2 6 = − 0
I 2. 3
I 1. 1 −3
4 −1
1
6
3
1
6
3
4 −1
8
8
1
1
2
3 2
3 1
3 −4 3 −4 1
5 −4 5 −4 = 0
5 0 −3 = 0 −6
I 4. 0
I 3. 2
5 −4
0
1 −5 3
7
4
5 −4
7
7
En los ejercicios 5 a 9 halla los determinantes indicados mediante reducción a forma escalonada.
1
3
0
4 1
5 −3 5 −6 −2 −5
1
7
4 3 4 I 7. I 6. 3 −3
I 5. −1 −4
3
5
2
1 2 13 −7 −2 −7
9
1 −1
2 −3 1 −1 −3
1
0 3
3 −4 0
0
1
5
4 1
2 −5 I 9. I 8. −
1
2
8
5 2
5
4
−
3
3 −1 −2
−3 −7 −5
3
2
Halla los determinantes de los ejercicios 10 a 17 mediante un desarrollo por cofactores en la
fila o columna según te parezca mejor. Indica claramente cuál es la fila o columna utilizada.
1
2 −4
0
3
3
5 3 5
1 0
4 1
1 3
2 0 1
2 I 12. 3
I 13. 2
I 11. 4 −3
I 10. 2
3
1
0
2
4
1
4 −1 4
2
5 −1 4
8
5 −2
2
1
6 3 −4 4 3
0 3 −5 5
2 0
3 0
5 I 17. 4
I 16. 6
I 14. 4
I 15. 0
3 −2
5
2 −4
9
7
3
5
7
1
6
Para cada una de las matrices dadas en los ejercicios 18 a 20 halla una forma escalonada
sin utilizar operaciones de reescalado. Realiza el menor número posible de operaciones e indica
claramente el número de intercambios de filas utilizado y el valor del determinante de la matriz
original.






6
0
0
5
1 −2
5
2
3
5 −8
4
1
0
0 −2
7
2 −5
0
3
0
3 −7



I 18. 
I 19. 
I 20. 
2




0
0
0
2 −6 −7
5
0
0
1
5
8
3
1
8
5
0
4
4
0
0
0
2
Halla los determinantes de los ejercicios 21 a 26 combinando los métodos de la traspuesta,
reducción por filas y desarrollo por cofactores, según sea más apropiado.
33
3.9. Ejercicios adicionales
2
3
I 21. −6
4
−3
1
I 24. −3
3
5
0
0
10
−3 −1 1 −3 −4
9 −4 −1 −2
1
3
0
4 −2
−4
0
−4 −3 8 4
3. Matrices y Determinantes
−1
3
I 22. 5
4
4
0
I 25. 7
5
0
2
4
4
2
0
0
3
0
0
2
5
4
4
7
6
I 23. 6
−
2
−
4
−6
7
7
6
3
2
−5 9
0
−
4
0 6
−8 I 26. 8 −5
3
0
0
−3 4
2
3
2
1 2 0 0
0 0 6 3
3
3
6
4
−7
3
2
0
−6
4
5
2
9 −1
4
1
7
0
2
0 0 1 0 0


a b c
Halla los determinantes de los ejercicios 27 a 32 sabiendo que det  d e f  = 7.
g h i
b c a
a b c
a 3b c a
b c e f I 30. h i g I 29. g h i I 28. d 3e f I 27. d
e f d
d e f g 3h i 5g 5h 5i a
a + d b + e 2( c + f ) b
c e
2f
I 31. 2d + a 2e + b 2 f + c I 32. d
g
g
h
i h
2i
I 33. Usa determinantes para establecer cuáles de las siguientes matrices son inversibles:
3 −9
4 −9
6 −9
(a)
(b)
(c)
2
6
0
5
−4
6
Utiliza determinantes para averiguar si las matrices dadas en los ejercicios 34 a 36 son inversibles.






2
0
0 8
2 3 0
5
0 −1
1 −7 −5 0

I 34. 1 3 4
I 35. 1 −3 −2
I 36. 
3
8
6 0
1 2 1
0
5
3
0
7
5 4
Utiliza determinantes para averiguar si los conjuntos de vectores dados en los ejercicios 37 a
39 son linealmente independientes.

    
4
−7
−3
I 37.  6,  0, −5
−7
2
6


3
 5

I 39. 
 −6,
4


2
 −6
 ,
 0
7

    
7
−8
7
I 38. −4,  5,  0
−6
7
−5

  
−2
0
 −1  0
 ,  
 3  0
0
−3
Halla los determinantes de las matrices elementales dadas en los ejercicios 40 a 45.
34
3. Matrices y Determinantes

0
1
k


0 0
1 0
0 1
1
I 40. 0
0
1
I 41. 0
k

0
0
1
3.9. Ejercicios adicionales


0 0
1 0
0 1

0
k
0
k
I 42. 0
0
1
I 43. 0
0

0
0
1

1
0
0

0
0
1

0
1
0

1
0
0
0
I 44. 1
0
0
I 45. 0
1
Usa los ejercicios 40 a 45 para contestar razonadamente las preguntas 46 y 47 siguientes.
I 46. ¿Cuál es el determinante de una matriz elemental de reemplazo de filas?
I 47. ¿Cuál es el determinante de una matriz elemental escalonada con el número k en la diagonal?
En los ejercicios 48 a 51, comprueba
que det( EA) = (det E)(det A), donde E es la matriz
a b
elemental que se muestra y A =
.
c d
0 1
1 0
1 k
1 0
I 48.
I 49.
I 50.
I 51.
1 0
0 k
0 1
k 1
I 52. Sea A =
I 53. Sean A =
1
. Calcula 5A. ¿Es det(5A) = 5 det A?
2
3
4
a
c
b
d
y k un número. Halla una fórmula que relacione det(kA) con k y det A.
I 54. Sea A es una matriz n × n. Para cada una de las siguientes afirmaciones indica si es verdadera
o falsa. Justifica tus respuestas.
(a) Un determinante de una matriz n × n está definido por determinantes de submatrices de
orden (n − 1) × (n − 1).
(b) El cofactor (i, j) de una matriz A es la matriz Aij que se obtiene al eliminar de A su i-ésima
fila y su j-ésima columna.
(c) El desarrollo por cofactores de A bajando por una columna da un resultado opuesto al
que se obtiene mediante el desarrollo por cofactores a lo largo de una fila.
(d) El determinante de una matriz triangular es la suma de los elementos de la diagonal
principal.
3
1
yv=
. Calcula el área del paralelogramo cuyos vértices son 0, u, v, y
0
2
3 1
u + v. Halla el determinante de la matriz [u v] =
. ¿Qué diferencia hay entre los dos
0 2
resultados? Reemplaza el primer elemento de v por un número indeterminado x, y repite el
problema. Haz un dibujo que explique tus resultados.
I 55. Sean u =
35
3.9. Ejercicios adicionales
3. Matrices y Determinantes
a
c
I 56. Sean u =
yv=
donde a, b y c son positivos (por simplificar). Calcula el área del
b
0
paralelogramo cuyos vértices son u, v, u + v, y 0. Halla los determinantes de las matrices [u v]
y [v u]. Haz un dibujo que explique tus resultados.
En los ejercicios 57 y 58, A y B son matrices n × n. Indica para cada enunciado si es verdadero
o falso. Justifica tus respuestas.
I 57.
(a) Una operación de reemplazo de filas no afecta al determinante de una matriz.
(b) El determinante de A es el producto de los pivotes presentes en cualquier forma escalonada U de A, multiplicado por (−1)r , donde r es el número de intercambios de fila
realizados durante la reducción por filas de A a U.
(c) Si las columnas de A son linealmente dependientes, entonces det A = 0.
(d) det( A + B) = det A + det B.
I 58.
(a) Si se realizan dos intercambios sucesivos de fila, entonces el nuevo determinante es igual
al determinante antiguo.
(b) El determinante de A es el producto de los elementos diagonales de A.
(c) Si det A es cero, entonces dos filas o dos columnas son iguales, o una fila o una columna
es cero.
(d) det A T = (−1) det A.

1
I 59. Halla det B5 , donde B = 1
1
0
1
2

1
2.
1
I 60. Utiliza las propiedades de los determinantes sobre cómo les afectan las operaciones elementales por filas para demostrar que si dos filas de una matriz cuadrada A son iguales, entonces
det A = 0. Esto se cumple también para dos columnas. ¿Por qué?
I 61. Demuestra que si A es inversible, entonces det A−1 =
1
.
det A
I 62. Halla una fórmula para det(rA) cuando A es una matriz n × n.
I 63. Sean A y B matrices cuadradas. Demuestra que aunque AB y BA no sean iguales, siempre
es cierto que det( AB) = det( BA).
I 64. Sean A y P matrices cuadradas, con P inversible. Demuestra que det( PAP−1 ) = det A.
I 65. Sea U una matriz cuadrada tal que U T U = I. Demuestra que det U = ±1.
I 66. Supongamos que A es una matriz cuadrada tal que det( A4 ) = 0. Explica por qué A no
puede ser inversible.
Comprueba que det( AB) = (det A)(det B) para las matrices de los ejercicios 67 y 68.
36
3. Matrices y Determinantes
I 67. A =
3
6
0
,
1
B=
2
5
3.9. Ejercicios adicionales
0
4
I 68. A =
3
−1
6
,
−2
B=
4
−1
2
−1
I 69. Sean A y B matrices 3 × 3, con det A = 4 y det B = −3. Utiliza las propiedades de los
determinantes para calcular:
(a) det( A5 )
(b) det(5A)
(c) det( B T )
(d) det( A−1 )
(e) det( A3 )
I 70. Sean A y B matrices de 4 × 4, con det A = −1 y det B = 2. Halla:
(b) det( B5 )
(a) det( AB)
(c) det(2A)
(d) det( A T A)
I 71. Comprueba que det A = det B + det C, donde
a+e b+ f
a b
A=
,
B=
,
c
d
c d
I 72. Demuestra que det A = det B + det C, donde



a b u1 + v1
a
A =  c d u2 + v2  ,
B = c
e f u3 + v3
e
b
d
f

u1
u2  ,
u3
C=
e
c

a
C = c
e
(e) det( B − AB)
f
d
b
d
f

v1
v2 
v3
Sin embargo, observa que A no es igual a B + C.
I 73. Sean A =
1
0
0
1
yB=
a
c
b
. Demuestra que det( A + B) = det A + det B si, y sólo si,
d
a + d = 0.
I 74. La multiplicación por la derecha por una matriz elemental afecta a las columnas de A en
la misma forma en que la multiplicación por la izquierda afecta a las filas. Sea E una matriz
elemental y A una matriz cuadrada del mismo orden. Utiliza el hecho evidente de que E T es
otra matriz elemental, para demostrar que
det( AE) = (det E)(det A)
37