京都大学 情報学研究科 知能情報学専攻 2021年8月実施 情報学基礎 F1-1
Author
Isidore
Description
設問1
以下で定義される実行列 A A A と B B B 、および実数ベクトル x x x に関して、以下の問いに答えよ。
ここで、 ∥ x ∥ \|x\| ∥ x ∥ は x x x の長さを表す。
A = ( 1 − 3 3 1 ) , B = ( r 0 0 r ) , x = ( x 1 x 2 ) , ∥ x ∥ > 0 A =
\begin{pmatrix}
1 & -\sqrt{3} \\
\sqrt{3} & 1
\end{pmatrix}, \quad
B =
\begin{pmatrix}
r & 0 \\
0 & r
\end{pmatrix}, \quad
x =
\begin{pmatrix}
x_1 \\
x_2
\end{pmatrix}, \quad
\|x\| > 0 A = ( 1 3 − 3 1 ) , B = ( r 0 0 r ) , x = ( x 1 x 2 ) , ∥ x ∥ > 0
(1) A 3 A^3 A 3 を求めよ。
(2) A − 1 A^{-1} A − 1 を求めよ。
(3) A 15 A^{15} A 15 を求めよ。
(4) lim n → ∞ ∥ ( A B ) n x ∥ \lim_{n \to \infty} \|(AB)^n x\| lim n → ∞ ∥ ( A B ) n x ∥ を求めよ。ここで、 n n n は自然数とする。
設問2
m × n m \times n m × n ( m > n m > n m > n ) の実行列 A A A と以下で定義される行列 B B B と C C C を考える。ここで、 T T T は転置を表すものとする。
また、内積 ( a , b ) (a, b) ( a , b ) を以下で定義する。
ここで、 a a a と b b b は l l l 次元ベクトル、 a k a_k a k は a a a の k k k 番目の要素、 b k b_k b k は b b b の k k k 番目の要素とする。
B = A T A , C = A A T , ( a , b ) = ∑ k = 1 l a k b k B = A^T A, \quad C = A A^T, \quad (a, b) = \sum_{k=1}^l a_k b_k B = A T A , C = A A T , ( a , b ) = k = 1 ∑ l a k b k
以下の問いに答えよ。必要であれば、実対称行列に関する以下の性質を使用せよ:
実対称行列の固有値は全て実数である。
実対称行列のどの固有値に対しても、実ベクトルからなる固有ベクトルをとることができる。
(1) B B B と C C C がともに実対称行列であることを示せ。
(2) B B B の全ての固有値が非負であることを示せ。
(3) 実対称行列は直交行列で対角化できる。 B B B を対角化する直交行列の列ベクトルは、いずれも B B B の固有ベクトルである。
これらのうち正の固有値 λ i , λ j \lambda_i, \lambda_j λ i , λ j に対応する固有ベクトルを p i , p j p_i, p_j p i , p j とし、 q i q_i q i と q j q_j q j を
q i = 1 λ i A p i , q j = 1 λ j A p j q_i = \frac{1}{\sqrt{\lambda_i}} A p_i, \ q_j = \frac{1}{\sqrt{\lambda_j}} A p_j q i = λ i 1 A p i , q j = λ j 1 A p j
で与えられるベクトルとする。
q i , q j q_i, q_j q i , q j が C C C の固有ベクトルであって対応する固有値はそれぞれ λ i , λ j \lambda_i, \lambda_j λ i , λ j であること、および
( q i , q j ) = { 0 ( i ≠ j ) 1 ( i = j ) (q_i, q_j) =
\begin{cases}
0 & (i \neq j) \\
1 & (i = j)
\end{cases} ( q i , q j ) = { 0 1 ( i = j ) ( i = j )
であることを示せ。また、
p i = 1 λ i A T q i p_i = \frac{1}{\sqrt{\lambda_i}} A^T q_i p i = λ i 1 A T q i
であることを示せ。
题目描述
回答以下两题。
给定实数 r r r 、非零实向量 x \boldsymbol{x} x 以及实矩阵
A = ( 1 − 3 3 1 ) , B = ( r 0 0 r ) , x = ( x 1 x 2 ) , ∥ x ∥ > 0 , A=
\begin{pmatrix}
1&-\sqrt{3}\\
\sqrt{3}&1
\end{pmatrix},
\qquad
B=
\begin{pmatrix}
r&0\\
0&r
\end{pmatrix},
\qquad
\boldsymbol{x}=
\begin{pmatrix}
x_1\\
x_2
\end{pmatrix},
\qquad
\|\boldsymbol{x}\|>0, A = ( 1 3 − 3 1 ) , B = ( r 0 0 r ) , x = ( x 1 x 2 ) , ∥ x ∥ > 0 ,
其中 ∥ x ∥ \|\boldsymbol{x}\| ∥ x ∥ 表示向量长度。完成:
求 A 3 A^3 A 3 。
求 A − 1 A^{-1} A − 1 。
求 A 15 A^{15} A 15 。
当自然数 n → ∞ n\to\infty n → ∞ 时,求
∥ ( A B ) n x ∥ \|(AB)^n\boldsymbol{x}\| ∥ ( A B ) n x ∥ 的极限;答案应按参数 r r r 的取值讨论,包括极限不存在或发散的情形。
设 A A A 是 m × n m\times n m × n 实矩阵,且 m > n m>n m > n 。定义
B = A ⊤ A , C = A A ⊤ , ( a , b ) = ∑ k = 1 l a k b k , B=A^\top A,\qquad
C=AA^\top,\qquad
(\boldsymbol{a},\boldsymbol{b})
=\sum_{k=1}^{l}a_kb_k, B = A ⊤ A , C = A A ⊤ , ( a , b ) = k = 1 ∑ l a k b k ,
其中 ⊤ \top ⊤ 表示转置,最后一个公式定义 l l l 维向量
a , b \boldsymbol{a},\boldsymbol{b} a , b 的内积。必要时可以使用以下事实:实对称矩阵的特征值全为实数;对它的任一特征值,都可以选取实特征向量。
证明 B B B 和 C C C 都是实对称矩阵。
证明 B B B 的所有特征值均非负。
实对称矩阵可以用正交矩阵对角化。设该正交矩阵中对应于 B B B 的正特征值
λ i , λ j \lambda_i,\lambda_j λ i , λ j 的列特征向量分别为
p i , p j \boldsymbol{p}_i,\boldsymbol{p}_j p i , p j ,并定义
q i = 1 λ i A p i , q j = 1 λ j A p j . \boldsymbol{q}_i
=\frac{1}{\sqrt{\lambda_i}}A\boldsymbol{p}_i,
\qquad
\boldsymbol{q}_j
=\frac{1}{\sqrt{\lambda_j}}A\boldsymbol{p}_j. q i = λ i 1 A p i , q j = λ j 1 A p j .
证明 q i , q j \boldsymbol{q}_i,\boldsymbol{q}_j q i , q j 是 C C C 的特征向量,对应特征值分别为
λ i , λ j \lambda_i,\lambda_j λ i , λ j ;证明
( q i , q j ) = { 0 , i ≠ j , 1 , i = j ; (\boldsymbol{q}_i,\boldsymbol{q}_j)
=
\begin{cases}
0,&i\ne j,\\
1,&i=j;
\end{cases} ( q i , q j ) = { 0 , 1 , i = j , i = j ;
并证明
p i = 1 λ i A ⊤ q i . \boldsymbol{p}_i
=\frac{1}{\sqrt{\lambda_i}}A^\top\boldsymbol{q}_i. p i = λ i 1 A ⊤ q i .
Kai
設問1
(1)
A 3 = ( − 8 0 0 − 8 ) A^{3} = \begin{pmatrix}
-8 & 0 \\
0 & -8 \\
\end{pmatrix} A 3 = ( − 8 0 0 − 8 )
(2)
A − 1 = 1 4 ( 1 3 − 3 1 ) A^{-1} = \frac{1}{4} \begin{pmatrix}
1 & \sqrt{3} \\
-\sqrt{3} & 1 \\
\end{pmatrix} A − 1 = 4 1 ( 1 − 3 3 1 )
(3)
A 15 = ( A 3 ) 5 = ( − 8 E ) 5 = − 2 15 E A^{15}=(A^3)^5=(-8E)^5=-2^{15}E A 15 = ( A 3 ) 5 = ( − 8 E ) 5 = − 2 15 E
(4)
Let the required value denote I I I , and we have
I = lim n → ∞ ∥ ( A B ) n x ∥ = lim n → ∞ ∥ ( − 8 r 3 E ) n 3 x ∥ I = \lim_{n \rightarrow \infty}\|(AB)^nx\| = \lim_{n \rightarrow \infty}\|(-8r^3E)^{\frac{n}{3}}x\| I = n → ∞ lim ∥ ( A B ) n x ∥ = n → ∞ lim ∥ ( − 8 r 3 E ) 3 n x ∥
Obviously, the limit converges if and only if
− 8 r 3 ∈ [ − 1 , 1 ] -8r^3 \in [-1, 1] − 8 r 3 ∈ [ − 1 , 1 ]
Hence, when r ∈ ( − 1 2 , 1 2 ) r \in (-\frac{1}{2}, \frac{1}{2}) r ∈ ( − 2 1 , 2 1 ) , we have
Particularly, when r = ± 1 2 r=\pm \frac{1}{2} r = ± 2 1 , we have
I = lim n → ∞ ∥ ( − 1 ) n 3 x ∥ = ∥ x ∥ = x 1 2 + x 2 2 I = \lim_{n \rightarrow \infty}\|(-1)^{\frac{n}{3}}x\| = \|x\| = \sqrt{x_1^2+x_2^2} I = n → ∞ lim ∥ ( − 1 ) 3 n x ∥ = ∥ x ∥ = x 1 2 + x 2 2
設問2
(1)
Obviously, as A A A is a real matrix, B B B and C C C are both real matrices. Also,
B T = ( A T A ) T = A T ( A T ) T = A T A = B C T = ( A A T ) T = ( A T ) T A T = A A T = C \begin{aligned}
B^T
&= (A^T A)^T
= A^T (A^T)^T
= A^T A = B\\
C^T
&= (A A^T)^T
= (A^T)^T A^T
= A A^T = C
\end{aligned} B T C T = ( A T A ) T = A T ( A T ) T = A T A = B = ( A A T ) T = ( A T ) T A T = A A T = C
hence B B B and C C C are real symmetric matrices.
(2)
Consider B's eigenvalue, we have
A T A x = λ x A^TAx=\lambda x A T A x = λ x
hence,
x T A T A x = λ x T x x^TA^TAx = \lambda x^Tx x T A T A x = λ x T x
hence,
λ = x T A T A x x T x = ∥ A x ∥ ∥ x ∥ ≥ 0 \lambda = \frac{x^TA^TAx}{x^Tx} = \frac{\|Ax\|}{\|x\|} \geq 0 λ = x T x x T A T A x = ∥ x ∥ ∥ A x ∥ ≥ 0
(3)
Consider the eigenvalue λ i \lambda_i λ i with its eigenvector p i p_i p i , we have
A T A p i = λ i p i A^TAp_i=\lambda_ip_i A T A p i = λ i p i
Left-multiply by 1 λ i A \frac{1}{\sqrt{\lambda_i}}A λ i 1 A , we have
1 λ i A A T A p i = λ i A p i A A T ( 1 λ i A p i ) = λ i ( 1 λ i A p i ) \frac{1}{\sqrt{\lambda_i}}AA^TAp_i = \sqrt{\lambda_i}Ap_i \\
AA^T(\frac{1}{\sqrt{\lambda_i}}Ap_i) = \lambda_i(\frac{1}{\sqrt{\lambda_i}}Ap_i) λ i 1 A A T A p i = λ i A p i A A T ( λ i 1 A p i ) = λ i ( λ i 1 A p i )
since q i = 1 λ i A p i q_i = \frac{1}{\sqrt{\lambda_i}}Ap_i q i = λ i 1 A p i , we have
C q i = λ i q i Cq_i = \lambda_i q_i C q i = λ i q i
Similarly, we have C q j = λ j q j Cq_j = \lambda_j q_j C q j = λ j q j for a different eigenvalue λ j \lambda_j λ j . Therefore, q i q_i q i and q j q_j q j are eigenvectors of C corresponding to eigenvalues λ i \lambda_i λ i and λ j \lambda_j λ j .
Now we consider the inner product of q i q_i q i and q j q_j q j
q i T q j = 1 λ i λ j ( A p i ) T ( A p j ) = λ j λ i p i T p j q_i^Tq_j = \frac{1}{\sqrt{\lambda_i \lambda_j}}(Ap_i)^T(Ap_j) = \frac{\sqrt{\lambda_j}}{\sqrt{\lambda_i}}p_i^Tp_j q i T q j = λ i λ j 1 ( A p i ) T ( A p j ) = λ i λ j p i T p j
Obviously, if i ≠ j i \neq j i = j , as B B B is a real symmetric matrix, p i T p j = 0 p_i^Tp_j = 0 p i T p j = 0 . If i = j i = j i = j , p i T p j = 1 p_i^Tp_j = 1 p i T p j = 1
Therefore, we have
q i T q j = { 0 if i ≠ j 1 if i = j q_i^Tq_j = \begin{cases}
0 & \text{if } i \neq j \\
1 & \text{if } i = j
\end{cases} q i T q j = { 0 1 if i = j if i = j
Finally, insert q i = 1 λ i A p i q_i = \frac{1}{\sqrt{\lambda_i}}Ap_i q i = λ i 1 A p i to 1 λ i A T q i \frac{1}{\sqrt{\lambda_i}}A^Tq_i λ i 1 A T q i , we immediately have
1 λ i A T q i = 1 λ i A T 1 λ i A p i = 1 λ i A T A p i = p i \frac{1}{\sqrt{\lambda_i}}A^Tq_i = \frac{1}{\sqrt{\lambda_i}}A^T\frac{1}{\sqrt{\lambda_i}}Ap_i = \frac{1}{\lambda_i}A^TAp_i = p_i λ i 1 A T q i = λ i 1 A T λ i 1 A p i = λ i 1 A T A p i = p i
Q.E.D