https://ocw.mit.edu/courses/18-06-linear-algebra-spring-2010/video_galleries/video-lectures/
이번 강의에서는 다음 내용을 다룬다.
- "Diagonalization a matrix S−1AS=Λ"
- "Powers of A"
- "equation uk+1=Auk"
우선 지난 강의에서 eigenvalues and eigenvectors를 배웠다.
- Ax=λx
- A−λI:singular (det(A−λI)=0)
먼저, 행렬 A의 n개의 independent한 egienvectors를 S라는 행렬의 칼럼으로 놔보자.
S=⎣⎢⎡∣x1∣∣x2∣.........∣xn∣⎦⎥⎤
그리고 다음과 같이 AS를 구해보자.
AS=A⎣⎢⎡∣x1∣∣x2∣.........∣xn∣⎦⎥⎤=⎣⎢⎡∣λ1x1∣∣λ2x2∣.........∣λnxn∣⎦⎥⎤
=⎣⎢⎡∣x1∣∣x2∣.........∣xn∣⎦⎥⎤diagonal eigenvalues matrix⎣⎢⎢⎢⎡λ10.00λ2..............0..λn⎦⎥⎥⎥⎤=SΛ
위 과정에서 AS=SΛ 라는 것을 알았다.
그러면 위로부터 다음과 같은 수식도 알 수 있을 것이다.
S−1AS=Λ
A=SΛS−1
이제 위 수식을 이용해서 더 알아보자.
만약 Ax=λx인 경우에, Ak은 어떻게 나올까? 다음 식을 보자.
A2x=λAx=λ2x
이전에 배웠다시피 A가 변해도 eigenvector에는 변화가 없고 eigenvalue만 변한다.
이어서 다음을 보자.
A=SΛS−1
A2=SΛS−1SΛS−1=SΛ2S−1
Ak=SΛkS−1
마찬가지로 A=SΛS−1 수식에서 A가 변해도 eigenvectors의 행렬은 변화가 없고 빅 람다 Λ 값만 변한다.
따라서 위로부터 다음과 같은 추정이 가능하다.
If all ∣λi∣<1,then Ak→0 when k→∞
왜냐하면 Λ<1이 되어, Λ∞→0가 될 것이기 때문이다. 따라서 Ak=SΛkS−1에 의해 Ak→0가 된다.
그리고 다음과 같은 추정도 가능할 것이다.
If all the λ′s are different, (no repeated λ′s)A is sure to have n independent eigenvectors. (and be diagonalizable)
이제 uk+1=Auk에 대해서 알아보자. uo부터 시작하여 아래 수식을 정리해보자.
uk+1=Auk
u2=Au1=AAu0=A2u0
uk=Aku0
즉 uk+1=Auk와 같은 규칙이 존재할 때, uk=Aku0 라는 것을 알 수 있다.
계속 이어서 보자.
u0=c1x1+c2x2+...+cnxn=cS=Sc (S=⎣⎢⎡∣x1∣∣x2∣.........∣xn∣⎦⎥⎤)
Au0=c1λ1x1+c2λ2x2+...+cnλnxn (Aixi=λixi)
A100u0=c1λ1100x1+c2λ2100x2+...+cnλn100xn=Λ100Sc
또한 uk+1=Auk와 같은 규칙이 존재할 때, Aku0=ΛkSc 라는 것을 알 수 있다.
한번 예시를 적용해보자.
다음은 Fibonacci 수열의 예시이다.
F0=0F1=1F2=1F3=2F4=3F5=5F6=8F7=13...Fk+2=Fk+1+Fk
여기에 약간의 변형을 가하여 uk+1=Auk 규칙을 적용해보자.
Fk+2=Fk+1+Fk
Fk+1=Fk+1
uk=[Fk+1Fk]
uk+1=A[1110]uk
위와 같이 변형을 통해 uk+1=Auk 형태로 만들었다.
그럼 계속해서 이어가 보자.
A=[1110]
∣A−λI∣=∣∣∣∣∣1−λ11−λ∣∣∣∣∣=λ2−λ−1=0
∴λ=21±5
λ1=21+5≈1.6...>1
λ2=21−5≈−0.6...<1
따라서 피보나치 수열 값이 커지면 λ1만 남을 것이다. (λ2는 1보다 작으므로 k→∞면 Ak→0이 됨.)
그러므로 F100은 다음과 같다.
F100≈c1(21+5)100+0
이제 x1,x2도 이어서 구해보자.
(A−λI)x=0
[1−λ11−λ][λ1]=[00]
∴x1=[λ11],x2=[λ21]
이제 u0=c1x1+c2x2에 x1,x2를 적용하면,
u0=[F1F0]=[10]
u0=c1x1+c2x2
c1,c2 값을 구할 수 있을 것이다.
따라서 피보나치 수열을 1부터 시작하지 않아도 바로 F100≈c1(21+5)100로 값을 구할 수가 있다.