Mas de 160 Problemas de Algebra Lineal
Mas de 160 Problemas de Algebra Lineal
problemas resueltos de
Álgebra Lineal
Prólogo vii
2. Matrices 25
3. Determinantes 53
Glosario 277
Bibliografía 279
v
Prólogo
Este manual recoge una colección de problemas resueltos de Álgebra Lineal; pues-
to que existen extraordinarios e imprescindibles textos (ver por ejemplo, [2], [3], [4],
[11] y [15]), nuestro principal objetivo en la elaboración del mismo ha sido acercarnos
al razonamiento que poseen los estudiantes en esta etapa del aprendizaje, lo que se
reeja en las explicaciones y detalles que aparecen en la resolución de los problemas
presentados. Se trata de un manual complementario que el alumno utilizará para re-
forzar sus conocimientos en Álgebra Lineal; por ello, a lo largo del texto aparecen
algunas deniciones y resultados que pueden facilitar al alumnado la resolución de
los problemas. Como novedad, queremos señalar que, en el último capítulo, aparecen
algunos ejercicios resueltos utilizando el programa de cálculo simbólico Mathematica.
Creemos que este manual puede ser bastante útil para los estudiantes de un primer
curso de cualquier Grado en las ramas de Ciencias o de Ingeniería.
viii
1 | Sistemas de ecuaciones
lineales
4x +3y +2z =1 x +3y +2z = 3
c) x −2y =4 d) y −3z = 1
−x +2y =2 y +2z = 4
x −2y +2z = 3 x −2y +3z =1
e) −2x +4y +z = −1
f) 2x −4y +6z =1
x +y = 3 x +y = 0
g) h)
2x −y = −9 x +z = 0
x −y = −2 y +z = 0
x −y = 2
2x +2y +2z = 1
2x +y = 1
i) j)
x +3y −z = 3
x −2y = 3
−x −3y +z = −7
x +2y = −1
x1 + x2 + x3 + x4 = 1 2x +y −z = 0
k) x1 − x2 + x3 − x4 = 0 l) 2x +3y −3z = −3
x1 + x2 − x3 − x4 = 2 2x −y +z = 3
9
10 Más de 160 problemas resueltos de Álgebra Lineal
x1 +x2 −3x3 + 5x4 =0
x2 −x3 + 5x4 =0
m)
−3x3 + 5x4 =0
x1 +3x3 − 5x4 =0
Solución:
2x −3y −4z = 3
a) 2y +z = 2
3z = −3
5x −3y +2z = 4
b)
y −2z = 3
De la segunda ecuación obtenemos y = 3 + 2z . Si sustituimos este valor en la
primera ecuación, se tiene:
5x − 3(3 + 2z) + 2z = 4.
4x +3y +2z = 1
c) x −2y = 4
−x +2y = 2
x −2y +2z = 3
e)
−2x +4y +z = −1
Sumando a la segunda ecuación el doble de la primera, se tiene 5z = 5. De donde
z = 1.
Sustituyendo el valor de z en el sistema, se tiene:
x −2y = 1
−2x +4y = −2
x = 1 + 2λ, y = λ, z = 1, λ ∈ R.
x −2y +3z = 1
f)
2x −4y +6z = 1
A la segunda ecuación le sumamos la primera multiplicada por (−2), entonces se
tiene:
x −2y +3z = 1
0 = −1
x +y = 0
h) x +z = 0
y +z = 0
x −y = 2
2x +y = 1
i)
x −2y = 3
x +2y = −1
x1 +x2 +x3 +x4 = 1
k) x1 −x2 +x3 −x4 = 0
x1 +x2 −x3 −x4 = 2
2x +y −z = 0
l) 2x +3y −3z = −3
2x −y +z = 3
3
De aquí se tiene x = . Si sustituimos este valor en las restantes ecuaciones,
4
y −z = −3/2
3y −3z = −9/2.
14 Más de 160 problemas resueltos de Álgebra Lineal
Observamos que la segunda ecuación es tres veces la primera. Por tanto, no nos
aporta ninguna información y podemos eliminarla. Así si y = λ, vemos que se
trata de un sistema compatible indeterminado, cuyo conjunto de soluciones viene
dado por:
3 3
x= , y = λ, z= + λ, λ ∈ R.
4 2
x1 +x2 −3x3 +5x4 = 0
x2 −x3 +5x4 = 0
m)
−3x3 +5x4 = 0
x1 +3x3 −5x4 = 0
Si ahora restamos las dos primeras ecuaciones, tenemos −2x3 = 0. Esto es, x3 = 0.
Si sustituimos este valor en la tercera ecuación obtenemos x4 = 0. Finalmente
sustituyendo estos valores en la primera ecuación, obtenemos x2 = 0. Así el sistema
homogéneo sólo tiene la solución trivial. Esto es, x1 = x2 = x3 = x4 = 0.
x −y +z −u +v = 0 3x +2y +z = 0
c) d)
x +y +2z +2u −v = 0 2x +3y +z = 0
−x +5y +z +7u −5v = 0 2x +y +3z = 0
2x −y +3z = 9 5x +2y −3z = −2
e) 3x −5y +z = −4 f) 2x −2y +5z = 0
4x −7y +z = 5 3x +4y +2z = −10
x +y +z +t = 6
x +y +z +t = 0
x −y +z −t = −2 x +z = 0
g) h)
3x −y +3z −t = 2
x +3y +z +3t = 0
7x −5y +7z −5t = −6 x −y +z +t = 0
x −y +t = 0 x +2y −3z +4t = 2
i) j)
−x +2y +z = 0 2x +5y −2z +t = 1
y −3z +2t = 0 5x +12y −7z +6t = 7
x +y −z +2t = 1
x −3y −2z +4t = 5
x −y +z −t = 2
k) l)
3x −8y −3z +8t = 18
2x +2y −z +5t = 0
2x −3y +5z −4t = 19
y +4z −3t = 5
Solución:
x1 −2x2 +x3 = 7
a) 2x1 −5x2 +2x3 = 6
3x1 +2x2 −x3 = 1
2x −3y +z +t = 0
x −y +z +t = 0
b)
3x −2y +z +t = 0
6x −6y +3z +3t = 0
x −y +z +t = 0
y +z +t = 0
z +t = 0
−3z −3t = 0.
x = 0, y = 0, z = −λ, t = λ, λ ∈ R.
x −y +z −u +v = 0
c) x +y +2z +2u −v = 0
−x +5y +z +7u −5v = 0
2x −y +3z = 9
e) 3x −5y +z = −4
4x −7y +z = 5
Multiplicamos por 1
2 la primera ecuación:
x − 12 y + 32 z = 9
2
3x −5y +z = −4
4x −7y +z = 5
x − 21 y + 32 z = 9
2
y +z = 5
13
y +z =
5
5x +2y −3z = −2
f) 2x −2y +5z = 0
3x +4y +2z = −10
Multiplicamos por 1
2 la nueva primera ecuación:
x −y + 52 z = 0
5x +2y −3z = −2
3x +4y +2z = −10
x = 2 − λ, y = 4 − µ, z = λ, t = µ, λ, µ ∈ R.
x +y +z +t = 0
x +z = 0
h)
x +3y
+z +3t = 0
x −y +z +t = 0
x = −λ, y = 0, z = λ, t = 0, λ ∈ R.
x −y +t = 0
i) −x +2y +z = 0
y −3z +2t = 0
x +2y −3z +4t = 2
j) 2x +5y −2z +t = 1
5x +12y −7z +6t = 7
22 Más de 160 problemas resueltos de Álgebra Lineal
x −3y −2z +4t = 5
k) 3x −8y −3z +8t = 18
2x −3y +5z −4t = 19
x +y −z +2t = 1
x −y +z −t = 2
l)
2x +2y −z +5t = 0
y +4z −3t = 5
0 1 0 1 −5 −6 −7 −8
1 3 −2 1 2
2 −2 4 2 5 1 −3 5
d) H = −1 3 4 , I=
−1 3 2 −3 −3
.
1 −2 −3
3 −4 −1 2 9
25
26 Más de 160 problemas resueltos de Álgebra Lineal
Solución: En este ejercicio utilizaremos el algoritmo que aparece en [9] para escalonar
una matriz:
1 2 3 −−−−−→ 1 2 3 −−−−−→ 1 2 3
−−−−−→
a) A = 4 5 6 F2 − 4F1 0 −3 −6 − 13 F2 0 1 2
1 F3 − F2
7 8 9 F3 − 7F1 0 −6 −12 − 6 F3 0 1 2
1 2 3
0 1 2 . El rango de A es 2.
0 0 0
0 1 2 −1 0 3 1 0 −3
−−−−−→ −−→
B = −1 0 3 0 1 2 0 1 2
F1 ↔ F2 −F1
−2 −3 0 −2 −3 0 −2 −3 0
1 0 −3 1 0 −3
−−−−−→ −−−−−→
0 1 2 0 1 2 . El rango de B es 2.
F3 + 2F1 F3 + 3F2
0 −3 −6 0 0 0
0 1 2 1 0 3 1 0 3
−−−−−→ −−−−−→ −−−−−→
C=1 0 3 0 1 2 0 1 2
F1 ↔ F2 F3 − 2F1 F3 − 3F2
2 3 0 2 3 0 0 3 −6
2. Matrices 27
1 0 3 1 0 3
−−−−−→
0 1 2 1 0 1 2 . El rango de C es 3.
− 12 F3
0 0 −12 0 0 1
2 1 3 −2 1 1 1 1 −−−−−→ 1 1 1 1
−−−−−→
b) D = 2 −1 5 2 2 −1 5 2 F2 − 2F1 0 −3 3 0
F1 ↔ F3
1 1 1 1 2 1 3 −2 F3 − 2F1 0 −1 1 −4
−−−→ 1 1 1 1 1 1 1 1
−−−−−→ 0 1 −1 4 −−−−−→
−F2 0 3 −3 0
F2 ↔ F3 F3 − 3F2
−F3 0 1 −1 4 0 3 −3 0
1 1 1 1 1 1 1 1
−−−−−→
0 1 −1 4 1 0 1 −1 4 . El rango de D es 3.
− 12 F3
0 0 0 −12 0 0 0 1
1 2 3 3 10 6 −−−−−→ 1 2 3 3 10 6
0 −3 −6 −6 −18 −9 −−−−−→
2 − 2F1
2 1 0 0 2 3 F
E=
0 −2 −4 −5 −15 −7 − 1 F2
2 2 2 1 5 5 F3 − 2F1 3
−1 1 3 2 5 2 F4 + F1 0 3 6 5 15 8
1 2 3 3 10 6 1 2 3 3 10 6
0 1 2 2 −−−−−→ −−−−−→
6 3 F3 + 2F2
0
1 2 2 6 3 −F3
0 −2 −4 −5 −15 −7 0 0 0 −1 −3 −1
F4 − 3F2
0 3 6 5 15 8 0 0 0 −1 −3 −1
1 2 3 3 10 6 1 2 3 3 10 6
−−−−−→
0 1 2 2 6 3 0 1 2 2 6 3
. El rango de E es 3.
0 0 0 1 3 1 F4 + F3 0 0 0 1 3 1
0 0 0 −1 −3 −1 0 0 0 0 0 0
1 0 0 1 1 0 0 1 1 0 0 1
−−−−−→ −−−−−→ −−−−−→
0 1 1 0 0 1 1 0 0 1 1 0
c) F =
1 0 1 0 F3 − F1 0 0 1 −1 F4 − F2 0 0 1 −1 F4 + F3
0 1 0 1 0 1 0 1 0 0 −1 1
1 0 0 1
0 1 1 0
. El rango de F es 3.
0 0 1 −1
0 0 0 0
1 2 3 4 −−−−−→ 1 2 3 4
0 −4 −8 −12 −−−1−−→
F2 − 5F1
5 6 7 8
G=
−1 −2 −3 −4 F3 + F1 0 0 0 0 − 4 F2
−5 −6 −7 −8 F4 + 5F1 0 4 8 12
28 Más de 160 problemas resueltos de Álgebra Lineal
3 41 2 1 2 3 4
−−−−−→
0
2 3 1 0 1 2 3
. El rango de G es 2.
0 0 F4 − 4F2 0 0 0 0
0 0
0
8 12 4 0 0 0 0
2 −2 4 1 −2 −3 −−−−−→ 1 −2 −3
−−−−−→
d) H = −1 3 4 −1 3 4 F2 + F1 0 1 1
F1 ↔ F3
1 −2 −3 2 −2 4 F3 − 2F1 0 2 10
1 −2 −3 1 −2 −3
−−−−−→ −−− −− → 0 1 1 . El rango de H es 3.
0 1 1 1
F3 − 2F2 8 F3
0 0 8 0 0 1
1 3 −2 1 2 −−−−−→ 1 3 −2 1 2
2 5 1 −3 5 F2 − 2F1 0 −1 5 −5 1 −−−−−→
I=
−1 3 2 −3 −3 F3 + F1 0 6 0 −2 −1 −F2
3 −4 −1 2 9 F4 − 3F1 0 −13 5 −1 3
1 3 −2 1 2 1 3 −2 1 2
−−−−−→
1 −5 5 −1 F3 − 6F2 0 1 −5 5 −1 −−− −−→
0
1
0 6 0 −2 −1 0 0 30 −32 5 30 F3
F4 + 13F2
0 −13 5 −1 3 0 0 −60 64 −10
1 3 −2 1 2 1 3 −2 1 2
0 1 −5 −1
5 −−−−−→ 0 1 −5 5 −1
.
1 − 16 1
− 16 1
0 0 15 6 F 4 + 60F 3 0 0 1 15 6
0 0 −60 64 −10 0 0 0 0 0
El rango de I es 3.
2. Resuelve:
−1 3 2 4 1 1
1 2 3 3 2 1
a) 4 5 6
+
2 3 4
, b) 4 0 1 + 3 2 −2 .
−2 1 5 1 2 −3
Solución:
1 2 3 3 2 1 1+3 2+2 3+1 4 4 4
a) + = = .
4 5 6 2 3 4 4+2 5+3 6+4 6 8 10
−1 3 2 4 1 1 −1 + 4 3 + 1 2 + 1 3 4 3
b) 4 0 1 + 3 2 −2 = 4 + 3 0 + 2 1 − 2 = 7 2 −1 .
−2 1 5 1 2 −3 −2 + 1 1 + 2 5 − 3 −1 3 2
2. Matrices 29
Solución:
(−6) · 1 (−6) · 3 (−6) · 5 (−6) · 7
(−6)A = (−6) · 2 (−6) · 6 (−6) · 10 (−6) · 14 =
(−6) · 1 (−6) · (−3) (−6) · 6 (−6) · (−7)
−6 −18 −30 −42
−12 −36 −60 −84 .
−6 18 −36 42
(3/5) · 1 (3/5) · 3 (3/5) · 5 (3/5) · 7
(3/5)A = (3/5) · 2 (3/5) · 6 (3/5) · 10 (3/5) · 14 =
(3/5) · 1 (3/5) · (−3) (3/5) · 6 (3/5) · (−7)
3/5 9/5 3 21/5
6/5 18/5 6 42/5 .
3/5 −9/5 18/5 −21/5
2 6 10 14
2A = 4 12 20 28 .
2 −6 12 −14
9 27 45 63
32 A = 9A = 18 54 90 126 .
9 −27 54 −63
1 0 1 0 0 1
1 2 1 0 −1
c) 3 4 0 1 −1
, d) 0 1 1 0 1 0 .
1 1 0 1 0 0
30 Más de 160 problemas resueltos de Álgebra Lineal
Solución:
2
3 9 −4 3 · 2 + 9 · 3 + (−4) · (−1) 37
a) 3 = = .
−9 7 11 (−9) · 2 + 7 · 3 + 11 · (−1) −8
−1
1 2 3 1 −1
b) 2 3 4 3 −2 =
4 5 6 −1 1
1 · 1 + 2 · 3 + 3 · (−1) 1 · (−1) + 2 · (−2) + 3 · 1 4 −2
2 · 1 + 3 · 3 + 4 · (−1) 2 · (−1) + 3 · (−2) + 4 · 1 = 7 −4 .
4 · 1 + 5 · 3 + 6 · (−1) 4 · (−1) + 5 · (−2) + 6 · 1 13 −8
1 2 1 0 −1 1 2 −3
c) = .
3 4 0 1 −1 3 4 −7
1 0 1 0 0 1 1 0 1
d) 0 1 1 0 1 0 = 1 1 0 .
1 1 0 1 0 0 0 1 1
calcular: A2 , B 2 , C 2 , A · B, B · A, A · C, B · C y (A − C)B.
Solución:
2 3 −2 2 3 −2 −7 −4 −17
A2 = −1 0 −3 −1 0 −3 = −14 −18 −4 .
4 5 2 4 5 2 11 22 −19
−2 3 −1 −2 3 −1 11 2 1
B2 = 4 2 1 4 2 1 = 5 14 2 .
5 −2 4 5 −2 4 2 3 9
0 6 2 0 6 2 22 8 34
C2 = 3 −1 4 3 −1 4 = 5 47 22 .
2 7 5 2 7 5 31 40 57
2. Matrices 31
2 3 −2 −2 3 −1 −2 16 −7
A · B = −1 0 −3 4 2 1 = −13 3 −11 .
4 5 2 5 −2 4 22 18 9
−2 3 −1 2 3 −2 −11 −11 −7
B·A= 4 2 1 −1 0 −3 = 10 17 −12 .
5 −2 4 4 5 2 28 35 4
−2 3 −1 0 6 2 7 −22 3
B·C = 4 2 1 3 −1 4 = 8 29 21 .
5 −2 4 2 7 5 2 60 22
2 3 −2 0 6 2 −2 3 −1
(A − C) · B = −1 0 −3 − 3 −1 4 4 2 1 =
4 5 2 2 7 5 5 −2 4
2 −3 −4 −2 3 −1 −36 8 −21
= −4 1 −7 4 2 1 = −23 4 −23 .
2 −2 −3 5 −2 4 −27 8 −16
b) (A − B)2 .
c) (A − B)(A + B).
Solución:
a) (A + B)2 = (A + B)(A + B) = A2 + AB + BA + B 2 .
b) (A − B)2 = (A − B)(A − B) = A2 − AB − BA + B 2 .
c) (A − B)(A + B) = A2 + AB − BA − B 2 .
32 Más de 160 problemas resueltos de Álgebra Lineal
1 1
7. Hallar todas las matrices que conmuten con A = 0 1
.
b11 b12
Solución: Sea B = tal que AB = BA. Si
b21 b22
1 1 b11 b12 b11 + b21 b12 + b22
AB = =
0 1 b21 b22 b21 b22
y
b11 b12 1 1 b11 b11 + b12
BA = = ,
b21 b22 0 1 b21 b21 + b22
Es decir,
b11 + b21 = b11
b12 + b22 = b11 + b12
b22 = b21 + b22
de donde
b21 = 0
b11 = b22 .
α β
Luego las matrices que conmutan con A tienen la forma , con α, β ∈ R.
0 α
de donde
1 1
B= 2 1 .
1 3
Esto es,
1 41 41 41 1 1 1
A= = .
41 82 82 123 2 2 3
34 Más de 160 problemas resueltos de Álgebra Lineal
Solución:
a b c x k
a) d e f y = g .
h i j z l
x
2 −1 −1 3
b) y = .
4 −1 2 8
z
1 −3 2 1 4 1 0 2 1 −1 −2
A= 2 1 −3 , B= 2 1 1 1 , C= 3 −2 −1 −1 .
4 −3 −1 1 −2 1 2 2 −5 −1 0
Calculemos ahora AC ,
1 −3 2 2 1 −1 −2 −3 −3 0 1
AC = 2 1 −3 3 −2 −1 −1 = 1 15 0 −5 .
4 −3 −1 2 −5 −1 0 −3 15 0 −5
12. Probar las siguientes propiedades que verican las matrices simétricas
y antisimétricas:
a) La suma de dos matrices simétricas del mismo orden es otra matriz
simétrica del mismo orden.
b) La traspuesta de una matriz simétrica es simétrica.
c) El producto de una matriz por su traspuesta es una matriz simétrica.
d) La traspuesta de una matriz antisimétrica es antisimétrica.
e) La suma de una matriz y de su traspuesta es simétrica.
f) Una matriz simétrica y antisimétrica es la matriz nula.
g) La suma de dos matrices antisimétricas es antisimétrica.
Matriz traspuesta
Dada una matriz A = (aij )i,j ∈ Mn×m (K), llamaremos matriz traspuesta
de A, a la matriz At = (aji )j,i ∈ Mm×n (K). Esto es, At será una matriz cuyas
las (columnas) son, respectivamente, las columnas (las) de la matriz A.
Matriz simétrica
Una matriz cuadrada A es simétrica si es igual a su traspuesta. Esto es,
A = At .
así como la denición de matriz antisimétrica que nos sera de utilidad en este ejercicio,
Matriz antisimétrica
Una matriz cuadrada A es antisimétrica si es igual a la opuesta de su
traspuesta. Esto es, A = −At .
b) Sea A ∈ Mn (K) una matriz simétrica. Esto es, A = At y por tanto At = (At )t , lo
que nos dice que At es matriz simétrica.
Solución:
14. Se pide:
a) Demostrar que toda matriz se puede descomponer de manera única
como suma de una matriz simétrica y otra antisimétrica.
b) Descomponer la matriz
3 5 1
H= 6 0 4
7 5 5
Solución:
Por tanto,
0 − 12 −3 3 11
2 4
1 1
A = (H − H t ) = 12
0 − 21
y S= (H + H t ) =
11
2 0 9
2 .
2 2
3 12 0 4 9
2 5
Matriz idempotente
I − 2A + A2 = 1 − 2A + A = I − A = B.
Y, en consecuencia, B es idempotente.
Veamos ahora que AB = BA = θ.
AB = A(I − A) = AI − AA = AI − A2 = A − A = θ.
BA = (I − A)A = IA − AA = IA − A2 = A − A = θ.
1 1 1 1 0 0
16. Calcula A2 , A3 , A4 , An , siendo A = 0 1 1 . Idem para B = 1 1 0 .
0 0 1 1 1 1
Solución:
1 1 1 1 1 1 1 2 3
A2 = 0 1 1 0 1 1 = 0 1 2 .
0 0 1 0 0 1 0 0 1
1 1 1 1 2 3 1 3 6
A3 = AA2 = 0 1 1 0 1 2 = 0 1 3 .
0 0 1 0 0 1 0 0 1
1 1 1 1 3 6 1 4 10
A4 = AA3 = 0 1 1 0 1 3 = 0 1 4 .
0 0 1 0 0 1 0 0 1
n
1 n 2 (n+ 1)
Se tiene que An = 0 1 n .
0 0 1
Solución: En este ejercicio utilizaremos el algoritmo, que aparece en [9], para obtener
la matriz inversa de una dada.
1 0 0 1
0 1 1 0
a) A =
0
.
0 1 0
0 1 0 1
2. Matrices 41
Por tanto:
1 1 −1 −1
0 1 −1 0
A−1 =
0
.
0 1 0
0 −1 1 1
1 2 3 4
5 6 7 8
b) B =
−1
.
−2 −3 −4
−5 −6 −7 −8
2 −2 −4
c) C = −1 3 4 .
1 −2 −3
42 Más de 160 problemas resueltos de Álgebra Lineal
1 1 2 1 1 2
18. Siendo: A = 3 4
,B=
1 1
,C=
1 3
.
Solución: En primer lugar calcularemos A−1 ya que nos será necesario para resolver
el ejercicio. Por el método de Gauss-Jordan, se tiene que:
−1 4 −1
A = .
−3 1
a b
d) Si XA + BX = C , supongamos que X = , entonces
c d
a b 1 1 2 1 a b 1 2
+ = .
c d 3 4 1 1 c d 1 3
Esto es,
a + 3b a + 4b 2a + c 2b + d 1 2
+ = .
c + 3d c + 4d a+c b+d 1 3
Así,
3a + 3b + c a + 6b + d 1 2
= .
a + 2c + 3d b + c + 5d 1 3
1 2 2
20. Calcula A2 , siendo A = 13 2 −2 1 . ¾Cuál es su inversa?
2 1 −2
2 1
21. Demuestra que la matriz A = 1 2
verica una ecuación del tipo
A2 + αA + βI = θ, determinando α, β . ¾Es A regular?
Luego:
5 4 2α α β 0
A2 + αA + βI = ++ =
4 5 α 2α 0 β
5 + 2α + β 4+α 0 0
= .
4+α 5 + 2α + β 0 0
Solución:
Denotemos por M , la matriz que representa los tornillos y soportes que lleva cada
estantería. Así, las columnas corresponden a tornillos y soportes, respectivamente,
y las las a los tamaños (grande, pequeño).
tornillos soportes
M = grande 16 6 .
pequeño 12 4
b) La matriz pedida, que denotaremos por P , será el producto de las dos matrices
anteriores. Esto es:
tornillos soportes
A 112000 38000
P =D·M =
B
.
200000 72000
C 136000 48000
Solución:
Administración Taller
C =P ·H = A 705 17650 .
B 459 11490
Teorema de Rouché-Frobenius
Dado un sistema de m ecuaciones lineales con n incógnitas, con matriz de
coecientes A y matriz ampliada (A|B), se tiene que:
1. El sistema es compatible si, y sólo si, rg(A) = rg(A|B).
2. El sistema es compatible determinado si, y sólo si, rg(A) = rg(A|B) = n.
Conviene observar que si B = θ, entonces el sistema siempre tiene solución
y, si la solución es única, ésta es la solución trivial.
x −2y +z = 5
a) 2x −y −2z = −1
x +3y +z = 0
48 Más de 160 problemas resueltos de Álgebra Lineal
Se tiene que rg(A) = 3 =rg(A|B), que coincide con el número de incógnitas. Así
el sistema es un sistema compatible determinado. Para resolverlo, consideramos el
sistema correspondiente a la matriz resultante de efectuar operaciones elementales
en la matriz (A|B). Este sistema es:
x −2y +z = 5
y = −1
z = 2
x = 1, y = −1, z = 2.
3x −4y +6z = 0
b) 5x +2y −4z = 0
x +3y −5z = 0
x +2y +2z = 0
3x +y +3z = 0
c)
2x +y −2z = 0
2x +y +z = 0
Se tiene que rg(A) = 3, que coincide con el número de incógnitas. Así el sistema
tendrá sólo la solución trivial. Esto es, x = y = z = 0.
2x −y +z = 5
−3x +2y +2z = 0
d) −x −2y +2z = 0
x +y +z = 0
5x −6y −2z = 0,
Se tiene que rg(A) = 3, mientras que rg(A|B) = 4. Por tanto el sistema es incom-
patible.
2. Matrices 51
x +y −z = 3
e) 2x −y +4z = 3
3x +2y −z = 8
1 1 −1 3 1 1 −1 3
A= 2 −1 4 3 , (A|B) = 2 −1 4 3 .
3 2 −1 8 3 2 −1 8
x +y −z = 3
y −2z = 1,
x = 2 − λ, y = 1 + 2λ, z = λ, λ ∈ R.
−2x +y +z = 1
f) x −2y +z = −2
x +y −2z = 4
Este capítulo está dedicado al estudio de los determinantes de una matriz cuadrada.
Así:
Calcularemos el determinante de una matriz cuadrada.
Aplicaremos las propiedades de los determinantes.
Veremos dos aplicaciones del cálculo de determinantes: cálculo del rango de una
matriz y cálculo de la inversa de una matriz regular.
Estudiaremos y resolveremos, aplicando en algunos casos la Regla de Cramer,
sistemas de ecuaciones lineales.
Discutiremos sistemas de ecuaciones dependientes de un parámetro.
1 3 0 5 −1 7 1 2 3
a) −1 2 −4 , b) 6 4 3 , c) 4 5 6 .
1 1 2 3 2 1 7 8 9
Regla de Sarrus
a11 a12 a13
Dada una matriz cuadrada A = a21 a22 a23 , se tiene que:
a31 a32 a33
|A| = (a11 a22 a33 +a12 a23 a31 +a13 a32 a21 )−(a13 a22 a31 +a12 a21 a33 +a11 a23 a32 ).
53
54 Más de 160 problemas resueltos de Álgebra Lineal
1 3 0
a) −1 2 −4 = (4 + 0 − 12) − (0 − 4 − 6) = 2.
1 1 2
5 −1 7
b) 6 4 3 = (20 + 84 − 9) − (84 + 30 − 6) = −13.
3 2 1
1 2 3
c) 4 5 6 = (45 + 96 + 84) − (105 + 48 + 72) = 0.
7 8 9
Solución: Para resolver este ejercicio, nos será útil recordar las siguientes deniciones
y propiedades que pueden ser consultadas en [9] y [7] entre otros.
Comenzamos recordando la denición de determinante de una matriz cuadrada.
Esta denición se hará de forma inductiva.
0 −2 1 9
−2 1 9
0 1 4 −9
= (1)(−1)4+1 1 4 −9 ,
0 1 −3 5
1 −3 5
1 0 1 −2
56 Más de 160 problemas resueltos de Álgebra Lineal
Podemos volver a descomponer este nuevo determinante como suma de dos. Esto
es,
Repitiendo el proceso,
1 2 3 1234 1 2 3 1230 1 2 3 4
2 3 4 2341 2 3 4 2340 2 3 4 1
= + .
3 4 1 3412 3 4 1 3410 3 4 1 2
4 1 2 4123 4 1 2 4120 4 1 2 3
1 2 3 1234 1 2 3 4
2 3 4 2341 2 3 4 1
= .
3 4 1 3412 3 4 1 2
4 1 2 4123 4 1 2 3
1 2 3 4 10 2 3 4 1 2 3 4
2 3 4 1 10 3 4 1 1 3 4 1
= = 10 .
3 4 1 2 10 4 1 2 1 4 1 2
4 1 2 3 10 1 2 3 1 1 2 3
1 2 3 4 1 2 3 4
1 3 4 1 0 1 1 −3
10 = 10 .
1 4 1 2 0 2 −2 −2
1 1 2 3 0 −1 −1 −1
1 2 3 4
1 1 −3
0 1 1 −3
10 = 10 · 1 · (−1)1+1 2 −2 −2 = 160.
0 2 −2 −2
−1 −1 −1
0 −1 −1 −1
58 Más de 160 problemas resueltos de Álgebra Lineal
Solución:
b) Observamos que la quinta columna tiene todos sus elementos iguales. Por ello,
aplicando las propiedades de los determinantes y desarrollando el determinante
por la última columna, se tiene
1 2 3 4 5 1 2 3 4 1 1 2 3 4 1
1 0 0 0
2 2 3 4 5 2 2 3 4 1 1 0 0 0 0
2 1 0 0
3 3 3 4 5 =5 3 3 3 4 1 =5 2 1 0 0 0 =5 = 5.
3 2 1 0
4 4 4 4 5 4 4 4 4 1 3 2 1 0 0
4 3 2 1
5 5 5 5 5 5 5 5 5 1 4 3 2 1 0
4. Calcula:
1 1 1 1
1 1 1
1 1 a b c d
a) a b
, b) a b c , c) a2 b2 c2 d2
.
a2 b2 c2
a3 b3 c3 d3
Solución:
1 1
a) = b − a.
a b
b) En este apartado y en los que siguen, escribiremos debajo del signo = la operación
elemental de las que se está realizando.
1 1 1 1 1 1 1 1 1
= =
a b c a b c 0 b−a c−a =
F3 − aF2 F2 − aF1
a2 b2 c2 0 b2 − ab c2 − ac 0 b2 − ab c2 − ac
b−a c−a 1 1
= = (b − a)(c − a) = (b − a)(c − a)(c − b).
b(b − a) c(c − a) b c
1 1 1 1 1 1 1 1
a b c d = a b c d =
c)
a2 b2 c2 d2 F4 − aF3 a2 b2 c2 d2 F3 − aF2
a3 b3 c3 d3 0 b3 − ab2 c3 − ac2 d3 − ad2
1 1 1 1 1 1 1 1
a b c d = 0 b−a c−a d−a
= =
0 b2 − ab c2 − ac d2 − ad F2 − aF1 0 b2 − ab c2 − ac d2 − ad
0 b3 − ab2 c3 − ac2 d3 − ad2 0 b3 − ab2 c3 − ac2 d3 − ad2
60 Más de 160 problemas resueltos de Álgebra Lineal
Solución:
a) Sabemos que si multiplicamos una la o columna de una matriz por un número,
el determinante queda multiplicado por ese número. Así,
|5A| = 53 |A| = 250.
Solución:
a) En primer lugar vemos que AB es una matriz invertible ya que |AB| = |A|·|B| =
6 0.
Veamos ahora que (AB)−1 = B −1 A−1 . Para ello basta ver que (B −1 A−1 )(AB) =
(AB)(B −1 A−1 ) = I . Vamos a verlo:
De la misma forma:
1 −1 1 1 1
A (λA) = (A−1 λ)A = (λA−1 )A = λ (A−1 A) = 1I = I.
λ λ λ λ
= A−1 −1 −1
n+1 An ·. . . ·A1 .
Solución:
a) =⇒ Tendremos que ver que A−1 = (A−1 )t . Esto es, (A−1 )t A = I = A(A−1 )t .
Utilizando las propiedades de la trasposición de matrices y que A es simétrica,
se tiene que,
A(A−1 )t = A(At )−1 = AA−1 = I.
Solución:
Solución:
= A3 − 4A2 + 5A − 2In = θn .
Por tanto,
A3 − 4A2 + 5A = 2In .
Así,
1 3
(A − 4A2 + 5A) = In .
2
Esto es,
1 2
A (A − 4A + 5In ) = In .
2
1
Luego A−1 = (A2 − 4A + 5In ).
2
Solución:
1 2 3 4
2 3 5 −2
a) A =
3
. Como |A| = 0, A no tiene inversa.
5 1 1
1 1 −5 −7
66 Más de 160 problemas resueltos de Álgebra Lineal
1 2 3 1
1 3 3 2
b) A =
2
. Como |A| = −1, A tiene inversa. Vamos a calcularla.
4 3 3
1 1 1 1
La expresión de la matriz inversa de una matriz regular viene dada por A−1 =
1
· (Adj(A))t , siendo Adj(A) = (αij )ij con αij = (−1)i+j |Aij |.
|A|
3 3 2 1 3 2
α11 = (−1)1+1 4 3 3 = −1, α12 = (−1)1+2 2 3 3 = −1,
1 1 1 1 1 1
1 3 2 1 3 3
α13 = (−1)1+3 2 4 3 = 0, α14 = (−1)1+4 2 4 3 = 2,
1 1 1 1 1 1
2 3 1 1 3 1
α21 = (−1)2+1 4 3 3 = 2, α22 = (−1)2+2 2 3 3 = 2,
1 1 1 1 1 1
1 2 1 1 2 3
α23 = (−1)2+3 2 4 3 = −1, α24 = (−1)2+4 2 4 3 = −3,
1 1 1 1 1 1
2 3 1 1 3 1
α31 = (−1)3+1 3 3 2 = −1, α32 = (−1)3+2 1 3 2 = −2,
1 1 1 1 1 1
1 2 1 1 2 3
α33 = (−1)3+3 1 3 2 = 1, α34 = (−1)3+4 1 3 3 = 2,
1 1 1 1 1 1
2 3 1 1 3 1
α41 = (−1)4+1 3 3 2 = 0, α42 = (−1)4+2 1 3 2 = 3,
4 3 3 2 3 3
1 2 1 1 2 3
α43 = (−1)4+3 1 3 2 = −1, α44 = (−1)4+4 1 3 3 = −3.
2 4 3 2 4 3
Luego,
−1 −1 0 2 −1 2 −1 0
2 2 −1 −3 −1 2 −2 3
Adj(A) = y Adj(A)t =
.
−1 −2 1 2 0 −1 1 −1
0 3 −1 −3 2 −3 2 −3
3. Determinantes 67
Por tanto,
1 −2 1 0
1 −2 2 −3
A−1 = .
0 1 −1 1
−2 3 −2 3
Solución:
Vemos que es un menor no nulo. Como no existe un menor de orden mayor cuyo
determinante no sea nulo, podemos armar que el rango de la matriz B es 3.
68 Más de 160 problemas resueltos de Álgebra Lineal
16. Estudiar, y resolver en su caso, cada uno de los sistemas lineales que
siguen:
2x +3y −4z = 4 4x +3y −z = −7
a) b) 5x
5x +y −7z = 0 +y +2z = −3
−3x +2y +3z = 4 8x +6y −2z = −14
x1 −x2 +x3 = 4
2x +y −z = −3
3x1 +2x2 −3x3 = 1
c)
−y + 21 z = −2
x
2x1 +3x2 −4x3 = −3 d)
−x +4y − 52 z = 1
x1 − 21 x2 + 23 x3 = 11
2
6x +6y −3z = 0
2x −3y = −4 x1 −3x2 +x3 = 2
f) 2x1
e) 12 x − 34 y = −1 +x2 +4x3 = −2
5x1 −x2 +9x3 = −2
1
5x +2y =
3
x1 −2x2 +3x3 +x4 = 1
x1 +x2 +x3 +x4 = 1
3x1 +2x2 −x3 = 1 x1 −x2 +x3 −x4 = 0
g) −3x +x h)
1 2 +2x3 −3x4 = 1 x1
+x2 −x3 −x4 = −1
x1 +x2 +4x3 −2x4 = 1 x1 +x2 +x3 −x4 = 2
x −2y +z −t = 1
x +y +z +t = −6
2x +y −3z +t = 0 −2x +2y −z = −2
i) j)
x −y +2z +3t = 1
−2y +3z −2t = 9
2x +y +3z −4t = 3 −x +y +3z −t = 0
Solución:
2x +3y −4z = 4
a) 5x +y −7z = 0
−3x +2y +3z = 4
Si llamamos z = λ, se tiene
2x +3y = 4 + 4λ
5x +y = 7λ
4x +3y −z = −7
b) 5x +y +2z = −3
8x +6y −2z = −14
x1 −2x2 +x3 = 4
3x1 +2x2 −3x3 = 1
c)
2x1 +3x2 −4x3 = −3
x1 − 12 x2 + 32 x3 = 11
2
2 1 −1 −3
1
1 −1 2 −2
|(A|B)| = = −12.
−1 4 − 52 1
6 6 −3 0
Esto implica que el rango de la matriz ampliada es 4. Este rango, podemos armar,
es distinto del rango de la matriz de los coecientes puesto que éste es como
máximo 3, ya que se trata de una matriz de orden 4 × 3. Así rg(A) 6= rg(A|B). En
consecuencia el sistema es incompatible.
2x −3y
= −4
e) 1 3
2x −4y = −1
1
5x +2y =
3
2 −3 −4
2 −3
1
A= − 34 , (A|B) = 1
− 43 −1 .
2 2
5 2 5 2 1
3
Así podemos armar que rg(A) = 2. Ahora calculamos el rango de (A|B). Para
ello calculamos su determinante. Esto es,
2 −3 −4
1
|(A|B)| = 2 − 34 −1 = 0.
1
5 2 3
Luego rg(A) = rg(A|B) = 2, que coincide con el número de incógnitas. Por ello
el sistema es compatible determinado. Para resolver el sistema consideraremos un
subsistema del anterior tal que la matriz de coecientes y matriz ampliada tengan
72 Más de 160 problemas resueltos de Álgebra Lineal
cuya solución es
7 62
x=− , y= .
19 57
x1 −3x2 +x3 = 2
f) 2x1 +x2 +4x3 = −2
5x1 −x2 +9x3 = −2
diferente de 0, se tiene que rg(A) = 2. Además, como todos los menores de orden
3 de la matriz ampliada tienen determinante nulo, se tiene que rg(A|B) = 2. Esto
es, rg(A) = 2 = rg(A|B), que es menor que el número de incógnitas. Por tanto, el
sistema es compatible indeterminado. Para resolverlo, consideremos un subsistema
cuya matriz de coecientes y matriz ampliada tengan rango 2. Esto es,
x1 −3x2 +x3 = 2
2x1 +x2 +4x3 = −2
Si llamamos x3 = λ, se tiene
x1 −3x2 = 2−λ
2x1 +x2 = −2 − 4λ
−6 − 2λ
Así x2 = . Por tanto, el conjunto de soluciones del sistema viene dado
7
por:
−4 − 13λ −6 − 2λ
x1 = , x2 = , λ ∈ R.
7 7
x1 −2x2 +3x3 +x4 =1
3x1 +2x2 −x3 =1
g)
−3x1 +x2 +2x3 −3x4 =1
x1 +x2 +4x3 −2x4 =1
Así podemos armar que la matriz de coecientes tiene rango menor que 4. Sin
embargo si calculamos el siguiente menor de la matriz ampliada, se tiene:
1 −2 3 1
3 2 −1 1
|(A|B)| = = −76.
−3 1 2 1
1 1 4 1
Esto implica que la matriz ampliada tiene rango 4, distinto del rango de la matriz
de coecientes. Por tanto el sistema es incompatible.
x1 +x2 +x3 +x4 = 1
x1 −x2 +x3 −x4 = 0
h)
x1 +x2 −x3 −x4 = −1
x1 +x2 +x3 −x4 = 2
Por tanto, rg(A) = 4. En este caso podemos armar que el rango de (A|B) es
también 4, ya que es una matriz de orden 4 × 5 y la matriz de coecientes es una
74 Más de 160 problemas resueltos de Álgebra Lineal
submatriz de ella con determinante no nulo. Así rg(A) = rg(A|B) = 4, que coincide
con el número de incógnitas. Por ello, el sistema es compatible determinado. Este
tipo de sistemas reciben el nombre de sistemas de Cramer. A continuación vamos
a recordar su denición que puede ser consultada en [1].
Sistemas de Cramer
Un sistema de ecuaciones lineales AX = B se llama sistema de Cramer
si su matriz de coecientes A es cuadrada y regular. Esto es, si el sistema
tiene tantas ecuaciones como incógnitas y el determinante de su matriz
de coecientes es no nulo. Así los sistemas de Cramer son de la forma:
a11 x1 + a12 x2 + · · · + a1n xn = b1
a21 x1 + a22 x2 + · · · + a2n xn = b2
..
.
an1 x1 + an2 x2 + · · · + ann xn = bn
Regla de Cramer
1 1 1 1 1 1 1 1
0 −1 1 −1 1 0 1 −1
−1 1 −1 −1 1 −1 −1 −1
2 1 1 −1 8 1 2 1 −1 −8
x1 = = = −1, x2 = = = 1,
|A| −8 −8 −8
3. Determinantes 75
1 1 1 1 1 1 1 1
1 −1 0 −1 1 −1 1 0
1 1 −1 −1 1 1 −1 −1
1 1 2 −1 −12 3 1 1 1 2 4 1
x3 = = = , x4 = = =− .
−8 −8 2 −8 −8 2
Por tanto, la solución del sistema es:
3 1
x1 = −1, x2 = 1, x3 = , x4 = − .
2 2
x −2y +z −t = 1
2x +y −3z +t = 0
i)
x −y +2z +3t = 1
2x +y +3z −4t = 3
1 −2 1 −1 1 −2 1 1
2 1 0 1 2 1 −3 0
1 −1 1 3 1 −1 2 1
2 1 3 −4 −61 61 2 1 3 3 18 9
z= = = , t= = =− .
−152 −152 152 −152 −152 76
76 Más de 160 problemas resueltos de Álgebra Lineal
x +y +z +t = −6
−2x +2y −z = −2
j)
−2y +3z −2t = 9
−x +y +3z −t = 0
Así podemos armar que rg(A) < 4. Si ahora consideramos una submatriz de la
matriz ampliada,
1 1 1 −6
−2 2 0 −2
0 −2 −2 9
−1 1 −1 0
se tiene que el determinante de esta submatriz es no nulo. Esto signica que el rango
de la matriz ampliada (A|B) es 4. Por tanto podemos armar que rg(A) 6= rg(A|B),
y por tanto el sistema es incompatible.
x −2y +3z = 0
a)
2x +y −2z = 0
1 −2 3 1 −2
La matriz de coecientes A = tiene rango 2 puesto que =
2 1 −2 2 1
5 6= 0.
Por tanto, existen soluciones diferentes de la trivial. Vamos a calcularlas. Para ello,
a la segunda ecuación le sumamos la primera multiplicada por (−2), quedando,
x −2y +3z = 0
5y −8z = 0.
2x −y +3z = 0
b)
2x −y +2z = 0
2 −1 3 −1 3
La matriz de coecientes A = tiene rango 2 puesto que =
2 −1 2 −1 2
1 6= 0.
x = λ, y = 2λ, z = 0, λ ∈ R.
2x −3y +z = 0
c) x +2y +z = 0
5x −4y +3z = 0
78 Más de 160 problemas resueltos de Álgebra Lineal
x +2y −3z = 0
d) 2x −y +2z = 0
3x +y +2z = 0
−2 −14 6
|A| = −1 −7 3 = 0.
2 −3 1
cuyo determinante es no nulo, armamos que rg(A) = 2, siendo, por tanto, menor
que el número de incógnitas. De ahí que el sistema tenga soluciones diferentes de
la trivial. Vamos a resolverlo. Para ello, consideramos un subsistema, cuya matriz
de coecientes tenga rango 2. Esto es,
−x −7y +3z = 0
2x −3y +z = 0.
2x1 +x2 −x3 +x4 = 0
f) x1 +2x2 +x3 −x4 = 0
3x1 −x2 −2x4 = 0
x1 +2x2 +x3 −x4 = 0
x2 +x3 −x4 = 0
−7x2 −3x3 +x4 = 0.
son nulos. Sin embargo, podemos encontrar un menor de orden 3, no nulo. Esto
es,
0 −1 2
1 0 0 = 1 6= 0.
0 1 −1
3. Determinantes 81
Llamando x4 = α, x5 = β, x6 = γ , se tiene
x1 = −β − γ
−x2 +2x3 = −α + β − 2γ
x3 = −2α + β − 2γ
x6 = γ, α, β, γ ∈ R.
Solución:
x −3y = 4
a) 2x +5y = 1
x −my = 0
Tenemos que imponer que rg(A|B) = 2. Para ello debe vericarse que |(A|B)| = 0.
Esto es
1 −3 4
|(A|B)| = 2 5 1 = −23 − 7m = 0.
1 −m 0
Por tanto, m = − 23
7 . Resolviendo el sistema para m = − 7 , se tiene que la solución
23
x +2y = 3
b) x −3y = 1
2x +y = m
Ahora imponemos a la matriz ampliada tener rango 2. Para ello debe vericarse
que |(A|B)| = 0. Es decir:
1 2 3
|(A|B)| = 1 −3 1 = −5m + 24 = 0.
2 1 m
Por tanto, m = 5 .
24
En este caso la solución del sistema es:
11 2
x= , y= .
5 5
x −2y +3z =0
3x −4y −z =0
x +y +mz =0
Si el sistema tiene solución diferente de la trivial, debe vericarse que rg(A) < 3, ya
que 3 es el número de incógnitas. Para ello |A| = 0. Como |A| = 24 + 2m, se tiene que
la solución será distinta de (0, 0, 0) si m = −12.
m = −1 − 2i, m = −1 + 2i, m = 2.
En cuyo caso, se tiene que rg(A|B) = 2. Así rg(A) 6= rg(A|B), y de ahí armamos
que el sistema es incompatible.
Si p = −2, la matriz de coecientes queda:
1 1 −2
A= 1 −2 1 ,
−2 1 1
1 1 1
1 −2 2 = −12 6= 0.
−2 1 1
Solución: Al ser el sistema homogéneo, siempre se tiene que rg(A) = rg(A|B). Por
ello, sólo estudiaremos rg(A). Para ello, comenzamos calculando
1 −5 0
|A| = 4 −6 2 = −42 + 14p.
5 −4 p
Entonces, |A| = 0 si, y sólo si, p = 3. Así, tendremos que considerar los siguientes
casos:
Si p 6= 3, entonces |A| =
6 0. Por tanto, el rango de A es 3, que coincide con el
número de incógnitas, y en este caso la solución trivial es la única solución del
sistema. Esto es, x = y = z = 0.
Si p = 3, entonces |A| = 0, pero podemos extraer de la matriz A, una submatriz
0 1 −1 0
A = ,
4 −6 2
que tiene rango 2. En este caso se vericará que el sistema tiene soluciones
diferente de la trivial. Vamos a resolverlo. Para ello, consideramos un sistema
cuya matriz de coecientes sea la matriz A0 , esto es,
x −5y = 0
4x −6y +2z = 0.
Multiplicamos por 14 ,
1
la segunda ecuación,
x −5y = 0
y + 71 z = 0.
Se tiene que |A| = 2 − 3a + a3 . Por tanto, |A| = 0 si, y sólo si, a = 1 ó a = −2.
Si a 6= 1 y a 6= −2, entonces rg(A) = rg(A|B) = 3. Así, en este caso, el sistema
es compatible determinado.
Si a = 1, entonces
1 1 1 3
(A|B) = 1 1 1 3 .
1 1 1 3
Observamos que todas las las de la matriz son iguales. Por tanto, rg(A) =
rg(A|B) = 1, que al ser menor que el número de incógnitas, se tiene que el
sistema es compatible indeterminado.
Si a = −2, entonces la matriz de coecientes,
−2 1 1
A= 1 −2 1 ,
1 1 −2
tiene determinante nulo. Así rg(A) < 3. Pero podemos considerar la submatriz
−2 1
1 −2
cuyo determinante es no nulo. Por ello, podemos armar que rg(A) = 2. Estu-
diemos ahora el rango de la matriz ampliada,
−2 1 1 3
(A|B) = 1 −2 1 3 .
1 1 −2 0
con determinante no nulo. Por ello, podemos armar que rg(A|B) = 3. Así
tenemos que rg(A) = 2 6= 3 = rg(A|B). Por tanto, el sistema es incompatible.
88 Más de 160 problemas resueltos de Álgebra Lineal
−1 −2
donde se tiene que |A| = 0, pero la submatriz , tiene determinante
1 1
no nulo. Por ello podemos armar que rg(A) = 2.
Consideremos ahora la matriz ampliada,
−1 0 −2 6
(A|B) = 1 0 1 2 .
−1 0 −2 1
cuyo determinante es no nulo. Por ello podemos armar que rg(A|B) = 3. Así
concluimos que rg(A|B) = 3 6= 2 = rg(A). En consecuencia el sistema es incom-
patible.
3. Determinantes 89
25. Discutir los siguientes sistemas en función de los parámetros que apa-
recen y resolver por Cramer si el sistema es compatible determinado:
3x −y = ax 2x −ay +bz = 4
a) 5x +y +2z = ay b) x +z = 2
4y +3z = az x +y +z = 2
ax +2z = 2 2x +ay +z = a − 2
c) 5x +2y = 1 d) ax −2ay +3z = 0
x −2y +bz = 3 2x +8y −z = 2a − 4
x +3y −az = 4
−ax +y +az = 0
e)
−x +2ay = a+2
2x −y −2z = 0
Solución:
3x −y = ax
a) 5x +y +2z = ay
4y +3z = az
2x −ay +bz = 4
b) x +z = 2
x +y +z = 2
ax +2z = 2
c) 5x +2y = 1
x −2y +bz = 3
2 0 2 a 2 2
1 2 0 5 1 0
3 −2 b 2b − 8 1 3 b ab − 10b + 28
x= = , y= = ,
2ab − 24 ab − 12 2ab − 24 2ab − 24
a 0 2
5 2 1
1 −2 3 4a − 12
z= = .
2ab − 24 ab − 12
Para b 6= 0 y a = 12
b , se tiene que rg(A) = 2. Vamos a calcular el rango de la
matriz ampliada. Para ello consideremos el menor,
0 2 2
2 0 1 = 4b − 16.
−2 b 3
2x +ay +z = a−2
d) ax −2ay +3z = 0
2x +8y −z = 2a − 4
1 2
Solución: La matriz A de coecientes tiene el menor = 5 6= 0. Por tanto,
−1 3
podemos armar que rg(A) = 2. Estudiemos ahora el rango de la matriz ampliada
(A|B). Le hemos de imponer que rg(A|B) = 2. Esto es,
a 2 1 1 2 3
1 2 3 = −1 3 2 = 0.
−1 3 2 2 b 0
2 4
Como la submatriz tiene determinante no nulo, entonces rg(A) = 2. Estu-
1 1
diemos ahora el rango de la matriz ampliada,
1 m 2
(A|B) = 2 4 1 .
1 1 n
Se tiene que |(A|B)| = −5+m+4n−2mn. Esto es, |(A|B)| = 0 si, y sólo si, m = 4n−5
2n−1
con n 6= 12 .
Así tendremos que distinguir los siguientes casos:
Si n 6= 12 , m 6= 4n−5
2n−1 , entonces rg(A) = 2 y rg(A|B) = 3. En consecuencia el
sistema es incompatible.
Si n 6= 21 , m = 4n−5
2n−1 , entonces rg(A|B) = 2 = rg(A). Luego el sistema es
compatible determinado.
Si n = 12 , entonces el determinante de la matriz ampliada será
1 m 2
|(A|B)| = 2 4 1 = −3 6= 0.
1
1 1 2
1
Luego |A| = 0, si, y sólo si, a = 0 ó a = −1 ó a = .
2
Vamos a distinguir los siguientes casos:
a 6= 0, a 6= −1 y a 6= 12 . En este caso, como |A| =
6 0, el rango de A es 3, que debe
coincidir necesariamente con el de la matriz ampliada, y además es el número
de incógnitas. Luego, el sistema es compatible determinado.
a = 0. Sustituyendo el valor de a en el sistema, la matriz ampliada es
1 −1 −1 2
(A|B) = 0 0 1 0 .
1 −1 −1 1
−2 0
Observamos que la submatriz de A, tiene determinante no nulo, y
−1 1
en consecuencia rg(A) = 2. La matriz ampliada es
0 −2 0 2
(A|B) = 0 −1 1 0 .
0 −3 −2 1
3. Determinantes 97
3
− 12
Análogamente a los casos anteriores, existe un menor, 2
1
= 43 , no nulo
0 2
de A, luego rg(A) = 2 y
3
2 − 12 2
1 3
0 2 0 =−
3
4
2 0 1
Solución:
√
a) Z[ 5] no tiene estructura de espacio vectorial sobre Q, ya que no es posible denir
la ley de composición externa,
√ √
Q × Z[ √ 5] −→ Z[ 5] √
(α, a + b 5) 7−→ αa + αb 5,
√
b) Denimos una operación interna en Q[ 5]:
√ √ √
Q[ 5] ×√Q[ 5] √ −→ Q[ 5] √
(a1 + b1 5, a2 + b2 5) 7−→ (a1 + a2 ) + (b1 + b2 ) 5.
√ √ √
Conmutativa: Dados a1 + b1 5, a2 + b2 5 ∈ Q[ 5], se tiene que:
√ √ √
(a1 + b1 5) + (a2 + b2 5) = (a1 + a2 ) + (b1 + b2 ) 5 =
√ √ √
(a2 + a1 ) + (b2 + b1 ) 5 = (a2 + b2 5) + (a1 + b1 5).
Existencia
√ √de elemento neutro: El elemento neutro es 0, ya que dados 0, a +
b 5 ∈ Q[ 5], se tiene que,
√ √ √
0 + (a + b 5) = (a + b 5) + 0 = a + b 5.
√ √ √
Existencia
√ de elemento simétrico: Dado a + b 5 ∈ Q[ 5], existe −a − b 5 ∈
Q[ 5] tal que:
√ √ √
a + b 5 + (−a − b 5) = (a − a) + (b − b) 5 = 0.
√ √
Por tanto, el elemento simétrico respecto de la suma de a + b 5 es −a − b 5.
Denimos ahora la siguiente operación externa:
√ √
Q × Q[ √ 5] −→ Q[ 5] √
(α, a + b 5) 7−→ αa + αb 5.
Demostremos que la suma usual de funciones es una operación interna en F . Esto es,
dadas f , g ∈ F , se tiene que f + g ∈ F . Para ello tendrá que vericarse:
2(f + g)(0) = (f + g)(1).
Veámoslo:
2(f + g)(0) = 2[f (0) + g(0)] = 2f (0) + 2g(0) = f (1) + g(1) = (f + g)(1).
Es fácil comprobar que (F, +) es un grupo abeliano. Esto es, la operación + verica
las propiedades: asociativa, conmutativa, existencia de elemento neutro (la función
nula) y existencia de elemento simétrico, ya que dado f ∈ F , existe −f ∈ F tal que
f + (−f ) = 0, donde (−f )(x) := −f (x).
Consideramos ahora, el producto usual de una función por un número real,
R×F −→ F
(α, f ) 7−→ (α · f )(x) := αf (x).
1. (α + β)f = αf + βf .
2. α(f + g) = αf + αg .
3. (αβ)f = α(βf ).
4. 1R f = f .
Por tanto, podemos concluir que (F, +, ·R) es un espacio vectorial.
Solución:
1 a
a) Consideremos el conjunto, V1 = : a, b ∈ R . Se tiene que la suma
b 1
1 1
usual de matrices no es una operación interna en V1 , ya que dadas ,
2 1
1 3
∈ V1 , se verica que
0 1
1 1 1 3 2 4
+ = ∈
/ V1 .
2 1 0 1 2 2
se tiene que
0 −a
∈ V2 ,
−b 0
Dígase, para cada uno de los casos anteriores, si se obtiene o no una es-
tructura de espacio vectorial en R2 .
Solución:
5. Se pide:
a) ¾Son los vectores (1, 2, 3) y (1, 1, 1) combinación lineal de los vectores
del sistema S = {(1, 0, 1), (0, 2, 2)}?
104 Más de 160 problemas resueltos de Álgebra Lineal
Solución:
−11 a
6. Calcular a y b para que la matriz
4 b
sea combinación lineal de
2 1 5 −3
las matrices 5 3 y 2 −1 .
4. Espacios vectoriales. Bases 105
−11 a 2 1
Solución: Si la matriz fuese combinación lineal de las matrices
4 b 5 3
5 −3
y , entonces existirían α, β ∈ R tales que
2 −1
2 1 5 −3 −11 a
α +β = .
5 3 2 −1 4 b
Esto es,
2α +5β = −11
α −3β = a
5α +2β = 4
3α −β = b
donde rg(A) = 3, que coincide con el número de incógnitas. Por tanto, la única
solución es la trivial. Esto es, α = β = γ = 0, y de ahí podemos armar que el
sistema S es libre.
b) Consideremos la combinación lineal:
α(1, 0, 0, −1) + β(2, 1, 1, −2) + γ(0, 1, 1, 0) + δ(1, −1, −1, −1) = (0, 0, 0, 0).
106 Más de 160 problemas resueltos de Álgebra Lineal
donde rg(A) = 3 que es menor que el número de incógnitas. Por tanto se tiene un
sistema compatible indeterminado y de ahí que el sistema sea ligado.
λ1 e1 + λ2 (e1 + e2 ) + · · · + λn (e1 + e2 + · · · + en ) = 0.
Esto es,
2α + 4λ = 0
−α − 3β + λ = 0
4α + 4β + 7λ = 0
5β − 5λ = 0.
Como rg(A) = 3, que coincide con el número de incógnitas, se tiene que el sistema
tiene únicamente la solución trivial. Esto es, α = β = λ = 0. Por tanto, el conjunto
es linealmente independiente.
Solución:
se tiene el sistema,
2α − β − λ + 6γ = 0
−2α − 2β − 3λ + γ = 0
5α + 3β + 3γ = 0
6α + 2β + 5λ = 0.
√ √
Se tiene que |A| = 0 si, y sólo si, a = −1, a = 12 1 − 17 ó a = 1
17 .
2 1+
Estudiemos los diferentes casos que pueden presentarse:
Si a = −1, entonces rg(A) = 2, que es menor que el número de incógnitas. Por
tanto el sistema es compatible indeterminado. Así, el sistema es linealmente
dependiente.
√
Si a = 12 1 − 17 , entonces rg(A) = 2, que es menor que el número de
Esto es:
x = 2α +β +γ
y = α +2β +γ
z = α +β +2γ
Como rg(A) = 3, que coincide con el número de incógnitas, el sistema tiene solución
única y, por tanto, {(2, 1, 1), (1, 2, 1), (1, 1, 2)} genera R3 .
c) Para que el conjunto dado genere el espacio vectorial R2 [x], dado ax2 + bx + c ∈
R2 [x], deben existir α, β, γ ∈ R tales que
Se tiene que rg(A) = 3, que coincide con el número de incógnitas. Por tanto el
sistema tiene solución única, siendo ésta:
a+b −a − b + 3c
α = a, β= , γ= ,
3 6
Solución: Para ver si el conjunto B = {p1 (x), p2 (x), p3 (x), p4 (x)} forma una base
del espacio vectorial de los polinomios de grado menor o igual que 3, al vericarse que
card(B) = 4 = dim(R3 [x]), bastará ver que son linealmente independientes. Para ello
consideramos la combinación lineal:
α1 p1 (x) + α2 p2 (x) + α3 p3 (x) + α4 p4 (x) = 0.
Esto es,
α1 + α2 (x + 1)2 + α3 x + α4 x3 = 0.
Solución: Recordamos un resultado que nos será útil en este ejercicio. Este resultado,
que puede consultarse en [9], nos presenta un método para construir una base de un
espacio vectorial,
es libre, y de ahí podemos armar que es una base de R4 con las condiciones
pedidas.
5α = 1
3α = 0
que es incompatible.
Añadimos dos elementos de la base canónica y veremos si estos vectores junto con
los ya existentes forman un sistema libre. Por tanto se trata de probar que B =
{(5, 3, −1, 2), (1, 0, 3, 2), (1, 0, 0, 0), (0, 1, 0, 0)} es libre. Para ello consideramos la ma-
triz A cuyas las con los vectores dados. Esto es,
5 3 −1 2
1 0 3 2
A=
.
1 0 0 0
0 1 0 0
Como rg(A) = 4, se tiene que el conjunto B forma un sistema libre. Al vericarse que
card(B) = 4 = dim(R4 ), se tiene que B es base.
4. Espacios vectoriales. Bases 113
Solución:
a) Bastará ver que es linealmente independiente ya que el cardinal del conjunto coin-
cide con la dimensión de R3 . Para ello construimos el sistema de ecuaciones corres-
pondiente a la combinación lineal α(2, 1, −3) + β(3, 2, −5) + γ(1, −1, 1) = (0, 0, 0).
Esto es:
2α + 3β + γ = 0
α + 2β − γ = 0
−3α − 5β + γ = 0.
b) Las coordenadas de (6, 2, −7) serían los escalares α, β y γ que permiten escribir
(6, 2, −7) como combinación lineal de ellos. Esto es:
cuya solución viene dada por α = β = γ = 1. Por tanto las coordenadas del vector
(6, 2, −7) en la base anterior (1, 1, 1).
obtenemos el sistema:
α1 +2α2 +α3 = 0
2α1 +3α2 +3α3 = 0
−α1 −α3 = 0
−2α1 −α2 +α4 = 0,
Las coordenadas del vector serán los escalares de la anterior combinación. Ahora
simplemente tenemos que resolver el siguiente sistema:
7 = λ1 +2λ2 +λ3
14 = 2λ1 +3λ2 +3λ3
−1 = −λ1 −λ3
2 = −2λ1 −λ2 +λ4
Solución:
a) En primer lugar vamos a ver qué vectores son linealmente independientes. Para
ello, consideremos la matriz cuyas las son los vectores dados.
1 0 −1 2
0 1 1 2
A=
1 1 0 4 .
1 1 −1 6
1 −1 −2 −4
4. Espacios vectoriales. Bases 115
b) Escribamos el vector (1, 1, 0, 4) como combinación lineal de los vectores que forman
la base. Esto es,
(1, 1, 0, 4) = α(1, 0, −1, 2) + β(0, 1, 1, 2) + γ(1, 1, −1, 6) + λ(1, −1, −2, −4).
Resolviendo el sistema resultante, obtenemos que las coordenados del vector (1, 1, 0, 4)
son: α = β = 1, γ = λ = 0. Es decir, (1, 1, 0, 4)B = (1, 1, 0, 0).
Solución: Probemos primero que A es una base. Para ello, como A es un conjunto
que tiene 4 elementos y este número coincide con la dimensión de R3 [x], será suciente
ver que A es un conjunto libre. Para ello, consideraremos la combinación lineal,
α · 1 + β · x + γ · x2 + δ · 2x3 = 0,
coordenadas de cada vector respecto de la otra base. Empecemos por 1. Para ello,
planteamos la ecuación:
1 = 2α + β(x − 2) + γ(x2 + 1) + δx3 .
Obteniendo el sistema:
1 = 2α − 2β + γ
0 = β
0 = γ
0 = δ
a) Demostrar que los polinomios p(x), p0 (x), p00 (x), p000 (x), forman una base
de R3 [x].
b) Si p(x) = 5x3 + 3x2 − 2x + 1, hallar las coordenadas del vector q(x) =
15x3 − 21x2 − 18x + 37 en la base B = {p(x), p0 (x), p00 (x), p000 (x)}.
Solución:
a) Veamos que B = {ax3 + bx2 + cx + d, 3ax2 + 2bx + c, 6ax + 2b, 6a} es libre. Para
ello consideramos la combinación lineal:
α(ax3 + bx2 + cx + d) + β(3ax2 + 2bx + c) + λ(6ax + 2b) + γ6a = 0.
15x3 − 21x2 − 18x + 37 = αp(x) + βp0 (x) + λp00 (x) + γp000 (x).
Solución:
1 −1 1
0 1 1 = 1 6= 0.
0 0 1
118 Más de 160 problemas resueltos de Álgebra Lineal
Es decir, la matriz cuyas las son las coordenadas de los vectores considerados,
tiene rango 3, que coincide con la dimensión de V . Así pues, la nueva base es
B 0 = {x, y, z} = {e1 − e2 + e3 , e2 + e3 , e3 }.
Por tanto,
e1 − e2 + e3 = y1 + y2 + y3 .
Como las coordenadas de un vector respecto de una base son únicas, se tiene que
1 = 1+a
−1 = b
1 = 1+c
Solución:
De aquí:
α1 + iα2 = β1 + iβ2 = 0.
Por tanto el conjunto {(1, 0), (0, 1)} es una base de C2 como C-espacio vectorial.
Por tanto:
α1 + α2 i = 2 + 3i, β1 + β2 i = −4 + 5i.
Así las coordenadas del vector (2+3i, −4+5i), respecto de la base {(1, 0), (0, 1)}
son (2 + 3i, −4 + 5i).
b.1) Veamos que el conjunto {(1, 1), (1, i), (i, 1), (i, −i)} es linealmente independiente
en C2 , considerado como R-espacio vectorial. Para ello consideramos la combi-
nación lineal α(1, 1) + β(1, i) + λ(i, 1) + γ(i, −i) = (0, 0), con α, β, λ, γ ∈ R.
Entonces:
α + β + λi + γi = 0
α + βi + λ − γi = 0.
120 Más de 160 problemas resueltos de Álgebra Lineal
Es decir:
α + β + λi + γi = a + bi
α + βi + λ − γi = c + di.
b.2) Razonando del mismo modo que en el apartado anterior, obtenemos el sistema,
α + β = 2
λ + γ = 3
α + λ = −4
β − γ = 5
Solución: En este ejercicio y en los que siguen nos será de bastante utilidad recor-
dar las ecuaciones de cambio de base en un espacio vectorial. Este resultado puede
consultarse en [9].
obteniendo que las coordenadas del vector (7, 4, 9) respecto de la base B 0 son
(4, 6, 11).
Solución: Sea B 0 una base donde v3 = ae1 +be2 +ce3 . Entonces la ecuación matricial
de cambio de base entre B 0 y B es
2 0 a
XB = −1 −1 b XB 0 .
0 0 c
Esto es:
1 1
a=− , b = 2, c= .
5 5
Así, la base pedida es B 0 = {v1 , v2 , v3 }, siendo
1 1
v1 = 2e1 − e2 , v2 = −e2 y v3 = − e1 + 2e2 + e3 .
5 5
Puesto que B = {u, v, w} es base, necesariamente se tiene que los vectores u, v , w son
linealmente independientes. Por tanto:
α −β +6γ = 0
α −β +γ = 0
α +5β −γ = 0.
Se pide:
a) Hallar la expresión de cambio de base entre B1 y B2 .
b) Determinar las coordenadas respecto de B1 del vector x cuyas coor-
denadas respecto de B2 son (2, 1, 0, −1).
124 Más de 160 problemas resueltos de Álgebra Lineal
Solución:
Así,
0 −1 −2 −1 2 1 3 0
−1 −1 0 −1 −2 1 0 −2
X = X .
0 −1 −1 B1 0 1 1 −1 B2
0
3 0 2 1 1 0 −1 1
Por tanto:
−1
0 −1 −2 −1 2 1 3 0
−1 −1 0 −1 −2 1 0 −2
XB1 = X .
1 −1 B2
0 0 −1 −1 0 1
3 0 2 1 1 0 −1 1
Esto es:
6/7 2/7 3/7 2/7
−3/7 −1/7 −5/7 −1/7
XB1 = X .
−11/7 1/7 −9/7 −6/7 B2
11/7 −8/7 2/7 13/7
Así obtenemos
12 6 15 1
XB1 = ,− ,− , .
7 7 7 7
Solución:
Solución: En este ejercicio consideraremos R4 como R-espacio vectorial con las ope-
raciones usuales de suma y producto por un número real.
127
128 Más de 160 problemas resueltos de Álgebra Lineal
= (3(αx2 + βy2 ), αx2 + βy2 , (αx4 + βy4 ) + (αx2 + βy2 ), αx4 + βy4 ),
ya que 1 · 1 6= 0.
1
c) Claramente C no es subespacio vectorial, puesto que tomando α = ∈ R y
2
1 1 1 1 1
x = (1, 1, 1, 1) ∈ C se tiene que αx = , , , / C ya que
∈ / Z.
∈
2 2 2 2 2
f) Supongamos que (3x, 2x, x, x), (3y, 2y, y, y) ∈ F y que α, β ∈ R. Tenemos que:
α(3x, 2x, x, x) + β(3y, 2y, y, y) = (3(αx + βy), 2(αx + βy), αx + βy, αx + βy) ∈ F.
Ahora,
n
X n
X n
X n
X n
X
(αxi + βyi ) = αxi + βyi = α xi + β yi = 0,
i=1 i=1 i=1 i=1 i=1
n n
ya que, por hipótesis, se verica que yi . De ahí que αx + βy ∈ F .
X X
xi = 0 =
i=1 i=1
Por tanto F es subespacio vectorial.
c) Sean A, B ∈ H3 y α, β ∈ R, entonces
a b d e αa αb βd βe
αA + βB = α +β = + =
−b a −e d −αb αa −βe βd
αa + βd αb + βe
= ∈ H3 .
−(αb + βe) αa + βd
d) Sean A, B ∈ H4 y α, β ∈ R, entonces
0 b 0 d 0 αb 0 βd
αA + βB = α +β = + =
a 0 c 0 αa 0 βc 0
0 αb + βd
= ∈ H4 .
αa + βc 0
Solución:
Así,
(αf +βg)(x) = αf (x)+βg(x) = αan xn +...+αa1 x+αa0 +βbn xn +...+βb1 x+βb0 =
= (αan + βbn )xn + (αan−1 + βbn−1 )xn−1 + ... + (αa1 + βb1 )x + (αa0 + βb0 ).
=⇒ 1) Como x, y ∈ L({x, y}), se tiene que x, y ∈ L({x, y}) = L({u, v}). De igual
forma obtenemos 2), observando que u, v ∈ L({u, v}).
⇐= Como x, y ∈ L({u, v}) entonces L({x, y}) ⊆ L({u, v}). De igual forma si u, v ∈
L({x, y}) entonces L({u, v}) ⊆ L({x, y}). Por lo que L({u, v}) = L({x, y}).
5. Subespacios vectoriales 133
Una vez demostrado el hecho anterior, veamos que x, y ∈ L({u, v}) y que u, v ∈
L({x, y}), y de ahí tendremos el resultado deseado.
Solución:
V3 = {(λ + µ, λ + γ, γ + δ, λ + δ) : λ, µ, γ, δ ∈ R}.
Solución:
es no nulo. Así la matriz ampliada tiene rango 3 y por tanto no coincidirá nunca
con el rango de la matriz de coecientes. Luego el sistema no tiene solución y en
consecuencia (1, 0, 1, −2) ∈
/ V2 .
Solución:
x1 = λ +2µ
x2 = −λ +3µ
a) x3 = 2λ +µ
x4 = λ −µ
x5 = µ
1 2
−1 3
es no nulo, el rango de la matriz A será 2 si, y sólo si, los menores de orden 3,
que obtenemos orlando el menor elegido con la tercera columna y cada una de las
restantes las, son nulos. Por tanto,
1 2 x1 1 2 x1 1 2 x1
−1 3 x2 = 0, −1 3 x2 = 0, −1 3 x2 = 0.
2 1 x3 1 −1 x4 0 1 x5
x1 = 2λ +µ
x2 = −3λ +2µ
b)
x3 = λ −µ
x4 = λ +2µ
5. Subespacios vectoriales 137
es distinto de cero, orlando con la tercera columna y cada una de las restantes
las, ha de vericarse
2 1 x1 2 1 x1
−3 2 x2 = 0, −3 2 x2 = 0.
1 −1 x3 1 2 x4
Por tanto, las ecuaciones implícitas del subespacio vienen dadas por:
x1 +3x2 +7x3 = 0
8x1 +3x2 −7x4 = 0.
x1 = λ +µ
c) x2 = λ
x3 = −λ −2µ
Solución:
x1 = λ +µ
x2 = λ −µ +δ
a) x3 = λ +µ +δ
x4 = λ −δ
x5 = 2λ +δ
es no nulo, el rango de la matriz A será 3 si, y sólo si, todos los menores de orden
4, que obtenemos orlando el menor elegido con la cuarta columna y cada una de
las demás las, son nulos. Por tanto,
1 1 0 x1 1 1 0 x1
1 −1 1 x2 1 −1 1 x2
= 0, = 0.
1 1 1 x3 1 1 1 x3
1 0 −1 x4 2 0 1 x5
es distinto de cero, orlando con la cuarta columna y cada una de las restantes las,
ha de vericarse que:
1 1 2 x1 1 1 2 x1
1 2 3 x2 1 2 3 x2
= 0, = 0.
1 1 0 x3 1 1 0 x3
0 0 0 x4 2 −1 0 x5
Solución:
F0
Si p(x) = ax3 + bx2 + cx + d ∈ F0 , tenemos que p(0) = 0. Luego se tiene que
p(0) = d = 0. Como F0 es un subespacio de R3 [x], dim(F0 ) ≤ dim(R3 [x]) = 4.
140 Más de 160 problemas resueltos de Álgebra Lineal
Calculemos una base de F0 . Si p(x) ∈ F0 entonces p(x) = ax3 + bx2 + cx. Luego
un sistema generador de F0 es {x3 , x2 , x}, que al ser libre, constituye una base
de F0 . Por tanto, dim(F0 ) = 3 y tiene una única ecuación implícita que es d = 0.
F1
Si p(x) = ax3 + bx2 + cx + d ∈ F1 , tenemos que p(1) = 0. Luego se tiene que
p(1) = a + b + c + d = 0. Como F1 es un subespacio de R3 [x], dim(F1 ) ≤
dim(R3 [x]) = 4. Calculemos una base de F1 . Si p(x) ∈ F1 entonces p(x) =
ax3 + bx2 + cx − (a + b + c) = a(x3 − 1) + b(x2 − 1) + c(x − 1). Luego un
sistema generador de F1 es {x3 − 1, x2 − 1, x − 1}. Como además es libre, es
base de F1 . Por tanto, dim(F1 ) = 3 y tiene una única ecuación implícita que es
a + b + c + d = 0.
F−1
Si p(x) = ax3 + bx2 + cx + d ∈ F−1 , tenemos que p(−1) = 0. Luego se tiene que
p(−1) = −a + b − c + d = 0. Como F−1 es un subespacio de R3 [x], dim(F−1 ) ≤
dim(R3 [x]) = 4. Calculemos una base de F−1 . Si p(x) ∈ F−1 entonces p(x) =
ax3 + bx2 + cx + (a − b + c) = a(x3 + 1) + b(x2 − 1) + c(x + 1). Luego un
sistema generador de F−1 es {x3 + 1, x2 − 1, x + 1}. Como además es libre, es
base de F−1 . Luego dim(F−1 ) = 3 y existe una única ecuación implícita que es
−a + b − c + d = 0.
Solución:
a) Sea B = {p(x), p0 (x), p00 (x), p000 (x)} = {x3 , 3x2 , 6x, 6}. Como card(B) = 4 =
dim(R3 [x]), basta ver que B es linealmente independiente. Veámoslo.
5. Subespacios vectoriales 141
Esto es,
(5α + 2β + γ)x2 + (3α + 3γ)x + (α − β + 3γ) = 0.
vemos que rg(A) = 2, que es menor que el número de incógnitas, es decir, el sis-
tema tiene solución distinta de la trivial y por tanto, los vectores p1 (x), p2 (x) y
p3 (x) son linealmente dependientes.
1 + x = α(1 + x2 ) + β(1 − x2 ).
Esto es,
1 + x = (α + β) + (α − β)x2 .
142 Más de 160 problemas resueltos de Álgebra Lineal
1 2
14. Sea A = 3 m
una matriz en el espacio vectorial de matrices cua-
dradas de orden 2.
Solución:
Solución: Para probar que F es un subespacio vectorial de R3 [x], tomemos p(x), q(x) ∈
F y α, β ∈ R. Entonces
Z 1 Z 1 Z 1 Z 1 Z 1
(αp(x)+βq(x))dx = αp(x)dx+ βq(x)dx = α p(x)dx+β q(x)dx = 0,
0 0 0 0 0
De aquí, si α, β, γ ∈ R, obtenemos:
a = 4α
b = 3β
c = 2γ
d = −α −β −γ
146 Más de 160 problemas resueltos de Álgebra Lineal
Por tanto, todo polinomio p(x) = ax3 + bx2 + cx + d ∈ F , se puede expresar como
ax3 +bx2 +cx+d = 4αx3 +3βx2 +2γx+(−α−β−γ) = α(4x3 −1)+β(3x2 −1)+γ(2x−1).
Esto es,
2ax3 + 2cx = 0.
Así, se tiene que p(x) = bx2 + d. Por tanto una base vendrá dada por {1, x2 }. Por
último, tenemos que las ecuaciones implícitas de P son a = c = 0.
Solución:
Base y dimensión de F .
Llamando a = α, b = β , c = γ , tenemos que d = β − γ . Así, dado (a, b, c, d) ∈ F ,
se tiene que:
(a, b, c, d) = (α, β, γ, β − γ) = α(1, 0, 0, 0) + β(0, 1, 0, 1) + γ(0, 0, 1, −1).
5. Subespacios vectoriales 147
Como estos vectores son libres, ya que el rango de la matriz cuyas las son
estos vectores es 3, se tiene que una base de F es BF = {(1, 0, 0, 0), (0, 1, 0, 1),
(0, 0, 1, −1)}. Así dim(F ) = 3.
Base y dimensión de G.
Llamando a = α, d = β , se tiene que b = −α y c = 2β . Así, dado (a, b, c, d) ∈ G,
entonces:
(a, b, c, d) = (α, −α, 2β, β) = α(1, −1, 0, 0) + β(0, 0, 2, 1).
Como estos vectores son libres, una base de G es BG = {(1, −1, 0, 0), (0, 0, 2, 1)}.
Luego dim(G) = 2.
Base y dimensión de F ∩ G.
Primero veamos qué conjunto es F ∩ G. Como ha de cumplir las condiciones de
los dos conjuntos, se tendrá que
F ∩ G = {(a, b, c, d) : b − c − d = 0, a + b = 0, c = 2d}.
Esto es,
a +b = 0
b −c −d = 0
c −2d = 0.
Base y dimensión de F + G.
Sabemos que un sistema generador de F + G es BF ∪ BG , es decir:
{(1, 0, 0, 0), (0, 1, 0, 1), (0, 0, 1, −1), (1, −1, 0, 0), (0, 0, 2, 1)}.
donde
a b a 0 0 b 1 0 0 1
= + =a +b .
a−b a+b a a −b b 1 1 −1 1
1 0 0 1
Luego B1 = , es sistema generador de V1 y, por ser libre,
1 1 −1 1
es también base.
a b
siendo α, β ∈ R. Esto es, si ∈ V2 , entonces
c d
a b α −3α α −3α 0 0
= = + =
c d β 2α − β 0 2α β −β
1 −3 0 0
α +β .
0 2 1 −1
1 −3 0 0
Luego una base de V2 es B2 = , .
0 2 1 −1
−a +b +c = 0
−a −b +d = 0
a +b +c +d = 0
2a −c −d = 0.
La solución de este sistema nos dará las ecuaciones paramétricas de V1 ∩ V2 . Esto es,
dado α ∈ R, se tiene que:
a = −α
b = 3α
c = −4α
d = 2α
−1 3
luego una base de V1 ∩ V2 viene dada por B = .
−4 2
Si además es libre, sería la base buscada. Vamos a verlo: sean α, β, γ, δ ∈ R tales que
1 0 0 1 1 −1 0 0 0 0
α +β +γ +δ = .
1 1 −1 1 0 0 1 −1 0 0
tiene rango 4. Así la única solución del sistema es la solución trivial y por tanto B 0
es libre, y de ahí podemos armar que B 0 es una base de V1 + V2 .
5. Subespacios vectoriales 151
Solución:
Recordamos ahora dos resultados que nos serán de bastante utilidad en los ejer-
cicios que vienen a continuación. Estos resultados pueden consultarse en [9]:
152 Más de 160 problemas resueltos de Álgebra Lineal
Fórmula de Grassmann
Dados dos subespacios vectoriales F y G de un espacio vectorial V de di-
mensión nita, se verica que:
dim(F ) + dim(G) = dim(F ∩ G) + dim(F + G).
21. Sea R3 [x] el espacio vectorial de los polinomios de grado menor o igual
que 3 con coecientes reales y consideremos la base B = {x3 , x2 , x, 1}. Dados
los subespacios
2 2 2 x2 + x3 = 0
U = L(x + 2x, −x + x, x + x), W ≡ ,
2x2 − x3 = 0
x1 = 0
x2 = −β
V ≡
x3 = 0
x4 = α+β
Calcula:
a) U ∩ W , U + W.
b) ¾Son U y W suplementarios?
c) Una base de W ∩ V.
d) Ecuaciones cartesianas de U + V.
Solución:
a) U ∩W
Calculamos una base de U a partir del sistema generador {x2 + 2x, −x2 +
x, x2 + x}. Comprobemos si es libre. Como vimos en ejercicios anteriores, esto
5. Subespacios vectoriales 153
U +W
Sabemos que un sistema generador de U +W viene dado por BU ∪BW , siendo
BU y BW bases de U y W respectivamente. Ya tenemos una base de U ,
calculemos ahora una base de W . Vemos que las ecuaciones cartesianas de W
son equivalentes a x2 = x3 = 0. Esto es, si p(x) = x1 x3 +x2 x2 +x3 x+x4 ∈ W ,
entonces p(x) = x1 x3 +x4 . Luego una base de W es BW = {x3 , 1}. Así U +W
está generado por el conjunto {x2 + 2x, x, x3 , 1} que, como además es libre,
constituye una base, luego BU +W = {x3 , x2 + x, x, 1} es una base de U + W .
b) Tenemos que ver que U ∩ W = {0} y que U + W = R3 [x]. Veamos primero que
su intersección es el vector nulo. Sea p(x) = x1 x3 + x2 x2 + x3 x + x4 ∈ U ∩ W .
Si p(x) ∈ U , entonces x1 = x4 = 0 y si p(x) ∈ W , entonces x2 = x3 = 0, luego
p(x) = 0. Así, U ∩ W = {0}. Veamos que su suma es el total. Por la fórmula de
Grassmann, se tiene que
dim(U ) + dim(W ) = dim(U ∩ W ) + dim(U + W ),
de donde
dim(U + W ) = 4 = dim(R3 [x]).
Por tanto, una base de U + V viene dada por BU +V = {x2 + 2x, x, 1}. Vamos a
calcular ahora las ecuaciones cartesianas de U + V . Si p(x) = x1 x3 + x2 x2 + x3 x +
x4 ∈ U + V, entonces p(x) = α(x2 + 2x) + βx + λ. Esto es,
x1 = 0
x2 = α
x3
= 2α + β
x4 = λ
En primer lugar calcularemos una base de F . Para ello, bastará ver si el conjunto
S es libre. Esto puede comprobarse calculando el rango de la matriz A, cuyas las
son los vectores dados. Esto es,
0 2 −1 1 0
0 3 0 0 1 .
0 1 1 −2 1
Como rg(A) = 3, se tiene que el conjunto S es libre y de ahí, que constituya una
base de F . En consecuencia, dim(F ) = 3. Por tanto, la dimensión de un subespacio
5. Subespacios vectoriales 155
suplementario tiene que ser 2. Como una base de R5 es S ∪{(1, 0, 0, 0, 0), (0, 0, 0, 0, 1)},
una base del subespacio suplementario vendrá dada por {(1, 0, 0, 0, 0), (0, 0, 0, 0, 1)}.
Solución:
a) De las ecuaciones paramétricas de U , podemos ver que {(1, 0, 1), (0, 1, 1), (1, 1, 2)}
es un sistema generador de U . Pero estos vectores no son linealmente inde-
pendientes puesto que (1, 0, 1) + (0, 1, 1) = (1, 1, 2). Luego podemos decir que
{(1, 0, 1), (0, 1, 1)} es base de U ya que ambos vectores no son proporcionales y
por tanto linealmente independientes.
luego {(1, 1, 0), (−2, 0, 1)} es sistema generador y, al ser libre, es base de W .
Luego
x = α + λ − 2δ
y = β + λ
z = α + β + δ
156 Más de 160 problemas resueltos de Álgebra Lineal
Por tanto {(1, 0, 1), (0, 1, 1), (1, 1, 0), (−2, 0, 1)} es sistema generador de U + W .
Podemos excluir (−2, 0, 1) puesto que (−2, 0, 1) = − 12 (1, 0, 1) + 32 (0, 1, 1) −
2 (1, 1, 0). Además {(1, 0, 1), (0, 1, 1), (1, 1, 0)} es un sistema libre, ya que el rango
3
de la matriz cuyas las son esos vectores es 3. En consecuencia {(1, 0, 1), (0, 1, 1),
(1, 1, 0)} es una base de U + V .
Esto es,
α + γ = 2
β + γ = 3
α + β = 5.
Por tanto, f (x) + f (−x) = 2p(x). Luego p(x) = 12 [f (x) + f (−x)] y f (x) − f (−x) =
2g(x). Es decir, g(x) = 21 [f (x) − f (−x)].
Solución: Observamos en primer lugar que (1, 0, 1, 1) + (1, 2, 1, 1) = (2, 2, 2, 2). Por
tanto podemos eliminar este último vector del conjunto, obteniendo la base de V , esto
es,
BV = {(1, 0, 1, 1), (1, 2, 1, 1)}.
A continuación, vamos a añadir vectores a esta base hasta obtener una base de R4 .
Puede comprobarse que el conjunto
constituye una base de R4 . Así una base de R4 /V es {[(0, 0, 1, 0)], [(0, 0, 0, 1)]}.
Solución: Primero vamos a calcular una base de U . Observamos que dados los
vectores (2, 0, −1), (1, 2, 0), (0, 4, 1), podemos construir la matriz cuyas las sean los
vectores dados:
2 0 −1
1 2 0 .
0 4 1
Así, podemos considerar como sistema generador a {(1, 2, 0), (0, 4, 1)} que, al ser libre,
será una base. Esto es, tomamos como base de U el conjunto BU = {(1, 2, 0), (0, 4, 1)}.
Así un sistema generador de W será {(0, 1, −1)}, que al ser libre es base. Por
tanto, una base de W será BW = {(0, 1, −1)}. Ahora, podemos completar BW con
los vectores (1, 0, 0), (0, 0, 1) para obtener una base de R3 y, así,
{[(1, 0, 0)], [(0, 0, 1)]}
((αx1 + βy1 ) − (αx2 + βy2 ), (αx1 + βy1 ) + (αx2 + βy2 ), αx1 + βy1 ) =
(α(x1 − x2 ) + β(y1 − y2 ), α(x1 + x2 ) + β(y1 + y2 ), αx1 + βy1 ) =
(α(x1 − x2 ), α(x1 + x2 ), αx1 ) + (β(y1 − y2 ), β(y1 + y2 ), βy1 ) =
α(x1 − x2 , x1 + x2 , x1 ) + β(y1 − y2 , y1 + y2 , y1 ) = αf (x1 , x2 ) + βf (y1 , y2 ).
Mostraremos algunas propiedades que verica una aplicación lineal, dependiendo del
núcleo y de la imagen,
Quedando pues
−c a−d 0 0
= .
0 c 0 0
a = 0
b = 0
c+d = 0,
Solución:
⇐= Demostraremos que ker(f ) = {0}. Supongamos que exista x ∈ ker(f ) tal que
x 6= 0. Así, el conjunto unitario {x} es linealmente independiente. Entonces, por
hipótesis, se tiene que {f (x)} es libre. Pero x ∈ ker(f ), luego {f (x)} = {0} es un
conjunto libre, y esto es una contradicción. Por tanto ker(f ) = {0} y de ahí que f
sea una aplicación inyectiva.
166 Más de 160 problemas resueltos de Álgebra Lineal
Solución:
c) Una base de R3 es {(1, 0, 0), (0, 1, 0), (0, 0, 1)}. Por tanto un sistema generador de
la imagen de f es
{f ((1, 0, 0)), f ((0, 1, 0)), f ((0, 0, 1))} = {(1, 0, 0), (0, 0, 0), (0, 0, 0)}.
Solución:
c) Tal y como vimos en el primer apartado, las imágenes de los elementos de la base
canónica de R2 son
f (1, 0) = (1, 1, 1), f (0, 1) = (−1, 1, 0),
luego B = {(1, 1, 1), (−1, 1, 0)} es sistema generador del conjunto Imagen y, puesto
que B es además linealmente independiente, armamos que B es base de Im(f ). El
rango de la aplicación es el rango de cualquier matriz asociada. Esto es, el rango
de f es el rg(A), que es en este caso 2.
En [9] puede consultarse el siguiente resultado que nos será de bastante utilidad en
los ejercicios que siguen:
Solución:
f inyectiva =⇒ f sobreyectiva.
Supongamos que f es inyectiva. Entonces ker(f ) = {0}. Esto es, dim(ker(f )) =
0. Utilizando la fórmula de las dimensiones se tiene dim(V ) = dim(ker(f )) +
dim(Im(f )) = 0 + dim(Im(f )). Como dim(V ) = dim(V 0 ), podemos armar que
dim(V 0 ) = dim(Im(f )) y de ahí que Im(f ) = V 0 . Por tanto f es sobreyectiva.
f sobreyectiva =⇒ f isomorsmo.
Supongamos que f es sobreyectiva. Entonces dim(Im(f )) = dim(V 0 ). Como
dim(V ) = dim(V 0 ), tenemos que dim(Im(f )) = dim(V ). Por tanto, utilizando
la fórmula de las dimensiones, podemos armar que dim(ker(f )) = 0, y de ahí
que ker(f ) = {0}. Así f es inyectiva, y al ser sobreyectiva, se tiene el resultado.
f isomorsmo =⇒ rg(f ) = n.
Si f es un isomorsmo, entonces f es sobreyectiva. Por tanto, rg(f ) = dim(Im(f ))
= dim(V 0 ) = n.
rg(f ) = n =⇒ f inyectiva.
Supongamos que rg(f ) = n. Entonces se tiene que dim(Im(f )) = n = dim(V ).
Utilizando la fórmula de las dimensiones se tiene que dim(ker(f ))+dim(Im(f )) =
dim(V ) = n. Por tanto, dim(ker(f )) = n − n = 0. Así f es inyectiva.
6. Aplicaciones lineales 169
Solución: Consideremos
1 0 0 1 0 0 0 0
B= , , ,
0 0 0 0 1 0 0 1
a b a b
Calculemos el núcleo de f . Si A = ∈ ker(f ), entonces f =
c d c d
0
0 . Esto es:
0
a+b+c+d 0
2a + 3b + 4c + 6d = 0 .
a + 3b + 5c + 9d 0
Esto es,
1 1 1 1
2 , 3 , 4 , 6
1 3 5 9
es sistema generador de Im(f ). Puede comprobarse que estos vectores son linealmente
dependientes, y por tanto, una base de Im(f ) viene dada por:
1 1
2 , 3 .
1 3
Solución:
donde f (ui )B 0 denota a las coordenadas del vector f (ui ) respecto de la base B 0 .
c) Sabemos que el rango de la aplicación f coincide con dim(Im(f )). Así, por el
apartado b), y utilizando la fórmula de las dimensiones, se tiene que el rango de
la aplicación es 3.
Solución:
Sea el vector u = αu1 + βu2 + γu3 ∈ V3 tal que f (u) = −3v1 + 2v3 . Entonces se
tiene que
−3 2 0 −1
α
0 1
= 1 0
β ,
2 −1 0 1
γ
0 1 −1 2
11. Sea f una aplicación lineal denida entre los espacios vectoriales V2 y
V4 tal que, respecto de las bases B = {u1 , u2 } y B1 = {v1 , v2 , v3 , v4 } tiene por
matriz asociada
1 −1
1 0
A= .
1 −1
1 0
Se eligen en V2 y V4 las nuevas bases B2 = {e1 , e2 } y B3 = {w1 , w2 , w3 , w4 }
donde:
w1 = v1 + v2 + v3 + v4
e1 = u1 + 2u2 w2 = v1 − v2
y
e2 = −u1 + u2
w3 = + v3 − v4
w4 = v4
6. Aplicaciones lineales 173
a) B y B3 , b) B2 y B1 , c) B2 y B3 .
Solución: Vamos a recordar la relación entre las matrices asociadas a una misma
aplicación lineal pero respecto a bases diferentes. Este resultado puede consultarse en
[9].
Entonces la relación que existe entre las matrices A y C viene dada por
C = Q−1 AP.
−1
1 1 0 0 1 −1/2
1 −1 0 0
A 1 0 0 −1/2
a) C = = .
1 0 1 0 0 1 0 −1/2
1 0 −1 1 0 0
−1
1 0 0 0 −1 −2
0 1 0 0
A 1 −1 1 −1
b) C = =
−1
.
0 0 1 0 2 1 −2
0 0 0 1 1 −1
−1
1 1 0 0 0 −3/2
1 −1 0 0
A 1 −1 −1 −1/2
c) C = =
−1
.
1 0 1 0 2 1 −1/2
1 0 −1 1 0 0
174 Más de 160 problemas resueltos de Álgebra Lineal
respectivamente. Hallar:
Solución: Recordamos a continuación los siguientes resultados que nos serán de gran
utilidad y que pueden consultarse en [9]:
a) Puesto que f y g están denidas en el mismo espacio vectorial y que sus matrices
asociadas están dadas respecto de las mismas bases, podemos armar que:
M (f + g) = M (f ) + M (g) = A + B =
−1 2 9 0 1 −1 −1 3 8
3 2 1 + 3 1 0 = 6 3 1 .
1 1 −2 −2 1 3 −1 2 1
6. Aplicaciones lineales 175
c) Se tiene que:
−1 3 8 1 17
(f + g)(u) = 6 3 1 −2 = 3 .
−1 2 1 3 −2
e) Como una aplicación lineal entre dos espacios de la misma dimensión es inyectiva
si, y sólo si, es sobreyectiva, como se ve en el ejercicio 7, bastará con ver una de las
dos propiedades para armar que es biyectiva. Estudiemos la inyectividad. Para
ello calculamos el núcleo de la aplicación. Esto es, sea (x, y, z) ∈ ker(f ), entonces
−1 2 9 x 0
3 2 1 y = 0 .
1 1 −2 z 0
y1 −5/28 13/28 −4/7 x1
y2 = 1/4 −1/4 1 x2 .
y3 1/28 3/28 −2/7 x3
6. Aplicaciones lineales 177
Por tanto:
f (1, 1) = f (e1 + e2 ) = f (e1 ) + f (e2 ) = 5v1 + 2v2 + 3v3 ,
Esto es:
f (e1 ) +f (e2 ) = 5v1 + 2v2 + 3v3
2f (e1 ) +3f (e2 ) = 2v1 + 4v3
De donde:
f (e1 ) = 13v1 + 6v2 + 5v3 ,
Por tanto,
13 −8
M [f : B2 , B3 ] = 6 −4 .
5 −2
Por tanto,
Esto es:
−2g(v1 ) +4g(v2 ) +2g(v3 ) = w1 + w2 + w3 + w4
g(v1 ) −g(v3 ) = 2w1 − w2 + 3w3 + 4w4
−g(v1 ) +2g(v2 ) = w2 + w4
De donde,
5 3 7 7
g(v1 ) = w1 − w2 + w3 + w4
2 2 2 2
5 1 7 9
g(v2 ) = w1 − w2 + w3 + w4
4 4 4 4
1 1 1 1
g(v3 ) = w1 − w2 + w3 − w4
2 2 2 2
Por tanto,
5/2 5/4 1/2
−3/2 −1/4 −1/2
M [g : B3 , B4 ] =
7/2
.
7/4 1/2
7/2 9/4 −1/2
Se pide:
a) Matriz de la aplicación lineal g ◦ f : A −→ C.
b) Hallar (g ◦ f )−1 (1, −1, 1).
c) Imagen, mediante g ◦ f , del subespacio intersección de los subespacios
de A siguientes:
x1 = 2λ + β
V1 ≡ x2 = λ − β V2 ≡ x1 − x2 + 2x3 = 0.
x3 = −λ
Solución:
Se tiene que
1 1
1 0 0
M [g : BB , BC ] = −1 −1 , M [f : BA , BB ] = .
−1 1 2
2 0
Por tanto:
1 1 0 1 2
1 0 0
M [g ◦ f : BA , BC ] = −1 −1 = 0 −1 −2 .
−1 1 2
2 0 2 0 0
b) Sea (x1 , x2 , x3 ) ∈ (g ◦ f )−1 (1, −1, 1). Entonces, se tendrá que vericar:
(g ◦ f )(x1 , x2 , x3 ) = (1, −1, 1).
Esto es,
0 1 2 x1 1
0 −1 −2 x2 = −1 .
2 0 0 x3 1
1
cuyas soluciones son x1 = , x2 = 1 − 2α, x3 = α con α ∈ R.
2
1
Así, (g ◦ f ) (1, −1, 1) =
−1
, 1 − 2α, α : α∈R .
2
c) Dadas las ecuaciones paramétricas de V1 ,
x1 = 2λ +β
V1 ≡ x2 = λ −β
x3 = −λ
su ecuación implícita viene dada por, x1 +x2 +3x3 = 0. Así, si (x1 , x2 , x3 ) ∈ V1 ∩V2 ,
se tiene que:
x1 −x2 +2x3 = 0
x1 +x2 +3x3 = 0.
Hallar:
a) Matriz asociada a (f + g).
b) Base y dimensión de ker(f + g) e Im(f + g).
c) ¾Qué tipo de aplicación es (f + g)?
d) Calcular f −1 (4, 6).
Solución:
se tiene que
1 −1 0
M [f ] = .
1 0 −1
Solución:
luego rg(A) = n.
Solución:
n
Así, si x ∈ V1 , entonces x = λi vi . Por tanto,
X
i=1
n
! n n
X X X
f (x) = f λi v i = λi f (vi ) = λi vi0 .
i=1 i=1 i=1
i=1 i=1
se tiene que
n
X n
X
αi vi0 = βi vi0 .
i=1 i=1
n
Es decir, (αi − βi )vi0 = 0. Como {v10 , . . . , vn0 } es una base, en particular
X
i=1
es libre. Por tanto, αi − βi = 0. Es decir, αi = βi con 1 ≤ i ≤ n. Así x = y .
184 Más de 160 problemas resueltos de Álgebra Lineal
Así,
dim(ker(f )) + 1 = dim(V ).
De donde
dim(ker(f )) = dim(V ) − 1.
Solución: Sea
a11 a12 a13
a21 a22 a23
A=
a31
a32 a33
a41 a42 a43
Obtenemos así que una base del ker(f ) viene dada por {(1, 0, 1)}. Por tanto
a11 a12 a13 0
a21 1
a22 a23 0
0 =
.
a31 a32 a33 0
1
a41 a42 a43 0
Se tiene que Im(f ) = {(y1 , y2 , y3 , y4 ) : y1 −y2 = 0, y2 −y3 = 0}. Como (a1i , a2i , a3i , a4i ) ∈
Im(f ) con i = 1, 2, 3, entonces a1i = a2i = a3i . Además, al vericarse que dim(Im(f )) =
2, se tendrá que rg(A) = 2. Por tanto, podemos armar que la matriz asociada a la
aplicación lineal f es una matriz de la forma
α β −α
α β −α
A= , α, β, γ, δ ∈ R.
α β −α
γ δ −γ
donde rg(A) = 2.
Solución:
Por tanto una base de U es {(2, 1, −2, 0), (0, 0, −1, 1)}. A continuación, comple-
taremos la base hasta obtener una base de R4 . Observemos que si añadimos los
vectores (0, 1, 0, 0) y (0, 0, 0, 1) la siguiente matriz tiene rango 4:
2 1 −2 0
0 1 0 0
.
0 0 −1 1
0 0 0 1
Quedándonos con esos vectores, obtenemos que una base de W viene dada por
b) Sea
a11 a12 a13 a14
a21 a22 a23 a24
A=
a31 a32 a33 a34
a41 a42 a43 a44
y
a11 a12 a13 a14 0 0
a21 a22 a23 a24 0 0
=
0 .
a31 a32 a33 a34 −1
a41 a42 a43 a44 1 0
6. Aplicaciones lineales 187
Solución: Consideremos el conjunto B1 = {(0, 0, 1), (1, −1, 0), (1, 1, 0)} y veamos que
se trata de una base de R3 . Para ello, bastará ver que es linealmente independiente, ya
que su cardinal coincide con la dimensión de R3 . Calculamos, pues, el determinante
de la matriz cuyas las son los vectores dados. Esto es,
0 0 1
1 −1 0 = 2 6= 0.
1 1 0
Luego
0 0 1
1
P = − 12 0 .
2
1 1
2 2 0
Pasemos ahora a calcular una base del núcleo de la aplicación. Sea (x, y, z) ∈ ker(f ),
entonces
1 1
2 −2 0 x 0
1 1
−2 0 y = 0 .
2
0 0 1 z 0
Esto es,
x = y
z = 0.
Luego una base del ker(f ) es {(1, 1, 0)} y por tanto dim(ker(f )) = 1. Para hallar el
subespacio imagen, calculamos las imágenes de los elementos de la base canónica, es
decir,
1 1 1 1
f (1, 0, 0) = ,− ,0 , f (0, 1, 0) = − , ,0 , f (0, 0, 1) = (0, 0, 1).
2 2 2 2
En consecuencia, una base de Im(f ) es {(1, −1, 0), (0, 0, 1)} y, por tanto, dim(Im(f )) =
2. Así:
{(1, 1, 0), (1, −1, 0), (0, 0, 1)},
constituye un sistema generador de ker(f ) + Im(f ) y, puesto que dicho conjunto, que
es en realidad B1 , es libre, se trata de una base. Así, dim(ker(f ) + Im(f )) = 3. Luego,
utilizando la fórmula de las dimensiones, se tiene que dim(ker(f ) ∩ Im(f )) = 0, y, por
tanto, ker(f ) ⊕ Im(f ) = R3 .
6. Aplicaciones lineales 189
22. Sea R2 [x] el espacio vectorial de los polinomios de grado menor o igual
que 2. Consideramos la aplicación lineal D : R2 [x] −→ R2 [x] que asocia a
cada polinomio su derivada. Se pide:
a) Calcula la matriz asociada a D en la base B = {1, x, x2 }.
b) Calcular la matriz asociada a D en la base B 0 = {x − 3, x2 − 3, 7 − x − x2 }.
Solución:
a) Si p(x) = 1, se tiene que p0 (x) = 0. Así las coordenadas de p0 (1) con respecto
a B son (0, 0, 0). Esto es, p0 (x)B = (0, 0, 0).
Si p(x) = x, se tiene que p0 (x) = 1. Así p0 (x)B = (1, 0, 0).
b) Sea ahora p(x) = x − 3, entonces p0 (x) = 1. Ahora tendremos que calcular las
coordenadas de este vector respecto de la base B 0 . Para ello calcularemos los
escalares α, β y γ tales que:
Sea ahora p(x) = (7 − x − x2 ). Se tiene que p0 (x) = −2x − 1. Para calcular sus
coordenadas respecto de la base B 0 , al igual que antes, calcularemos los escalares
α, β y γ tales que:
Solución:
cuyas soluciones dependen del parámetro α. Por ello vamos a discutir el sistema
dependiendo de los valores de α. Se tiene que |A| = 1 − 2α2 . Por tanto:
Si α 6= √12 , α 6= − √12 , entonces |A| 6= 0. Así la única solución del sistema es
la trivial. Esto es, x = y = z = 0. Luego ker(fα ) = {(0, 0, 0)}.
Si α = √1 ,
2
puede comprobarse que una base del ker(f √1 ) es {(− √12 , 1, − √12 )}.
2
Si α = − √12 , puede comprobarse que una base del ker(f− √1 ) es {( √12 , 1, √12 )}.
2
Sabemos que
1 α 0
!
A = M [fα : B] = α 1 α .
0 α 1
Así, la expresión de cambio de base, vendrá dada por M [fα : B 0 ] = Q−1 AP , siendo
1 0 0
!
P =Q= −1 1 0 .
0 −1 1
192 Más de 160 problemas resueltos de Álgebra Lineal
Por tanto,
1−α α 0
!
0
M [fα : B ] = 0 1 α .
−α α α+1
Solución:
Para que ambas matrices sean iguales debe vericarse M = Q−1 M P . Esto es,
P M = M P . Resolviendo el sistema resultante de la anterior igualdad, nos queda
que M debe de ser de la forma:
e d i−e−d
M = d e i − e − d .
0 0 i
6. Aplicaciones lineales 193
Observamos que un sistema generador de Im(f ) viene dado por {(e, d, 0), (d, e, 0), (i−
e − d, i − e − d)}. Por tanto, dim(Im(f )) ≤ 3.
b) En el apartado anterior hemos visto que una base de Im(f ) es {(1, −1, 0)}. Halle-
mos una base para ker(f ). Sea (x, y, z) ∈ ker(f ), entonces se tiene
1
2 − 12 0 x 0
− 12 1
0 y = 0 .
2
0 0 0 z 0
194 Más de 160 problemas resueltos de Álgebra Lineal
Esto es,
x = α
y = α
z = β,
donde α, β ∈ R. Por consiguiente, una base de ker(f ) viene dada por {(1, 1, 0), (0, 0, 1)}.
En primer lugar necesitamos calcular una base de R3 / ker(f ). Para ello calcularemos
una base de ker(f ). Sea (x1 , x2 , x3 ) ∈ R3 , entonces:
1 −1 0 x1 0
−1 1 0 x2 = 0 .
0 0 1 x3 0
verlo:
π(1, 0, 0) = [(1, 0, 0)]. Las coordenadas de [(1, 0, 0)] respecto de la base B serán los
escalares α, β ∈ R tales que:
[(1, 0, 0)] = α[(0, 1, 0)] + β[(0, 0, 1)].
Esto es, [(1, 0, 0)] = [(0, α, β)]. Por tanto se tendrá que vericar que (1, −α, −β) ∈
ker(f ). Así 1 = −α, β = 0. Por tanto, α = −1, β = 0. Esto es, π(1, 0, 0)B = (−1, 0).
Razonando de forma similar, se tiene que π(0, 1, 0)B = (1, 0), π(0, 0, 1)B = (0, 1).
f : R3 / ker(f ) −→ Im(f )
Puede comprobarse que una base de Im(f ) viene dada por B 0 = {(1, −1, 0), (0, 0, 1)}.
Así para calcular la matriz asociada a la aplicación f tendremos que calcular las
coordenadas de f ([(0, 1, 0)]) y f ([(0, 0, 1)]) respecto de la base B 0 . Puede comprobarse
que:
f ([(0, 1, 0)])B 0 = (−1, 0), f ([(0, 1, 0)])B 0 = (0, 1).
Así la matriz asociada a f en las bases indicadas es
−1 0
.
0 1
Solución:
a b
a) Sea A = c d la matriz pedida. Puede comprobarse que una base de ker(f )
e g
viene dada por {(2, −1)}. Además, por hipótesis se tiene que f (1, 1) = (0, 0, 1).
Así se tiene que vericar que:
a b 0 a b 0
c 2 c d 1
d = 0 , = 0 ,
−1 1
e g 0 e g 1
Éstos constituyen un sistema generador de Im(f ). Por tanto, una base de Im(f )
vendrá dada por {(0, 0, 1)} .
Calculamos ahora una base de R2 / ker(f ). Como una base de ker(f ) viene dada
por {(2, −1)}. Se tiene que una base de R2 / ker(f ) vendrá dada por B = {[(1, 0)]}.
c) Debemos calcular las matrices asociadas a las aplicaciones lineales que aparecen
en el siguiente diagrama conmutativo:
f
R2 −→ Rx
3
πy ]
i
R2 / ker(f ) −→ Im(f )
f
Matriz de π : R2 −→ R2 / ker(f ):
Tomamos la base canónica de R2 , {(1, 0), (0, 1)}, y la base que habíamos
calculado anteriormente de R2 / ker(f ), B = {[(1, 0)]}.
Tenemos que calcular las coordenadas de π((1, 0)) y de π((0, 1)) respecto de
la base B . Se tiene que π((1, 0)) = [(1, 0)]. Esto es, tendremos que calcular el
escalar α tal que
[(1, 0)] = α[(1, 0)].
Ahora tendremos que calcular las coordenadas de este vector respecto de la
base B . Esto es, tendremos que calcular el escalar α que verica 0, 0, 13 =
α(0, 0, 1). Por tanto, α = 13 . Así f ([(1, 0)])B = 13 . De donde podemos armar
que la matriz asociada, respecto a las bases indicadas es,
1
.
3
198 Más de 160 problemas resueltos de Álgebra Lineal
Matriz de i : Im(f ) −→ R3 :
Sea {(0, 0, 1)} la base de Im(f ) y consideremos en R3 la base canónica, B =
{(1, 0, 0), (0, 1, 0), (0, 0, 1)}.
Como i(0, 0, 1)B = (0, 0, 1)B = (0, 0, 1). La matriz asociada, respecto a las
bases indicadas, es
0
0 .
1
7 | Espacio dual
En este capítulo, dada una base B de un espacio vectorial V , calcularemos una base
B ∗ del espacio vectorial dual V ∗ y, recíprocamente, dada una base B ∗ del espacio
vectorial dual V ∗ , calcularemos una base B de V de la cual B ∗ es dual.
Solución: Vamos a recordar el concepto de base dual que será necesario en este
capítulo y que puede consultarse en [9].
Base dual
Sean B = {v1 , v2 , . . . , vn } una base de un espacio vectorial V y
B ∗ = {ϕ1 , ϕ2 , . . . , ϕn } una base de V ∗ , diremos que B ∗ es base dual
de B si para cada i = 1, 2, . . . , n, se verica:
Esto es,
1 0 0
(a11 a12 a13 ) −1 = 1, (a11 a12 a13 ) 1 = 0, (a11 a12 a13 ) 3 = 0.
3 −1 −2
199
200 Más de 160 problemas resueltos de Álgebra Lineal
Si (a31 a32 a33 ) es la matriz asociada a la forma lineal ϕ3 , también se tiene que:
1 0 0
(a31 a32 a33 ) −1 1 3 = 0 0 1 .
3 −1 −2
Esto es:
a11 a12 a13 1 0 0 1 0 0
a21 a22 a23 −1 1 3 = 0 1 0 .
a31 a32 a33 3 −1 −2 0 0 1
Vemos que se trata de calcular la inversa de la matriz cuyas columnas son las coor-
denadas de los vectores dados. Es decir:
−1
a11 a12 a13 1 0 0 1 0 0
a21 a22 a23 = −1 1 3 = 7 −2 −3 .
a31 a32 a33 3 −1 −2 −2 1 1
Se pide:
a) Probar que B ∗ = {f1 , f2 , f3 } es una base del espacio dual de R3 .
b) Hallar una base de R3 de la cual {f1 , f2 , f3 } es dual.
Solución:
α1 f1 + α2 f2 + α3 f3 = 0.
Tomamos la base canónica de R3 , B = {(1, 0, 0), (0, 1, 0), (0, 0, 1)}, y observemos
que:
= α1 · 2 + α2 · 0 + α3 · 0 = 2α1 .
Esto es,
2α1 + 0 + 0 = 0.
(α1 f1 +α2 f2 +α3 f3 )(0, 1, 0) = α1 f1 (0, 1, 0)+α2 f2 (0, 1, 0)+α3 f3 (0, 1, 0) = α2 −α1 .
Esto es,
−α1 + α2 + 0 = 0.
Y análogamente:
(α1 f1 +α2 f2 +α3 f3 )(0, 0, 1) = α1 f1 (0, 0, 1)+α2 f2 (0, 0, 1)+α3 f3 (0, 0, 1) = α1 +α2 +α3 .
Esto es,
α1 + α2 + α3 = 0.
Se pide:
a) Comprueba que B ∗ = {φ1 , φ2 , φ3 , φ4 } es una base de R3 [x]∗ .
b) Calcula una base de R3 [x] de la cual B ∗ es dual.
Solución:
Comprueba que {ϕ1 , ϕ2 } es una base de V∗ y halla la base {p1 (x), p2 (x)} de
V que tiene a {ϕ1 , ϕ2 } por dual.
Solución: Comenzamos viendo que {ϕ1 , ϕ2 } es base de V ∗ . Para ello, veremos que
es sistema generador del espacio y que es un conjunto linealmente independiente.
204 Más de 160 problemas resueltos de Álgebra Lineal
Observamos que:
Z 1
a
ϕ1 (ax + b) = (ax + b)dx = + b,
0 2
Z 2
ϕ2 (ax + b) = (ax + b)dx = 2a + 2b,
0
de donde:
1
a = −2ϕ1 (p(x)) + ϕ2 (p(x)), b = 2ϕ1 (p(x)) − ϕ2 (p(x)).
2
Es decir:
1
ϕ = [−2ϕ(x) + 2ϕ(1)]ϕ1 + ϕ(x) − ϕ(1) ϕ2 .
2
Nótese que, al ser card(B ∗ ) = 2, habría bastado con ver que el conjunto era
sistema generador o libre para poder armar que se trata de una base de V ∗ .
Ahora, para la base {1, x} de V se tiene que las matrices asociadas a las aplica-
ciones lineales respecto de dicha base son:
1
M [ϕ1 ] = 1 2 , M [ϕ2 ] = 2 2 .
La base {a11 + a21 x, a12 + a22 x} tal que {ϕ1 , ϕ2 } es su dual, deberá vericar:
1
1 2 a11 a12 1 0
= ,
2 2 a21 a22 0 1
de donde:
− 12
a11 a12 2
= .
a21 a22 −2 1
1
Por tanto, la base buscada es −2x + 2, x − .
2
Por tanto, f (e1 ) = −2, f (e2 ) = 3 y f (e3 ) = 1, así las coordenadas de f en la base B ∗
son (−2, 3, 1).
Solución:
Solución:
Por tanto,
ϕ1 (x1 , x2 , x3 ) = x1 + ax2 , (a ∈ R).
Por tanto,
ϕ2 (x1 , x2 , x3 ) = ϕ2 (x1 e1 + x2 e2 + x3 e3 ) = 3x2 + x3 .
Por tanto, {(−b, 0, 1), (−1, 1, 0)} es sistema generador de ker(ϕ3 ) y, puesto que ade-
más el conjunto es libre, es base. En consecuencia, ϕ3 (−b, 0, 1) = 0 y ϕ3 (−1, 1, 0) =
0. Además, por hipótesis se verica que ϕ3 (e1 ) = 2. Así ϕ3 verica las siguientes
condiciones:
ϕ3 (−be1 + e3 ) = 0,
ϕ3 (−e1 + e2 ) = 0,
ϕ3 (e1 ) = 2.
Esto es,
−bϕ3 (e1 ) + ϕ3 (e3 ) = 0,
−ϕ3 (e1 ) + ϕ3 (e2 ) = 0,
ϕ3 (e1 ) = 2.
Por tanto,
ϕ3 (x1 , x2 , x3 ) = 2x1 + 2x2 + 2bx3 , (b ∈ R dado).
c) Para que los vectores {ϕ1 , ϕ2 , ϕ3 } sean linealmente dependientes, el siguiente de-
terminante
1 a 0
0 3 1 = 2a + 6b − 2,
2 2 2b
debe de ser nulo. Por tanto, a y b deben ser tales que 2a + 6b − 2 = 0. Esto es,
a + 3b = 1.
8 | Espacio vectorial euclídeo
Se pide:
a) Matriz de Gram en la base {1, x}.
b) ¾Qué ángulo forman los polinomios x+3 y 2x + 4.?
c) Calcular una base ortonormal a partir de la base {1, x}.
Solución: Vamos a recordar la denición de matriz de Gram que nos será útil en
este ejercicio y en los que siguen. Esta denición, así como las propiedades de dicha
matriz pueden consultarse en [9].
Matriz de Gram
Sea (V, h , i un espacio vectorial ecuclídeo y sea B = {e1 , e2 , . . . , en } una
base de V . Llamaremos matriz de Gram respecto a la base B , a la matriz
A = (aij )1≤i,j≤n , donde aij = hei , ej i.
b) Vamos a recordar la denición de ángulo entre dos vectores que nos será útil en
este ejercicio y en los que siguen. Esta denición pueden consultarse en [9].
Vamos a calcular los distintos elementos que necesitamos para resolver el ejercicio,
Z 1
53
hx + 3, 2x + 4i = (x + 3)(2x + 4)dx = .
0 3
s r
Z 1
p 37
||x + 3|| = hx + 3, x + 3i = (x + 3)2 dx = .
0 3
s r
Z 1
p 19
||2x + 4|| = h2x + 4, 2x + 4i = (2x + 4)2 dx = 2 .
0 3
Por tanto:
53/3 53
cos(α) = q q = √ .
37
2 319 2 703
3
53
Así, α = arc cos √ .
2 703
c) Para obtener una base ortogonal a partir de una dada, utilizaremos el algoritmo
de Gram-Schmidt que puede consultarse en [9] y que recordamos a continuación:
8. Espacio vectorial euclídeo 211
w 1 = u1 ,
hu2 , w1 i
w 2 = u2 − w1 ,
||w1 ||2
hu3 , w1 i hu3 , w2 i
w 3 = u3 − w1 − w2 ,
||w1 ||2 ||w2 ||2
..
.
hun , w1 i hun , w2 i hun , wn−1 i
wn = un − w1 − w2 − · · · − wn−1 .
||w1 ||2 ||w2 ||2 ||wn−1 ||2
1
Observamos que x − no es un vector unitario, así que normalizamos este vector,
2
obteniendo: 1
x− 2 1
x− 2 √ √
1 = qR = 2 3x − 3.
||x − 2 || 1
(x − 1 )2 dx
0 2
√ √
Por tanto, la base ortonormal pedida es B = {1, 2 3x − 3}.
Solución:
h(1, 0, 0), (1, 0, 0)i = 2, h(1, 0, 0), (0, 1, 0)i = 1, h(1, 0, 0), (0, 0, 1)i = 1,
h(0, 1, 0), (1, 0, 0)i = 1, h(0, 1, 0), (0, 1, 0)i = 2, h(0, 1, 0), (0, 0, 1)i = 1,
h(0, 0, 1), (1, 0, 0)i = 1, h(0, 0, 1), (0, 1, 0)i = 1, h(0, 0, 1), (0, 0, 1)i = 2.
w1 = (1, 0, 0),
h(0, 1, 0), (1, 0, 0)i 1 1
w2 = (0, 1, 0) − (1, 0, 0) = (0, 1, 0) − (1, 0, 0) = − , 1, 0 ,
||(1, 0, 0)||2 2 2
h(0, 0, 1), (− 21 , 1, 0)i
h(0, 0, 1), (1, 0, 0)i 1
w3 = (0, 0, 1) − (1, 0, 0) − − , 1, 0 =
||(1, 0, 0)||2 ||(− 21 , 1, 0)||2 2
1 1 1 1 1
(0, 0, 1) − (1, 0, 0) − (− , 1, 0) = − , − , 1 .
2 3 2 3 3
Ahora, como u y v son unitarios, se tiene que ||u||2 = ||v||2 = 1. Así ||2u + v||2 =
4hu, vi + 5. Por otro lado, y puesto que los vectores u y v forman un ángulo de 60o se
tiene:
1
hu, vi = ||u||||v|| cos(60o ) = .
2
√
Luego ||2u + v||2 = 4hu, vi + 5 = 4
2 + 5 = 7 y, en consecuencia, ||2u + v|| = 7.
Solución:
a) Veamos que cumple las propiedades de un producto escalar, utilizando para ello
las propiedades del producto escalar usual h, i:
[x, y] = [y, x], siendo x = (x1 , x2 , x3 ), y = (y1 , y2 , y3 ) ∈ R3 ,
[x + z, y] = h(x1 + x3 + z1 + z3 , x1 + x2 + z1 + z2 , x2 + x3 + z2 + z3 ),
(y1 +y3 , y1 +y2 , y2 +y3 )i = h(x1 +x3 , x1 +x2 , x2 +x3 ), (y1 +y3 , y1 +y2 , y2 +y3 )i+
+h(z1 + z3 , z1 + z2 , z2 + z3 ), (y1 + y3 , y1 + y2 , y2 + y3 )i = [x + y] + [z + y].
[x, x] = 0 ⇔ x = 0.
⇒ Sea x ∈ R tal que [x, x] = 0, entonces h(x1 + x3 , x1 + x2 , x2 + x3 ), (x1 +
x3 , x1 + x2 , x2 + x3 )i = 0, de donde (x1 + x3 , x1 + x2 , x2 + x3 ) = (0, 0, 0) y
resolviendo el sistema resulta: x1 = 0, x2 = 0 y x3 = 0, es decir, x = 0.
⇐ Si x = 0, entonces (x1 + x3 , x1 + x2 , x2 + x3 ) = (0, 0, 0) y de ahí se tiene
que [x, x] = h(x1 + x3 , x1 + x2 , x2 + x3 ), (x1 + x3 , x1 + x2 , x2 + x3 )i = 0.
Solución:
b) En primer lugar, observamos que los vectores dados son linealmente independien-
tes. Por tanto forman una base de F . Vamos a calcular una base ortogonal y
luego, multiplicando cada vector por el inverso de su norma, obtendremos una
base ortonormal.
w1 = (1, 2, 1, 0),
h(0, −1, 1, 0), (1, 2, 1, 0)i 1 2 7
w2 = (0, −1, 1, 0) − (1, 2, 1, 0) = ,− , ,0 ,
||(1, 2, 1, 0)||2 6 3 6
h(1, 1, −2, 1), (1, 2, 1, 0)i
w3 = (1, 1, −2, 1) − (1, 2, 1, 0)−
||(1, 2, 1, 0)||2
216 Más de 160 problemas resueltos de Álgebra Lineal
hA, Bi = hB, Ai, A, B ∈ M2 (R). Sabiendo que tr(A) = tr(At ), se tiene que
hA, αBi.
a11 a12
hA, Ai ≥ 0, A ∈ M2 (R). Consideremos A = ∈ M2 (R). Así
a21 a22
pues
hA, Ai = tr(AAt ) = a211 + a212 + a221 + a222 ≥ 0.
Por tanto, existe una base ortonormal para el producto escalar dado, siendo esta
base B .
218 Más de 160 problemas resueltos de Álgebra Lineal
w1 = (1, 1, 1, 1),
h(1, 1, 0, 0), (1, 1, 1, 1)i 1 1 1 1
w2 = (1, 1, 0, 0) − (1, 1, 1, 1) = , ,− ,− ,
||(1, 1, 1, 1)||2 2 2 2 2
h(1, 1, 1, 0), (1, 1, 1, 1)i
w3 = (1, 1, 1, 0) − (1, 1, 1, 1)−
||(1, 1, 1, 1)||2
h(1, 1, 1, 0), ( 12 , 12 , − 12 , − 12 )i 1 1 1 1
1 1
, ,− ,− = 0, 0, , − ,
||( 12 , 12 , − 12 , − 12 )||2 2 2 2 2 2 2
h(1, 0, 0, 0), (1, 1, 1, 1)i
w4 = (1, 0, 0, 0) − (1, 1, 1, 1)−
||(1, 1, 1, 1)||2
h(1, 0, 0, 0), ( 12 , 12 , − 12 , − 12 )i 1 1 1 1
, ,− ,− −
||( 12 , 12 , − 12 , − 12 )||2 2 2 2 2
h(1, 0, 0, 0), (0, 0, 12 , − 12 )i
1 1 1 1
0, 0, , − = , − , 0, 0 .
||(0, 0, 12 , − 12 )||2 2 2 2 2
Solución:
Luego, un sistema generador de F es B = {(1, 0, 2, 0), (0, 1, 4, 0), (0, 0, 3, 1)}, que
al ser libre, es una base.
Solución: Se observa que los vectores (2, 1, 5, −1) y (−1, 1, 2, 0) son linealmente
independientes. Luego BU = {(2, 1, 5, −1), (−1, 1, 2, 0)} es base de U . Si (x, y, z, t) ∈
U ⊥ , entonces se debe vericar que:
h(x, y, z, t), (2, 1, 5, −1)i = 0, h(x, y, z, t), (−1, 1, 2, 0)i = 0.
Por tanto se tiene el sistema,
2x +y +5z −t = 0
−x +y +2z = 0,
cuyas soluciones son:
x = α + 2β, y = α, z = β, t = 3α + 9β, α, β ∈ R.
Así pues, una base de U ⊥ es BU ⊥ = {(1, 1, 0, 3), (2, 0, 1, 9)}.
220 Más de 160 problemas resueltos de Álgebra Lineal
Solución: Sea F el subespacio generado por {(2, 1, 0, 1), (0, 3, 1, 1)}. Como este con-
junto es libre, constituye una base de F . Calculemos ahora una base para F ⊥ . Sea
(x, y, z, t) ∈ F ⊥ entonces:
h(x, y, z, t), (2, 1, 0, 1)i = 0
h(x, y, z, t), (0, 3, 1, 1)i = 0.
Por tanto, una base de F ⊥ es BF ⊥ = {(1, −2, 6, 0), (0, −1, 2, 1)}. Escribamos ahora el
vector (1, 3, −1, 4) como suma de dos vectores, uno de F y otro de F ⊥ . Esto es,
(1, 3, −1, 4) = a(2, 1, 0, 1) + b(0, 3, 1, 1) + c(1, −2, 6, 0) + d(0, −1, 2, 1).
51 18 26 113
Obtenemos un sistema cuya solución es a = ,b = ,c = − yd= .
50 25 25 50
Por tanto,
51 159 18 87 26 9 43 113
(1, 3, −1, 4) = , , , + − ,− ,− , ,
25 50 25 50 25 50 25 50
51 159 18 87 26 9 43 113
donde , , , ∈ F y − ,− ,− , ∈ F ⊥.
25 50 25 50 25 50 25 50
Solución:
Diagonalización
1 0 0
1. Dada la matriz 2 1 0 , se pide:
3 0 2
a) Estudiar si 3 es o no un autovalor de A.
b) ¾Son los vectores u = (1, 1, 1) y v = (0, 0, 1) autovectores de A?. En caso
armativo buscar a qué autovalor están asociados.
c) Encontrar un polinomio que se anule para la matriz A.
Solución:
a) Para resolver este ejercicio, utilizaremos la denición de autovalor. Esto es, bus-
caremos un vector de la forma (x, y, z) ∈ R3 \ {(0, 0, 0)} tal que:
1 0 0 x x
2 1 0 y = 3 y .
3 0 2 z z
223
224 Más de 160 problemas resueltos de Álgebra Lineal
Esto es,
1 λ
3 = λ .
5 λ
Teorema de Cayley-Hamilton
Toda matriz verica su polinomio característico.
a 1 p
3. Hallar los autovalores de A = b 2 q , sabiendo que admite como
c −1 r
autovectores a los vectores (1, 1, 0), (−1, 0, 2), (0, 1, −1).
226 Más de 160 problemas resueltos de Álgebra Lineal
De aquí obtenemos, a − b = 1 y c = 1.
Por ser (−1, 0, 2) un autovector, existe λ ∈ R tal que
a 1 p −1 −1
b 2 q 0 = λ 0 .
1 −1 r 2 2
de donde obtenemos, 2r − 1 = 2a − 4p y b = 2q .
Por ser (0, 1, −1) un autovector, existirá λ ∈ R tal que
a 1 p 0 0
2q 2 q 1 = λ 1 .
1 −1 r −1 −1
de aquí, 2 − q = r + 1 y p = 1.
Si resolvemos el sistema formado por todas las ecuaciones obtenidas, se tiene que
1 1
a = 2, b = 1, c = 1, p = 1, q= , r= .
2 2
Calculemos ahora los autovalores de A. Para ello calculamos las raíces de su polinomio
característico. Esto es
2−λ 1 1
1
PA (λ) = 1 2−λ 2
= 0.
1
1 −1 2 −λ
Es decir,
9λ2 9λ
−λ3 + − = 0.
2 2
3
De aquí podemos armar que los autovalores son: λ1 = 0, λ2 = y λ3 = 3.
2
Solución: Calcularemos los autovalores de la matriz. Para ello calculamos las raíces
del polinomio característico:
α−λ α−1 1−α
α−2 α−λ 2−α = 2α−(2+3α)λ+(3+α)λ2 −λ3 = (λ−1)(λ−2)(λ−α) = 0.
α−2 α−1 3−α−λ
Solución:
b) Se tiene que αλ es autovalor de αA si, y sólo si, |αA − αλI| = 0. Veamos que se
cumple esta condición,
|αA − αλI| = |α(A − λI)| = αn |A − λI| = 0,
y de ahí el resultado, ya que λ es autovalor de A.
c) Se tiene que λ es autovalor de A. Por tanto
|(A − αI) − (λ − α)I| = |A − αI − λI + αI| = |A − λI| = 0.
7. Se considera la matriz:
1 1 1 1
1 1 −1 −1
A= .
1 −1 1 −1
1 −1 −1 1
Solución: Vamos a recordar dos resultados que nos serán útiles en los siguientes
ejercicios y que pueden consultarse en [10]:
Endomorsmo diagonalizable
Sean V un K-espacio vectorial de dimensión n 6= 0 y f : V −→ V un
endomorsmo. Entonces f es diagonalizable si y solo si se verican las dos
propiedades siguientes:
1. el polinomio característico de f posee n raíces en K, cuyos órdenes
de multiplicidad suman n (propiedad que siempre se verica cuando
K = C),
2. para cada raíz λ ∈ K del polinomio caracteristico de f , se verica
dim(V (λ)) = d, siendo d el orden de multiplicidad de λ.
Como consecuencia inmediata, podemos armar que si A es una matriz cua-
drada de orden n que tiene n autovalores distintos en K, entonces A es
diagonalizable.
230 Más de 160 problemas resueltos de Álgebra Lineal
1−λ 1 1 1
1 1−λ −1 −1
= −16 + 16λ − 4λ3 + 4λ4 = (λ − 2)3 (λ + 2).
1 −1 1−λ −1
1 −1 −1 1 − λ
Así, una base de V (2) vendrá dada por BV (2) = {(1, 0, 0, 1), (0, 1, 0, −1), (0, 0, 1, −1)}.
En este caso se tiene que la multiplicidad de λ = 2 coincide con la dimensión de su
subespacio propio, V (2). Por tanto podemos armar que la matriz A es diagonalizable.
Por tanto, tenemos que:
−1
−2 0 0 0 1 1 0 0 1 1 1 1 1 1 0 0
0 2 0 0 −1 0 1 0 1 1 −1 −1 −1 0 1 0
= ,
0 0 2 0 −1 0 0 1 1 −1 1 −1 −1 0 0 1
0 0 0 2 −1 1 −1 −1 1 −1 −1 1 −1 1 −1 −1
donde
1 1 0 0
−1 0 1 0
P =
−1
0 0 1
−1 1 −1 −1
Diagonalicemos ahora A−1 . Esto es, hallemos una matriz diagonal semejante a A−1 .
Por las anteriores igualdades, se tiene que
A−1 = (P DP −1 )−1 = P D−1 P −1 .
Solución:
A
Empezamos por calcular los valores propios de A. Para ello calcularemos su
polinomio característico.
1 −1 0
PA (λ) = |A − λI| = −1 1 0 = λ(λ − 2)(5 − λ).
0 0 5
Así, una base de V (2) será BV (2) = {(1, −1, 0)}. Finalmente, para V (5), resol-
vemos la ecuación matricial (A − 5I)X = 0. Esto es,
−4 −1 0 x1 0
−1 −4 0 x2 = 0 ,
0 0 0 x3 0
B
En primer lugar calculamos el polinomio característico de B . Esto es,
2−λ −1 0 0
−1 2−λ 0 0
PB (λ) = = λ(λ − 2)(λ − 3)(λ − 1).
0 0 1−λ −1
0 0 −1 1−λ
Por último, calculemos una base del subespacio V (3). Para ello considere-
mos el sistema matricial (A − 3I)X = 0. Esto es,
−1 −1 0 0 x1 0
−1 −1 0 0 x2 0
=
0 .
0 0 −2 −1 x3
0 0 −1 −2 x4 0
9. Valores y vectores propios. Diagonalización 235
C
Hallemos su polinomio característico
1−λ 2 1 √ √
PC (λ) = |C −λI| = 0 1−λ 1 = −(λ−1)(λ−(1+ 2))(λ−(1− 2)).
0 2 1−λ
√ √
Por tanto, los valores propios de C son λ = 1, λ = 1 + 2 y λ = 1 − 2, con
multiplicidad uno, por lo que la dimensión de los correspondientes subespacios
propios asociados es uno, luego C es diagonalizable. Calculamos una base de
V (1). Para ello, comenzamos resolviendo el sistema matricial (C − 1I)X = 0.
Esto es,
0 2 1 x1 0
0 0 1 x2 = 0 ,
0 2 0 x3 0
Esto es, √
− 2 √2 1 x1 0
0 − 2 √1 x2 = 0 ,
0 2 − 2 x3 0
es decir √
− 2x1 +2x2
√ +x3 = 0
− 2x2 +x3 = 0,
√ √
= (1 + 2)α, x2 = α y x3√= 2α
cuya solución es x1 √ √ , siendo α ∈ R. Por tanto
una base de V (1 + 2) es BV (1+√2) = {(1 + 2, 1, 2)}.
√
Por último, para V (1 − 2), repetimos el mismo procedimiento. Esto es, se
considera el sistema matricial,
√
(C − (1 − 2)I)X = 0.
Es decir, √
2 √2 1 x1 0
0 2 1 x2 = 0 ,
√
0 2 2 x3 0
quedándonos el sistema
√
2x1 +2x2
√ +x3 = 0
2x2 +x3 = 0,
√ √
cuya solución es x1 = (1 − 2)α√, x2 = √α y x3 = − 2α, siendo α √ ∈ R. En
consecuencia BV (1−√2) = {(1 − 2, 1, − 2)} es una base de V (1 − 2). Así
pues, podemos armar que la matriz C es diagonalizable. Esto es, se tiene que
D = P −1 CP . Es decir,
se tiene que
1 0√ 0
0 1+ 2 0√ =
0 0 1− 2
√ √ −1 √ √
1 1+ 2 1− 2 1 2 1 1 1+ 2 1− 2
0
√1 1 0 1 1 0
√1 1 .
√ √
0 2 − 2 0 2 1 0 2 − 2
1
están asociados respectivamente a los autovalores λ=1 y λ= . Se pide:
2
a) Matriz diagonal semejante a la matriz A.
b) Calcular la matriz 2A4 − 7A3 + 9A2 − 5A + I.
c) Calcular la matriz A−3 − 4A−2 + 5A−1 + 4I.
Solución:
Solución: Observamos que A y B son matrices simétricas reales. Estas matrices son
siempre diagonalizables como arma el siguiente resultado que puede consultarse en
[9]:
Así, B = BV (2) ∪ BV (5) = {(−1, 1, 0), (−1, 0, 1), (1, 1, 1)} es una base de R3 .
Los vectores de BV (2) son ortogonales al vector de BV (5) , pero los vectores de
BV (2) no son ortogonales entre sí. Por ello, utilizaremos el procedimiento de
Gram Schmidt para calcular una base ortogonal de V (2).
w1 = (−1, 1, 0),
h(−1, 0, 1), (−1, 1, 0)i 1 1
w2 = (−1, 0, 1) − (−1, 1, 0) = − ,− ,1 .
||(−1, 1, 0)||2 2 2
Y normalizando los vectores, nos queda la base
1 1 1 1 2 1 1 1
B0 = −√ , √ , 0 , −√ , −√ , √ , √ ,√ ,√ .
2 2 6 6 6 3 3 3
Así, B = BV (2) ∪ BV (−1) = {(1, 0, −1), (0, 1, −1), (1, 1, 1)} es una base de R3 .
Los vectores de BV (2) son ortogonales al vector de BV (−1) , pero los vectores
de BV (2) no lo son entre sí. Para ello, utilizaremos el procedimiento de Gram-
Schmidt para obtener una base ortogonal de V (2):
w1 = (1, 0, −1),
h(0, 1, −1), (1, 0, −1)i 1 1
w2 = (0, 1, −1) − (1, 0, −1) = − , 1, − .
||(0, 1, −1)||2 2 2
El último de los eslabones, ker(f − λI)k , para el que ker(f − λI)k = ker(f −
λI)k+1 pero ker(f − λI)k−1 6= ker(f − λI)k recibe el nombre de subespacio
máximo del autovalor λ.
A continuación vamos ver cómo podemos obtener una base del subespacio máxi-
mo del autovalor λ.
Con estos datos calculados, vamos a ver cómo podemos interpretarlos para obtener
la matriz de Jordan.
−4 0 −4 11
−6 3 −6 9
A=
5 0 5 −7
−1 0 −1 5
Calculemos primero el polinomio característico:
2α + γ
con soluciones x1 = α, x2 = β, x3 = γ, x4 = con α, β, γ ∈ R.
3
Así, B3 = {(3, 0, 0, 2), (0, 1, 0, 0), (0, 0, 3, 1)} es una base para ker(f − 3I)3 . De
donde, la dimensión de ker(f − 3I)3 es 3. Se puede comprobar que dim(ker(f −
3I)n ) = 3, para n ≥ 3. Por tanto nos detenemos aquí.
244 Más de 160 problemas resueltos de Álgebra Lineal
vericándose que J = P −1 AP .
1 3 2 0
0 1 0 0
B= .
0 2 −2 0
6 6 −3 4
Calculemos primero el polinomio característico. Esto es,
3α 7α
cuyas soluciones son x1 = α, x2 = 0, x3 = − , x4 = − con α ∈ R.
2 4
Luego una base de V (−2) es BV (−2) = {(4, 0, −6, −7)}. Observamos que la
dimensión de V (−2) coincide con la multiplicidad de λ = −2.
Calculemos ahora una base de V (1). Para ello consideramos el sistema matricial
(B − I)X = 0. Esto es,
0 3 2 0 x1 0
0 0 0 0 x2 0
= .
0 2 −3 0 x3 0
6 6 −3 3 x4 0
Esto es,
4x2 −6x3 = 0
−6x2 +9x3 = 0
18x1 +30x2 +12x3 +9x4 = 0,
2β 38β
cuyas soluciones son x1 = α, x2 = β, x3 = , x4 = −2α − con α, β ∈ R.
3 9
Así, una base es B2 = {(1, 0, 0, −2), (0, 9, 6, −38)}. Por tanto, dim ker(f − I)2 =
246 Más de 160 problemas resueltos de Álgebra Lineal
vericándose que J = P −1 BP .
2 −1 0 0
−1 2 0 0
C= .
0 0 1 −1
0 0 −1 1
9. Valores y vectores propios. Diagonalización 247
Por tanto, la matriz es diagonalizable, ya que todos sus autovalores son diferen-
tes, siendo éstos λ = 3, λ = 2, λ = 1 y λ = 0. Estudiemos ahora los subespa-
cios propios. Empezaremos por V (3). Para ello planteamos el sistema matricial
(C − 3I)X = 0, que resolviéndolo igual que en apartados anteriores, nos queda
x1 = α, x2 = −α, x3 = 0, x4 = 0, siendo α ∈ R. Por tanto, una base de V (3)
viene dada por BV (3) = {(1, −1, 0, 0)}.
cumpliéndose que D = P −1 CP .
Se pide:
a) Calcula la forma canónica de Jordan J de la matriz de f en la base
canónica.
b) Calcula la base de R3 en la que la matriz de f es J .
248 Más de 160 problemas resueltos de Álgebra Lineal
c) Calcula A347 .
Solución:
(A − 0I)2 v1 −→ (A − 0I)v1 −→ v1
donde v1 es un vector que pertenece a ker(f − 0I)3 y no pertenece a ker(f − 0I)2 .
x −y +t = 0 x +2y −3z +4t = 2
i) j)
−x +2y +z = 0 2x +5y −2z +t = 1
y −3z +2t = 0 5x +12y −7z +6t = 7
x +y −z +2t = 1
x −3y −2z +4t = 5
x −y +z −t = 2
k) l)
3x −8y −3z +8t = 18
2x +2y −z +5t = 0
2x −3y +5z −4t = 19
y +4z −3t = 5
resultado.
h) In[2]:= Solve[{x + y + z + t == 0, x + z == 0, x + 3 y + z + 3 t == 0,
x - y + z + t == 0}, {x, y, z, t}]
Solve::svars: Equations may not give solutions for all "solve"variables.
Out[2]= {{y → 0, z → −x, t → 0}}
Igual que en el apartado anterior observamos que el sistema es compatible inde-
terminado.
i ) In[3]:= Solve[{x - y + t == 0, -x + 2 y + z == 0,
y - 3 z + 2 t == 0}, {x, y, z, t}]
Solve::svars: Equations may not give solutions for all "solve"variables.
5x x 4x
Out[3]= y→ ,z → − ,t → −
9 9 9
Nuevamente el sistema es compatible indeterminado.
j ) In[4]:= Solve[{x + 2 y - 3 z + 4 t == 2, 2 x + 5 y - 2 z + t == 1,
5 x + 12 y - 7 z + 6 t == 7}, {x, y, z, t}]
Out[4]= {}
Con el símbolo {}, Mathematica nos informa que el sistema es incompatible.
k) In[5]:= Solve[{x - 3 y - 2 z + 4 t == 5,
3 x - 8 y - 3 z + 8 t == 18, 2 x - 3 y + 5 z - 4 t == 19},
{x, y, z, t}]
Solve::svars: Equations may not give solutions for all "solve"variables.
21 3x 7y
Out[5]= z → 8 − x + 2y, t → − +
4 4 4
Nuevamente el sistema es compatible indeterminado.
l ) In[6]:= Solve[{x + y - z + 2 t == 1, x - y + z - t == 2,
2 x + 2 y - z + 5 t == 0, y + 4 z - 3 t == 5}, {x, y, z, t}]
44 30 7 31
Out[6]= x→ ,y → ,z → − ,t → −
19 19 19 19
En este caso, el sistema es compatible determinado.
Solución: En este ejercicio vamos a ver cómo podemos podemos obtener la forma
escalonada de una matriz utilizando la función RowReduce, que nos calcular la forma
escalonada reducida de la matriz. También mostraremos cómo la función MatrixRank
proporciona directamente el rango de la matriz dada.
Comenzamos introduciendo la matriz A,
In[1]:= a={{1, 2, 3}, {4, 5, 6}, {7, 8, 9}}
Out[1]== {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}}
Ahora aplicamos la orden RowReduce a la matriz,
In[2]:= RowReduce[a]
Out[4]= 2
Para la matriz B , se tiene:
In[5]:= RowReduce[{{0, 1, 2}, {-1, 0, 3}, {-2, -3, 0}}] // MatrixForm
1 0 −3
Out[5]//MatrixForm= 0 1 2
0 0 0
El rango de B es por tanto 2. Con la función MatrixRank, obtenemos el mismo
resultado:
In[6]:= MatrixRank[{{0, 1, 2}, {-1, 0, 3}, {-2, -3, 0}}]
Out[6]= 2
Para la matriz C , se tiene:
In[7]:= RowReduce[{{0, 1, 2}, {1, 0, 3}, {2, 3, 0}}] // MatrixForm
1 0 0
Out[7]//MatrixForm= 0 1 0
0 0 1
254 Más de 160 problemas resueltos de Álgebra Lineal
Out[8]= 3
[C2P2] Resuelve:
−1 3 2 4 1 1
1 2 3 3 2 1
a) 4 5 6
+
2 3 4
, b) 4 0 1 + 3 2 −2 .
−2 1 5 1 2 −3
Solución: En este ejercicio veremos cómo Mathematica multiplica un escalar por una
matriz. Si no queremos que nos devuelva los datos que hemos introducido, bastará
escribir al nal de la línea el símbolo ;
In[1]:= a = {{1, 3, 5, 7}, {2, 6, 10, 14}, {1, -3, 6, -7}};
In[4]:= 2 a // MatrixForm
2 6 10 14
Out[4]//MatrixForm= 4 12 20 28
2 −6 12 −14
In[5]:= 3^2 a // MatrixForm
9 27 45 63
Out[5]//MatrixForm= 18 54 90 126
9 −27 54 −63
calcular A2 , B 2 , C 2 , A · B , B · A, A · C , B · C y (A − C)B .
Solución:
mx −3y +z =0
x +my +z =0
3x −y +mz =0
Calculamos su determinante:
In[2]:= Det[a]
Out[2]= −10 + m + m3
Aquí vemos que la única raíz real es m = 2. Luego si m 6= 2, se tiene que rg(A) = 3.
En este caso la única solución del sistema sera la trivial. Si m = 2, entonces rg(A) = 2
y el sistema tendrá soluciones diferentes de la trivial.
10. Trabajando con Mathematica 259
Se pide:
a) Hallar la expresión de cambio de base entre B1 y B2 .
b) Determinar las coordenadas respecto de B1 del vector x cuyas coor-
denadas respecto de B2 son (2, 1, 0, −1).
Solución:
Solución:
In[2]:= D[p[x], x]
Out[2]= 3x2
Out[10]= 2a
y
In[11]:= Det[{{1,-1,b},{0,0,c},{1,1,d}}]
Out[11]= −2c
tenemos que las ecuaciones de F vendrán dadas por a = c = 0.
Vamos ahora a calcular las ecuaciones cartesianas de W . Así, si ax3 + bx2 +
cx + d ∈ W , entonces
ax3 + bx2 + cx + d = α(1 + 3x + 5x2 ) + β(−1 + 2x2 ) = (5α + 2β)x2 + 3αx + (α − β).
262 Más de 160 problemas resueltos de Álgebra Lineal
Out[13]= −6a
y
In[14]:= Det[{{5,2,b},{3,0,c},{1,-1,d}}]
Out[14]= −3b + 7c − 6d
tendremos que las ecuaciones cartesianas de W vendrán dadas por:
a = 0
3b − 7c + 6d = 0.
Out[15]= 3
Este resultado nos arma que estos 4 vectores son ligados, pero que existen 3 que
son linealmente independientes. Como
In[16]:= MatrixRank[{{1, 0, -1}, {1, 3, 5}, {-1, 0, 2}}]
Out[16]= 3
se tiene que una base de F +W viene dada por B = {1−x2 , 1+3x+5x2 , −1+2x2 }.
10. Trabajando con Mathematica 263
Solución:
Out[1]= 3
Como el rango vale 3, son libres y de ahí que forme una base de F . Por tanto
dim(F ) = 3.
b) De las ecuaciones cartesianas de G, obtenemos que un sistema generador viene
dado por {(1, −1, 0, 0), (0, 0, 2, 1)}, y procediendo como en el apartado anterior, se
tiene que:
In[2]:= MatrixRank[{{1, -1, 0, 0}, {0, 0, 2, 1}}]
Out[2]= 2
De ahí obtenemos que una base de G es {(1, −1, 0, 0), (0, 0, 2, 1)} y dim(G) = 2.
c) Bastará resolver el sistema formado por las ecuaciones cartesianas de F y las
ecuaciones cartesianas de G,
In[3]:= Solve[{a + b == 0, b - c - d == 0, c - 2 d == 0}]
Solve::svars: Equations may not give solutions for all "solve"variables.
Out[3]= {{a → −3d, b → 3d, c → 2d}}
Esto es, un sistema generador de F ∩ G viene dado por {(−3, 3, 2, 1)}, que al ser
libre, forma una base. Así dim(F ∩ G) = 1.
d) Se tratará de calcular el rango de la matriz cuyas las son los vectores de la base
de F y la base de G. Como
In[4]:= MatrixRank[{{1, 0, 0, 0}, {0, 1, 0, 1}, {0, 0, 1, -1},
{1, -1, 0, 0}, {0, 0, 2, 1}}]
Out[4]= 4
y
In[5]:= MatrixRank[{{1, 0, 0, 0}, {0, 1, 0, 1}, {0, 0, 1, -1},
{1, -1, 0, 0}}]
264 Más de 160 problemas resueltos de Álgebra Lineal
Out[5]= 4
se tiene que una base de F + G viene dada por {(1, 0, 0, 0), (0, 1, 0, 1), (0, 0, 1, −1),
(1, −1, 0, 0)} y de ahí que dim(F + G) = 4.
Out[1]= 2
En este caso, se tiene que los vectores dados son ligados. Como
In[2]:= MatrixRank[{{1, 0, 1, 1}, {1, 2, 1, 1}}]
Out[2]= 2
una base de V viene dada por {(1, 0, 1, 1), (1, 2, 1, 1)}. Ampliamos esta base hasta
obtener una base de R4 . Al vericarse que
In[3]:= MatrixRank[{{1, 0, 1, 1}, {1, 2, 1, 1}, {0, 0, 1, 0},
{0, 0, 0, 1}}]
Out[3]= 4
se tiene que una base de R4 /V será {[(0, 0, 1, 0)], [(0, 0, 0, 1)]}.
Hallar:
a) Matriz asociada a (f + g).
b) Base y dimensión de ker(f + g) e Im(f + g).
c) ¾Qué tipo de aplicación es (f + g)?
d) Calcular f −1 (4, 6).
Solución:
Por tanto, una base de ker(f + g) es {(1, 0, 1), (0, 1, 0)} y dim(ker(f + g)) = 2.
Un sistema generador Im(f +g) viene dado por {(f +g)(1, 0, 0), (f +g)(0, 1, 0), (f +
g)(0, 0, 1)}. Esto es,
Comprueba que {ϕ1 , ϕ2 } es una base de V∗ y halla la base {p1 (x), p2 (x)} de
V que tiene a {ϕ1 , ϕ2 } por dual.
[C8P2] Calcula la matriz de Gram, respecto de la base B = {(1, 0, 0), (0, 1, 0),
(0, 0, 1)}, del producto escalar dado por:
Solución:
Solución:
Solución: Comenzamos calculando una base del espacio ortogonal al dado. Para ello
utilizaremos,
In[1]:= Solve[{{x, y, z, t}.{2, 1, 0, 1} == 0,
{x, y, z, t}.{0, 3, 1, 1} == 0
}, {x,y,z,t}]
Solve::svars: Equations may not give solutions for all "solve"variables.
Out[1]= {{x → −t/3 + z/6, y → −t/3 − z/3}}
10. Trabajando con Mathematica 269
Así una base del espacio ortogonal al dado es B = {(1, −2, 6, 0), (0, −1, 2, 1)}.
Ahora se trata de expresar el vector (1, 3, −1, 4) como,
(1, 3, −1, 4) = a(2, 1, 0, 1) + b(0, 3, 1, 1) + c(1, −2, 6, 0) + d(0, −1, 2, 1).
Solución:
En este ejercio veremos cómo Mathematica puede calcular directamente los elementos
necesarios para ver si una matriz es diagonalizable.
Introducimos la matriz A
In[1]:= a = {{1, -1, 0}, {-1, 1, 0}, {0, 0, 5}};
Calculamos su polinomio característico,
In[2]:= CharacteristicPolynomial[a, x]
Vemos que los autovalores de la matriz A son: {5, 2, 0}. Éstos se podrían haber
calculado directamente mediante la orden Eigenvalues:
270 Más de 160 problemas resueltos de Álgebra Lineal
In[4]:= Eigenvalues[a]
Out[4]= {5, 2, 0}
Luego la matriz es diagonalizable ya que es de rango 3 y tiene tres autovalores
diferentes.
Vamos a calcular los autovectores correspondientes a cada valor propio.
Para el valor propio λ = 5, tendremos que resolver el sistema matricial
(A − 5I)X = 0.
In[5]:= Solve[(a-5 IdentityMatrix[3]).{x,y,z}==0,{x,y,z}]
Solve::svars: Equations may not give solutions for all "solve"variables.
Out[5]= {{x → 0, y → 0}}
Así una base de V (5) vendrá dada por {(0, 0, 1)}.
Para el valor propio λ = 2, tendremos que resolver el sistema matricial
(A − 5I)X = 0.
In[9]:= Eigensystem[a]
Out[11]= {{3, 2, 1, 0}, {{−1, 1, 0, 0}, {0, 0, −1, 1}, {1, 1, 0, 0}, {0, 0, 1, 1}}}
Por lo tanto, podemos armar que la matriz B es diagonalizable, siendo la matriz
diagonal D, semejante a la matriz A, y la matriz de paso P ,
3 0 0 0 −1 0 1 0
0 2 0 0 1 0 1 0
D=
, P =
.
0 0 1 0 0 −1 0 1
0 0 0 0 0 1 0 1
Para la matriz C ,
In[12]:= c = {{1, 2, 1}, {0, 1, 1}, {0, 2, 1}};
In[13]:= Eigensystem[c]
√ √ √ √
Out[13]= {{1 + 2, 1, 1 − 2}, {{1/2
√ + 1/2(1 + 2), −1/2
√ + 1/2(1 + 2), 1},
{1, 0, 0}, {1/2 + 1/2(1 − 2), −1/2 + 1/2(1 − 2), 1}}}
Así, podemos armar que la matriz C es diagonalizable, siendo la matriz dia-
gonal D, semejante a la matriz C , y la matriz de paso P ,
√ 1
√ 1
√
1+ 2 0 0 2 2+ 2 1 2 2− 2
D= 0 1 0√ , P = √1 0 − √12 .
2
0 0 1− 2 1 0 1
Solución:
In[2]:= Eigensystem[a]
Para la matriz B ,
In[4]:= b = {{1, -1, -1}, {-1, 1, -1}, {-1, -1, 1}};
In[5]:= Eigensystem[b]
Todas las órdenes, así como la sintaxis utilizada en los ejercicios anteriores, se
describen a continuación,
Orden Signicado
% El último resultado generado.
t /. x → a Reemplaza x por a en la expresión de t.
x.y Producto escalar entre los vectores x e y .
x.m Producto entre el vector x y la matriz m.
m .x Producto entre la matriz m y el vector x.
m + n Suma las matrices m y n.
x m Producto entre la matriz m un número x.
m .n Producto entre las matrices m y n.
CharacteristicPo\ Polinomio característico de la matriz m.
lynomial[m,x]
Det[m] Determinante de la matriz m.
D[f ,x] Derivada de la función f respecto de la variable x.
Eigenvectors[m] Autovectores de la matriz m.
Eigenvalues[m] Autovalores de la matriz m.
Eigensystem[m] Autovalores y autovectores de la matriz m.
Factor[poly ] Factoriza el polinomio poly sobre los números enteros.
FullSimplify[expr] Simplica la expresión expr.
IdentityMatrix[n] Genera una matriz identidad de orden n.
Inverse[m] Matriz inversa de la matriz m.
Integrate[f (x), Integra la función f en el rango indicado.
{x,xmin ,xmax }]
JordanDecom\ Devuelve una lista {s,j }, donde s es la matriz de paso
position[m,x] y j es la forma canónica de Jordan de la matriz m.
MatrixForm[m] Muestra la matriz m en forma matricial.
MatrixPower[m,n] n-ésima potencia la matriz m.
MatrixRank[m] Rango de la matriz m.
Normalize[v ] Normaliza el vector v .
Orthogonalize[{e1 , Base ortonormal a partir de los vectores ei y el
e2 ,...},f ] producto escalar dado por f . (En caso de no
especicar f se utiliza el producto escalar usual)
RowReduce[m] Matriz reducida por las equivalente a la matriz m.
Solve[expr,vars] Resuelve el sistema expr de ecuaciones respecto de las
variables vars.
Transpose[m] Matriz traspuesta de la matriz m.
Glosario
A M
Algoritmo Matriz
de Gram-Schmidt . . . . . . . . . . . 211 antisimétrica . . . . . . . . . . . . . . . . . 36
para calcular la matriz inversa 40 de Gram . . . . . . . . . . . . . . . . . . . . 209
para obtener un sistema de Jordan. . . . . . . . . . . . . . . . . . . 242
escalonado por las . . . . . . 14 idempotente . . . . . . . . . . . . . . . . . . 38
para obtener una matriz simétrica . . . . . . . . . . . . . . . . . . . . . 36
escalonada por las . . . . . . 26 traspuesta . . . . . . . . . . . . . . . . . . . . 35
Ángulo entre dos vectores . . . . . . . . 210
Aplicación lineal P
caracterización . . . . . . . . . . . . . . 161 Propiedades
imagen . . . . . . . . . . . . . . . . . . . . . . 163 del núcleo e imagen . . . . . . . . . 163
núcleo . . . . . . . . . . . . . . . . . . . . . . .163 determinante de una matriz . . 54
operaciones . . . . . . . . . . . . . . . . . 174 trasposición de matrices . . . . . . 35
B R
Base Rango de una matriz . . . . . . . . . . . . . .26
caracterización . . . . . . . . . . . . . . 111 Regla
del subespacio máximo . . . . . . 241 de Cramer . . . . . . . . . . . . . . . . . . . 74
dual . . . . . . . . . . . . . . . . . . . . . . . . .199 de Sarrus . . . . . . . . . . . . . . . . . . . . .53
C S
Cambio de base Sintaxis y funciones con Mathematica
en aplicaciones lineales . . . . . . 173 274
en espacios vectoriales . . . . . . . 121 Sistema generador del subespacio
imagen . . . . . . . . . . . . . . . . . 163
D Sistemas de Cramer . . . . . . . . . . . . . . . 74
Determinante de una matriz . . . . . . 54 Subespacio
Diagonalización por semejanza caracterización . . . . . . . . . . . . . . 127
ortogonal . . . . . . . . . . . . . . . 238 generado por un conjunto de
vectores . . . . . . . . . . . . . . . . 132
E suplementario . . . . . . . . . . . . . . . 152
Subespacios propios generalizados.
Endomorsmo diagonalizable . . . . 229 Subespacio máximo . . . . . 240
F T
Fórmula Teorema
de Grassmann . . . . . . . . . . . . . . . 151 de Cayley-Hamilton . . . . . . . . . 224
de las dimensiones . . . . . . . . . . .168 de Rouché-Frobenius . . . . . . . . . 47
277 Traza de una matriz cuadrada . . . 216
Bibliografía
[6] D. C. Lay , Álgebra Lineal y sus aplicaciones (3a Edición.), PRENTICE HALL
MEXICO, 2007.
279