Loading [MathJax]/jax/output/CommonHTML/jax.js
본문 바로가기

행렬(Matrix) 10. 행렬의 고유값 문제

행렬이 고유값 문제는 다음 형태의 연립방정식과 관련된다.

        Ax=λx        (1)

                A는 정행방렬(2×2,3×3,n×n),    x, λ는 구하려는 미지의 벡터, 스칼라

        x=0, 0=0  λ와 관계없이 식(1)을 만족

        식(1)의 x0인 해, 행렬 A의 고유 벡터

         x0인 해, 어떤 특정한 λ에 대해서만 존재, 이 값을 행렬 A의 고유값

식(1)을 푸는 것은 벡터 x에 행렬 A를 곱하는 것이 벡터 x에 스칼라 λ를 곱하는 것과 같은 영향을 갖는 x를 찾는 것

얻어지는 스칼라 곱 λxx의 성분에 비례하는 성분을 가지며 λ는 비례상수

 

■ 고유값, 고유벡터

        주어진 n×n 행렬 A[ajk]에 대하여 다음의 벡터 방정식을 고려

                Ax=λx,        x : 벡터,    λ : 스칼라

                0벡터 : 모든 x의 해(실용적으로 무의미)

                      x0 인 해,  λ : 고유값, x : 고유벡터

 

■ 고유값과 고유벡터 계산법

    ex 1)

        A=[5222]

    a) 고유값 : 우선적으로 결정되어야 한다.

        Ax=[5222][x1x2]=λ[x1x2]

        성분별로 풀어쓰면

                5x1+2x2=λx1

                2x12x2=λx2

        우변을 좌변으로 옮기면

                (5λ)x1+2x2=0        (2*)

                2x1+(2λ)x2=0

        행렬 표기

                ([5λ222λ]λ[1001])[x1x2]=[00]

                (AλI)x=0

                Crammer 정리에 의해 자명하지 않은 해 x0을 가진 필요충분조건은 행렬식이 0

                D(λ)=det(AλI)=0

                    =|5λ222λ|=(5λ)(2λ)2×2

                    =λ2+7λ+6=0

                 λ=7±724×1×62×1=7±49242

                 ∴ λ1=1,    λ2=6    (A의 고유값)

 

    b1) λ1에 대응하는 A의 고유벡터

        식(2*),        λ=λ1=1

                4x1+2x2=02x1x2=0}x2=2x1

                이 식에서 λ1=1에 대응하는 A 의 고유벡터가 상수배까지 결정

        만일 x1=1을 택하면 x2=2 따라서 λ=1에 대응하는 의 고유벡터는

                x1=[12]        Ax=λx

        검증

                                       A         x           λx

                Ax1=[5222][12]=[12]=(1)x1=λ1x1

 

    b2) λ2에 대응하는 A의 고유벡터

                λ=λ2=6        식(2*)는

                (5x)x12x2=0

                2x1+(2λ)x2=0에 대입

                x1+2x2=02x1+4x2=0}x2=x1/2

        만일 x1=2를 택하면 x2=1

        λ2=6에 대응하는 A의 고유벡터는

                x2=[21]

 

    일반화

                a11x1+a12x2++a1nxn=λx1a21x1+a22x2++a2nxn=λx2an1x1+an2x2++annxn=λxn

 

                (a11λ)x1+a12x2++a1nxn=0a21x1+(a22λ)x2++a2nxn=0an1x1+an2x2++(annλ)xn=0

 

    행렬식 표현

                (AλI)x=0

        Crammer  정리에서 x0detA=0 (필요충분 조건)

                D(λ)=det(AλI)=|a11λa12a1na21a21λa2nan1an2annλ|    특성방정식

                D(λ) 전개 → A의 특성 방정식

 

정리 1 (고유값)

        정방행렬 A의 고유값은 A의 특성방정식의 근이다. 따라서, n×n 행렬은 적어도 하나 이상, 만아야 n개의 서로 다른 고유값을 갖는다.

 

정리 2 (고유벡터)

        만일 x가 행렬 A의 고유값 λ에 대응하는 고유벡터인 경우, 임의의 k0에 대하여 kx도 고유벡터가 된다.

        Ax=λx

        k(Ax)=A(kx)=λ(kx)

 

○ 조립제법

        x34x2+3x+1=(x2)+    ( : 몪,    : 나머지)