[선형대수학] 19. 조르단 표준형
앞에서 다루었던 행렬B=(t100t100t)(t∈R)은 대각화가 가능하지 않은 행렬이다. 이러한 행렬은 다음의 조르단 표준형이라고 하는 형태이다.
n×n행렬 A가 s개의 선형독립인 고유벡터들을 가지면, 조르단 표준형(Jordan canonical form)이라고 불리는 다음의 행렬과 닮음이다.J=Q−1AQ=(J1J2⋱Js)여기서 Ji는 조르단 블럭(Jordan block)으로 다음과 같은 삼각행렬이고Ji=(λi1⋱⋱⋱1λi)λi는 행렬 A의 단일 고유값, s는 A의 선형독립인 고유벡터의 수이다.
주의할 점은 선형독립인 고유벡터의 수가 1개 이상이면, 동일한 고유값 λi는 여러 블럭에 나타날 수 있다. A가 n개의 선형독립인 고유벡터를 가지면, n개의 조르단 블럭을 갖고, 그 조르단 블럭은 1×1행렬이므로 A는 대각행렬이고, 대각행렬은 조르단 표준형의 특별한 경우이다.
5×5행렬 A가 중복도가 5인 하나의 고유값 λ를 갖는다고 하자. 그러면 가능한 조르단 표준형의 개수는 다음의 7개이다.
(1) A가 하나의 선형독립인 고유벡터를 갖는다고 하면 가능한 A의 조르단 표준형은 다음과 같다.J=(λ10000λ10000λ10000λ10000λ)이 조르단 표준형은 1개의 조르단 블럭으로 되어있고, 대각선에 고유값 λ가 있다.
(2) A가 두 개의 선형독립인 고유벡터를 갖는다고 하면 가능한 A의 조르단 표준형은 다음과 같다.J=(λ10λλ100λ100λ),J=(λλ1000λ1000λ1000λ)이 조르단 표준형은 2개의 조르단 블럭으로 되어있고, 대각선에 고유값 λ가 있다.
(3) A가 세 개의 선형독립인 고유벡터를 갖는다고 하면 가능한 A의 조르단 표준형은 다음과 같다.J=(λλ10λλ10λ),J=(λλλ100λ100λ)이 조르단 표준형은 3개의 조르단 블럭으로 되어있고, 대각선에 고유값 λ가 있다.
(4) A가 네 개의 선형독립인 고유벡터를 갖는다고 하면 가능한 A의 조르단 표준형은 다음과 같다.J=(λλλλ10λ)이 조르단 표준형은 4개의 조르단 블럭으로 되어있고, 대각선에 고유값 λ가 있다.
(5) A가 다섯 개의 선형독립인 고유벡터를 갖는다고 하면 가능한 A의 조르단 표준형은 다음과 같다.J=(λλλλλ)대각행렬과 동일하다(정확히는 5개의 조르단 블럭).
행렬 J가 다음과 같은 조르단 표준형이라고 하자.J=((6106)(0100)(0))=(J1J2J3)이 행렬의 고유값은 6과 0이고 각각 중복도가 2와 3이다. 고유값 λ=6에 대한 고유벡터를 x1(계산해보면 1개가 나온다), 고유값 λ=0에 대한 고유벡터를 x3, x5(계산해보면 2개가 나온다)라고 하면 다음과 같다.x1=(10000),x3=(00100),x5=(00001)Q=[x1x2⋯x5]라 하고 Q−1AQ=J를 AQ=QJ로 나타내면A(x1x2⋯x5)=(x1x2⋯x5)(610601000)이므로 따라서[Ax1Ax2Ax3Ax4Ax5]=[6x1x1+6x20x3x30x5]orAx1=6x1,Ax2=6x2+x1,Ax3=0x3,Ax4=0x4+x3,Ax5=0x5여기서 x1, x3, x5는 A의 고유벡터이지만 x2, x4는 아니다. 그러나 이 두 벡터는 다음의 식이 성립하고(A−6I)x2=x1,(A−0I)x4=x3따라서 다음이 성립한다.(A−6I)2x2=0,(A−0I)2x4=0여기서의 두 벡터 x2와 x4를 각각 고유값 λ=6과 λ=0에 대한 일반 고유벡터라고 한다. 참고로 다음이 성립한다.(A−6I)x2≠0,(A−0I)x4≠0x가 계수(rank)가 k이고 고유값이 λ인 A의 일반 고유벡터(generalized eigenvector)라는 것은 다음이 성립하는 것이다.(A−λI)kx=0,(A−λI)k−1x≠0k=1이면 고유벡터의 정의이다. x를 고유값이 λ이고 계수가 k인 일반 고유벡터라고 하고 x1,x2,...,xk를 다음과 같이 정의하자xk=xxk−1=(A−λI)x=(A−λI)xkxk−2=(A−λI)2x=(A−λI)xk−1⋮x2=(A−λI)k−2x=(A−λI)x3x1=(A−λI)k−1x=(A−λI)x2위의 벡터들의 집합 {x1,x2,...,xk}를 고유값이 λ인 일반 고유벡터들의 사슬(chain of generalized eigenvectors)이라고 한다.
x가 고유값이 λ이고 계수가 k(>1)인 A의 일반 고유벡터이면, 1<l≤k인 l에 대하여 (A−λI)lxl=(A−λI)kx=0이고 (A−λI)l−1xl=(A−λI)k−1x≠0이다. 따라서 벡터 xl=(A−λI)k−lx는 고유값이 λ이고 차원이 l인 일반 고유벡터이다. x1=(A−λI)k−1x는 항상 고유값이 λ인 고유벡터이므로 이 사슬의 초기벡터(initial vector)라고 한다. l≥i에 대하여 (A−λI)lxi=0이다.
다음은 일반 고유벡터의 사슬에 관한 성질인데 증명은 하지 않겠다.
1. 고유값이 λ인 일반 고유벡터들의 사슬 S={x1,x2,...,xk}는 선형독립이다.
2. 정방행렬 A의 서로 다른 고유값들에 대한 일반 고유벡터의 사슬의 합집합은 선형독립이다.
3. S={x1,x2,...,xk}와 T={y1,y2,...,yl}을 고유값이 λ인 행렬 A의 일반 고유벡터들의 사슬이라고 하자. 초기벡터 x1과 y1이 서로 독립이면, S∪T는 선형독립이다.
다음의 행렬에 대한 조르단 표준형을 구하자.A=(21402−1003)행렬 A의 고유값은 λ1=λ2=2, λ3=3이고 고유값 2에 대한 고유벡터를 u1, 고유값 3에 대한 고유벡터를 u3이라고 하면 다음이 성립한다.u1=(100),u3=(3−11)E(2)=2(고유값 2에 대한 고유공간의 차원이 2)이므로 u2는 고유값 2에 대한 일반 고유벡터이다.(A−2I)x=(01400−1001)u2≠0(A−2I)2=(00300−1001)u2=0이므로u2=(ab0)의 형태이나u1=(100)이므로u2=(010)이어야 한다. 그러면 다음이 성립한다.(A−2I)u2=u1(A−2I)2u2=(A−2I)u1=0또한 일반 고유벡터들의 사슬 {u1,u2,u3}은 선형독립이므로 이 벡터들로 행렬 Q와 그 역행렬 Q−1을 다음과 같이 구할 수 있다.Q=(10301−1001),Q−1=(10−3011001)그러면 다음과 같이 조르단 표준형을 구할 수 있다.Q−1AQ=(210020003)=(J100J2)(J1=(2102),J2=(3))
다음의 행렬에 대해 조르단 표준형을 구하자.A=(010000100001−14−64)고유값을 구하면det(A−λI)=λ4−4λ3+6λ2−4λ+1=(λ−1)4이므로 A의 고유값은 λ=1이고 중복도가 4이다.A−I=(−11000−11000−11−14−63)이고 rank(A−I)=3이므로 E(1)=1이다. 고유벡터를 구하면x=(1111)이므로 계수가 4인 일반 고유벡터를 찾아야 한다. 직접 계산해보면 (A−I)4=O(영행렬)이므로 식 (A−I)4x=0을 만족하고 계수가 4인 일반 고유벡터 x를 임의로 다음과 같이 설정할 수 있다.x=(−1000)그러면x3=(A−I)x4=(−11000−11000−11−14−63)(−1000)=(1001)x2=(A−I)x3=(−1012),x1=(A−I)x2=(1111)이고 일반 고유벡터들의 사슬 {x1,x2,x3,x4}은 선형독립이므로Q=(1−11−1100011001210),Q−1=(01000−11001−21−13−31)이고Q−1AQ=(1100011000110001)=J이 구하려는 조르단 표준형이다.
참고자료:
Linear Algebra, jinho Kwak, sungpyo Hong, Birkhauser
'대수학 > 선형대수학' 카테고리의 다른 글
[선형대수학] 21. 지수행렬의 계산 (0) | 2020.02.26 |
---|---|
[선형대수학] 20. 지수행렬 (0) | 2020.02.25 |
[선형대수학] 18. 행렬의 대각화 (0) | 2018.12.13 |
[선형대수학] 17. 고유값과 고유벡터 (0) | 2018.12.12 |
[선형대수학] 16. 내적공간 (4: 그람-슈미츠 직교화) (0) | 2017.11.06 |