요르단 표준형의 장점
고윳값, 고유 벡터의 모양이 보입니다.
거듭제곱을 구체적으로 계산할 수 있습니다.
이 두 가지는 폭주 판정과 관련이 있습니다.
요르단 표준형은 블록대각이므로 고윳값이나 거듭제곱 계산은 블록별로 보면 알 수 있습니다.
다음 요르단 셀을 이용하여 학습하겠습니다.
$B = \begin{pmatrix}
7 & 1 & 0 & 0\\
0 & 7 & 1 & 0\\
0 & 0 & 7 & 1\\
0 & 0 & 0 & 7\\
\end{pmatrix}$
요르단 표준형의 고윳값
요르단 셀 $B$의 고윳값은 7밖에 없습니다. (상삼각 행렬)
$\mathbf {p} = (\alpha, 0, 0, 0)^{T}$가 고윳값 7의 고유 벡터입니다. ($\alpha \neq 0$) 이 외의 고유 벡터는 없습니다.
요르단 표준형 $J$의 고윳값, 고유벡터는 한눈에 알 수 있습니다.
대각 성분이 고윳값 $\lambda$
대각 성분의 $\lambda$ 개수가 고윳값 $\lambda$가 몇 중 해인지(대수적 중복도)에 대응
대각 성분이 $\lambda$인 요르단 셀의 개수가 고윳값 $\lambda$에 선형 독립인 고유 벡터 개수(기하적 중복도)에 대응
위의 성질에서 $A$가 요르단 표준형으로 변환되면 $A$의 고윳값, 고유 벡터가 어떻게 되어 있는지 알 수 있습니다.
특히 고윳값에 중해가 없을 경우 요르단 표준형은 대각 행렬이 될 수밖에 없습니다. 즉,
고윳값에 중해가 없을 경우 대각 화할 수 있습니다.
요르단 표준형의 거듭제곱
$B = 7I + Z$
$Z = \begin {pmatrix}
0 & 1 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1\\
0 & 0 & 0 & 0\\
\end {pmatrix}$
위와 같이 분해해 두면 계산할 때 예측하기 쉽습니다. 특징은 $Z$가 왼쪽에 곱하면 1행을 밀고 오른쪽에 곱하면 1열 밀기 작용을 합니다.
$\begin {pmatrix}
0 & 1 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1\\
0 & 0 & 0 & 0\\
\end {pmatrix} \begin {pmatrix}
a\\
b\\
c\\
d\\
\end{pmatrix} = \begin{pmatrix}
b\\
c\\
d\\
0\\
\end {pmatrix}$
$\begin {pmatrix}
a & e & i & m\\
b & f & j & n\\
c & g & k & o\\
d & h & l & p\\
\end {pmatrix} \begin {pmatrix}
0 & 1 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1\\
0 & 0 & 0 & 0\\
\end{pmatrix} = \begin{pmatrix}
0 & a & e & i\\
0 & b & f & j\\
0 & c & g & k\\
0 & d & h & l\\
\end {pmatrix}$
모양입니다.
이 성질을 이용하면 $Z$의 거듭제곱은 간단합니다.
$Z^{2} = \begin {pmatrix}
0 & 0 & 1 & 0\\
0 & 0 & 0 & 1\\
0 & 0 & 0 & 0\\
0 & 0 & 0 & 0\\
\end {pmatrix}$, $Z^{3} = \begin{pmatrix}
0 & 0 & 0 & 1\\
0 & 0 & 0 & 0\\
0 & 0 & 0 & 0\\
0 & 0 & 0 & 0\\
\end{pmatrix}$, $Z^{4} = Z^{5} = \cdots = O$
위와 같은 식으로 2 제곱, 3 제곱...으로 늘릴 때마다 1의 위치가 오른쪽 위로 밀리는 것을 알 수 있습니다.
이제 $B^{2} = (7I + Z)^{2} = 7^{2} I + 2 \cdot 7Z + Z^{2}$에서
$B^{2} = \begin {pmatrix}
7^{2} & 2 \cdot 7 & 1 & 0\\
0 & 7^{2} & 2 \cdot 7 & 1\\
0 & 0 & 7^{2} & 2 \cdot 7\\
0 & 0 & 0 & 7^{2}\\
\end {pmatrix}$
$B^{3} = (7I + Z)^{3} = 7^{3} I + 3 \cdot 7^{2} Z + 3 \cdot 7Z^{2} + Z^{3}$에서
$B^{3} = \begin {pmatrix}
7^{3} & 3 \cdot 7^{2} & 3 \cdot 7 & 1\\
0 & 7^{3} & 3 \cdot 7^{2} & 3 \cdot 7\\
0 & 0 & 7^{3} & 3 \cdot 7^{2}\\
0 & 0 & 0 & 7^{3}\\
\end {pmatrix}$
$B^{4} = (7I + Z)^{4} = 7^{4} I + 4 \cdot 7^{3} Z + 6 \cdot 7^{2} Z^{2} + 4 \cdot 7Z^{3} + Z^{4}$ ($Z^{4} = O$)
$B^{4} = \begin {pmatrix}
7^{4} & 4 \cdot 7^{3} & 6 \cdot 7^{2} & 4 \cdot 7\\
0 & 7^{4} & 4 \cdot 7^{3} & 6 \cdot 7^{2}\\
0 & 0 & 7^{4} & 4 \cdot 7^{3}\\
0 & 0 & 0 & 7^{4}\\
\end {pmatrix}$
더 큰 t 제곱에서도
$B^{t} = (7I + Z)^{t} = 7^{t} I + t \cdot 7^{t-1} Z + _{2}^{t}\textrm {C} \cdot 7^{t-2} Z^{2} + _{3}^{t}\textrm {C} \cdot 7^{t-3} Z^{3} + _{4}^{t}\textrm {C} \cdot 7^{t-4} Z^{4} + \cdots + _{t}^{t-2}\textrm {C} \cdot 7^{2} Z^{t-2} + t \cdot 7Z^{t-1} + Z^{t}$에서 ($Z^{4} = Z^{5} = Z^{6} = \cdots = O$)
$B^{t} = \begin {pmatrix}
7^{t} & t \cdot 7^{t-1} & _{2}^{t}\textrm {C} \cdot 7^{t-2} & _{3}^{t}\textrm{C} \cdot 7^{t-3}\\
0 & 7^{t} & t \cdot 7^{t-1} & _{2}^{t}\textrm{C} \cdot 7^{t-2}\\
0 & 0 & 7^{t} & t \cdot 7^{t-1}\\
0 & 0 & 0 & 7^{t}\\
\end {pmatrix}$
도 구해집니다(t = 1, 2,...).
대각 성분이 $\lambda$인 요르단 셀 $B$에 대해 $B^{t}$를 계산하면 $_{t}^{s}\textrm {C}\lambda^{t-s}$ 항이 나타납니다(단, s > t인 경우에는 $_{t}^{s}\textrm{C}\lambda^{t-s} = 0$으로 간주하고 $\lambda^{0} = 1$이라 간주합니다).
다음과 같이도 말할 수 있습니다.
$f(\lambda) = \lambda^{t}$로 두고, $f$를 $\lambda$로 $s$회 미분한 식을
$f^{(s)}(\lambda) = \frac {d^{s}}{d\lambda^{s}}f(\lambda)$
로 두면 크기 m의 요르단 셀
$B = \begin {pmatrix}
\lambda & 1 & & \\
& \ddots & \ddots & \\
& & \ddots & 1\\
& & & \lambda
\end {pmatrix}$
의 t제곱은
$\begin {pmatrix}
f(\lambda) & f^{(1)}(\lambda) & \frac {1}{2} f^{(2)}(\lambda) & \frac{1}{3!}f^{(3)}(\lambda) & \cdots & \frac{1}{(m-1)!}f^{(m-1)}(\lambda)\\
& f(\lambda) & f^{(1)}(\lambda) & \frac{1}{2}f^{(2)}(\lambda) & \cdots & \frac {1}{(m-2)!} f^{(m-2)}(\lambda)\\
& & \ddots & \ddots & \ddots&\vdots \\
& & & \ddots &\ddots &\frac {1}{2} f^{(2)}(\lambda)\\
& & & &\ddots &f^{(1)}(\lambda)\\
& & & & &f(\lambda)\\
\end {pmatrix}$ (4.21)
이유는 $f^{s}(\lambda) = t(t-1) \cdots (t - s + 1)\lambda^{t-s}$이기 때문입니다.
$\frac {1}{s!} f^{s}(\lambda) = _{s}^{t}\textrm {C}\lambda^{t-s}$
식 (4.20)의 $J$라면 각 요르단 셀
$B_{1} = \begin {pmatrix}
3 & 1& 0 & 0 &0 \\
0 & 3 & 1 & 0& 0\\
0 & 0& 3& 1 &0 \\
0 & 0& 0 & 3&1 \\
0 & 0 & 0&0 & 3
\end {pmatrix}$, $B_{2} = \begin{pmatrix}
3 & 1 \\
0 & 3 \\
\end{pmatrix}$, $B_{3} = (4)$, $B_{4} = \begin{pmatrix}
5 & 1& 0\\
0& 5 &1 \\
0& 0 &5
\end{pmatrix}$
에 대해 각각 t제곱을 구해두면
$J^{t} = \begin{pmatrix}
B_{1}^{t} & O & O &O \\
O& B_{2}^{t}&O & O\\
O& O& B_{3}^{t}& O\\
O& O & O& B_{4}^{t}
\end{pmatrix}$
로 구해집니다.
요르단 표준형이 아닌 정방 행렬 $A$의 거듭제곱 $A^{t}$는 계산 가능할까요?
$A$를 요르단 표준형 $J$로 변환하면 전과 같이 $A^{t}$를 계산할 수 있습니다.
$P^{-1} AP = J$라는 것은 $A = PJP^{-1}$이고, 그 t 제곱은
$A^{t} = (PJP^{-1})^{t} = PJ^{t} P^{-1}$입니다.
내용이 결고 쉽지 않습니다. 중간중간 막힌 부분은 다시 복습해야 합니다.
2019/09/04 - [선형대수/고윳값, 대각화, 요르단 표준형 - 폭주 위험] - 4.5.2 고윳값, 고유 벡터의 성질
감사합니다.
'프로그래머를 위한 선형대수 > 고윳값, 대각화, 요르단 표준형 - 폭주 위험' 카테고리의 다른 글
4.7.2 대각까지는 못하더라도 - 요르단 표준형 (0) | 2019.09.24 |
---|---|
4.7 대각화할 수 없는 경우 (0) | 2019.09.20 |
4.6.4 대각화할 수 있는 경우, 4.6.4 결론: 고윳값(실수부)의 부호 (0) | 2019.09.19 |
4.6.3 대각행렬일 때 (0) | 2019.09.18 |
4.6.2 1차원일 때 (0) | 2019.09.17 |
댓글