요르단 표준형의 장점
고윳값, 고유 벡터의 모양이 보입니다.
거듭제곱을 구체적으로 계산할 수 있습니다.
이 두 가지는 폭주 판정과 관련이 있습니다.
요르단 표준형은 블록대각이므로 고윳값이나 거듭제곱 계산은 블록별로 보면 알 수 있습니다.
다음 요르단 셀을 이용하여 학습하겠습니다.
B=(7100071000710007)
요르단 표준형의 고윳값
요르단 셀 B의 고윳값은 7밖에 없습니다. (상삼각 행렬)
p=(α,0,0,0)T가 고윳값 7의 고유 벡터입니다. (α≠0) 이 외의 고유 벡터는 없습니다.
요르단 표준형 J의 고윳값, 고유벡터는 한눈에 알 수 있습니다.
대각 성분이 고윳값 λ
대각 성분의 λ 개수가 고윳값 λ가 몇 중 해인지(대수적 중복도)에 대응
대각 성분이 λ인 요르단 셀의 개수가 고윳값 λ에 선형 독립인 고유 벡터 개수(기하적 중복도)에 대응
위의 성질에서 A가 요르단 표준형으로 변환되면 A의 고윳값, 고유 벡터가 어떻게 되어 있는지 알 수 있습니다.
특히 고윳값에 중해가 없을 경우 요르단 표준형은 대각 행렬이 될 수밖에 없습니다. 즉,
고윳값에 중해가 없을 경우 대각 화할 수 있습니다.
요르단 표준형의 거듭제곱
B=7I+Z
Z=(0100001000010000)
위와 같이 분해해 두면 계산할 때 예측하기 쉽습니다. 특징은 Z가 왼쪽에 곱하면 1행을 밀고 오른쪽에 곱하면 1열 밀기 작용을 합니다.
(0100001000010000)(abcd)=(bcd0)
(aeimbfjncgkodhlp)(0100001000010000)=(0aei0bfj0cgk0dhl)
모양입니다.
이 성질을 이용하면 Z의 거듭제곱은 간단합니다.
Z2=(0010000100000000), Z3=(0001000000000000), Z4=Z5=⋯=O
위와 같은 식으로 2 제곱, 3 제곱...으로 늘릴 때마다 1의 위치가 오른쪽 위로 밀리는 것을 알 수 있습니다.
이제 B2=(7I+Z)2=72I+2⋅7Z+Z2에서
B2=(722⋅7100722⋅7100722⋅700072)
B3=(7I+Z)3=73I+3⋅72Z+3⋅7Z2+Z3에서
B3=(733⋅723⋅710733⋅723⋅700733⋅7200073)
B4=(7I+Z)4=74I+4⋅73Z+6⋅72Z2+4⋅7Z3+Z4 (Z4=O)
B4=(744⋅736⋅724⋅70744⋅736⋅7200744⋅7300074)
더 큰 t 제곱에서도
Bt=(7I+Z)t=7tI+t⋅7t−1Z+t2C⋅7t−2Z2+t3C⋅7t−3Z3+t4C⋅7t−4Z4+⋯+t−2tC⋅72Zt−2+t⋅7Zt−1+Zt에서 (Z4=Z5=Z6=⋯=O)
Bt=(7tt⋅7t−1t2C⋅7t−2t3C⋅7t−307tt⋅7t−1t2C⋅7t−2007tt⋅7t−10007t)
도 구해집니다(t = 1, 2,...).
대각 성분이 λ인 요르단 셀 B에 대해 Bt를 계산하면 stCλt−s 항이 나타납니다(단, s > t인 경우에는 stCλt−s=0으로 간주하고 λ0=1이라 간주합니다).
다음과 같이도 말할 수 있습니다.
f(λ)=λt로 두고, f를 λ로 s회 미분한 식을
f(s)(λ)=dsdλsf(λ)
로 두면 크기 m의 요르단 셀
B=(λ1⋱⋱⋱1λ)
의 t제곱은
(f(λ)f(1)(λ)12f(2)(λ)13!f(3)(λ)⋯1(m−1)!f(m−1)(λ)f(λ)f(1)(λ)12f(2)(λ)⋯1(m−2)!f(m−2)(λ)⋱⋱⋱⋮⋱⋱12f(2)(λ)⋱f(1)(λ)f(λ)) (4.21)
이유는 fs(λ)=t(t−1)⋯(t−s+1)λt−s이기 때문입니다.
1s!fs(λ)=tsCλt−s
식 (4.20)의 J라면 각 요르단 셀
B1=(3100003100003100003100003), B2=(3103), B3=(4), B4=(510051005)
에 대해 각각 t제곱을 구해두면
Jt=(Bt1OOOOBt2OOOOBt3OOOOBt4)
로 구해집니다.
요르단 표준형이 아닌 정방 행렬 A의 거듭제곱 At는 계산 가능할까요?
A를 요르단 표준형 J로 변환하면 전과 같이 At를 계산할 수 있습니다.
P−1AP=J라는 것은 A=PJP−1이고, 그 t 제곱은
At=(PJP−1)t=PJtP−1입니다.
내용이 결고 쉽지 않습니다. 중간중간 막힌 부분은 다시 복습해야 합니다.
2019/09/04 - [선형대수/고윳값, 대각화, 요르단 표준형 - 폭주 위험] - 4.5.2 고윳값, 고유 벡터의 성질
감사합니다.
'프로그래머를 위한 선형대수 > 고윳값, 대각화, 요르단 표준형 - 폭주 위험' 카테고리의 다른 글
4.7.2 대각까지는 못하더라도 - 요르단 표준형 (0) | 2019.09.24 |
---|---|
4.7 대각화할 수 없는 경우 (0) | 2019.09.20 |
4.6.4 대각화할 수 있는 경우, 4.6.4 결론: 고윳값(실수부)의 부호 (0) | 2019.09.19 |
4.6.3 대각행렬일 때 (0) | 2019.09.18 |
4.6.2 1차원일 때 (0) | 2019.09.17 |
댓글