当記事は「統計学のための数学入門$30$講(朝倉書店)」の読解サポートにあたってChapter.$22$の「固有値と固有ベクトル」の章末問題の解答の作成を行いました。
基本的には書籍の購入者向けの解説なので、まだ入手されていない方は購入の上ご確認ください。また、解説はあくまでサイト運営者が独自に作成したものであり、書籍の公式ページではないことにご注意ください。
・書籍解答まとめ
https://www.hello-statisticians.com/answer_textbook_math#math_stat
Contents
本章のまとめ
固有値・固有ベクトルの性質
問題$22.1$で使用するので、$22.1$節、P.$140$の「固有値・固有ベクトルの性質」の内容を以下にまとめる。
$p$次の正方行列$A$の固有値$\lambda_i$に関して下記が成立する。
$(1) \,$ 固有方程式は$p$次方程式であり、固有値は$p$個ある。ただし、重根の場合には重複度を含めて$p$個と考える。
$(2) \,$ $\lambda_{1} + \lambda_{2} + \cdots \lambda_{p} = \mathrm{Tr}(A)$
$(3) \,$ $\lambda_{1} \lambda_{2} \cdots \lambda_{p} = |A|$
$(4) \,$ $|A| \, \iff \,$ 少なくとも$1$つの固有値が$0$
$(5) \,$ $A$の逆行列が存在する $\, \iff \,$ 全ての固有値が$0$ではない
$(6) \,$ $A$の逆行列$A^{-1}$が存在するとき、$A$の固有値が$\lambda_{1}, \lambda_{2}, \cdots \lambda_{p}$であれば$A^{-1}$の固有値は$1/\lambda_{1}, 1/\lambda_{2}, \cdots 1/\lambda_{p}$である。また$\lambda_i$に対応する$A$の固有ベクトルが$\mathbf{x}_{i}$であれば、$1/\lambda_i$に対応する$A^{-1}$の固有ベクトルは$\mathbf{x}_{i}$である。
$(7) \,$ $A$の固有値と$A^{\mathrm{T}}$の固有値は同じである。
演習問題解答
問題$22.1$
$$
\large
\begin{align}
A = \left( \begin{array}{cc} 1 & 1 \\ 0 & 2 \end{array} \right)
\end{align}
$$
上記の$A$の逆行列$A^{-1}$は下記のように表せる。
$$
\large
\begin{align}
A^{-1} &= \frac{1}{2-0} \left( \begin{array}{cc} 2 & -1 \\ 0 & 1 \end{array} \right) \\
&= \frac{1}{2} \left( \begin{array}{cc} 2 & -1 \\ 0 & 1 \end{array} \right) \\
&= \left( \begin{array}{cc} 1 & -1/2 \\ 0 & 1/2 \end{array} \right)
\end{align}
$$
$A$の固有値・固有ベクトル
固有方程式$\det(A – \lambda I_2)=0$は下記のように解ける。
$$
\large
\begin{align}
\det(A – \lambda I_2) &= \left| \begin{array}{cc} 1-\lambda & 1 \\ 0 & 2-\lambda \end{array} \right| \\
&= (1-\lambda)(2-\lambda) = 0 \\
\lambda &= 1, \, 2
\end{align}
$$
上記の長さ$1$の固有ベクトルを$\displaystyle \mathbf{u} = \left( \begin{array}{c} x \\ y \end{array} \right)$とおくと、$\lambda = 1, \, 2$に対応する$\mathbf{u}$はそれぞれ下記のように得られる。
・$\lambda = 1$
$$
\large
\begin{align}
\left( \begin{array}{cc} 1 & 1 \\ 0 & 2 \end{array} \right)\left( \begin{array}{c} x \\ y \end{array} \right) &= \left( \begin{array}{c} x \\ y \end{array} \right) \\
\left( \begin{array}{c} x+y \\ 2y \end{array} \right) &= \left( \begin{array}{c} x \\ y \end{array} \right) \\
\mathbf{u} &= \left( \begin{array}{c} x \\ y \end{array} \right) = \left( \begin{array}{c} 1 \\ 0 \end{array} \right)
\end{align}
$$
・$\lambda = 2$
$$
\large
\begin{align}
\left( \begin{array}{cc} 1 & 1 \\ 0 & 2 \end{array} \right)\left( \begin{array}{c} x \\ y \end{array} \right) &= 2 \left( \begin{array}{c} x \\ y \end{array} \right) \\
\left( \begin{array}{c} x+y \\ 2y \end{array} \right) &= \left( \begin{array}{c} 2x \\ 2y \end{array} \right) \\
\mathbf{u} &= \left( \begin{array}{c} x \\ y \end{array} \right) = \frac{1}{\sqrt{2}} \left( \begin{array}{c} 1 \\ 1 \end{array} \right)
\end{align}
$$
$A^{-1}$の固有値・固有ベクトル
固有方程式$\det(A^{-1} – \lambda I_2)=0$は下記のように解ける。
$$
\large
\begin{align}
\det(A^{-1} – \lambda I_2) &= \left| \begin{array}{cc} 1-\lambda & -1/2 \\ 0 & 1/2-\lambda \end{array} \right| \\
&= \left( 1-\lambda \right) \left( \frac{1}{2} – \lambda \right) = 0 \\
\lambda &= 1, \, \frac{1}{2}
\end{align}
$$
上記の長さ$1$の固有ベクトルを$\displaystyle \mathbf{u} = \left( \begin{array}{c} x \\ y \end{array} \right)$とおくと、$\lambda = 1, \, 1/2$に対応する$\mathbf{u}$はそれぞれ下記のように得られる。
・$\lambda = 1$
$$
\large
\begin{align}
\left( \begin{array}{cc} 1 & -1/2 \\ 0 & 1/2 \end{array} \right)\left( \begin{array}{c} x \\ y \end{array} \right) &= \left( \begin{array}{c} x \\ y \end{array} \right) \\
\left( \begin{array}{c} x-y/2 \\ y/2 \end{array} \right) &= \left( \begin{array}{c} x \\ y \end{array} \right) \\
\mathbf{u} &= \left( \begin{array}{c} x \\ y \end{array} \right) = \left( \begin{array}{c} 1 \\ 0 \end{array} \right)
\end{align}
$$
・$\lambda = 1/2$
$$
\large
\begin{align}
\left( \begin{array}{cc} 1 & -1/2 \\ 0 & 1/2 \end{array} \right)\left( \begin{array}{c} x \\ y \end{array} \right) &= \frac{1}{2} \left( \begin{array}{c} x \\ y \end{array} \right) \\
\left( \begin{array}{c} x-y/2 \\ y/2 \end{array} \right) &= \left( \begin{array}{c} x/2 \\ y/2 \end{array} \right) \\
\mathbf{u} &= \left( \begin{array}{c} x \\ y \end{array} \right) = \frac{1}{\sqrt{2}} \left( \begin{array}{c} 1 \\ 1 \end{array} \right)
\end{align}
$$
上記に対し、$(1)$〜$(3)$と$(5)$、$(6)$が成立することを以下確認を行う。
$(1) \,$ $2$次の正方行列$A, A^{-1}$に対し、固有方程式$\det(A – \lambda I_2)=0, \det(A^{-1} – \lambda I_2)$はそれぞれ$2$次方程式であり、それぞれ$2$個の方程式の解に対応する固有値を持つ。
$(2) \,$ $A$に関して$\lambda_1+\lambda_2 = 1+2 = 3 = \mathrm{Tr}(A)$が成立し、同様に$A^{-1}$に関して$\lambda_1+\lambda_2 = 1+1/2 = 3/2 = \mathrm{Tr}(A^{-1})$が成立する。
$(3) \,$ $A$に関して$\lambda_1 \times \lambda_2 = 2 = |A|$が成立し、同様に$A^{-1}$に関して$\lambda_1 \times \lambda_2 = 1/2 = |A^{-1}|$が成立する。
$(5) \,$ $A$の固有値は$1 \neq 0, 2 \neq 0$である。また、$A^{-1}$も存在する。
$(6) \,$ $A, A^{-1}$の固有値・固有ベクトルをそれぞれ確認することで、確認できる。
$A^{n}$の計算
$A$に関する固有ベクトルを縦に並べた行列を$U$、固有値を対角に並べた行列を$\Lambda$とおくと、$U, \Lambda$は下記のように表せる。
$$
\large
\begin{align}
U &= \left( \begin{array}{cc} 1 & 1/\sqrt{2} \\ 0 & 1/\sqrt{2} \end{array} \right) \\
\Lambda &= \left( \begin{array}{cc} 1 & 0 \\ 0 & 2 \end{array} \right)
\end{align}
$$
このとき、$A, U, \Lambda$に関して下記が成立する。
$$
\large
\begin{align}
AU &= U \Lambda \\
U^{-1}AU &= \Lambda \\
(U^{-1}AU)^{n} &= \Lambda^{n} \\
U^{-1} A^{n} U &= \Lambda^{n} \\
A^{n} &= U \Lambda^{n} U^{-1} \\
A^{n} &= \left( \begin{array}{cc} 1 & 1/\sqrt{2} \\ 0 & 1/\sqrt{2} \end{array} \right) \left( \begin{array}{cc} 1 & 0 \\ 0 & 2^{n} \end{array} \right) \left( \begin{array}{cc} 1 & -1 \\ 0 & \sqrt{2} \end{array} \right) \\
&= \left( \begin{array}{cc} 1 & 2^{n}/\sqrt{2} \\ 0 & 2^{n}/\sqrt{2} \end{array} \right) \left( \begin{array}{cc} 1 & -1 \\ 0 & \sqrt{2} \end{array} \right) \\
&= \left( \begin{array}{cc} 1 & 2^{n}-1 \\ 0 & 2^{n} \end{array} \right)
\end{align}
$$