跳到主要内容

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

Author

Isidore

Description

設問1

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

A=(1331),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) A1A^{-1} を求めよ。

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

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

設問2

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

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) BBCC がともに実対称行列であることを示せ。

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

(3) 実対称行列は直交行列で対角化できる。 BB を対角化する直交行列の列ベクトルは、いずれも BB の固有ベクトルである。 これらのうち正の固有値 λi,λj\lambda_i, \lambda_j に対応する固有ベクトルを pi,pjp_i, p_j とし、 qiq_iqjq_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_jCC の固有ベクトルであって対応する固有値はそれぞれ λi,λj\lambda_i, \lambda_j であること、および

(qi,qj)={0(ij)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=(1331),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. A1A^{-1}

  3. A15A^{15}

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

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

B=AA,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. 证明 BBCC 都是实对称矩阵。
  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}_jCC 的特征向量,对应特征值分别为 λi,λj\lambda_i,\lambda_j;证明

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

并证明

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

Kai

設問1

(1)

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

(2)

A1=14(1331)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)

Let the required value denote II , and we have

I=limn(AB)nx=limn(8r3E)n3xI = \lim_{n \rightarrow \infty}\|(AB)^nx\| = \lim_{n \rightarrow \infty}\|(-8r^3E)^{\frac{n}{3}}x\|

Obviously, the limit converges if and only if

8r3[1,1]-8r^3 \in [-1, 1]

Hence, when r(12,12)r \in (-\frac{1}{2}, \frac{1}{2}) , we have

I=0I = 0

Particularly, when r=±12r=\pm \frac{1}{2} , we have

I=limn(1)n3x=x=x12+x22I = \lim_{n \rightarrow \infty}\|(-1)^{\frac{n}{3}}x\| = \|x\| = \sqrt{x_1^2+x_2^2}

設問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=Axx0\lambda = \frac{x^TA^TAx}{x^Tx} = \frac{\|Ax\|}{\|x\|} \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 for a different eigenvalue λj\lambda_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

Obviously, if iji \neq j , as BB is a real symmetric matrix, piTpj=0p_i^Tp_j = 0 . If i=ji = j , piTpj=1p_i^Tp_j = 1 Therefore, we have

qiTqj={0if ij1if 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