고유값과 고유 벡터 - Eigenvalue and Eigenvector
Definition - Eigenvalue
n×n 행렬 A가 있다. 0이 아닌 열 벡터 v∈Rn에 대해 Av=λv를 만족하는 스칼라 λ∈R를 고윳값(eigenvalue)이라고 한다.
Eigenvalue는 Characteristic value라고도 부른다
Definition - Eigenvector
Av=λv를 만족하는 영이 아닌 열 벡터 v∈Rn를 고윳값 λ에 대응하는 A의 고유 벡터(eigenvector)라고 한다.
Eigenvector는 Characteristic vector라고도 부른다.
영 벡터는 결코 고유 벡터가 될 수 없다.
예시 1.
A=(2003)이 있다.
(2003)(10)=(20)=2(10) 이므로,
(1,0)T은 고윳값 2에 대응하는 A의 고유 벡터이다. (λ=2, v=(10))
(2003)(0−2)=(0−6)=3(0−2) 이므로,
(0,−2)T은 고윳값 3에 대응하는 A의 고유 벡터이다. (λ=3, v=(0−2))
예시 2.
A=(0110)이 있다.
(0110)(11)=(11),
(0110)(1−1)=(−11)=(−1)(1−1) 이므로,
(1,1)T는 고윳값 1에 대응하는 A의 고유 벡터이고, (1,−1)T는 고윳값 −1에 대응하는 A의 고유 벡터이다.
벡터 v1=(1,1)T, v2=(1,−1)T는 R2의 기저를 형성한다.
L(x)=Ax로 주어지는 선형 연산자 L:R2→R2를 생각하라.
기저 v1,v2에 대한 L의 행렬 B는 다음과 같다.
B=(100−1)
왜 그런지 기억이 안난다면 L(v1),L(v2)를 생각해보자.
Theorem
n×n 행렬 A가 있다. L(x)=Ax로 주어지는 선형 연산자 L:Rn→Rn를 생각하라.
v1,v2,...,vn는 Rn의 표준 기저가 아닌 기저이고 이 기저에 대한 L의 행렬을 B라고 하자.
벡터 v1,v2,...,vn가 A의 고유 벡터라면 행렬 B는 대각 행렬이다.
또한 행렬 B가 대각 행렬이라면 벡터 v1,v2,...,vn는 A의 고유 벡터이다.
이 경우 B의 주대각 성분은 A의 고윳값과 같다.
고유 공간 - Eigenspace
n×n 행렬 A와 고윳값 λ, 고유 벡터 v가 있다고 하자.
Av=λv이므로, Av=(λI)v⇒(A−λI)v=0 이다.
λ가 A의 고윳값이라면 (A−λI)v=0은 0이 아닌 해(nontrivial value)를 가진다.
그러므로 벡터 v는 행렬 (A−λI)의 영 공간이다.
v∈N(A−λI)
즉, λ가 A의 고윳값이라면 N(A−λI)={0}이고, N(A−λI) 안의 영이 아닌 벡터는 λ에 대응하는 고유 벡터이다.
반대로, 만약 벡터 x가 N(A−λI) 안의 벡터라면 Ax=λx이다.
그러므로 고윳값 λ에 대응하는 A의 고유 벡터들은 N(A−λI)의 영이 아닌 벡터들이다.
Definition
N(A−λI)={0}이면, N(A−λI)는 고윳값 λ에 대응하는 행렬 A의 고유 공간(eigenspace)이다.
고윳값과 고유벡터 구하기
정방행렬 A와 스칼라 λ가 주어졌을 때 다음 문장들은 동등한 의미를 가진다.
- λ는 A의 고윳값이다.
- N(A−λI)={0}
- 행렬 (A−λI)는 특이(singular) 행렬이다.
- det(A−λI)=0
Definition
det(A−λI)=0은 행렬 A의 특성 방정식(characteristic equation)이다.
A의 고윳값 λ들은 특성 방정식의 근이다.
A의 고유벡터들은 방정식 (A−λI)x=0의 0이 아닌 해다.
고윳값과 관련된 고유벡터를 구하는 절차
p(λ)=det(λIn−A)
이는 A의 고윳값들이다.
- Step 2. 각 고윳값 λ마다, 동차 방정식 (λIn−A)x=0의 0이 아닌 모든 해를 찾는다.
이는 고윳값 λ에 대응하는 A의 고유벡터이다.
예시 1.
A=⎝⎜⎛a11a12a13a21a22a23a31a32a33⎠⎟⎞
det(A−λI)=∣∣∣∣∣∣∣a11−λa21a31a12a22−λa32a13a23a33−λ∣∣∣∣∣∣∣=−λ3+c1λ2−c2λ+c3,
c1=a11+a22+a33 (trace of A) ,
c2=∣∣∣∣∣a11a21a12a22∣∣∣∣∣+∣∣∣∣∣a11a31a13a32∣∣∣∣∣+∣∣∣∣∣a22a32a23a33∣∣∣∣∣,
c3=detA
예시 2.
A=(2112)
Characteristic equation :
det(A−λI)=∣∣∣∣∣2−λ112−λ∣∣∣∣∣=0
(2−λ)2−1=0⇒λ1=1, λ2=3
(A−λ1I)x=(A−I)x=0⇔(1111)(xy)=(00)
⇔(1010)(xy)=(00)⇔x+y=0
위 방정식의 일반 해는 (−t,t)=t(−1,1), t∈R이다.
그러므로 벡터 v1=(−1,1)은 고윳값 1에 대한 고유벡터이다.
이에 대응하는 고유공간은 v1에 의해 생성되는 선이다.
(A−λ2I)x=(A−3I)x=0⇔(−111−1)(xy)=(00)
⇔(10−10)(xy)=(00)⇔x−y=0
위 방정식의 일반 해는 (t,t)=t(1,1), t∈R이다.
그러므로 벡터 v2=(1,1)은 고윳값 3에 대한 고유벡터이다.
이에 대응하는 고유공간은 v2에 의해 생성되는 선이다.
특성 다항식 - Characteristic polynomial
다항식 p(λ)=det(A−λI)를 행렬 A의 특성 다항식(characteristic polynomial)이라고 한다.
det(A−λI)=(−1)nλn+c1λn−1+⋯+cn
(−1)n−1c1=a11+a22+⋯+ann=trace(A),
cn=detA
p(λ)=det(A−λI)=∣∣∣∣∣∣∣∣∣∣a11−λa21⋮an1a12a22−λ⋮an2⋯⋯⋮⋯a1na2n⋮ann−λ∣∣∣∣∣∣∣∣∣∣=(a11−λ)(a22−λ)⋯(ann−λ)=(−1)nλn+c1λn−1+⋯+cn
이 식은 다음과 같이 인수분해 될 수 있다.
p(λ)=(−1)n(λ−λ1)(λ−λ2)⋯(λ−λn)=(λ−λ1)(λ−λ2)⋯(λn−λ)
여기서 우리는 다음과 같은 사실을 볼 수 있다.
p(0)=λ1⋅λ2⋯λn=det(A)
∑i=1nλi=∑i=1naii=tr(A)
즉, 고윳값들의 곱은 A의 행렬식과 같고 고윳값들의 합은 A의 대각합(trace)과 같다.
고윳값의 합과 곱 예시
A=(51−18−1)
det(A)=−5+18=13, tr(A)=5−1=4
p(λ)=∣∣∣∣∣5−λ1−18−1−λ∣∣∣∣∣=λ2−4λ+13
λ1=2+3i, λ2=2−3i
λ1+λ2=4=tr(A)
λ1⋅λ2=13=det(A)
연산자의 고윳값과 고유벡터
선형 연산자 L:V→V가 0이 아닌 벡터 v∈V에 대해서 L(v)=λv일 때,
λ는 연산자 L의 고윳값이고 v는 그에 대한 고유벡터이다.
만약 V가 함수 공간이라면 고유벡터는 고유 함수(eigenfunction)이라고도 불린다.
L(x)=Ax라면(A는 정방행렬) L의 고윳값과 고유벡터는 A의 고윳값과 고유벡터이다.
명제
V는 벡터 공간이고 L:V→V는 선형 연산자라고 하자.
- 벡터 v∈V가 L의 고유벡터라면 이와 관련된 고윳값은 유일하다.
- 벡터 v1, v2가 서로 다른 고윳값 λ1, λ2에 대응하는 고유벡터라면 v1과 v2는 선형 독립이다.
⇒ Theorem
벡터 v1,v2,⋯,vk가 고윳값 λ1,λ2,⋯,λk에 대응하는 선형 연산자 L의 고유벡터라면 v1,v2,⋯,vk는 선형 독립이다.
행렬의 닮음 - Similarity
정방 행렬 A와 B가 있다.
B가 A의 유사행렬이라면 두 행렬은 같은 특성 다항식과 고윳값을 가진다.
<참고자료>
Steve J. Leon, Linear Algebra with Application, 10th edition, 2021.