跳到主要内容

京都大学 情報学研究科 システム科学専攻 2026年8月実施 数学【I】

Author​

AnonTokyo #83056

Description​

以下の問題において,R\mathbb{R} は実数全体からなる集合を,行列 AA に対し ATA^\mathsf{T} は AA の転置 (transpose) を表す。

問題 1​

正方行列

A=(10−110−10−1−1)A=\begin{pmatrix} 1&0&-1\\ 1&0&-1\\ 0&-1&-1 \end{pmatrix}

に関する以下の設問に答えよ。

(1) AA の階数 (rank) を求めよ。

(2) AA を表現行列とする線形写像 (linear mapping) の核 (kernel) ker⁡A\ker A を求めよ。

(3) x∈R3x\in\mathbb{R}^3 を変数とする線形方程式 (linear equations)

Ax=(−1−1−1)Ax=\begin{pmatrix}-1\\-1\\-1\end{pmatrix}

の一般解 (general solution) を求めよ。解が存在しない場合は「存在しない」と答えよ。

問題 2​

d,nd,n を n>d>1n>d>1 なる整数とする。i=1,2,…,ni=1,2,\ldots,n に対して xi=(xi,1,xi,2,…,xi,d)Tx_i=(x_{i,1},x_{i,2},\ldots,x_{i,d})^\mathsf{T} とし,d×nd\times n 行列 XX を X=(x1,x2,…,xn)X=(x_1,x_2,\ldots,x_n) により表す。また,A=1n∑i=1nxixiTA=\frac{1}{n}\sum_{i=1}^{n}x_ix_i^\mathsf{T} とする。AA は相異なるゼロでない固有値 (eigenvalue) をもつとし,これらを降順に並べたものを λ1>λ2>⋯>λd\lambda_1>\lambda_2>\cdots>\lambda_d とする。また対応する固有ベクトル (eigenvector) を j=1,2,…,dj=1,2,\ldots,d に対してそれぞれ vj=(vj,1,vj,2,…,vj,d)Tv_j=(v_{j,1},v_{j,2},\ldots,v_{j,d})^\mathsf{T} とし,これらは ∥vj∥=1\|v_j\|=1 を満たすとする。ただし u=(u1,u2,…,ud)T∈Rdu=(u_1,u_2,\ldots,u_d)^\mathsf{T}\in\mathbb{R}^d に対して ∥u∥=uTu\|u\|=\sqrt{u^\mathsf{T}u} は uu のユークリッドノルムを表す。さらに m<dm<d なる正の整数 mm に対して d×md\times m 行列 VV を V=(v1,v2,…,vm)V=(v_1,v_2,\ldots,v_m) と定義する。以下の設問に答えよ。

(1) AA を n,Xn,X を用いて表せ。

(2) v1,v2,…,vmv_1,v_2,\ldots,v_m を基底 (basis) とする Rd\mathbb{R}^d の線形部分空間 (linear subspace) への xix_i の正射影 (orthogonal projection) を zi∈Rdz_i\in\mathbb{R}^d とする。ziz_i を V,xiV,x_i を用いて表せ。ただし,v1,v2,…,vdv_1,v_2,\ldots,v_d が互いに直交する (mutually orthogonal) ことは説明なしに用いてよい。

(3) 設問 (2) で定義した z1,z2,…,znz_1,z_2,\ldots,z_n について,射影誤差 (projection error)

E=1n∑i=1n∥zi−xi∥2E=\frac{1}{n}\sum_{i=1}^{n}\|z_i-x_i\|^2

を A,VA,V を用いて表せ。必要であれば単位行列 (identity matrix) II,行列式 (determinant) det⁡\det,トレース (trace) tr⁡\operatorname{tr} を用いてよい。

(4) 設問 (3) の射影誤差 EE を AA の固有値を用いて表せ。

Kai​

問題 1​

(1) 階数​

第 2 行は第 1 行と等しく,第 1 行と第 3 行は一次独立である。行基本変形により

A∼(10−10−1−1000)A\sim\begin{pmatrix}1&0&-1\\0&-1&-1\\0&0&0\end{pmatrix}

を得る。したがって,rank⁡(A)=2\boxed{\operatorname{rank}(A)=2} である。

(2) 核​

x=(x1,x2,x3)Tx=(x_1,x_2,x_3)^\mathsf{T} とおくと,Ax=0Ax=0 は

x1−x3=0,−x2−x3=0x_1-x_3=0,\qquad -x_2-x_3=0

に等しい。x3=tx_3=t とおけば x1=t, x2=−tx_1=t,\ x_2=-t であるから,

ker⁡A={t(1,−1,1)T∣t∈R}.\boxed{\ker A=\{t(1,-1,1)^\mathsf{T}\mid t\in\mathbb{R}\}}.

(3) 線形方程式の一般解​

Ax=(−1,−1,−1)TAx=(-1,-1,-1)^\mathsf{T} は

x1−x3=−1,−x2−x3=−1x_1-x_3=-1,\qquad -x_2-x_3=-1

に帰着する。x3=t∈Rx_3=t\in\mathbb{R} として,x1=t−1, x2=1−tx_1=t-1,\ x_2=1-t である。よって解は存在し,一般解は

x=(−1,1,0)T+t(1,−1,1)T(t∈R).\boxed{x=(-1,1,0)^\mathsf{T}+t(1,-1,1)^\mathsf{T}\quad(t\in\mathbb{R})}.

問題 2​

(1) 行列表示​

XX の各列は xix_i なので,XXT=∑i=1nxixiTXX^\mathsf{T}=\sum_{i=1}^n x_ix_i^\mathsf{T} となる。ゆえに

A=1nXXT.\boxed{A=\frac1n XX^\mathsf{T}}.

(2) 正射影​

AA は実対称行列であり,v1,…,vmv_1,\ldots,v_m は正規直交系をなす。従って VTV=ImV^\mathsf{T}V=I_m であり,射影行列は P=VVTP=VV^\mathsf{T} である。したがって

zi=VVTxi=∑j=1m(vjTxi)vj.\boxed{z_i=VV^\mathsf{T}x_i=\sum_{j=1}^m (v_j^\mathsf{T}x_i)v_j}.

(3) 射影誤差​

P2=PP^2=P,PT=PP^\mathsf{T}=P より,(I−P)T(I−P)=I−P(I-P)^\mathsf{T}(I-P)=I-P である。このこととトレースの性質を用いると

E=1n∑i=1nxiT(I−P)xi=tr⁡ ⁣[(I−P)1n∑i=1nxixiT]=tr⁡[(I−VVT)A].\begin{aligned} E&=\frac1n\sum_{i=1}^n x_i^\mathsf{T}(I-P)x_i\\ &=\operatorname{tr}\!\left[(I-P)\frac1n\sum_{i=1}^n x_ix_i^\mathsf{T}\right]\\ &=\operatorname{tr}[(I-VV^\mathsf{T})A]. \end{aligned}

従って

E=tr⁡(A)−tr⁡(VTAV).\boxed{E=\operatorname{tr}(A)-\operatorname{tr}(V^\mathsf{T}AV)}.

(4) 固有値による表示​

Avj=λjvjAv_j=\lambda_jv_j であるから,VTAV=diag⁡(λ1,…,λm)V^\mathsf{T}AV=\operatorname{diag}(\lambda_1,\ldots,\lambda_m) である。また,tr⁡(A)=∑j=1dλj\operatorname{tr}(A)=\sum_{j=1}^d\lambda_j なので

E=∑j=m+1dλj.\boxed{E=\sum_{j=m+1}^{d}\lambda_j}.