跳到主要内容

東京工業大学 工学院 システム制御系 2019年8月実施 筆答専門試験問題2

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

問1​

(1) aa を実数として以下の行列 AA を考える。

A=(01−a110a01)A=\begin{pmatrix}0&1&-a\\1&1&0\\a&0&1\end{pmatrix}

AA が正則にならない条件を述べよ。また、このときの AA の核空間(Ax=0A\boldsymbol x=\boldsymbol0 を満たすベクトル x\boldsymbol x のなす空間)の次元を求め、その理由を述べよ。

(2) 以下の行列 AA とベクトル x\boldsymbol x を考える。

A=(−111−1),x=(02)A=\begin{pmatrix}-1&1\\1&-1\end{pmatrix},\qquad \boldsymbol x=\begin{pmatrix}0\\2\end{pmatrix}

実数 t≥0t\ge0 に対してベクトル y(t)=exp⁡(At)x\boldsymbol y(t)=\exp(At)\boldsymbol x と lim⁡t→∞y(t)\lim_{t\to\infty}\boldsymbol y(t) を求めよ。

問2​

以下の行列 AA を考える。

A=(100011)A=\begin{pmatrix}1&0&0\\0&1&1\end{pmatrix}

(1) 行列 B=AATB=AA^{\mathsf T} の固有値と正規化した固有ベクトルの組をすべて求めよ。ここで T\mathsf T は転置を表す。

(2) 行列 C=ATAC=A^{\mathsf T}A の固有値と正規化した固有ベクトルの組をすべて求めよ。

(3) 行列 BB の固有ベクトルを適切に並べて構成した直交行列を UU、行列 CC の固有ベクトルを適切に並べて構成した直交行列を VV とすると、行列 AA を A=UWVTA=UWV^{\mathsf T} という形で表せる。このとき行列 WW を求めよ。

题目描述​

问1

(1) 设 aa 为实数,考虑矩阵

A=(01−a110a01).A=\begin{pmatrix}0&1&-a\\1&1&0\\a&0&1\end{pmatrix}.

写出 AA 不可逆的条件。求此时 AA 的核空间,即满足 Ax=0A\boldsymbol x=\boldsymbol0 的向量 x\boldsymbol x 所构成空间的维数,并说明理由。

(2) 考虑矩阵与向量

A=(−111−1),x=(02).A=\begin{pmatrix}-1&1\\1&-1\end{pmatrix},\qquad \boldsymbol x=\begin{pmatrix}0\\2\end{pmatrix}.

对实数 t≥0t\ge0,求向量 y(t)=exp⁡(At)x\boldsymbol y(t)=\exp(At)\boldsymbol x 及 lim⁡t→∞y(t)\lim_{t\to\infty}\boldsymbol y(t)。

问2 考虑矩阵

A=(100011).A=\begin{pmatrix}1&0&0\\0&1&1\end{pmatrix}.

(1) 求矩阵 B=AATB=AA^{\mathsf T} 的全部特征值及对应的归一化特征向量组。其中 T\mathsf T 表示转置。

(2) 求矩阵 C=ATAC=A^{\mathsf T}A 的全部特征值及对应的归一化特征向量组。

(3) 将 BB 的特征向量适当排列,构成正交矩阵 UU;将 CC 的特征向量适当排列,构成正交矩阵 VV,则 AA 可写为 A=UWVTA=UWV^{\mathsf T}。求此时的矩阵 WW。

Kai​

問1​

(1) det⁡A=a2−1\det A=a^2-1 より、正則でない条件は

a=±1.\boxed{a=\pm1}.

どちらの場合も左上の 2×22\times2 小行列式は −1≠0-1\ne0 であり、det⁡A=0\det A=0 なので rank⁡A=2\operatorname{rank}A=2。階数・退化次数の定理より、

dim⁡ker⁡A=3−2=1.\boxed{\dim\ker A=3-2=1}.

(2) AA の固有値は 0,−20,-2 であり、正規化固有ベクトルを列に並べた直交行列

P=12(111−1)P=\frac1{\sqrt2}\begin{pmatrix}1&1\\1&-1\end{pmatrix}

を用いると PTAP=diag⁡(0,−2)P^{\mathsf T}AP=\operatorname{diag}(0,-2) である。したがって、

eAt=P(100e−2t)PT=12(1+e−2t1−e−2t1−e−2t1+e−2t).e^{At}=P\begin{pmatrix}1&0\\0&e^{-2t}\end{pmatrix}P^{\mathsf T} =\frac12\begin{pmatrix}1+e^{-2t}&1-e^{-2t}\\1-e^{-2t}&1+e^{-2t}\end{pmatrix}.

ゆえに、

y(t)=(1−e−2t1+e−2t),lim⁡t→∞y(t)=(11).\boxed{\boldsymbol y(t)=\begin{pmatrix}1-e^{-2t}\\1+e^{-2t}\end{pmatrix}}, \qquad \boxed{\lim_{t\to\infty}\boldsymbol y(t)=\begin{pmatrix}1\\1\end{pmatrix}}.

問2​

(1) B=diag⁡(1,2)B=\operatorname{diag}(1,2) であるから、固有値と正規化固有ベクトルの組は

(1,(10)),(2,(01)).\boxed{\left(1,\begin{pmatrix}1\\0\end{pmatrix}\right),\qquad \left(2,\begin{pmatrix}0\\1\end{pmatrix}\right)}.

(2)

C=(100011011)C=\begin{pmatrix}1&0&0\\0&1&1\\0&1&1\end{pmatrix}

より、求める組は

(1,(100)),(2,12(011)),(0,12(01−1)).\boxed{ \left(1,\begin{pmatrix}1\\0\\0\end{pmatrix}\right),\quad \left(2,\frac1{\sqrt2}\begin{pmatrix}0\\1\\1\end{pmatrix}\right),\quad \left(0,\frac1{\sqrt2}\begin{pmatrix}0\\1\\-1\end{pmatrix}\right)}.

固有ベクトルの符号を反転してもよい。

(3) 上の順序に合わせて

U=I2,V=(10001/21/201/2−1/2)U=I_2,\qquad V=\begin{pmatrix}1&0&0\\0&1/\sqrt2&1/\sqrt2\\0&1/\sqrt2&-1/\sqrt2\end{pmatrix}

とすれば、

W=UTAV=(100020).\boxed{W=U^{\mathsf T}AV=\begin{pmatrix}1&0&0\\0&\sqrt2&0\end{pmatrix}}.

実際に UWVT=AUWV^{\mathsf T}=A が成立する。