Álgebra Lineal

Curso 2014/15
Universidade de Vigo
Departamento de Matemática Aplicada II 13 de enero de 2015
E. T. S. E. de Minas

Soluciones a los ejercicios del examen final

1 Se consideran    
1 1 1 1
A =  1 −1 1  , b = 0 .
−1 1 −1 1
a) Usar operaciones elementales por filas para calcular una matriz inversible F tal
que T = F A sea triangular superior.
b) Probar que el sistema de ecuaciones lineales Ax = b es incompatible.
c) Calcular una solución del sistema Ax = b en el sentido de mı́nimos cuadrados
cuya primera coordenada sea cero.

Solución:
a) Realizando operaciones elementales sobre las filas de la matriz A, tenemos:
     
1 1 1 F21 (−1) 1 1 1 1 11
F32 (1)
A =  1 −1 1  −→  0 −2 0  −→  0 −2 0  = T.
−1 1 −1 F31 (1) 0 20 0 00
La matriz triangular T se obtiene como T = F A, donde F = F32 (1)F31 (1)F21 (−1). La
matriz F resulta de hacer las mismas operaciones elementales a la matriz identidad:
     
1 0 0 F21 (−1) 100 100
F32 (1)
I =  0 1 0  −→  −1 1 0  −→  −1 1 0  = F.
001 F31 (1) 101 011

b) El sistema Ax = b es equivalente a F Ax = F b, donde

      
1 11 100 1 1
F A = T =  0 −2 0  , F b = b0 =  −1 1 0   0  =  −1  .
0 00 011 1 1

Es claro que rg(T ) = 2 6= rg(T |b0 ) = 3 y por tanto el sistema es incompatible.

c) Las soluciones de Ax = b en el sentido de mı́nimos cuadrados son las soluciones del

sistema At Ax = At b, que resulta
    
3 −1 3 x 0 (
3x − y + 3z = 0
 −1 3 −1   y  =  2  ⇐⇒
3 −1 3 z 0 −x + 3y − z = 2
 Como buscamos la solución con x = 0, se tiene:

x=0

 
x = 0


  
 
 
 3
−y + 3z = 0 ⇐⇒ y = 4
  
z = 1
  
3y − z = 2
  

4
La solución buscada es (0, 3/4, 1/4).

2) Se considera la forma cuadrática ω : R3 → R definida por

ω(x, y, z) = 2x2 + ay 2 + 2z 2 + 2axz,

donde a es un parámetro real.

a) Calcular la matriz M asociada a ω.
b) Hallar el espectro de M en función de los valores de a.
c) Determinar los valores de a para los que se cumplen cada una de las siguientes
propiedades:
c.1) ω es no degenerada.
c.2) ω es definida positiva.
c.3) ω es semidefinida positiva.
d ) Para a = 1, calcular la mejor aproximación de rango 1 de M usando su descom-
posición espectral.

Solución:
a) La expresión matricial de la forma cuadrática es
  
20a x
ω(x, y, z) = (x, y, z)  0 a 0   y  = v t M v.
a02 z

b) El polinomio caracterı́stico de M es

2−x 0 a
2−x a F12 (1) (2 + a − x) (2 + a − x)
qM (x) = 0 a−x 0 = (a − x) = (a − x) =
a 2−x a 2−x
a 0 2−x

1 1
= (a − x)(2 + a − x) = (a − x)(2 + a − x)(2 − a − x).
a 2−x

Por tanto, Sp(M ) = {a, 2 + a, 2 − a}.

c.1) Como |M | = a(4 − a2 ) = a(2 − a)(2 + a):

ω es no degenerada ⇐⇒ a(2 + a)(2 − a) 6= 0 ⇐⇒ a 6= −2, 0, 2.

c.2) Los menores principales de M son ∆1 = 2 > 0, ∆2 = 2a, ∆3 = |M | = a(4 − a2 ). Por

tanto:
( ) ( )
2a > 0 a>0
ω es definida positiva ⇐⇒ ⇐⇒ ⇐⇒ 0 < a < 2.
a(4 − a2 ) > 0 a2 < 4

c.3) La forma cuadrática ω sólo puede ser semidefinida positiva en un caso degenerado, es
decir, para a ∈ {−2, 0, 2}.
Para a = −2, Sp(M ) = {−2, 0, 4} y por tanto ω es indefinida.
Para a = 0, Sp(M ) = {0, 2, 2} y por tanto ω es semidefinida positiva.
Para a = 2, Sp(M ) = {2, 4, 0} y por tanto ω es semidefinida positiva.
En resumen, ω es semidefinida positiva ⇐⇒ a = 0 o a = 2.

d) Para a = 1,  
201
M = 0 1 0 , Sp(M ) = {3, 1, 1}.
102
Como M es simétrica, es ortogonalmente diagonalizable y admite una descomposición
espectral
M = 3u1 ut1 + u2 ut2 + u3 ut3 .
Como todos los autovalores de M son positivos, su aproximación de rango 1 es

M1 = 3u1 ut1 ,

donde {u1 } es una base ortonormal de V (3) = Ker(M − 3I).

 
−1 0 1  
3 z =x
V (3) = Ker  0 −2 0 = (x, y, z) ∈ R /
 =< {(1, 0, 1)} > .
y=0
1 0 −1

Por tanto,
 √   √  3 3
1/ 2 1/ 2  √ √  2 0 2
u1 =  0√  =⇒ M1 = 3  0√  1/ 2 , 0 , 1/ 2 =  0 0 0  .
3 3
1/ 2 1/ 2 2 0 2
3) Se considera el subespacio vectorial

U = {X ∈ M2×2 (R) / u y v son autovectores de X} ,

donde u = (1, 1), v = (1, −1).

a) Hallar una base de U .
b) Calcular las matrices X ∈ U tales que Sp(X) = {2, 4}.

Se proponen dos soluciones:

Primera solución:
 
xy
a) La condición para que una matriz X = pertenezca al subespacio U es que existan
z t
valores λ, µ ∈ R tales que
        
xy 1 1 x+y λ
=λ ⇐⇒ = ⇐⇒ x + y = z + t ⇐⇒ x + y − z − t = 0
z t 1 1 z+t λ
        
xy 1 1 x−y µ
=µ ⇐⇒ = ⇐⇒ x − y = −(z − t) ⇐⇒ x − y + z − t = 0.
z t −1 −1 z−t −µ

Por tanto,
        
xy x+y−z−t=0 xy x−t=0 xy t=x
U= / = / = / =
z t x−y+z−t=0 z t y−z =0 z t z=y
      
xy 10 01
= / x, y ∈ R =< , >.
yx 01 10
   
10 01
El conjunto B = , es una base de U .
01 10
 
xy
b) Del apartado anterior obtenemos que los autovalores de una matriz X = ∈ U son
yx
λ = x + y, µ = x − y.
Por tanto, hay dos posibilidades para que Sp(X) = {2, 4}:
     
x+y =2 x=3 3 −1
λ = 2, µ = 4 =⇒ ⇐⇒ ⇐⇒ X = .
x−y =4 y = −1 −1 3
     
x+y =4 x=3 31
λ = 4, µ = 2 =⇒ ⇐⇒ ⇐⇒ X = .
x−y =2 y=1 13

Segunda solución:
a) Si (1, 1) y (1, −1) son autovectores de X entonces X es diagonalizable y existen dos
autovalores (no necesariamente distintos) λ, µ tales que X = P DP −1 , con
   
1 1 λ0
P = , D= .
1 −1 0µ
 Por tanto,  
−1 1 λ+µ λ−µ
X ∈ U ⇐⇒ X = P DP = .
2 λ−µ λ+µ
Si denotamos x = λ/2, y = µ/2, se tiene:
        
x+y x−y 11 1 −1
U= / x, y ∈ R = x +y / x, y ∈ R =
x−y x+y 11 −1 1
   
11 1 −1
=< , >.
11 −1 1
   
11 1 −1
El conjunto B = , es una base de U .
11 −1 1

b) Del apartado anterior obtenemos que si λ y µ son los autovalores de una matriz X ∈ U
entonces  
1 λ+µ λ−µ
X= .
2 λ−µ λ+µ
Por tanto, hay dos posibilidades para que Sp(X) = {2, 4}:
 
3 −1
λ = 2, µ = 4 =⇒ X = .
−1 3
 
31
λ = 4, µ = 2 =⇒ X = .
13

4) a) Sea A ∈ Mn×n (R) una matriz de rango 1. Probar que det(cos(A)) = cos(tr(A)).
b) Hallar la matriz cos(A) para
 
0 0 0
A = 0 0 0.
πππ

Solución:
a) Como rg(A) = 1, 0 es un autovalor de A y m.g.(0) = n − rg(A) = n − 1. Como la suma de
los autovalores es la traza de A, se tiene que Sp(A) = {0, 0, . . . , 0, tr(A)}. Por otra parte,
| {z }
n−1

Sp(cos(A)) = {cos(0), cos(0), . . . , cos(0), cos(tr(A))} = {1, 1, . . . , 1, cos(tr(A))}.

| {z } | {z }
n−1 n−1

Utilizando que el determinante de una matriz es el producto de sus autovalores, se obtiene:

det(cos(A)) = 1 · 1 · · · 1 · cos(tr(A)) = cos(tr(A)).

b) Definimos la función f (x) = cos(x), de modo que cos(A) = f (A).
Como rg(A) = 1 y tr(A) = π, se obtiene inmediatamente que Sp(A) = {0, 0, π}. El
conjunto de valores de f sobre A es

Vf,A = {f (0), f 0 (0), f (π)} .

Dado que f 0 (x) = − sen(x), se tiene que f (0) = cos(0) = 1, f 0 (0) = − sen(0) = 0,
f (π) = cos(π) = −1.
Tenemos que determinar un polinomio r(x) = a + bx + cx2 tal que


 r(0) = a = f (0) = 1

 
 a=1

 
 

r0 (0) = b = f 0 (0) = 0 ⇐⇒ b = 0
c = − 2

 
 

r(π) = a + bπ + cπ 2 = f (π) = −1
  
π2
2x2
Por tanto r(x) = a + bx + cx2 = 1 − y se tiene:
π2
     
100 0 0 0 1 0 0
2 2  2
cos(A) = f (A) = r(A) = I − 2 A = 0 1 0  − 2  0 0 0  =  0 1 0  .
π π
001 π2 π2 π2 −2 −2 −1