Loading [MathJax]/jax/output/CommonHTML/jax.js
본문 바로가기
프로그래머를 위한 선형대수/고윳값, 대각화, 요르단 표준형 - 폭주 위험

4.7.3 요르단 표준형의 성질

by 지식광부키우기 2019. 10. 1.

 

요르단 표준형의 장점

 

고윳값, 고유 벡터의 모양이 보입니다.

 

거듭제곱을 구체적으로 계산할 수 있습니다.

 

이 두 가지는 폭주 판정과 관련이 있습니다.

 

요르단 표준형은 블록대각이므로 고윳값이나 거듭제곱 계산은 블록별로 보면 알 수 있습니다.

 

다음 요르단 셀을 이용하여 학습하겠습니다.

 

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+27Z+Z2에서

 

B2=(72271007227100722700072)

 

B3=(7I+Z)3=73I+372Z+37Z2+Z3에서

 

B3=(7337237107337237007337200073)

 

B4=(7I+Z)4=74I+473Z+672Z2+47Z3+Z4 (Z4=O)

 

B4=(7447367247074473672007447300074)

 

더 큰 t 제곱에서도

 

Bt=(7I+Z)t=7tI+t7t1Z+t2C7t2Z2+t3C7t3Z3+t4C7t4Z4++t2tC72Zt2+t7Zt1+Zt에서 (Z4=Z5=Z6==O)

 

Bt=(7tt7t1t2C7t2t3C7t307tt7t1t2C7t2007tt7t10007t)

 

도 구해집니다(t = 1, 2,...).

 

대각 성분이 λ인 요르단 셀 B에 대해 Bt를 계산하면 stCλts 항이 나타납니다(단, s > t인 경우에는 stCλts=0으로 간주하고 λ0=1이라 간주합니다).

 

다음과 같이도 말할 수 있습니다.

 

f(λ)=λt로 두고, fλs회 미분한 식을

 

f(s)(λ)=dsdλsf(λ)

 

로 두면 크기 m의 요르단 셀

 

B=(λ11λ)

 

의 t제곱은 

 

(f(λ)f(1)(λ)12f(2)(λ)13!f(3)(λ)1(m1)!f(m1)(λ)f(λ)f(1)(λ)12f(2)(λ)1(m2)!f(m2)(λ)12f(2)(λ)f(1)(λ)f(λ)) (4.21)

 

이유는 fs(λ)=t(t1)(ts+1)λts이기 때문입니다. 

 

1s!fs(λ)=tsCλts

 

식 (4.20)의 J라면 각 요르단 셀

 

B1=(3100003100003100003100003)B2=(3103)B3=(4)B4=(510051005) 

 

에 대해 각각 t제곱을 구해두면 

 

Jt=(Bt1OOOOBt2OOOOBt3OOOOBt4)

 

로 구해집니다.

 

요르단 표준형이 아닌 정방 행렬 A의 거듭제곱 At는 계산 가능할까요? 

 

A를 요르단 표준형 J로 변환하면 전과 같이 At를 계산할 수 있습니다.

 

P1AP=J라는 것은 A=PJP1이고, 그 t 제곱은

 

At=(PJP1)t=PJtP1입니다. 

 

 

내용이 결고 쉽지 않습니다. 중간중간 막힌 부분은 다시 복습해야 합니다.

 

2019/09/04 - [선형대수/고윳값, 대각화, 요르단 표준형 - 폭주 위험] - 4.5.2 고윳값, 고유 벡터의 성질

 

감사합니다.

댓글