跳到主要内容

京都大学 情報学研究科 知能情報学専攻 2021年7月実施 情報学基礎 F1-1

Author​

Isidore, 祭音Myyura

Description​

大学公表の原題

設問1​

以下で定義される実行列 AA と BB 、および実数ベクトル xx に関して、以下の問いに答えよ。 ここで、 ∥x∥\|x\| は xx の長さを表す。

A=(1−331),B=(r00r),x=(x1x2),∥x∥>0A = \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

(1) A3A^3 を求めよ。

(2) A−1A^{-1} を求めよ。

(3) A15A^{15} を求めよ。

(4) lim⁡n→∞∥(AB)nx∥\lim_{n \to \infty} \|(AB)^n x\| を求めよ。ここで、 nn は自然数とする。

設問2​

m×nm \times n ( m>nm > n ) の実行列 AA と以下で定義される行列 BB と CC を考える。ここで、 TT は転置を表すものとする。 また、内積 (a,b)(a, b) を以下で定義する。 ここで、 aa と bb は ll 次元ベクトル、 aka_k は aa の kk 番目の要素、 bkb_k は bb の kk 番目の要素とする。

B=ATA,C=AAT,(a,b)=∑k=1lakbkB = A^T A, \quad C = A A^T, \quad (a, b) = \sum_{k=1}^l a_k b_k

以下の問いに答えよ。必要であれば、実対称行列に関する以下の性質を使用せよ:

  • 実対称行列の固有値は全て実数である。
  • 実対称行列のどの固有値に対しても、実ベクトルからなる固有ベクトルをとることができる。

(1) BB と CC がともに実対称行列であることを示せ。

(2) BB の全ての固有値が非負であることを示せ。

(3) 実対称行列は直交行列で対角化できる。 BB を対角化する直交行列の列ベクトルは、いずれも BB の固有ベクトルである。 これらのうち正の固有値 λi,λj\lambda_i, \lambda_j に対応する固有ベクトルを pi,pjp_i, p_j とし、 qiq_i と qjq_j を

qi=1λiApi, qj=1λjApjq_i = \frac{1}{\sqrt{\lambda_i}} A p_i, \ q_j = \frac{1}{\sqrt{\lambda_j}} A p_j

で与えられるベクトルとする。 qi,qjq_i, q_j が CC の固有ベクトルであって対応する固有値はそれぞれ λi,λj\lambda_i, \lambda_j であること、および

(qi,qj)={0(i≠j)1(i=j) (q_i, q_j) = \begin{cases} 0 & (i \neq j) \\ 1 & (i = j) \end{cases}

であることを示せ。また、

pi=1λiATqip_i = \frac{1}{\sqrt{\lambda_i}} A^T q_i

であることを示せ。

题目描述​

回答以下两题。

  1. 给定实数 rr、非零实向量 x\boldsymbol{x} 以及实矩阵
A=(1−331),B=(r00r),x=(x1x2),∥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,

其中 ∥x∥\|\boldsymbol{x}\| 表示向量长度。完成:

  1. 求 A3A^3。

  2. 求 A−1A^{-1}。

  3. 求 A15A^{15}。

  4. 当自然数 n→∞n\to\infty 时,求 ∥(AB)nx∥\|(AB)^n\boldsymbol{x}\| 的极限;答案应按参数 rr 的取值讨论,包括极限不存在或发散的情形。

  5. 设 AA 是 m×nm\times n 实矩阵,且 m>nm>n。定义

B=A⊤A,C=AA⊤,(a,b)=∑k=1lakbk,B=A^\top A,\qquad C=AA^\top,\qquad (\boldsymbol{a},\boldsymbol{b}) =\sum_{k=1}^{l}a_kb_k,

其中 ⊤\top 表示转置,最后一个公式定义 ll 维向量 a,b\boldsymbol{a},\boldsymbol{b} 的内积。必要时可以使用以下事实:实对称矩阵的特征值全为实数;对它的任一特征值,都可以选取实特征向量。

  1. 证明 BB 和 CC 都是实对称矩阵。
  2. 证明 BB 的所有特征值均非负。
  3. 实对称矩阵可以用正交矩阵对角化。设该正交矩阵中对应于 BB 的正特征值 λi,λj\lambda_i,\lambda_j 的列特征向量分别为 pi,pj\boldsymbol{p}_i,\boldsymbol{p}_j,并定义
qi=1λiApi,qj=1λjApj.\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.

证明 qi,qj\boldsymbol{q}_i,\boldsymbol{q}_j 是 CC 的特征向量,对应特征值分别为 λi,λj\lambda_i,\lambda_j;证明

(qi,qj)={0,i≠j,1,i=j;(\boldsymbol{q}_i,\boldsymbol{q}_j) = \begin{cases} 0,&i\ne j,\\ 1,&i=j; \end{cases}

并证明

pi=1λiA⊤qi.\boldsymbol{p}_i =\frac{1}{\sqrt{\lambda_i}}A^\top\boldsymbol{q}_i.

Kai​

設問1​

(1)​

A3=(−800−8)A^{3} = \begin{pmatrix} -8 & 0 \\ 0 & -8 \\ \end{pmatrix}

(2)​

A−1=14(13−31)A^{-1} = \frac{1}{4} \begin{pmatrix} 1 & \sqrt{3} \\ -\sqrt{3} & 1 \\ \end{pmatrix}

(3)​

A15=(A3)5=(−8E)5=−215EA^{15}=(A^3)^5=(-8E)^5=-2^{15}E

(4)​

Since A=2RA=2R, where RR is an orthogonal rotation matrix, and B=rEB=rE,

∥(AB)nx∥=∣2r∣n∥x∥.\|(AB)^n x\|=|2r|^n\|x\|.

Therefore,

lim⁡n→∞∥(AB)nx∥={0,∣r∣<12,∥x∥=x12+x22,∣r∣=12,+∞,∣r∣>12.\lim_{n\to\infty}\|(AB)^n x\| = \begin{cases} 0,& |r|<\dfrac12,\\ \|x\|=\sqrt{x_1^2+x_2^2},& |r|=\dfrac12,\\ +\infty,& |r|>\dfrac12. \end{cases}

設問2​

(1)​

Obviously, as AA is a real matrix, BB and CC are both real matrices. Also,

BT=(ATA)T=AT(AT)T=ATA=BCT=(AAT)T=(AT)TAT=AAT=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}

hence BB and CC are real symmetric matrices.

(2)​

Consider B's eigenvalue, we have

ATAx=λxA^TAx=\lambda x

hence,

xTATAx=λxTxx^TA^TAx = \lambda x^Tx

hence,

λ=xTATAxxTx=∥Ax∥2∥x∥2≥0\lambda = \frac{x^TA^TAx}{x^Tx} = \frac{\|Ax\|^2}{\|x\|^2} \geq 0

(3)​

Consider the eigenvalue λi\lambda_i with its eigenvector pip_i , we have

ATApi=λipiA^TAp_i=\lambda_ip_i

Left-multiply by 1λiA\frac{1}{\sqrt{\lambda_i}}A , we have

1λiAATApi=λiApiAAT(1λiApi)=λi(1λiApi)\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)

since qi=1λiApiq_i = \frac{1}{\sqrt{\lambda_i}}Ap_i , we have

Cqi=λiqiCq_i = \lambda_i q_i

Similarly, we have Cqj=λjqjCq_j = \lambda_j q_j . Therefore, qiq_i and qjq_j are eigenvectors of C corresponding to eigenvalues λi\lambda_i and λj\lambda_j .

Now we consider the inner product of qiq_i and qjq_j

qiTqj=1λiλj(Api)T(Apj)=λjλipiTpjq_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

If i≠ji \neq j , then piTpj=0p_i^Tp_j = 0 because pi,pjp_i,p_j are distinct columns of the orthogonal diagonalizing matrix. If i=ji = j , piTpj=1p_i^Tp_j = 1 Therefore, we have

qiTqj={0if i≠j1if i=jq_i^Tq_j = \begin{cases} 0 & \text{if } i \neq j \\ 1 & \text{if } i = j \end{cases}

Finally, insert qi=1λiApiq_i = \frac{1}{\sqrt{\lambda_i}}Ap_i to 1λiATqi\frac{1}{\sqrt{\lambda_i}}A^Tq_i , we immediately have

1λiATqi=1λiAT1λiApi=1λiATApi=pi\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

Q.E.D