跳到主要内容

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

Author

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

Description

問1

(i) 複素行列 H=(hij)H=(h_{ij}) がエルミート行列、すなわち hij=hji\overline{h_{ij}}=h_{ji} なら、固有値はすべて実数であることを証明せよ。横線は複素共役を表す。

(ii) 次の行列の nnAn,BnA^n,B^n を求め、その結果を証明せよ。nn は自然数、a,ba,b は実数である。

A=(1a01),B=(1b0001010000001b0000010001001b000001).A=\begin{pmatrix}1&a\\0&1\end{pmatrix},\qquad B=\begin{pmatrix}1&b&0&0&0&1\\0&1&0&0&0&0\\0&0&1&b&0&0\\0&0&0&1&0&0\\0&1&0&0&1&b\\0&0&0&0&0&1\end{pmatrix}.

問2

ある都市の8月の気温を長年調べたところ、翌日が真夏日になるかは当日だけに依存し、当日が真夏日なら翌日も真夏日となる確率は 2/32/3、当日が真夏日でなければ翌日も真夏日でない確率は 3/43/4 であった。8月 kk 日が真夏日である確率を hkh_k、そうでない確率を ckc_k とすると

(hk+1ck+1)=A(hkck)\begin{pmatrix}h_{k+1}\\c_{k+1}\end{pmatrix}=A\begin{pmatrix}h_k\\c_k\end{pmatrix}

と書ける。

(i) AA とその固有値を求めよ。

(ii) 正則行列 PP を用いて D=P1APD=P^{-1}AP と対角化できることを示せ。

(iii) 8月1日が真夏日であったとして、8月 (d+1)(d+1) 日(d=0,,30d=0,\ldots,30)が真夏日となる確率を求めよ。

题目描述

问1 (i) 若复矩阵 H=(hij)H=(h_{ij}) 为 Hermite 矩阵,即 hij=hji\overline{h_{ij}}=h_{ji},证明其全部特征值都是实数。(ii) 求并证明下列矩阵的 nn 次幂;nn 为自然数,a,ba,b 为实数。

A=(1a01),B=(1b0001010000001b0000010001001b000001).A=\begin{pmatrix}1&a\\0&1\end{pmatrix},\qquad B=\begin{pmatrix}1&b&0&0&0&1\\0&1&0&0&0&0\\0&0&1&b&0&0\\0&0&0&1&0&0\\0&1&0&0&1&b\\0&0&0&0&0&1\end{pmatrix}.

问2 某城市8月每日是否为盛夏日仅依赖前一天。若当天是盛夏日,次日仍为盛夏日的概率为 2/32/3;若当天不是,次日也不是的概率为 3/43/4。以 hk,ckh_k,c_k 表示8月 kk 日为、非盛夏日的概率,满足 (hk+1,ck+1)T=A(hk,ck)T(h_{k+1},c_{k+1})^{\mathrm T}=A(h_k,c_k)^{\mathrm T}。(i) 求 AA 及特征值。(ii) 给出可逆矩阵 PP 使 P1APP^{-1}AP 为对角阵。(iii) 若8月1日为盛夏日,求8月 d+1d+1 日(d=0,,30d=0,\ldots,30)为盛夏日的概率。

Kai

問1

(i) Hv=λvHv=\lambda vv0v\ne0 とする。H=HH^*=H より vHvv^*Hv は実数である。従って λ=vHv/(vv)R\lambda=v^*Hv/(v^*v)\in\mathbb R

(ii) A=I+NA,B=I+NBA=I+N_A,B=I+N_B とすると、直接乗算により NA2=NB2=0N_A^2=N_B^2=0。二項定理から

An=(1na01),Bn=(1nb000n010000001nb000001000n001nb000001).\boxed{A^n=\begin{pmatrix}1&na\\0&1\end{pmatrix}},\qquad \boxed{B^n=\begin{pmatrix}1&nb&0&0&0&n\\0&1&0&0&0&0\\0&0&1&nb&0&0\\0&0&0&1&0&0\\0&n&0&0&1&nb\\0&0&0&0&0&1\end{pmatrix}}.

問2

(i) 条件付き確率より

A=(2/31/41/33/4),det(λIA)=(λ1)(λ5/12).\boxed{A=\begin{pmatrix}2/3&1/4\\1/3&3/4\end{pmatrix}},\qquad \det(\lambda I-A)=(\lambda-1)(\lambda-5/12).

従って固有値は 1,5/12\boxed{1,5/12}

(ii) 対応する固有ベクトルとして (3,4)T,(1,1)T(3,4)^{\mathrm T},(1,-1)^{\mathrm T} をとれる。よって

P=(3141),P1AP=diag(1,5/12).\boxed{P=\begin{pmatrix}3&1\\4&-1\end{pmatrix}},\qquad \boxed{P^{-1}AP=\operatorname{diag}(1,5/12)}.

(iii) 初期値は (h1,c1)=(1,0)(h_1,c_1)=(1,0)Ad=Pdiag(1,(5/12)d)P1A^d=P\operatorname{diag}(1,(5/12)^d)P^{-1} より

hd+1=37+47(512)d.\boxed{h_{d+1}=\frac37+\frac47\left(\frac5{12}\right)^d}.