跳到主要内容

東京工業大学 工学院 システム制御系 2023年度 数学 問題2(2022年8月実施)

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

問1​

行列

A=(2112)A=\begin{pmatrix}2&1\\1&2\end{pmatrix}

について、次の問いに答えよ。

  1. 固有値と対応する固有ベクトルをすべて求めよ。固有ベクトルの長さは 11 に正規化すること。
  2. 自然数 nn に対して AnA^n を求めよ。

問2​

図1の点は、原点 OO を通る軸の周りを角速度 ω>0\omega>0 で等速円運動する。軸の向きを単位ベクトル p=(p1,p2,p3)T\boldsymbol p=(p_1,p_2,p_3)^T、時刻 tt の点の位置を x(t)\boldsymbol x(t) と表す。

図1 原点を通る軸の周りの等速円運動

原図:大学公式問題(4ページ)

速度を dx(t)/dt=Cx(t)d\boldsymbol x(t)/dt=C\boldsymbol x(t) と書くとき、3×33\times3 行列 CC を p1,p2,p3,ωp_1,p_2,p_3,\omega で表せ。上付きの TT は転置を意味する。

問3​

実数成分の正方行列 XX が XT=−XX^T=-X を満たすとき、XX を実交代行列という。行列指数関数を

Y=exp⁡X=∑n=0∞Xnn!Y=\exp X=\sum_{n=0}^{\infty}\frac{X^n}{n!}

で定める。以下の結果を導出とともに示せ。

  1. det⁡Y\det Y を求めよ。実交代行列は対角化でき、非零固有値は純虚数で、±λ1,±λ2,…\pm\lambda_1,\pm\lambda_2,\ldots の対をなすことを使ってよい。
  2. YYTYY^T と YTYY^TY を求めよ。

题目描述​

问1:对 A=(2112)A=\begin{pmatrix}2&1\\1&2\end{pmatrix},求全部特征值及对应的单位特征向量,并求自然数 nn 次幂 AnA^n。

问2:如图1,质点以角速度 ω>0\omega>0 绕过原点的轴做匀速圆周运动。转轴方向的单位向量为 p=(p1,p2,p3)T\boldsymbol p=(p_1,p_2,p_3)^T,位置为 x(t)\boldsymbol x(t)。将速度写成 dx/dt=Cxd\boldsymbol x/dt=C\boldsymbol x,用 p1,p2,p3,ωp_1,p_2,p_3,\omega 表示三阶矩阵 CC。TT 表示转置。

问3:设实反对称矩阵 XX 满足 XT=−XX^T=-X,定义 Y=exp⁡X=∑n=0∞Xn/n!Y=\exp X=\sum_{n=0}^{\infty}X^n/n!。推导 det⁡Y\det Y、YYTYY^T 和 YTYY^TY。可使用实反对称矩阵可对角化,且非零特征值为成对的纯虚数 ±λ1,±λ2,…\pm\lambda_1,\pm\lambda_2,\ldots 这一事实。

Kai​

問1​

(1)​

特性方程式は

det⁡(A−λI)=(2−λ)2−1=(λ−1)(λ−3)=0.\det(A-\lambda I)=(2-\lambda)^2-1=(\lambda-1)(\lambda-3)=0.

したがって、固有値と正規化された実固有ベクトルは

λ=1:±12(1−1),λ=3:±12(11).\boxed{\lambda=1:\quad\pm\frac1{\sqrt2}\begin{pmatrix}1\\-1\end{pmatrix}, \qquad \lambda=3:\quad\pm\frac1{\sqrt2}\begin{pmatrix}1\\1\end{pmatrix}}.

(2)​

固有ベクトルを列に並べた直交行列

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

を用いると PTAP=diag⁡(1,3)P^TAP=\operatorname{diag}(1,3) である。よって、

An=P(1003n)PT=12(1+3n3n−13n−11+3n).\boxed{A^n=P\begin{pmatrix}1&0\\0&3^n\end{pmatrix}P^T =\frac12\begin{pmatrix}1+3^n&3^n-1\\3^n-1&1+3^n\end{pmatrix}}.

問2​

図の回転方向に対する角速度ベクトルは ωp\omega\boldsymbol p である。したがって、

dxdt=ωp×x=ω(p2x3−p3x2p3x1−p1x3p1x2−p2x1).\frac{d\boldsymbol x}{dt}=\omega\boldsymbol p\times\boldsymbol x =\omega\begin{pmatrix} p_2x_3-p_3x_2\\p_3x_1-p_1x_3\\p_1x_2-p_2x_1 \end{pmatrix}.

各成分の係数を並べれば、

C=ω(0−p3p2p30−p1−p2p10).\boxed{C=\omega\begin{pmatrix} 0&-p_3&p_2\\p_3&0&-p_1\\-p_2&p_1&0 \end{pmatrix}}.

問3​

(1)​

X=PDP−1X=PDP^{-1} と対角化すると、べき級数の定義から Y=Pexp⁡(D)P−1Y=P\exp(D)P^{-1} となる。零固有値がある場合、その指数は 11 であり、非零固有値は対をなすので、

det⁡Y=det⁡(exp⁡D)=∏jeλje−λj=1.\boxed{\det Y=\det(\exp D) =\prod_j e^{\lambda_j}e^{-\lambda_j}=1}.

(2)​

べき級数の各項を転置して、

YT=∑n=0∞(XT)nn!=exp⁡(−X).Y^T=\sum_{n=0}^{\infty}\frac{(X^T)^n}{n!} =\exp(-X).

XX と −X-X は可換であるから、

YYT=eXe−X=I,YTY=e−XeX=I.\boxed{YY^T=e^Xe^{-X}=I,\qquad Y^TY=e^{-X}e^X=I}.