https://ocw.mit.edu/courses/18-06-linear-algebra-spring-2010/video_galleries/video-lectures/
이번 강의에서는 다음과 같은 내용을 다룬다.
- Singular Value Decomposition ( SVD )
- A=uΣvT
- u,v : orthogonal
- Σ : diagonal
지난 강의에서 A=QΛQT 를 배운 적이 있다.
- 대칭(Symmetric) 행렬 A 는 항상 직교(Orthogonal) 행렬 Q 와 대각 행렬 Λ 로 분해될 수 있다.
- 만약 A 가 positive definite matrix 라면, Λ 의 모든 대각 성분(고유값)이 양수가 된다.
A=uΣvT 를 확인하기 위해 다음 예시를 보자.
A=[4−343]
v1,v2∈R2(rowspace)
u1,u2∈R2(columnspace)
u,v:orthonormal vector
Av1=σ1u1
Av2=σ2u2
따라서 다음과 같이 표현할 수가 있다.
Av=uΣ (Σ=[σ100σ2])
→A=uΣv−1=uΣvT (v:orthonormal→v−1=vT)
이제 A=uΣvT 식이 어떻게 나왔는지 알았다. 이어서 보자.
ATA=vΣTuTuΣvT=vΣTΣvT=v[σ1200σ22]vT
위에서 ATA 수식을 봤을 때, eigenvalue σ2 와 eigenvector 를 구할 수가 있다.
ATA=[44−33][4−343]=[257725]
- ATA=32[11]
- ATA=18[−11]
위에서 구한 걸 A=uΣvT 에 적용해보자.
A=uΣvT
[4−343]=u[320018][1/2−1/21/21/2]
이제 AAT를 통해 u 를 구해보자.
AAT=uΣvTvΣTuT=uΣΣTuT
AAT=[4−343][44−33]=[320018]=ΣΣT
∴u=[1001]
따라서 최종적으로 A=uΣvT를 증명할 수가 있다.
A=uΣvT
[4−343]=[1001][320018][1/2−1/21/21/2]
새로운 예시를 들어보자.
A=[4836]
보다시피 A 는 singular matrix 이다.
그리고 다음과 같은 해석이 가능하다.
- row space in [43]
- v1=[4/53/5],v2=[−3/54/5]
- column space in [12]
- u1=[1/52/5],u2=[−2/51/5]
이어서 과정을 진행해보자.
ATA=[4386][4836]=[80606045]
ATA는 singular matrix 이므로 λ=0,125 가 나올 것이다.
따라서 σ1=125 임을 알 수 있다.
이제 이어서 진행해보자.
A=uΣvT
[4386]=[1/52/5−2/51/5][125000][4/5−3/53/54/5]
이처럼 singular matrix 에 대해서도 A=uΣvT 를 적용할 수 있다.
즉, SVD 를 최종적으로 정리하면 다음과 같다.
v1,...,vr:orthonormal basis for rowspace
u1,...,ur:orthonormal basis for columnspace
vr+1,...,vn:orthonormal basis for nullspace
ur+1,...,un:orthonormal basis for N(AT)
Avi=σiui