30.15 대각화의 절차: 고유벡터 행렬과 고유값 대각 행렬의 구성

30.15 대각화의 절차: 고유벡터 행렬과 고유값 대각 행렬의 구성

1. 대각화 절차의 개요

n \times n 행렬 A의 대각화 A = PDP^{-1}을 구성하는 절차는 다음의 네 단계로 이루어진다.

  1. 특성 다항식의 계산과 고유값 결정
  2. 각 고유값에 대한 고유 공간의 기저 계산
  3. 대각화 가능성 판별
  4. 고유벡터 행렬 P와 대각 행렬 D의 구성

2. 단계 1: 고유값의 결정

특성 다항식 p_A(\lambda) = \det(\lambda I - A)를 계산하고, p_A(\lambda) = 0의 근을 구한다.

p_A(\lambda) = (\lambda - \lambda_1)^{m_1}(\lambda - \lambda_2)^{m_2} \cdots (\lambda - \lambda_k)^{m_k}

서로 다른 고유값 \lambda_1, \ldots, \lambda_k와 각각의 대수적 중복도 m_1, \ldots, m_k를 기록한다. \sum_{i=1}^k m_i = n이어야 한다.

3. 단계 2: 고유 공간의 기저 계산

각 고유값 \lambda_i에 대하여 동차 선형 계

(A - \lambda_i I)v = 0

을 풀어 고유 공간 E_{\lambda_i} = \ker(A - \lambda_i I)의 기저를 구한다.

구체적 절차:

  1. A - \lambda_i I를 구성한다.
  2. 가우스 소거법으로 기약 행 사다리꼴(reduced row echelon form)로 환원한다.
  3. 자유 변수를 식별하고, 자유 변수에 단위값을 순차적으로 대입하여 기저 벡터를 구한다.

기하학적 중복도 d_i = \dim(E_{\lambda_i})를 확인한다.

4. 단계 3: 대각화 가능성 판별

모든 고유값에 대하여 d_i = m_i인지 확인한다.

  • 모든 i에서 d_i = m_i이면 대각화 가능하다. 단계 4로 진행한다.
  • 어떤 i에서 d_i < m_i이면 대각화 불가능이다. 절차를 종료한다.

5. 단계 4: PD의 구성

5.1 고유벡터 행렬 P

각 고유 공간의 기저 벡터를 열로 배치하여 n \times n 행렬 P를 구성한다.

P = (v_1 \mid v_2 \mid \cdots \mid v_n)

열의 순서는 자유롭게 선택할 수 있으나, D의 대각 성분의 순서와 일치시켜야 한다.

5.2 고유값 대각 행렬 D

D = \begin{pmatrix} \lambda_{(1)} & 0 & \cdots & 0 \\ 0 & \lambda_{(2)} & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda_{(n)} \end{pmatrix}

Dj번째 대각 성분 \lambda_{(j)}Pj번째 열 v_j에 대응하는 고유값이다.

5.3 대응 관계의 핵심

Pj번째 열과 Dj번째 대각 성분은 반드시 동일한 고유쌍에 속해야 한다. 즉, Av_j = \lambda_{(j)} v_j.

6. 완전 계산 예시

6.1 예시 1: 3 \times 3 행렬의 대각화

A = \begin{pmatrix} 1 & 0 & 0 \\ -2 & 3 & 0 \\ 0 & 2 & 3 \end{pmatrix}

단계 1: 하삼각 행렬이므로 고유값은 대각 성분이다.

p_A(\lambda) = (\lambda - 1)(\lambda - 3)^2

\lambda_1 = 1 (m_1 = 1), \lambda_2 = 3 (m_2 = 2).

단계 2:

\lambda_1 = 1:

A - I = \begin{pmatrix} 0 & 0 & 0 \\ -2 & 2 & 0 \\ 0 & 2 & 2 \end{pmatrix} \xrightarrow{\text{RREF}} \begin{pmatrix} 1 & 0 & -1 \\ 0 & 1 & 1 \\ 0 & 0 & 0 \end{pmatrix}

x_3 자유. x_1 = x_3, x_2 = -x_3. v_1 = \begin{pmatrix} 1 \\ -1 \\ 1 \end{pmatrix}. d_1 = 1 = m_1. \checkmark

\lambda_2 = 3:

A - 3I = \begin{pmatrix} -2 & 0 & 0 \\ -2 & 0 & 0 \\ 0 & 2 & 0 \end{pmatrix} \xrightarrow{\text{RREF}} \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 0 \end{pmatrix}

x_3 자유. x_1 = 0, x_2 = 0. v_2 = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix}. d_2 = 1.

d_2 = 1 < m_2 = 2이므로 대각화 불가능이다.

6.2 예시 2: 대각화 가능한 3 \times 3 행렬

A = \begin{pmatrix} 2 & 0 & 0 \\ 0 & 3 & 4 \\ 0 & 4 & 9 \end{pmatrix}

단계 1:

p_A(\lambda) = (\lambda - 2)[(\lambda - 3)(\lambda - 9) - 16] = (\lambda - 2)(\lambda^2 - 12\lambda + 11) = (\lambda - 2)(\lambda - 1)(\lambda - 11)

\lambda_1 = 1, \lambda_2 = 2, \lambda_3 = 11. 모두 단순 고유값이므로 대각화 가능하다.

단계 2:

\lambda_1 = 1: (A - I)v = 0 \Rightarrow v_1 = \begin{pmatrix} 0 \\ -2 \\ 1 \end{pmatrix}

\lambda_2 = 2: (A - 2I)v = 0 \Rightarrow v_2 = \begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix}

\lambda_3 = 11: (A - 11I)v = 0 \Rightarrow v_3 = \begin{pmatrix} 0 \\ 1 \\ 2 \end{pmatrix}

단계 4:

P = \begin{pmatrix} 0 & 1 & 0 \\ -2 & 0 & 1 \\ 1 & 0 & 2 \end{pmatrix}, \quad D = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 11 \end{pmatrix}

검증: AP의 각 열이 \lambda_i v_i임을 확인한다.

AP의 1열: A\begin{pmatrix}0\\-2\\1\end{pmatrix} = \begin{pmatrix}0\\-2\\1\end{pmatrix} = 1 \cdot v_1. \checkmark

AP의 2열: A\begin{pmatrix}1\\0\\0\end{pmatrix} = \begin{pmatrix}2\\0\\0\end{pmatrix} = 2 \cdot v_2. \checkmark

AP의 3열: A\begin{pmatrix}0\\1\\2\end{pmatrix} = \begin{pmatrix}0\\11\\22\end{pmatrix} = 11 \cdot v_3. \checkmark

7. 중복 고유값의 처리

중복 고유값이 존재하고 d_i = m_i인 경우, 해당 고유 공간의 d_i개 기저 벡터를 P의 연속된 열로 배치한다. D의 대응하는 대각 성분은 모두 \lambda_i가 된다.

P = (\cdots \mid v_{i,1} \mid v_{i,2} \mid \cdots \mid v_{i,d_i} \mid \cdots), \quad D = \text{diag}(\cdots, \underbrace{\lambda_i, \lambda_i, \ldots, \lambda_i}_{d_i\text{개}}, \cdots)

8. 실용적 고려 사항

8.1 고유벡터의 스케일링

고유벡터는 비영 스칼라 배의 자유도를 가진다. 실용적으로는 정수 성분을 갖도록 스케일링하거나, 단위 벡터(\|v\| = 1)로 정규화하는 것이 일반적이다.

8.2 열 순서의 선택

P의 열 순서(따라서 D의 대각 성분 순서)는 자유롭게 선택할 수 있다. 고유값을 크기 순으로 정렬하거나 절대값 순으로 정렬하는 것이 통상적이다.

8.3 P^{-1}의 계산

A = PDP^{-1} 관계에서 P^{-1}이 명시적으로 필요한 경우, 가우스-조르단 소거법 또는 수반 행렬 공식으로 계산한다. 대칭 행렬의 경우 P를 직교 행렬로 취하면 P^{-1} = P^T이므로 역행렬 계산이 불필요하다.