跳到主要内容

京都大学 情報学研究科 知能情報学専攻 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

であることを示せ。

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