Tema 3 - Operaciones Con Matrices
Tema 3 - Operaciones Con Matrices
Tema 3 - Operaciones Con Matrices
Ideas clave . . . . . . . . . . . . . . . . . . . . . . 3
3.8 A fondo . . . . . . . . . . . . . . . . . . . . 30
3.9 Test . . . . . . . . . . . . . . . . . . . . . 33
Esquema
El lenguaje matricial es fundamental en el estudio del Álgebra Lineal, donde las ideas
principales se expresan por medio de operaciones entre matrices. En este tema se ex-
ponen los principales tipos de matrices, sus características y las operaciones y trans-
formaciones que se pueden hacer sobre ellas. Veremos que las matrices se pueden
sumar y multiplicar, y por tanto conforman un sistema algebraico de alguna forma
análogo al conjunto de los números reales R, propiedad estudiada por primera vez de
forma sistemática en 1858 a manos del matemático Arthur Cayley.
I Conocer el concepto de vector, matriz y las operaciones sobre ellos de suma, pro-
ducto y producto por un escalar.
I Matrices cero-uno.
Las operaciones con matrices son importantes porque permiten realizar cálculos ma-
temáticos y resolver problemas que de otra manera serían difíciles o imposibles de
abordar. Además, tienen aplicaciones prácticas en diversas áreas de la ciencia y la in-
geniería. Las operaciones fundamentales con matrices son la suma, la resta y la mul-
tiplicación.
Por otro lado, la multiplicación de matrices es una operación más compleja que re-
quiere el uso de reglas específicas. En general, para que la multiplicación de matrices
esté definida, el número de columnas de la matriz izquierda debe ser igual al número
de filas de la matriz derecha. La matriz resultante de la multiplicación tiene un número
de filas igual al número de filas de la matriz izquierda y un número de columnas igual
al número de columnas de la matriz derecha.
Definición 1
a11 a12 · · · a1n
a21 a22 · · · a2n
A= .. .. . . ..
.
. . .
an1 an2 · · · amn
Las entradas o coeficientes de la matriz A se denotan por aij , indicando que son el
elemento de la matriz A en la fila i y columna j. Generalmente los coeficientes aij
Diremos que dos matrices A y B son matrices iguales si tienen el mismo orden aij =
bij para todo i, j.
Definición 2
Ejemplo 1.
La matriz:
3 4 1
A=
0 2 −1
Las matrices:
4 1
B= y C= 3 4
2 −1
Una matriz fila es aquella que solo tiene una fila y una matriz columna es aque-
lla que solamente tiene una columna. Son lo que usualmente se conoce como
vectores.
Ejemplo 2.
3
La matriz A = 0 −1 0 2 es una matriz (o vector) fila, y B = 4
−1
es una matriz (o vector) columna.
Definición 4
Una matriz cuadrada es aquella matriz A ∈ Mn×n que tiene el mismo número
de filas que de columnas. Al conjunto de matrices cuadradas de orden n × n lo
denotaremos por Mn .
Ejemplo 3.
2 −1 0
La matriz A = 0 8 1 es una matriz cuadrada de orden 3.
2 −1 5
Definición 5
Una matriz diagonal es una matriz cuadrada que únicamente tiene coeficientes
no nulos en la diagonal, es decir, aquella matriz A tal que aij = 0, ∀i 6= j.
Decimos que una matriz A estriangular superior si todos los coeficientes que se
encuentran por debajo de la diagonal son cero, es decir, si aij = 0, ∀i > j.
Ejemplo 4.
De las siguientes matrices A es diagonal, B es triangular superior y C triangular
inferior.
a11 0 0 0 a11 a12 a13 a14
0 a22 0 0 0 a22 a23 a24
A=
;
B=
0 0 a33 0 0 0 a33 a34
0 0 0 a44 0 0 0 a44
a11 0 0 0
a21 a22 0 0
=
a31 a32 a33 0
a41 a42 a43 a44
Como caso particular de matriz diagonal está la matriz identidad In , la cual se carac-
teriza por tener todos los coeficientes de la diagonal aij igual a 1. Por ejemplo:
1 0 0
A= 0 1 0
0 0 1
Definición 6
Una matriz simétrica es una matriz A tal que aij = aji para todo i, j. Es decir,
es una matriz en la que los coeficientes son simétricos con respecto a la diagonal
principal.
La matriz:
1 4 0
A = 4 −3 6
0 6 9
En esta sección vamos a describir las distintas operaciones aritméticas que podemos
realizar con matrices.
Suma de matrices
a + b11 a12 + b12 ··· a1n + b1n
11
a21 + b21 a22 + b22 ··· a2n + b2n
A + B := (aij + bij ) = .. .. .. ..
.
. . .
an1 + bn1 an2 + bn2 · · · amn + bmn
Se observa por tanto que para que la suma esté definida es necesario que las matrices
Ejemplo 6.
1 4 0 2 0 3 3 4 3
4 −3 6 + 1 3 1 = 5 0 7
0 6 9 1 2 0 1 8 9
I Elemento neutro: Existe una matriz 0 ∈ Mm×n tal que A + 0 = A para cualquier
A ∈ Mm×n .
I Elemento opuesto: Para toda A ∈ Mm×n existe la matriz −A ∈ Mm×n tal que
A + (−A) = 0.
Estas propiedades hacen que el conjunto de matrices junto con la operación suma,
denotado por (Mm×n , +), forman un grupo Abeliano (ver Castellet y Llerena (1996)
para más información sobre teoría de grupos).
Ejemplo 7.
2 −1 0 4 −2 0
Dada la matriz A = 0 8 1 , se tiene que 2A = 0 16 2
2 −1 5 4 −2 10
El producto AB de dos matrices se define de manera algo más elaborada ya que para
que la operación esté definida, se requiere que A ∈ Mm×p y B ∈ Mp×n , es decir:
p
X
cij = aik bkj = ai1 b1j + ai2 b2j + ...aip bpj
k=1
Para n = 3, m = 4, p = 2:
a a
11 12 b11 b12 b13 b14
· =
a21 a22
b21 b22 b23 b24
a31 a32 p×m
n×p
a b + a12 b21 a11 b12 + a12 b22 a11 b13 + a12 b23 a11 b14 + a12 b24
11 11
= a21 b11 + a22 b21 a21 b12 + a22 b22 a21 b13 + a22 b23 a21 b14 + a22 b24
a31 b11 + a32 b21 a31 b12 + a32 b22 a31 b13 + a32 b23 a31 b14 + a32 b24
n×m
Ejemplo 8.
1 2 0 1 1·0+2·2 1·1+2·3 4 7
· = =
0 1 2 3 0·0+1·2 0·1+1·3 2 3
I k(AB) = (kA)B
tenemos que:
2 1 1 2 1 5
AB = =
0 −1 −1 1 1 −1
1 2 2 1 2 −1
BA = =
−1 1 0 −1 −2 −2
Este hecho junto con las operaciones que de suma y producto que acabamos de ver,
hacen que (Mn , +, ·) sea un anillo no conmutativo (Castellet y Llerena, 1996).
En muchas ocasiones, en las que las matrices que se desean multiplicar son grandes,
resulta práctico realizar la multiplicación pro bloques: Sean A ∈ Mn×p y B ∈ Mp×m
a11 a12 a13 a14
b11 b12 b13
a21 a22 a23 a24
b21 b22 b23
· =
a31 a32 a33 a34
b31 b32 b33
a41 a42 a43 a44
b41 b42 b43
a51 a52 a53 a54 1×3|3×3
5×1|5×3
a11 a12 a13 a14
a21 a a23 a24 b b b
22 21 22 23
= a31 · b11 b12 b13 + a32 · b31 b32 b33
a33 a34
a41 a42 a43 a44 b41 b42 b43
a51 a52 a53 a54
Ejemplo 9.
0 0
0 0 2 −1
0 0
Sean A = 0 0 0 0 y B = . Realicemos el producto por bloques.
5 0
0 0 0 0
1 0
0 0
0 0 2 −1
0 0
· =
0 0 0 0
5 0
0 0 0 0
3×2|3×2 1 0
2×2|2×2
0 0 2 −1
0 0 5 0
= 0 0 · + 0 · =
0
0 0 1 0
0 0 0 0
Antes de terminar la sección de la multiplicación de matrices, cabe decir que el algoritmo para
obtener la multiplicación de matrices puede ser escrito y ejecutado de forma sencilla, aunque
no de forma única. A continuación se presenta el pseudocódigo del algoritmo basado en la
definición del producto de dos matrices.
para i := 1, 2, . . . , n
paraj := 1, 2, . . . , m
hacer
cij := 0
para q := 1, 2, ..., k
cij := cij + aiq bqj
fin
fin
C = [cij ] es el producto de A por B
Definición 10
Sea A una matriz cuadrada de orden n y sea k un número entero no negativo, podemos
definir la potencia de A como Ak = A
| · .{z
. . · A} ·In
k veces
Existen muchas otras operaciones que se realizan con matrices, o trasformaciones de matrices
con sus aplicaciones en distintos ámbitos. A continuación destacamos tres de las mas impor-
tantes, la trasposcición de matrices, la inversa de una matriz y el determinante de una matriz
cuadrada
AB = BA = In
En caso de existir esta pareja de matrices, decimos que A es una matriz invertible y escri-
bimos B = A−1 .
Ejemplo 10.
1 1 2 −1
Sea A = es inversible y su matriz inversa es A−1 =
1 2 −1 1
1 1 2 −1 1 0
AA−1 = = y
1 2 −1 1 0 1
2 −1 1 1 1 0
A−1 A =
−1 1 1 2 0 1
I (AB)−1 = B −1 A−1 .
Demostración
3. AL igual que en la demostración (2), AT (A−1 )T = (A−1 A)T = InT = In y por otra
parte, (A−1 )T AT = In .
En el siguiente vídeo puedes ver un ejemplo del cálculo de la inversa de una matriz de tamaño
2 × 2.
Ejemplo 11.
2 0
2 1 −1
La traspuesta de la matriz A = es AT =
1 −1
0 −1 0
−1 0
I (A + B)T = AT + B T
I (AT )T = A
I (kA)T = kAT , ∀k ∈ K
I (AB)T = B T AT
a11 a12
det(A) = = a11 a22 − a12 a21
a21 a22
a11 a12 a12
I Si A es de orden 3, es decir A = a21 a22 a23 , entonces:
a31 a32 a33
I det(A) = 3 − (−2) = 5
a11 a12 · · · a1m
a21 a22 · · · a2m
A=
.. .. .. ..
. . . .
an1 an2 · · · anm
Y el determinante de A como:
Esta fórmula también es conocida como el desarrollo de Laplace del determinante de A por
su primera columna.
El cálculo del determinante de una matriz de orden n se realiza a partir de los menores adjuntos
de orden n − 1. La operación puede expresarse en términos de determinantes de orden 2 o
3, que pueden ser resueltos por la regla de Sarrus.
1 0 −1 2 1 0 −1 2
0 0 0 1 0 0 0 1
det
=
=
2 1 0 1 2 1 0 1
1 −1 1 −1 1 −1 1 −1
0 0 1 1 −1 2 0 −1 2 0 −1 2
= 1· 1 0 1 −0· 1 0 1 +2· 0 0 1 −1· 0 0 1 =
−1 1 −1 −1 1 −1 −1 1 −1 1 0 1
= 1 · 1 + 2 · 1 − 1(−1) = 4
I Si tres matrices cuadradas A, A0 y A00 son idénticas salvo en que la i-ésima fila (o columna)
de A es la suma de las filas (o columnas) correspondientes de A0 y A00 , entonces:
I Si se multiplican los elementos de una fila o columna de la matriz A por un escalar k, en-
tonces el determinante queda multiplicado por k.
I Si a una fila o columna de una matriz cuadrada A se le suma otra multiplicada por un escalar
k, entonces su determinante no varía.
EL calculo de los determinantes de una matriz, resulta de gran utilidad para obtener la inversa
de esta. Para ello, en primer lugar definimos la matriz adjunta de una matriz cuadrada:
A∗ := (αij )
Se cumple que si A es una matriz regular (se puede obtener la matriz inversa de A), entonces;
(A+ )T
A−1 =
det(A)
En el siguiente vídeo puedes ver un ejemplo del cálculo de la inversa de una matriz utilizando
el cálculo de los determinantes.
Para finalizar este tema, se profundiza en el estudio de las matrices cero-uno o matrices bolea-
nas, introducido en el Tema 4 y que son de gran utilidad para la representación de estructuras
discretas, tal como se verá en el tema 9, en el que se utilizan las matrices cero-uno para repre-
sentar grados.
Se llama matriz cero-uno o matriz boleana a aquella matriz en la que todos los elemen-
tos que la forman son 0 ó 1 y los algoritmos que trabajan con ellas, están basados en la
aritmética boleana sobre matrices de ceros y unos. Esta aritmética se construye con las
operaciones boleanas ∨, ∧e = sobre pares de bits que se vieron en el Tema 1:
1 si b1 = b2 = 1
b1 ∧ b2 =
0 otro caso
1 si b1 = 1 ó b2 = 1
b1 ∨ b2 =
0 otro caso
Es habitual encontrar las operaciones de unión e intersección bajo terminología inglesa como
meet para denotar a la matriz unión de A y B y join para la matriz intersección de A con B.
1 0 1 0 1 0
A= y B=
0 1 0 1 1 0
solución:
1∨0 0∨1 1∨0 1 1 1
A∨B = =
0∨1 1∨1 0∨0 1 1 0
1∧0 0∧1 1∧0 0 0 0
A∧B = =
0∧1 1∧1 0∧0 0 1 0
n
_
cij = aik ∧ bkj
k=1
Solución:
(1 ∧ 1) ∨ (0 ∧ 0) (1 ∧ 1) ∨ (0 ∧ 1) (1 ∧ 0) ∨ (0 ∧ 1)
A B = (0 ∧ 1) ∨ (1 ∧ 0) (0 ∧ 1) ∨ (1 ∧ 1) (0 ∧ 0) ∨ (1 ∧ 1)
(1 ∧ 1) ∨ (0 ∧ 0) (1 ∧ 1) ∨ (0 ∧ 1) (1 ∧ 0) ∨ (0 ∧ 1)
(1 ∨ 0) (1 ∨ 0) (0 ∨ 0)
= (0 ∨ 0) (0 ∨ 1) (0 ∨ 1)
(1 ∨ 0) (1 ∨ 0) (0 ∨ 0)
1 1 0
= 0 1 1
1 1 0
Al igual que en el caso del producto usual, se puede implementar el algoritmo para el producto
boleano en base al siguiente pseudocódigo:
para i := 1, 2, . . . , n
paraj := 1, 2, . . . , m
hacer
cij := 0
para q := 1, 2, ..., k
cij := cij ∨ (aiq ∧ bqj )
fin
fin
C = [cij ] es el producto boleano de A por B
Definición 17
Sea A una matriz cuadrada cero-uno de orden n y sea k un número entero positivo, po-
demos definir la potencia boleana de A como Ak = A
| .{z
.. A} ·In
k veces
Este producto está bien definido ya que el producto boleano es asociativo. Por definición
A0 = In .
Grimaldi, R. P. (2006). Discrete and Combinatorial Mathematics, 5/e. Pearson Education India.
Merino, L. M., & Santos Aláez, E. (2006). Álgebra lineal con métodos elementales. Madrid: Ed.
Thomson.
Castellet, M., & Llerena, I. (1996). Álgebra Lineal y Geometría. Barcelona: Reverte.
1
A = −5 3 4 1
2
Ejercicio 2.
¿Para qué valores de α ∈ R la matriz no es invertible?.
Ejercicio 3.
Calcular AB si:
1 1 0 4
a) A = y B=
3 2 3 1
1 −1
3 −2 −1
b) A = 0 1 y B =
1 0 2
2 3
Ejercicio 4.
1 0 1 0 1 1
Sean A = 1 1 0 y B = 1 0 1 .
0 0 1 1 0 1
Obtener: a) A ∨ B; b) B ∧ A y c)A B
Ejercicio 6.
Sean A y B matrices boleanas de orden m × n. Demostrar que:
a) A ∨ B = B ∨ A; b) A ∧ B = B ∧ A
3 4 −1
A = −15 −20 5
6 8 −2
Solución 2.
Tuna matriz es no invertible si su determinante es igual a cero. Por tanto, como |A| = 0
para cualquier valor de α ∈ R, la matriz A no es invertible para todo valor de α ∈ R
Solución 3.
2 −2 −3
1 11
a= ; b= 1 0
2
2 18
9 −4 4
Solución 4.
1 1 1 0 0 1 1 1 1
a) 1 1 1 ; b) 1 0 0 ; c) 1 1 1
1 0 1 0 0 1 1 0 1
Solución 5.
1 0 0 1 0 0 1 0 1
a) 1 1 0 ; b) 1 0 1 ; c) 1 1 1
1 0 1 1 1 0 1 1 1
Considerado durante mucho tiempo como un clásico en su género, este fue el primer libro
en inglés para incluir tres campos básicos del análisis de matrices, matrices simétricas y
formas cuadráticas, matrices y ecuaciones diferenciales, y las matrices positivas y su uso en
la teoría de probabilidades y matemáticas economía. Escrito en términos lúcidos y concisos,
este volumen abarca todos los aspectos clave de la matriz de análisis y presenta una gran
variedad de métodos fundamentales.
Publicado originalmente en 1960, este libro se sustituye la primera edición publicada
previamente por SIAM en la serie Classics. Aquí encontrarás una guía básica para las
operaciones con matrices y la teoría de matrices simétricas, además de la comprensión de
las matrices cuadradas generales, orígenes de matrices de Markov y matrices no negativas,
en general, la caracterización, la mínima-máxima de raíces características, productos de
Kronecker, funciones de matrices, y mucho más. Estas ideas y métodos servirán como
herramientas analíticas.
http://unir.summon.serialssolutions.com /
MatLab
MatLab es una herramienta muy utilizada para el cálculo de operaciones con matrices. En esta
herramienta las variables son matrices y el lenguaje proporciona ayuda para las operaciones
básicas con matrices. Además proporciona multitud de funciones para operar con estas
matrices. MatLab es de pago pero puedes utilizar la versión de prueba.
http://www.mathworks.es/programs/nrd/matlab-trial-request.html?region=es&s_eid=ppc_1817
Manual: http://www.mathworks.es/academia/student_center/tutorials/
Octave
http://www.gnu.org/software/octave/
Manual: http://en.wikibooks.org/wiki/Octave_Programming_Tutorial
6. El determinante representa:
A. El valor más importante de la matriz.
B. El valor menos importante de la matriz.
C. Cuando es cero significa que la matriz es invertible.
D. Cuando es cero significa que la matriz no es invertible.