Algebra Matricial
Algebra Matricial
Algebra Matricial
Álgebra matricial
1.1. Matrices
A partir de ahora K denotará un cuerpo.
Definición 1.1. Una matriz A = (aij ) de orden m × n sobre un cuerpo K es una tabla
de doble entrada con m · n elementos del cuerpo K, dispuestos en m filas y n columnas.
Al elemento que ocupa la fila i y la columna j de la matriz A lo denotaremos por aij o
por A(i, j). Así pues:
a11 a12 . . . a1n
a21 a22 . . . a2n
A= .
.. .. ..
.. . . .
am1 am2 . . . amn
Dos matrices A = (aij ) y B = (bij ) son iguales si tienen el mismo orden y para cada
par de índices i y j se tiene que aij = bij .
Al conjunto de las matrices de orden m × n sobre un cuerpo K se le denota por
Mm×n (K). Si m = n, se denota por Mn (K).
Si A ∈ M1×n (K) se dice que A es una matriz fila.
Si A ∈ Mm×1 (K) se dice que A es una matriz columna.
A las matrices de orden n × n se les llama matrices cuadradas de orden n.
Denotaremos por In la matriz diagonal de orden n, (δij ), con δii = 1 para i = 1, . . . , n
y δij = 0 para i 6= j. A esta matriz se le llamará matriz identidad n-ésima o matriz
identidad de orden n.
1
2 CAPÍTULO 1. ÁLGEBRA MATRICIAL
a11 + b11 a12 + b12 ... a1n + b1n
a21 + b21 a22 + b22 ... a2n + b2n
A+B = .
.. .. .. ..
. . . .
am1 + bm1 am2 + bm2 . . . amn + bmn
2. (α + β)A = αA + βA.
3. α(A + B) = αA + αB.
4. α(βA) = (αβ)A.
5. 1K A = A.
El elemento neutro del grupo (Mm×n (K), +) es la matriz con todas sus entradas 0
y la denotaremos por (0).
Definición 1.5. Sean A = (aij ) ∈ Mm×n (K) y B = (bij ) ∈ Mn×s (K) se define el
producto de A por B como una matriz AB en Mm×s (K) en donde:
n
P
(AB)(i, j) = aik bkj , para i = 1, . . . , m y j = 1, . . . , s.
k=1
Observaciones 1.6.
1. Nótese que, para hacer el producto AB es necesario que el número de columnas de
A coincida con el número de filas de B.
2. En general, el producto de matrices no es conmutativo. Por ejemplo,
0 0 1 0 0 0 0 0 1 0 0 0
= 6= = .
1 1 1 0 2 0 0 0 1 0 1 1
1.2. OPERACIONES CON MATRICES 3
(AB)C = A(BC).
A (B + C) = AB + AC.
(B + C)A = BA + CA.
1. Denotemos por F1 (A) = (a11 a12 · · · a1n ) ∈ M1×n (K) la matriz fila correspon-
diente a la fila 1 de la matriz A, entonces
b11 b12 ··· b1s
b21 b22 ··· b2s
F1 (A) B = (a11 a12 · · · a1n ) . .. =
.. ..
.. . . .
bn1 bn2 · · · bns
n
X n
X n
X
= a1k bk1 a1k bk2 ··· a1k bks =
k=1 k=1
| {z } | {z } |k=1 {z }
11 12 1s
= a11 (b11 b12 · · · b1s ) +a12 (b21 b22 · · · b2s ) + · · · + a1n (bn1 bn2 · · · bns ) =
| {z } | {z } | {z }
F1 (B) F2 (B) Fn (B)
Ejemplo 1.13.
1 2 0 7 −1
1 −1
1 2 −1
3
6 −3
0 =
1 1 0 4 −1
1 2
0 4 3 15 6
Nótese que, como hemos visto en 2 y 3 de las observaciones 1.12, se tiene que:
(7, −1) = 1(1, −1) + 2(3, 0) + 0(1, 2),
(6, −3) = 1(1, −1) + 2(3, 0) + (−1)(1, 2),
7 1 2 0
= 1 1 + 3 2 + 1 −1
6
4 1 1 0
15 0 4 3
y análogamente para las restantes filas y columnas.
6 CAPÍTULO 1. ÁLGEBRA MATRICIAL
Im −−−→ EαFi
αF i
EαFi (k, k) = 1 si k 6= i,
EαFi (i, i) = α,
y las demás entradas son cero.
Fi (EFi ↔Fj A) = 0F1 (A) + · · · + 0Fi (A) + · · · + 1Fj (A) + · · · + 0Fn (A) = Fj (A),
| {z } | {z }
i j
Fj (EFi ↔Fj A) = 0F1 (A) + · · · + 1Fi (A) + · · · + 0Fj (A) + · · · + 0Fn (A) = Fi (A),
| {z } | {z }
i j
y , si k 6= i, j, entonces Fk (EFi ↔Fj A) = Fk (A).
Fi (EFi +αFj A) = 0F1 (A) + · · · + 1Fi (A) + · · · + αFj (A) + · · · + 0Fn (A),
| {z } | {z }
i j
es decir, Fi (EFi +α Fj A) = Fi (A) + α Fj (A),
y, si k 6= i, entonces Fk (EFi +αFj A) = Fk (A).
Corolario 1.16. Toda matriz elemental es invertible y su inversa también es una matriz
elemental. Concretamente,
(EFi ↔Fj )−1 = EFi ↔Fj , (EαFi )−1 = Eα−1 Fi y (EFi +αFj )−1 = EFi −αFj .
Definición 1.17. Dos matrices A, B ∈ Mm×n (K) son equivalentes por filas si existen
matrices elementales E1 , . . . , Es tales que B = Es · · · E1 A.
Según esta observación la matriz EFi ↔Fj es la matriz obtenida de la identidad inter-
cambiando la columna i-ésima con la j-ésima, EαFi es la matriz obtenida de la identidad
multiplicando la columna i-ésima por un escalar α 6= 0 y EFi +αFj es la matriz obtenida
de la identidad sumándole a la columna j-ésima la i-ésima multiplicada por α ∈ K.
Esto justifica la notación siguiente:
EFi ↔Fj = ECi ↔Cj , EαFi = EαCi , EFi +αFj = ECj +αCi .
Si además todos los pivotes son 1 y los elementos que aparecen en la misma columna
que el pivote de una fila son todos cero, se dice que es escalonada reducida por filas.
Nótese que si una matriz es escalonada reducida por filas, el número de pivotes
coincide con el número de filas no nulas.
Demostración.
“⇒” Si A es invertible ninguna de sus filas es nula y, como es escalonada, el número
de pivotes coincide con el número de filas y, por tanto, con el número de columnas.
Además, por ser A reducida, se tiene que A = In .
“⇐” Trivial.
1.4. FORMA ESCALONADA Y RANGO POR FILAS DE UNA MATRIZ 9
Teorema 1.22. Toda matriz es equivalente por filas a una matriz escalonada y a una
única matriz escalonada reducida por filas.
Demostración. Ver álgebra Lineal de L. Merino y E. Santos, pág. 21, o Álgebra Lineal
y aplicaciones de J. Arvesú, R. Álvarez y F. Marcellán, pág. 43.
Ejercicio 1.23. Calcular una matriz escalonada equivalente por filas a la matriz A y
una escalonada reducida equivalente por filas a la matriz B, siendo
1 −2 3 9 1 −2 3 9
A = −1 3 0 −4 , B = 0 1 3 5 ∈ M3×4 (R).
2 −5 5 17 0 0 1 2
Solución.
Para la matriz A se tiene que:
1 −2 3 9 1 −2 3 9
−1 3 0 −4 −−−−−→ 0 1 3 5 −−−−−−→
F2 + F1 F3 −2 F1
2 −5 5 17 2 −5 5 17
1 −2 3 9 1 −2 3 9
−→ 0 1 3 5 −−−−−→ 0 1 3 5 .
F3 + F2
0 −1 −1 −1 0 0 2 4
Análogamente para la matriz B,
1 −2 3 9 1 0 9 19
0 1 3 5 −−−−−−→ 0 1 3 5 −−−−−−→
F1 +2 F2 F1 −9 F3
0 0 1 2 0 0 1 2
1 0 0 1 1 0 0 1
−→ 0 1 3 5 −−−−−−→ 0 1 0 −1
F2 −3 F3
0 0 1 2 0 0 1 2
Corolario 1.24. Sea A ∈ Mn (K).
nos indica como calcular la inversa de una matriz no singular, usando transformaciones
elementales por filas.
10 CAPÍTULO 1. ÁLGEBRA MATRICIAL
Ejemplo 1.25.
−1 −2
Sea A = ∈ M2 (R). Para calcular la inversa de la matriz A procedemos
3 8
del modo siguiente:
−1 −2 1 0 1 2 −1 0
(A|I2 ) = −−−→ −−−−−→
3 8 0 1 −F1 3 8 0 1 F2 −3F1
1 2 −1 0 1 0 −4 −1
! !
1 2 −1 0
−→ −−−→ 3 1 −−−−−→ 3 1
0 2 3 1 1
F
2 2
0 1 F1 −2F2 0 1
2 2 2 2
−4 −1
!
y, así, A−1 = 3 1 .
2 2
Definición 1.26. Sea A ∈ Mm×n (K), llamaremos rango por filas de A, y lo denotare-
mos por rf (A), al número de pivotes (o el número de filas no nulas) de cualquier matriz
escalonada por filas que sea equivalente por filas a A.
Observación 1.27. El rango por filas está bien definido, ya que cualquiera de las matrices
escalonadas que sean equivalentes por filas a A tiene el mismo número de pivotes. En
efecto, si B es una matriz escalonada por filas, su número de pivotes coincide con el
número de pivotes de la única escalonada reducida a la que es equivalente por filas, ya
que:
1.5. Determinantes
En este apartado, las matrices que se consideran son matrices cuadradas.
Sea n > 1. Supuesto conocido el valor del determinante de una matriz de orden
n − 1, se define el determinante de la matriz A como:
|A| = det(A) = a1j α1j + · · · + anj αnj , para cualquier j ∈ {1, . . . , n}.
También se puede calcular cambiando el papel de las columnas por filas, es decir,
|A| = det(A) = ai1 αi1 + · · · + ain αin , para cualquier i ∈ {1, . . . , n}.
Ejemplos 1.28.
3.
1 2 3
1 3
4 6 = 6 − 12 = −6 y 0 1 2 = 0 + 0 + 8 − 6 − 2 − 0 = 0.
2 1 0
3.
F1 (A) F1 (A) F1 (A) F1 (A) F1 (A)
· · · ··· ··· ··· ···
Fi (A) + Fj (A) Fi (A) Fi (A) Fj (A) Fj (A)
0 = ··· = ··· + ··· + ··· + ···
Fj (A) + Fi (A) Fj (A) Fi (A) Fj (A) Fi (A)
· · · ··· ··· ··· ···
Fn (A) F (A) F (A) F (A) F (A)
n n n n
F1 (A) F1 (A)
··· ···
Fi (A) Fj (A)
= · · · + · · ·
Fj (A) Fi (A)
··· ···
F (A) F (A)
n n
4. Sea A = (aij ) una matriz de orden n, usaremos inducción en n para probar esta
propiedad. Si n = 1 es evidente. Supuesto n > 1 y que el resultado es cierto para
matrices de orden n − 1, si
a11 . . . a1n
··· ... ···
0
A = βai1 . . . βain
... ... ...
an1 . . . ann
A es singular si, y solo si, A es equivalente por filas a una matriz escalonada por
filas A0 que tiene una fila de ceros si, y solo si, existen E1 , . . . , Es matrices elementales
tal que A = E1 · · · Es A0 , de lo que se deduce que
Demostración.
Caso 1) Si A es no singular, hemos visto que A es producto de matrices elementales, es
decir A = E1 · · · Ek .
Utilizando reiteradamente las propiedades 3, 4 y 5 de la proposición 1.29 y que el
producto de matrices es asociativo, se tiene que:
|AB| = |E1 (E2 · · · Ek B)| = |E1 | |E2 · · · Ek B| = · · · = |E1 | |E2 | · · · |Ek | |B|
= |E1 | · · · |Ek−1 Ek | |B| = · · · = |A| |B| .
puesto que la matriz A0 B tiene también alguna fila de ceros, y |AB| = 0 = |A| |B|.
Pretendemos ahora ver que el desarrollo de Laplace para el cálculo del determinante
puede hacerse a lo largo de cualquier columna o de cualquier fila y, para ello, veremos
que el determinante de una matriz coincide con el de su traspuesta.
Lema 1.32. La traspuesta de una matriz elemental, E, es una matriz elemental que
tiene mismo determinante que E.
t
A = (Ek )t · · · (E1 )t = |Ek | · · · |E1 | = |A| .
Corolario 1.34. Las proposición 1.29 es válida si se cambia la palabra fila por columna.
se tiene que:
= a11 (−1)1+1 det(A11 ) + · · · + a1n (−1)n+1 det(A1n ) = a11 α11 + · · · + a1n α1n .
|A| = a21 (−1)2+1 det(A21 ) + · · · + a2n (−1)2+n det(A2n ) = a21 α21 + · · · + a2n α2n .
|A| = a31 (−1)3+1 det(A31 ) · · · + a3n (−1)3+n det(A3n ) = a31 α31 + · · · + a3n α3n .
1
A−1 = (adj(A))t , en donde adj(A) = (αij ) ∈ Mn (K).
|A|
1
Demostración. Como |A| =
6 0, podemos definir la matriz B := (adj(A))t . La matriz
|A|
B es de orden n y vamos a comprobar que es inversa de A. En efecto,
n
P n
P 1 1
(AB)(i, i) = aik B(k, i) = aik αik = |A| = 1
k=1 k=1 |A| |A|
y, si i 6= j, se tiene que
n
P n
P 1 1 P n
(AB)(i, j) = aik B(k, j) = aik αjk = ( aik αjk ) = 0
k=1 k=1 |A| |A| k=1
n
P
ya que aik αjk = 0, porque es el desarrollo de Lapace, a lo largo de la fila j-ésima, del
k=1
determinante de una matriz que tiene todas sus filas como las de A, excepto la j-ésima
que coincide con la i-ésima de A, es decir, es el determinante de una matriz con dos
filas iguales.