跳到主要内容

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

Author​

犬 (finalized by 祭音Myyura with assistance from GPT 6 Astra)

Description​

問1​

n×nn\times n 実行列 An=(aij)A_n=(a_{ij}) の成分が aii=1a_{ii}=1、ai,i+1=−ba_{i,i+1}=-b、その他は 00 であるとき、An−1A_n^{-1} を求めよ。

問2​

A=I+abTA=I+ab^T とする。II は nn 次単位行列、a,ba,b は nn 次元実列ベクトルである。

(i) スカラー数列 {ck}\{c_k\} によって Ak=I+ckabTA^k=I+c_kab^T(k=1,2,…k=1,2,\ldots)と表せることを数学的帰納法で示せ。

(ii) bTa≠0b^Ta\ne0 とする。ckc_k を求め、実数 tt に対し exp⁡(At)=∑k=0∞Aktk/k!\exp(At)=\sum_{k=0}^\infty A^kt^k/k! を求めよ。

問3​

n×nn\times n 実対称行列 A=(aij),B=(bij)A=(a_{ij}),B=(b_{ij}) を考える。

(i) A=UDUTA=UDU^T と表せることを示せ。UU は直交行列、DD は対角行列である。すなわち aij=∑p=1nuipdpujpa_{ij}=\sum_{p=1}^n u_{ip}d_pu_{jp} と表される。

(ii) cij=aijbijc_{ij}=a_{ij}b_{ij} によって行列 CC を定義する。A,BA,B が正定値ならば CC も正定値であることを、xTCx=∑i,jcijxixjx^TCx=\sum_{i,j}c_{ij}x_ix_j を調べて証明せよ。

题目描述​

问1​

实矩阵 AnA_n 的主对角线为 11,上方紧邻的副对角线为 −b-b,其余为 00。求其逆矩阵。

问2​

设 A=I+abTA=I+ab^T,a,ba,b 为 nn 维实列向量。

(i) 用数学归纳法证明存在标量序列 ckc_k,使 Ak=I+ckabTA^k=I+c_kab^T。

(ii) 若 bTa≠0b^Ta\ne0,求 ckc_k 及 exp⁡(At)=∑k=0∞Aktk/k!\exp(At)=\sum_{k=0}^\infty A^kt^k/k!,其中 tt 为实数。

问3​

A,BA,B 为 nn 阶实对称矩阵。

(i) 证明 A=UDUTA=UDU^T,其中 UU 正交、DD 对角,即 aij=∑puipdpujpa_{ij}=\sum_pu_{ip}d_pu_{jp}。

(ii) 定义 cij=aijbijc_{ij}=a_{ij}b_{ij}。通过分析二次型 xTCxx^TCx,证明 A,BA,B 正定时 CC 也正定。

Kai​

問1​

上方の副対角成分だけが 11 の行列を SS とすると、Sn=0S^n=0、An=I−bSA_n=I-bS。したがって

An−1=∑k=0n−1bkSk,(An−1)ij={bj−ij≥i,0j<i.\boxed{A_n^{-1}=\sum_{k=0}^{n-1}b^kS^k},\qquad (A_n^{-1})_{ij}=\begin{cases}b^{j-i}&j\ge i,\\0&j

実際、(I−bS)∑k=0n−1bkSk=I−bnSn=I(I-bS)\sum_{k=0}^{n-1}b^kS^k=I-b^nS^n=I。

問2​

(i) c1=1c_1=1 とする。β=bTa\beta=b^Ta とおけば (abT)2=βabT(ab^T)^2=\beta ab^T であるから

(I+ckabT)(I+abT)=I+{1+(1+β)ck}abT.(I+c_kab^T)(I+ab^T)=I+\{1+(1+\beta)c_k\}ab^T.

したがって ck+1=1+(1+β)ckc_{k+1}=1+(1+\beta)c_k と定めれば帰納法が成立する。

(ii) 漸化式を解くと ck=((1+β)k−1)/βc_k=((1+\beta)^k-1)/\beta。c0=0c_0=0 として級数に代入すれば

eAt=etI+e(1+β)t−etβabT=et(I+eβt−1βabT).\boxed{e^{At}=e^tI+\frac{e^{(1+\beta)t}-e^t}{\beta}ab^T =e^t\left(I+\frac{e^{\beta t}-1}{\beta}ab^T\right)}.

問3​

(i) 単位球面上で xTAxx^TAx を最大にする単位ベクトル uu を選ぶ。ラグランジュの未定乗数法から Au=duAu=du。対称性より u⊥u^\perp は AA の不変部分空間である。同じ議論をその部分空間に帰納的に適用すれば、正規直交固有ベクトル u1,…,unu_1,\ldots,u_n が得られる。これらを列とする UU と D=diag⁡(d1,…,dn)D=\operatorname{diag}(d_1,\ldots,d_n) により A=UDUTA=UDU^T。

(ii) yp=(u1px1,…,unpxn)Ty_p=(u_{1p}x_1,\ldots,u_{np}x_n)^T とおくと

xTCx=∑p=1ndp ypTByp.x^TCx=\sum_{p=1}^nd_p\,y_p^TBy_p.

すべての dp>0d_p>0 であり、各項は非負。また x≠0x\ne0 ならば

∑p∥yp∥2=∑ixi2∑puip2=∥x∥2>0\sum_p\|y_p\|^2=\sum_i x_i^2\sum_pu_{ip}^2=\|x\|^2>0

なので少なくとも一つの yp≠0y_p\ne0。その項は正であり、xTCx>0x^TCx>0。よって CC は正定値である。