早稲田大学 創造理工学研究科 経営システム工学専攻 2013年8月実施 线性代数
Author
思齐塾 , 祭音Myyura
Description
次の行列の固有値と固有ベクトルをすべて求めよ.
( 1 1 0 0 0 2 1 0 0 0 2 1 0 0 0 1 ) \begin{pmatrix} 1 & 1 & 0 & 0 \\ 0 & 2 & 1 & 0 \\ 0 & 0 & 2 & 1 \\ 0 & 0 & 0 & 1 \end{pmatrix} 1 0 0 0 1 2 0 0 0 1 2 0 0 0 1 1
题目描述
求矩阵
( 1 1 0 0 0 2 1 0 0 0 2 1 0 0 0 1 ) \begin{pmatrix}
1&1&0&0\\
0&2&1&0\\
0&0&2&1\\
0&0&0&1
\end{pmatrix} 1 0 0 0 1 2 0 0 0 1 2 0 0 0 1 1
的全部特征值及其对应的全部特征向量。
Kai
Let A = ( 1 1 0 0 0 2 1 0 0 0 2 1 0 0 0 1 ) A = \begin{pmatrix} 1 & 1 & 0 & 0 \\ 0 & 2 & 1 & 0 \\ 0 & 0 & 2 & 1 \\ 0 & 0 & 0 & 1 \end{pmatrix} A = 1 0 0 0 1 2 0 0 0 1 2 0 0 0 1 1 .
To find the eigenvalues, we solve the characteristic equation det ( A − λ I ) = 0 \det(A - \lambda I) = 0 det ( A − λ I ) = 0 .
A − λ I = ( 1 − λ 1 0 0 0 2 − λ 1 0 0 0 2 − λ 1 0 0 0 1 − λ ) A - \lambda I = \begin{pmatrix} 1-\lambda & 1 & 0 & 0 \\ 0 & 2-\lambda & 1 & 0 \\ 0 & 0 & 2-\lambda & 1 \\ 0 & 0 & 0 & 1-\lambda \end{pmatrix} A − λ I = 1 − λ 0 0 0 1 2 − λ 0 0 0 1 2 − λ 0 0 0 1 1 − λ
det ( A − λ I ) = ( 1 − λ ) ( 2 − λ ) ( 2 − λ ) ( 1 − λ ) = ( 1 − λ ) 2 ( 2 − λ ) 2 = 0 \det(A - \lambda I) = (1-\lambda)(2-\lambda)(2-\lambda)(1-\lambda) = (1-\lambda)^2(2-\lambda)^2 = 0 det ( A − λ I ) = ( 1 − λ ) ( 2 − λ ) ( 2 − λ ) ( 1 − λ ) = ( 1 − λ ) 2 ( 2 − λ ) 2 = 0
So the eigenvalues are λ 1 = 1 \lambda_1 = 1 λ 1 = 1 (with multiplicity 2) and λ 2 = 2 \lambda_2 = 2 λ 2 = 2 (with multiplicity 2).
For λ 1 = 1 \lambda_1 = 1 λ 1 = 1 :
A − I = ( 0 1 0 0 0 1 1 0 0 0 1 1 0 0 0 0 ) A - I = \begin{pmatrix} 0 & 1 & 0 & 0 \\ 0 & 1 & 1 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 \end{pmatrix} A − I = 0 0 0 0 1 1 0 0 0 1 1 0 0 0 1 0
Solving ( A − I ) v = 0 (A - I)v = 0 ( A − I ) v = 0 , we get:
v 2 = 0 v_2 = 0 v 2 = 0
v 2 + v 3 = 0 v_2 + v_3 = 0 v 2 + v 3 = 0 , so v 3 = 0 v_3 = 0 v 3 = 0
v 3 + v 4 = 0 v_3 + v_4 = 0 v 3 + v 4 = 0 , so v 4 = 0 v_4 = 0 v 4 = 0
Thus v = ( v 1 0 0 0 ) = v 1 ( 1 0 0 0 ) v = \begin{pmatrix} v_1 \\ 0 \\ 0 \\ 0 \end{pmatrix} = v_1 \begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix} v = v 1 0 0 0 = v 1 1 0 0 0 .
Eigenvector for λ 1 = 1 \lambda_1 = 1 λ 1 = 1 is v 1 = ( 1 0 0 0 ) v_1 = \begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix} v 1 = 1 0 0 0 .
Let's find a generalized eigenvector. We need to solve ( A − I ) 2 w = 0 (A-I)^2 w = 0 ( A − I ) 2 w = 0 . We want to find w w w such that ( A − I ) w = v 1 (A-I)w = v_1 ( A − I ) w = v 1 .
( A − I ) w = ( 0 1 0 0 0 1 1 0 0 0 1 1 0 0 0 0 ) ( w 1 w 2 w 3 w 4 ) = ( 1 0 0 0 ) (A - I) w = \begin{pmatrix} 0 & 1 & 0 & 0 \\ 0 & 1 & 1 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} w_1 \\ w_2 \\ w_3 \\ w_4 \end{pmatrix} = \begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix} ( A − I ) w = 0 0 0 0 1 1 0 0 0 1 1 0 0 0 1 0 w 1 w 2 w 3 w 4 = 1 0 0 0
w 2 = 1 w_2 = 1 w 2 = 1
w 2 + w 3 = 0 w_2 + w_3 = 0 w 2 + w 3 = 0 , so w 3 = − 1 w_3 = -1 w 3 = − 1
w 3 + w 4 = 0 w_3 + w_4 = 0 w 3 + w 4 = 0 , so w 4 = 1 w_4 = 1 w 4 = 1
Thus w = ( w 1 1 − 1 1 ) w = \begin{pmatrix} w_1 \\ 1 \\ -1 \\ 1 \end{pmatrix} w = w 1 1 − 1 1 . Choosing w 1 = 0 w_1 = 0 w 1 = 0 , we get w = ( 0 1 − 1 1 ) w = \begin{pmatrix} 0 \\ 1 \\ -1 \\ 1 \end{pmatrix} w = 0 1 − 1 1 . This is a generalized eigenvector.
For λ 2 = 2 \lambda_2 = 2 λ 2 = 2 :
A − 2 I = ( − 1 1 0 0 0 0 1 0 0 0 0 1 0 0 0 − 1 ) A - 2I = \begin{pmatrix} -1 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & -1 \end{pmatrix} A − 2 I = − 1 0 0 0 1 0 0 0 0 1 0 0 0 0 1 − 1
Solving ( A − 2 I ) v = 0 (A - 2I)v = 0 ( A − 2 I ) v = 0 , we get:
− v 1 + v 2 = 0 -v_1 + v_2 = 0 − v 1 + v 2 = 0 , so v 1 = v 2 v_1 = v_2 v 1 = v 2
v 3 = 0 v_3 = 0 v 3 = 0
v 4 = 0 v_4 = 0 v 4 = 0
− v 4 = 0 -v_4 = 0 − v 4 = 0
Thus v = ( v 2 v 2 0 0 ) = v 2 ( 1 1 0 0 ) v = \begin{pmatrix} v_2 \\ v_2 \\ 0 \\ 0 \end{pmatrix} = v_2 \begin{pmatrix} 1 \\ 1 \\ 0 \\ 0 \end{pmatrix} v = v 2 v 2 0 0 = v 2 1 1 0 0 .
Eigenvector for λ 2 = 2 \lambda_2 = 2 λ 2 = 2 is v 2 = ( 1 1 0 0 ) v_2 = \begin{pmatrix} 1 \\ 1 \\ 0 \\ 0 \end{pmatrix} v 2 = 1 1 0 0 .
Let's find a generalized eigenvector. We need to solve ( A − 2 I ) 2 w = 0 (A-2I)^2 w = 0 ( A − 2 I ) 2 w = 0 . We want to find w w w such that ( A − 2 I ) w = v 2 (A-2I)w = v_2 ( A − 2 I ) w = v 2 .
( A − 2 I ) w = ( − 1 1 0 0 0 0 1 0 0 0 0 1 0 0 0 − 1 ) ( w 1 w 2 w 3 w 4 ) = ( 1 1 0 0 ) (A - 2I) w = \begin{pmatrix} -1 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & -1 \end{pmatrix} \begin{pmatrix} w_1 \\ w_2 \\ w_3 \\ w_4 \end{pmatrix} = \begin{pmatrix} 1 \\ 1 \\ 0 \\ 0 \end{pmatrix} ( A − 2 I ) w = − 1 0 0 0 1 0 0 0 0 1 0 0 0 0 1 − 1 w 1 w 2 w 3 w 4 = 1 1 0 0
− w 1 + w 2 = 1 -w_1 + w_2 = 1 − w 1 + w 2 = 1
w 3 = 1 w_3 = 1 w 3 = 1
w 4 = 0 w_4 = 0 w 4 = 0
− w 4 = 0 -w_4 = 0 − w 4 = 0
Thus w = ( w 1 w 1 + 1 1 0 ) w = \begin{pmatrix} w_1 \\ w_1 + 1 \\ 1 \\ 0 \end{pmatrix} w = w 1 w 1 + 1 1 0 . Choosing w 1 = 0 w_1 = 0 w 1 = 0 , we have w = ( 0 1 1 0 ) w = \begin{pmatrix} 0 \\ 1 \\ 1 \\ 0 \end{pmatrix} w = 0 1 1 0 .
Therefore, the eigenvalues are 1 1 1 and 2 2 2 , each with algebraic multiplicity 2 2 2 , and all corresponding eigenvectors are
ker ( A − I ) = span { ( 1 0 0 0 ) } , ker ( A − 2 I ) = span { ( 1 1 0 0 ) } . \ker(A-I)=\operatorname{span}\left\{
\begin{pmatrix}1\\0\\0\\0\end{pmatrix}
\right\},
\qquad
\ker(A-2I)=\operatorname{span}\left\{
\begin{pmatrix}1\\1\\0\\0\end{pmatrix}
\right\}. ker ( A − I ) = span ⎩ ⎨ ⎧ 1 0 0 0 ⎭ ⎬ ⎫ , ker ( A − 2 I ) = span ⎩ ⎨ ⎧ 1 1 0 0 ⎭ ⎬ ⎫ .
All eigenvectors are the nonzero vectors in these two eigenspaces.
The vectors ( 0 1 − 1 1 ) \begin{pmatrix}0\\1\\-1\\1\end{pmatrix} 0 1 − 1 1 and
( 0 1 1 0 ) \begin{pmatrix}0\\1\\1\\0\end{pmatrix} 0 1 1 0 found above are generalized eigenvectors, not eigenvectors.