Introducción a Matrices y Operaciones
Introducción a Matrices y Operaciones
El álgebra matricial es una rama de las matemáticas en la que los números se tratan colec-
tivamente (como arreglos rectangulares de números llamados matrices) en lugar de individual-
mente, como en el álgebra ”ordinaria”. El término matriz se define formalmente más adelante.
a11 a12 a13 ··· a1j ··· a1n
a21 a22 a23 ··· a2j ··· a2n
.. .. ..
. . .
A=
ai1 ai2 ai3 · · · aij ··· a3n
. .. .. .. .. ..
.. . . . . .
am1 am2 am3 · · · amj · · · amn
Las matrices se denotan por letras mayúsculas A, B, C, ...,etc.
Los elementos o entradas, que en nuestro caso son números reales (llamados escalares) se
denotan por aij .
{
i − ésima fila
aij es el elemento situado en la
j − ésima columna
Definición 1.2 Una matriz de orden m × n es una matriz que tiene m filas y n columnas. Si
m = n,se dice que la matriz es cuadrada.
En lo que sigue una matriz de m filas y n columnas será denotada por Am×n = [aij ] o
simplemente por A.
1
es un vector columna de orden 4. Similarmente, una matriz que tiene solo una fila es un vector
fila [ ]
Y = 4 6 −7 1
es un vector fila de orden 4.
Definición 1.6 (Igualdad de matrices) Dos matrices A y B del mismo orden son iguales
si todos sus elementos correspondientes son iguales. En tal caso se usa la notación A = B.
1.3.2 Propiedades
Para A, B y C matrices de orden m × n cualesquiera y para r, s ∈ R, se cumplen
1. A + B = B + A (propiedad conmutativa)
2. A + (B + C) = (A + B) + C (propiedad asociativa)
2
3. Existe una única matriz (nula) 0 tal que para toda matriz A, A + 0 = A
4. Para toda matriz A, existe una única matriz −A (opuesta de A) tal que
A + (−A) = 0.
5. r(A + B) = rA + rB
6. (r + s)A = rA + sA
7. r(sA) = (rs)A
8. 1A = A.
en donde λi ∈ R.
Es decir
λ1 0 0 ··· 0
0 λ2 0 ··· 0
0 0 λ3 ··· 0
A=
.. .. .. ..
. . . .
0 0 0 ··· λn
Propiedad fundamental
A·I =I ·A=A
Ejemplo 1.8
a11 a12 a13 a14 a15
0 a22 a23 a24 a25
A=
0 0 a33 a34 a35
0 0 0 a44 a45
0 0 0 0 a55
es una matriz triangular superior
3
4. Matriz triangular inferior.- Es la matriz cuadrada que verifica aij = 0, ∀i < j
Ejemplo 1.9
a11 0 0 0
a21 a22 0 0
A=
a31 a32 a33 0
4
et1 = (1, 0) , et2 = (0, 1)
(a, b) = (a, 0) + (0, b) = a (1, 0) + b (0, 1)
= ae1 + be2
Obsérvese que en general, Eij = ei etj es la matriz nula, excepto para el elemento (i, j) que es
∑n
igual a uno. También In = Eii .
i=1
El producto de dos matrices no está definida de manera obvia; esto es, el producto de dos
matrices no se obtiene multiplicando sus componentes correspondientes.
Antes de definir la multiplicación matricial, se requiere una definición previa.
[ ]
Definición 1.12 Sea A = a11 a12 · · · a1n una matriz o vector fila n−dimensional y sea
b11
b21
B = ..
.
bn1
una matriz o vector columna n−dimensional. Entonces el producto, AB, de A y B está dado
por
b11
b12
[ ] .
AB = a11 a12 . . . a1n
.
.
bn1
= a11 b11 + a12 b21 + ... + a1n bn1
Obsérvese que el producto de un vector fila y un vector columna no se puede definir a menos
que sean de tamaños compatibles. Además, el vector fila debe escribirse a la izquierda.
[ ] 4
Ejemplo 1.13 El producto de A = 3 2 5 y B = 3 es
2
5
[ ] 4
AB = 3 2 5 3 = 3 (4) + 2 (3) + 5 (2) = 12 + 6 + 10 = 28.
2
Usando la definición de un vector fila por un vector columna se puede definir la multiplica-
ción matricial.
A = [aij ] de orden m × p
y
B = [bij ] de orden p × n;
el producto A × B, en ese orden, es la matriz C = [cij ] de orden m × n cuya componente cij es
el producto de la fila i de A y la columna j de B.
c11 = (fila 1 de
A)×(columna 1 de B) = 2 (3) + 1 (2) + (−3)1 = 5
c12 = (fila 1 de
A)×(columna 2 de B) = 2(−1) + 1 (4) + (−3)5 = −13
c21 = (fila 2 de
A)×(columna 1 de B) = 4 (3) + (−5)2 + 1 (1) = 3
c22 = (fila 2 de
A)×(columna 2 de B) = 4(−1) + (−5)4 + 1 (5) = −19
[ ] [ ]
c11 c12 5 −13
De lo anterior se concluye que C = = .
c21 c22 3 −19
Observación 1.17 Dadas dos matrices, por ejemplo, A2×2 y B2×4 es posible efectuar (AB) =
C2×4 ,debido a que sus tamaños son compatibles, es decir, el número de columnas de A es igual
al número de filas de B; sin embargo, no es posible calcular BA.
6
Propiedades de la multiplicación
Suponiendo que A, B y C son matrices compatibles con respecto la suma y el producto,se
cumplen
1. A(B + C) = AB + AC 1a propiedad distributiva
2. (A + B)C = AC + BC 2a propiedad distributiva
3. A(BC) = (AB)C propiedad asociativa
4. (kA)B = k(AB) = A(kB), k es un escalar.
Además se debe tener en cuenta
1. Por lo general AB ̸= BA
2. AB = 0 no implica necesariamente que A = 0 o B = 0
3. AB = AC no implica necesariamente que B = C
Matriz simétrica y matriz antisimétrica.
Una matriz cuadrada A se llama simétrica si At = A y antisimétrica si At = −A.
2 1 5
Ejemplo 1.18 Dada la matriz A = 1 4 6 , su transpuesta es la matriz
5 6 9
2 1 5
At = 1 4 6
5 6 9
es decir At = A por lo que A es una matriz simétrica. Observe que los elementos opuestos por
la diagonal principal son iguales.
0 1 −2 3
−1 0 6 7
Ejemplo 1.19 Dada la matriz A = 2
, su transpuesta es la matriz
−6 0 −5
−3 −7 5 0
0 −1 2 −3
1 0 −6 −7
At =
−2
,
6 0 5
3 7 −5 0
es decir At = −A, por lo que A es una matriz antisimétrica.
Observe que todos los elementos de la diagonal principal son nulos y los elementos de la
matriz A opuestos a ella son iguales pero de signo contrario.
Propiedades
1. Si A es una matriz cuadrada de orden n entonces A + At es una matriz simétrica y A − At
es una matriz antisimétrica
La justificación es inmediata puesto que
(A + At )t = At + (At )t = At + A = A + At
(A − At )t = At − (At )t = At − A = −(A − At ).
7
2. Toda matriz cuadrada A se puede descomponer en la suma de dos matrices, una simétrica
y otra antisimétrica. Es decir
1 1
A = (A + At ) + (A − At ),
2 2
en donde B = 21 (A + At ) es simétrica y C = 21 (A − At ) es antisimétrica.
Producto interno de un vector sumante con si mismo es un escalar, por ejemplo si tomamos el
vector sumante de orden n entonces
1tn · 1n = n
y otros productos de diferente orden son matrices con todos sus elementos iguales a uno, así
1 [ ] 1 1
13 · 1t2 = 1 1 1 = 1 1 = J 3×2 .
1 1 1
8
Ejemplo 1.22 Supogamos que un experimentado taxista tiene la certeza que cuando está en
la ciudad 1 tiene una probabilidad de 0.2 que el próximo destino será dentro de la ciudad
1 y una probabilidad de 0.8 que ésta será en la ciudad 2. Pero cuando está en la ciudad 2
sus probabilidades son de 0.4 de ir a la ciudad 1 y 0.6 de quedarse en la ciudad 2. Estas
probabilidades pueden ser juntadas en una matriz
[ ]
0.2 0.8
P =
0.4 0.6
luego pij es la probabilidad cuando está en la ciudad 1 y que la próxima parada es la ciudad j.
y P 2 1 = P (P · 1) = P · 1 = 1, y extendiendo
P n 1 = 1,
demostrando que la suma de filas de una potencia de una matriz de probabilidad de transición
son también unitarios. En general, una matriz de probabilidad de transición es de la forma
p11 p12 · · · p1j · · · p1m
p21 p22 · · · p2j · · · p2m
. ..
. .. ..
.
.. ..
.
. . . .
P = = (pij )m×m , i, j = 1, 2, · · · , m,
pi1 pi2 · · · pij · · · pim
. .. .. .. .. ..
.. . . . . .
pm1 pm2 · · · pmj · · · pmm
donde
∑
m
pij = pi1 + pi2 + · · · + pim = 1,
j=1
9
Teorema 1.24 La matriz centralizadora satisface
ii.) Cn 1 = 0.
iii.) Cn J = J Cn = 0.
es el vector dato con cada observación expresada como una desviación de x. Post multiplicando
xt Cn por x se tiene
( ) ( )
xt Cn x = xt − x1t x = xt x − x 1t .x = xt x − nx2 .
1 t
Así tenemos que xt es un vector dato, x .1 es la media, xt Cn es el vector de desviaciones de la
n
media, y xt Cn x es la suma de cuadrados alrededor de la media. La expresión xt Cn x es un caso
especial de la forma xt Ax, conocida como una forma cuadrática posteriormente estudiada.
Definición 1.25 La matriz A = (aij )n×n es inversible o regular o no singular, si existe una
matriz B tal que
AB = BA = In .
A la inversa de la matriz A, si existe, se le denota por A−1 y es igual a la matriz B.
Definición 1.26 Una matriz cuadrada no singular es ortogonal si y sólo si su inversa es igual
a su transpuesta. Esto es
A es ortogonal ⇔ A−1 = At
Teorema 1.27 Dadas las matrices A = (aij )n×n y B = (bij )n×n entonces
10
1.7 Matrices particionadas
Consideremos la matriz
1 6 8 9 1 −4
2 0 −1 0 4 5
A=
4 5 3 2 6 9 .
9 −2 1 −3 7 2
6 0 7 6 8 1
Supongamos que podemos trazar líneas entre ciertas filas y columnas como en
..
1 6 8 9 . 1 −4
..
2 0 −1 0 . 4 5
..
. 9
A= 4 5 3 2 6
··· ··· ··· ··· ··· ··· ···
.
9 −2 1 −3 .. 7 2
..
6 0 7 6 . 8 1
Cada uno de los arreglos de números en las 4 secciones de A producidas por el trazo de líneas
es una matriz
1 6 8 9 1 −4
A11 = 2 0 −1 0 , A12 = 4 5 ,
[ 4 5 3 2 ] 6[ 9 ]
9 −2 1 −3 7 2
A21 = y A22 = .
6 0 7 6 8 1
De esta manera, la matriz A puede ser escrita como una matriz de matrices
[ ]
A11 A12
A= (1.1)
A21 A22
Esta especificación de A es llamada una partición de A, y las matrices A11 , A12 , A21 y A22 se
dicen que son submatrices de A; además A de (1.1) es llamada una matriz particionada.
Observaciones.
1. A11 y A21 tienen el mismo número de columnas, como A12 y A22 . De igual modo A11 y
A12 tienen el mismo número de filas, como A21 y A22 .
11
donde A21 y A22 es como en el caso anterior, y
[ ] [ ]
1 6 8 9 1 −4
A01 = , A02 = ,
[ 2 0 −1 0] [ 4 ]5
A03 = 4 5 3 2 , A04 = 6 9 .
donde K, L, M y N son las submatrices con sus ordenes mostrados como subí[Link]
podemos decir que la matriz A = (aij )m×n puede ser particionada en p filas y q columnas
de submatrices como
A11 A12 · · · A1q
A21 A22 · · · A2q
A = .. .. . . ..
. . . .
A1p A2p · · · Apq
donde Aij es la submatriz en la i−ésima fila y la j−ésima columna de submatrices. Si
la la i−ésima fila de submatrices tiene mi filas de elementos y la j−ésima columna de
submatrices tiene nj columnas, entonces Aij tiene orden mi × nj , donde
∑
p
∑
p
mi = m y nj = n.
i=1 j=1
12
[ ] 3 2
3 7 5
Por ejemplo, si A = entonces At = 7 0 . Por otro lado, si particionamos la
2 0 6
5 [6 ] [ ]
[ ] 3 7 5
matriz A de la forma A = A11 ..
. A12 , con A11 = 2 y A12 = 0 6 tenemos que
[ . ]t [ ] 3 2
3 .. 7 5 [ ]t At11 ··· ···
t
A = = .. = = .
. A11 . A12 At12 7 0
2 .. 0 6
5 6
13
[ ][ ] [ ]
A11 A12 B1 A11 B1 + A12 B2
AB = =
A21 A22 B2 A21 B1 + A22 B2
[ ] 6 4 [ ]
A21 B1 = 0 −4 −2 −2 1 = 14 −18
−3 7
[ ]
[ ] −1 3 [ ]
A22 B2 = 7 −1 = −12 19
5 2
−5 4
−6 2
=
··· ···
2 1
Por tanto
−5 4
−6 2
AB =
··· ···
2 1
es importante que cada producto por bloques Aij Bi , i, j = 1, 2 sea escrito en ese orden, ya que
la multiplicación de matrices no es conmutativa.
particionar adecuadamente A y B, según sea el caso, para calcular por bloques AB.
1.8 Ejercicios
1. Determinar el valor de verdad de las siguientes afirmaciones, justificando debidamente su
respuesta.
(a) ¿Es det (αA) = α det (A) para toda matriz A = (aij )n×n ?
(b) Sean A = (aij )n×n y B = (bij )n×n . Entonces det (A · B) = det (A) · det (B).
(a) Si AB = I = CA entonces B = C
(b) Sean A = (aij )n×n y B = (bij )n×n . Entonces det (A + B) = det (A) + det (B).
14
4 − (i + j) , i < j
4. Si A = (aij ) donde aij = i, j = 1, 2, 3 . Hallar tr (A − I).
j − i − 1, i ≥ j
a p q
5. Si A = r a s , es una matriz diagonal, calcular Ak , para k ∈ Z+ .
u v a
6. Probar que si A y B son antisimétricas, entonces A + B y αA son antisimétricas, para
todo escalar α.
8. Probar que toda matriz cuadrada A puede ser escrita como la suma de uma matriz
simétrica y uma anti-simétrica. (Sugerencia. Observe el resultado de la suma de A + At
con A − At .)
10. (a) Probar que si A es una matriz m × n tal que AX = 0, para toda matriz X, n × 1,
entonces A = 0.
(b) Sean B y C matrices m × n, tales que BX = CX, para todo X, n × 1. Probar que
B = C.
(Sugerencia: Use el item anterior.)
11. Probar que la matriz identidad In es la única matriz tal que AIn = In A = A para cualquier
matriz A, n × n. (Sugerencia: Sea Jn una matriz tal que AJn = Jn A = A. Mostrar que
Jn = In .)
15
( )
1 q
(d) Simplificar (pI + qJ) 1− J .
p p + qn
(e) Explicar porqué JXJ = xJ.
(f) Probar que I − J n es idempotente, pero I − J no.
(g) Para at de orden n mostrar que |λI + 1at | = λn−1 (λ + at 1).
(h) Probar |aIn + bJn | = an−1 (a + bn) .
18. Probar que P 2 es una matriz de probabilidad de transición, si P es una matriz de proba-
bilidad de transición.
(a) |J n | = 0.
(b) J 2n = nJ n .
∑
n
(c) La matriz centralizante Cn satisface X t Cn X = (xi − x)2 .
i=1
16