Matemática Aplicada a la Computación

Profesor: Dr. Daniel Alexis Gutierrez Pachas

dgutierrezp@[Link]
Departamento de Ciencia de la Computación - UCSP

23 de octubre de 2024

Ejercicios resueltos sobre SVD

1. Encontrar los valores singulares de las siguientes matrices:
" # " #
1 0 −5 0
(a) A = (b) A = .
0 −3 0 0
√
Solución: Recordar que los valores singulares son σi = λi , donde los λ’s son los autovalores de AT A. En
cada caso, la matriz es simétrica significa que
(a) " #
T 2 1 0 T
A A = A = AA = ⇒ σ1 = 3 , σ2 = 1.
0 9

(b) " #
T 2 1 0 T
√
A A = A = AA = ⇒ σ1 = 5 , σ2 = 0.
0 9

Notar que los valores singulares se ordenan de mayor a menor.

2. Construir el SVD para la matriz: " #
2 −1
A= .
2 2
Solución: Vamos a seguir paso a paso, para representar A = U ΣV T .
(a) Determinando V :
" #
8 2
AT A = , luego P (λ) = λ2 − 13λ + 36 = 0. Entonces λ1 = 9 y λ2 = 4.
2 5
" #
−1 2
• Para λ1 = 9, se tiene (AT A − 9I) = . Luego,
2 −4
" #
n 2 o
T
N (A A − 9I) = Spam .
1
√ #
"
2/√5
Entonces v1 =
1/ 5
" #
4 2
• Para λ2 = 4, se tiene (AT A − 4I) = . Luego,
2 1
" #
n −1 o
T
N (A A − 4I) = Spam .
2
√ #
"
−1/√ 5
Entonces v2 = .
2/ 5

1
 Finalmente √" √ #
2/√5 −1/√ 5
V = .
1/ 5 2/ 5

(b) Determinando Σ:
"√ # " #
λ1 √0 3 0
Σ= =
0 λ2 0 2

(c) Determinando U :
" #" √ # " √ #
1 1 2 −1 −2/√ 5 1/√5
u1 = Av1 = = .
σ1 3 2 2 1/ 5 2/ 5
" #" √ # " √ #
1 1 2 −1 −1/√ 5 −2/√ 5
u2 = Av2 = = .
σ2 2 2 2 2/ 5 1/ 5

Finalmente, √ " √ #
1/√5 −2/√ 5
U= .
2/ 5 1/ 5

3. Construir el SVD para la matriz: " #

3 2 2
A= .
2 3 −2

Solución: Vamos a seguir paso a paso, para representar A = U ΣV T .

(a) Determinando V :
 
13 12 2
Calculando AT A = 12 13 −2, luego P (λ) = λ3 − 34λ2 + 225λ = 0. Entonces λ1 = 25, λ2 = 9
 
2 −2 8
y λ3 = 0.
 
−12 12 2
• Para λ1 = 25, se tiene (AT A − 25I) =  12 −12 −2 . Luego
 
2 −2 −17
 √ 
n 1/√2 o
N (AT A − 25I) = Spam 1/ 2 .
 
0
 √ 
1/√2
Entonces v1 = 1/ 2.
 
0
 
4 12 2
• Para λ2 = 9, se tiene (AT A − 9I) = 12 4 −2. Luego
 
2 −2 −1
√  
1/ √18 on
N (AT A − 9I) = Spam −1/√ 18 .
 
4/ 18
 √ 
1/ √18
Entonces v2 = −1/√ 18.
 
4/ 18

2
  
13 12 2
• Para λ3 = 0, se tiene (A A) = 12
T 13 −2. Luego
 
2 −2 8
 
2/3 o n
N (AT A) = Spam −2/3 .
 
−1/3
 
2/3
Entonces v3 = −2/3.
 
−1/3
Finalmente √ √ 
1/√2 1/ √18 2/3
V = 1/ 2 −1/√ 18 −2/3 .
 
0 4/ 18 −1/3

(b) Determinando Σ:
"√ # " #
λ1 √0 0 5 0 0
Σ= =
0 λ2 0 0 3 0

(c) Determinando U :
 √ 
" # 1/ 2 " √ #
1 1 3 2 2  √  1/√2
u1 = Av1 = 1/ 2 = .
σ1 5 2 3 −2 1/ 2
0
√  "
√ #

" # 1/ √18
1 1 3 2 2  1/ √2
u2 = Av2 = −1/√ 18 = .

σ2 3 2 3 −2 −1/ 2
4/ 18

Finalmente √ " √ #
1/√2 1/ √2
U= .
1/ 2 −1/ 2

4. Construir el SVD para la matriz:  

7 1
A = 0 0 .
 
5 5
Solución: Vamos a seguir paso a paso, para representar A = U ΣV T .
(a) Determinando V :
" #
74 32
AT A = luego P (λ) = λ2 − 100λ + 900 = 0. Entonces λ1 = 90 y λ2 = 10.
2 5
" #
−16 32
• Para λ1 = 90, se tiene (AT A − 90I) = . Luego
32 −64
" #
n 2 o
T
N (A A − 90I) = Spam .
1
√ #
"
2/√5
Entonces v1 = .
1/ 5

3
 " #
64 32
• Para λ = 10, se tiene (AT A − 10I) = . Luego
32 16
" #
n −1 o
N (AT A − 10I) = Spam .
2
"√ #
−1/√ 5
Entonces v2 = .
2/ 5
Finalmente √ " √ #
2/√5 −1/√ 5
V = .
1/ 5 2/ 5

(b) Determinando Σ:
√   √ 
λ1 √0 3 10 √0
Σ= 0 λ2  =  0 10
   
0 0 0 0

(c) Determinando U :
 √ 
7 1 " √ #
 
1/ 2
1 1   2/√5
u1 = Av1 = √ 0 0 =  0√  .
 
σ1 3 10 5 5 1/ 5 1/ 2
√ 
√ #
  
7 1 " −1/ 2
1 1   −1/√ 5
u2 = Av2 = √ 0 0 =  0√  .
 
σ2 10 5 5 2/ 5
1/ 2
En este problema, para el SVD completo, se necesita el autoespacio para AAT , donde λ = 0:
   
50 0 40 0
T
AA =  0 0 0  ⇒ u3 = 1 .
   
40 0 50 0

Finalmente √ √ 
1/ 2 −1/ 2 0
U =  0√ 0√ 1 .
 
1/ 2 1/ 2 0