跳到主要内容

東京大学 理学系研究科 天文学専攻 2022年8月実施 数学

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

以下、行列・ベクトルは実数とし、ATA^T は転置を表す。

  1. 実対称行列の異なる固有値に属する固有ベクトルが直交することを示せ。
  2. 対称行列 AA が相異なる固有値 λ1,,λn\lambda_1,\ldots,\lambda_n と単位固有ベクトル u1,,unu_1,\ldots,u_n をもつ。U=(u1,,un)U=(u_1,\ldots,u_n) として次を示せ。
    1. UT=U1U^T=U^{-1}
    2. A=Udiag(λ1,,λn)UTA=U\operatorname{diag}(\lambda_1,\ldots,\lambda_n)U^T
    3. xTAx>0x^TAx>0 が全ての x0x\ne0 で成り立つなら、全固有値は正である。
  3. PRn×nP\in\mathbb R^{n\times n}M=PPTM=PP^TN=PTPN=P^TP とする。
    1. M,NM,N の固有値が一致することを示せ。
    2. PP が零行列でなければ固有値が正であることを示せ。
    3. 同じ固有値 λ\lambda に属する M,NM,N の単位固有ベクトル u,vu,v の関係を示せ。
    4. P=Udiag(λ1,,λn)VTP=U\operatorname{diag}(\sqrt{\lambda_1},\ldots,\sqrt{\lambda_n})V^T を示せ。
  4. 512×512512\times512 画像行列を上の分解で階数 rr に近似する。各実数を4バイトとするとき、r=100,10r=100,10 のファイルサイズと元画像に対する割合を求めよ。

Kai

問1

Au=λu,Av=μvAu=\lambda u,Av=\mu v とする。A=ATA=A^T より

λuTv=(Au)Tv=uTAv=μuTv.\lambda u^Tv=(Au)^Tv=u^TAv=\mu u^Tv.

λμ\lambda\ne\mu なら (λμ)uTv=0(\lambda-\mu)u^Tv=0 なので、uTv=0\boxed{u^Tv=0}

問2

(a)

問1より uiTuj=δiju_i^Tu_j=\delta_{ij}。したがって

UTU=I,U1=UT.U^TU=I,\qquad\boxed{U^{-1}=U^T}.

(b)

AU=UΛAU=U\LambdaΛ=diag(λ1,,λn)\Lambda=\operatorname{diag}(\lambda_1,\ldots,\lambda_n) より

A=UΛUT.\boxed{A=U\Lambda U^T}.

(c)

ii について

0<uiTAui=λiuiTui=λi.0<u_i^TAu_i=\lambda_i u_i^Tu_i=\lambda_i.

よって λi>0 (i=1,,n)\boxed{\lambda_i>0\ (i=1,\ldots,n)}

問3

(a)

Mu=λuMu=\lambda uλ0\lambda\ne0 なら

N(PTu)=PTMu=λPTu,N(P^Tu)=P^TMu=\lambda P^Tu,

かつ PTu0P^Tu\ne0。逆も同様で、非零固有値とその重複度は一致する。また

dimkerM=dimkerPT=nrankP=dimkerN.\dim\ker M=\dim\ker P^T=n-\operatorname{rank}P =\dim\ker N.

したがって零固有値の重複度も一致する。

(b)

任意の xx に対し

xTMx=PTx20,xTNx=Px20.x^TMx=\|P^Tx\|^2\ge0,\qquad x^TNx=\|Px\|^2\ge0.

よって固有値は非負だが、この命題は一般には成り立たない。例えば

P=(1000)P=\begin{pmatrix}1&0\\0&0\end{pmatrix}

は零行列でないが、M=N=diag(1,0)M=N=\operatorname{diag}(1,0) である。PP が正則なら全固有値は正となる。

(c)

Nv=λvNv=\lambda vλ>0\lambda>0 に対して

u=Pvλu=\frac{Pv}{\sqrt\lambda}

とおけば u=1\|u\|=1Mu=λuMu=\lambda u であり、

u=Pvλ,v=PTuλ.\boxed{u=\frac{Pv}{\sqrt\lambda},\qquad v=\frac{P^Tu}{\sqrt\lambda}}.

固有空間が1次元なら、単位固有ベクトルの向きにより両式に同じ符号 ±\pm が付く。

(d)

λi>0\lambda_i>0 では問(c)の対応を用い、λi=0\lambda_i=0 では kerP,kerPT\ker P,\ker P^T の正規直交基底を選ぶ。このとき Pvi=λiuiPv_i=\sqrt{\lambda_i}u_i だから

PV=UΣ,Σ=diag(λ1,,λn).PV=U\Sigma,\qquad \Sigma=\operatorname{diag}(\sqrt{\lambda_1},\ldots,\sqrt{\lambda_n}).

V1=VTV^{-1}=V^T より

P=UΣVT.\boxed{P=U\Sigma V^T}.

問4

Ur,VrU_r,V_r は各 512r512r 個、対角成分は rr 個の実数をもつ。したがって

Sr=4(512r+r+512r)=4100r bytes.S_r=4(512r+r+512r)=4100r\ \text{bytes}.

元画像は 51224=1,048,576512^2\cdot4=1{,}048{,}576 bytes なので、

rSrSr/1,048,576100410,000 bytes39.1%40%1041,000 bytes3.91%4%\begin{array}{c|c|c} r&S_r&S_r/1{,}048{,}576\\ \hline 100&410{,}000\ \text{bytes}&39.1\%\simeq\boxed{40\%}\\ 10&41{,}000\ \text{bytes}&3.91\%\simeq\boxed{4\%} \end{array}

Reference