広島大学 先進理工系科学研究科 情報科学プログラム 2017年8月実施 専門科目I 問題1
Author
samparker, 祭音Myyura
Description
A = [ 1 3 3 1 2 − 4 − 1 0 1 ] A = \begin{bmatrix} 1 & 3 & 3 \\ 1 & 2 & -4 \\ -1 & 0 & 1 \end{bmatrix} A = 1 1 − 1 3 2 0 3 − 4 1 とする。
(1) A A A を対称行列 S ( S T = S ) S \ (S^T = S) S ( S T = S ) と交代行列 T ( T T = − T ) T \ (T^T = -T) T ( T T = − T ) の和 ( A = S + T ) (A = S + T) ( A = S + T ) に分解せよ。ただし、A T A^T A T は、行列 A A A の転置を表す。
(2) S S S のすべての固有値と対応する固有空間を求めよ。
(3) T T T のすべての固有値と対応する固有空間を求めよ。
(4) 一般に、実交代行列の固有値は 0 0 0 または純虚数であることを示せ。
Let A = [ 1 3 3 1 2 − 4 − 1 0 1 ] A = \begin{bmatrix} 1 & 3 & 3 \\ 1 & 2 & -4 \\ -1 & 0 & 1 \end{bmatrix} A = 1 1 − 1 3 2 0 3 − 4 1 .
(1) Decompose A A A into the sum A = S + T A = S + T A = S + T of the symmetric matrix S S S (S T = S S^T = S S T = S ) and the alternative matrix T T T (T T = − T T^T = -T T T = − T ). Here A T A^T A T denotes the transpose of the matrix A A A .
(2) Find all the eigenvalues of the symmetric matrix S S S and a basis of the corresponding eigenspaces.
(3) Find all the eigenvalues of the alternative matrix T T T and a basis of the corresponding eigenspaces.
(4) Show that the eigenvalues of the real alternative matrix are 0 0 0 or purely imaginary numbers.
题目描述
给定
A = [ 1 3 3 1 2 − 4 − 1 0 1 ] . A=\begin{bmatrix}1&3&3\\1&2&-4\\-1&0&1\end{bmatrix}. A = 1 1 − 1 3 2 0 3 − 4 1 .
将 A A A 分解为对称矩阵 S S S (S T = S S^T=S S T = S )与实反对称矩阵 T T T (T T = − T T^T=-T T T = − T )之和 A = S + T A=S+T A = S + T ,其中 A T A^T A T 表示 A A A 的转置。
求 S S S 的全部特征值及各特征值对应的特征空间的一组基。
求 T T T 的全部特征值及各特征值对应的特征空间的一组基。
一般地,证明实反对称矩阵的特征值只能是 0 0 0 或纯虚数。
特征值与特征向量:先由矩阵的对称与反对称分解得到 S , T S,T S , T ,再分别求其谱与特征空间,并利用反对称性证明实反对称矩阵的特征值类型。
Kai
(1)
Let S = [ a b c b d e c e f ] S = \begin{bmatrix} a & b & c \\ b & d & e \\ c & e & f \end{bmatrix} S = a b c b d e c e f and T = [ 0 g h − g 0 i − h − i 0 ] T = \begin{bmatrix} 0 & g & h \\ -g & 0 & i \\ -h & -i & 0 \end{bmatrix} T = 0 − g − h g 0 − i h i 0 . Then we have
[ 1 3 3 1 2 − 4 − 1 0 1 ] = [ a b + g c + h b − g d e + i c − h e − i f ] \begin{bmatrix} 1 & 3 & 3 \\ 1 & 2 & -4 \\ -1 & 0 & 1 \end{bmatrix} = \begin{bmatrix} a & b+g & c+h \\ b-g & d & e+i \\ c-h & e-i & f \end{bmatrix} 1 1 − 1 3 2 0 3 − 4 1 = a b − g c − h b + g d e − i c + h e + i f
Solving the equations we have
a = 1 , d = 2 , f = 1 , h = 2 , c = 1 , g = 1 , b = 2 , i = − 2 , e = − 2 a = 1, d=2, f = 1, h=2, c=1, g=1, b=2, i=-2, e=-2 a = 1 , d = 2 , f = 1 , h = 2 , c = 1 , g = 1 , b = 2 , i = − 2 , e = − 2
Hence
S = [ 1 2 1 2 2 − 2 1 − 2 1 ] , T = [ 0 1 2 − 1 0 − 2 − 2 2 0 ] S = \begin{bmatrix} 1 & 2 & 1 \\ 2 & 2 & -2 \\ 1 & -2 & 1 \end{bmatrix}, T = \begin{bmatrix} 0 & 1 & 2 \\ -1 & 0 & -2 \\ -2 & 2 & 0 \end{bmatrix} S = 1 2 1 2 2 − 2 1 − 2 1 , T = 0 − 1 − 2 1 0 2 2 − 2 0
(2)
Eigenvalues
λ 1 = 2 , λ 2 = 4 , λ 3 = − 2 \lambda_1 = 2, \lambda_2 = 4, \lambda_3 = -2 λ 1 = 2 , λ 2 = 4 , λ 3 = − 2
A basis of the corresponding eigenspaces
v 1 = ( 1 , 0 , 1 ) , v 2 = ( − 1 , − 2 , 1 ) , v 3 = ( − 1 , 1 , 1 ) v_1=(1,0,1), v_2=(-1,-2,1), v_3=(-1,1,1) v 1 = ( 1 , 0 , 1 ) , v 2 = ( − 1 , − 2 , 1 ) , v 3 = ( − 1 , 1 , 1 )
(3)
Eigenvalues
λ 1 = 3 i , λ 2 = − 3 i , λ 3 = 0 \lambda_1 = 3i, \lambda_2 = -3i, \lambda_3 = 0 λ 1 = 3 i , λ 2 = − 3 i , λ 3 = 0
A basis of the corresponding eigenspaces
v 1 = ( 1 − 3 i , 1 + 3 i , 4 ) , v 2 = ( 1 + 3 i , 1 − 3 i , 4 ) , v 3 = ( − 2 , − 2 , 1 ) v_1=(1-3i,1+3i,4), v_2=(1+3i,1-3i,4), v_3=(-2,-2,1) v 1 = ( 1 − 3 i , 1 + 3 i , 4 ) , v 2 = ( 1 + 3 i , 1 − 3 i , 4 ) , v 3 = ( − 2 , − 2 , 1 )
(4)
Let λ \lambda λ be an eigenvalue of A A A and let x \mathbf{x} x be an eigenvector corresponding to the eigenvalue λ \lambda λ . That is, we have
A x = λ x . A\mathbf{x}=\lambda \mathbf{x}. A x = λ x .
Multiplying by x ˉ T \bar{\mathbf{x}}^{T} x ˉ T from the left, we have
x ˉ T A x = λ x ˉ T x = λ ∣ ∣ x ∣ ∣ 2 . \begin{align}
\bar{\mathbf{x}}^{T}A\mathbf{x}=\lambda \bar{\mathbf{x}}^{T} \mathbf{x}=\lambda ||\mathbf{x}||^2. \tag{*}
\end{align} x ˉ T A x = λ x ˉ T x = λ ∣∣ x ∣ ∣ 2 . ( * )
Note that the left hand side x ˉ T A x \bar{\mathbf{x}}^{T}A\mathbf{x} x ˉ T A x is the dot (inner) product of x ˉ \bar{\mathbf{x}} x ˉ and A x A\mathbf{x} A x . Since the dot product is commutative, we have
The left hand side of (*) = x ˉ T A x = ( A x ) T x ˉ = x T A T x ˉ . \begin{aligned}
&\text{The left hand side of (*)}\\
&=\bar{\mathbf{x}}^{T}A\mathbf{x}=(A\mathbf{x})^{T}\bar{\mathbf{x}}\\
&=x^{T}A^{T}\bar{\mathbf{x}}.
\end{aligned} The left hand side of (*) = x ˉ T A x = ( A x ) T x ˉ = x T A T x ˉ .
Since A A A is skew-symmetric, we have A T = − A A^{T}=-A A T = − A . Substituting this into the above equality, we have
The left hand side of (*) = x T A T x ˉ = − x T A x ˉ \begin{aligned}
&\text{The left hand side of (*)}\\
&=x^{T}A^{T}\bar{\mathbf{x}}=-\mathbf{x}^{T}A\bar{\mathbf{x}}
\end{aligned} The left hand side of (*) = x T A T x ˉ = − x T A x ˉ
Taking conjugate of A x = λ x A\mathbf{x}=\lambda\mathbf{x} A x = λ x and use the fact that A A A is real, we have
A x ˉ = λ ˉ x ˉ . A\bar{\mathbf{x}}=\bar{\lambda}\bar{\mathbf{x}}. A x ˉ = λ ˉ x ˉ .
Thus, we have
The left hand side of (*) = − x T A x ˉ = − x T λ ˉ x ˉ = − λ ˉ ∣ ∣ x ∣ ∣ 2 . \begin{aligned}
&\text{The left hand side of (*)}\\
&=-\mathbf{x}^{T}A\bar{\mathbf{x}}\\
&=-\mathbf{x}^{T}\bar{\lambda}\bar{\mathbf{x}}=-\bar{\lambda}||\mathbf{x}||^2.
\end{aligned} The left hand side of (*) = − x T A x ˉ = − x T λ ˉ x ˉ = − λ ˉ ∣∣ x ∣ ∣ 2 .
Therefore comparing the left and right hand sides of (*) yields
− λ ˉ ∣ ∣ x ∣ ∣ 2 = λ ∣ ∣ x ∣ ∣ 2 . -\bar{\lambda}||\mathbf{x}||^2=\lambda ||\mathbf{x}||^2. − λ ˉ ∣∣ x ∣ ∣ 2 = λ ∣∣ x ∣ ∣ 2 .
Since x \mathbf{x} x is an eigenvector, it is nonzero by definition. Thus ∣ ∣ x ∣ ∣ ≠ 0 ||\mathbf{x}||\neq 0 ∣∣ x ∣∣ = 0 .
Hence we have
− λ ˉ = λ , -\bar{\lambda}=\lambda, − λ ˉ = λ ,
and this implies that λ \lambda λ is either 0 0 0 or purely imaginary number.