跳到主要内容

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

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

問1​

2 次正方行列

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

について、以下の問に答えよ。

(1) AA の固有値とこれに対応する固有ベクトルをすべて求めよ。固有ベクトルは大きさが 11 になるように正規化せよ。

(2) AA について、AnA^n(nn は自然数)を求めよ。

問2​

図1のように、ある点が角速度 ω\omega(ω>0\omega>0)で3次元空間内を等速円運動している。回転軸は原点 OO を通るとし、その方向を単位ベクトル p\boldsymbol p で表す。この回転する点の時刻 tt における位置を x(t)\boldsymbol x(t) とおく。x(t)\boldsymbol x(t) における点の速度を、3次正方行列 CC を用いて

dx(t)dt=Cx(t)\frac{d\boldsymbol x(t)}{dt}=C\boldsymbol x(t)

と表す。CC を p=(p1,p2,p3)T\boldsymbol p=(p_1,p_2,p_3)^T の要素 pip_i(i=1,2,3i=1,2,3)と角速度 ω\omega を用いて表せ。ただし、行列 AA の転置を ATA^T と表す。

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

問3​

XX は実交代行列であるとする。実交代行列とは実数からなる正方行列で、XT=−XX^T=-X を満たすものをいう。Y=exp⁡(X)Y=\exp(X) について、以下の問に答えよ。答えだけでなく導出も示すこと。ただし、exp⁡(X)\exp(X) は次式で表される。

exp⁡(X)=∑n=0∞1n!Xn\exp(X)=\sum_{n=0}^{\infty}\frac1{n!}X^n

(1) YY の行列式 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},

回答以下问题。

(1) 求 AA 的全部特征值及对应的特征向量。特征向量应归一化,使其长度为 11。

(2) 求 AnA^n,其中 nn 为自然数。

问2 如图1,一个点在三维空间内以角速度 ω\omega(ω>0\omega>0)做匀速圆周运动。旋转轴通过原点 OO,其方向用单位向量 p\boldsymbol p 表示。设该旋转点在时刻 tt 的位置为 x(t)\boldsymbol x(t)。使用三阶方阵 CC,将位于 x(t)\boldsymbol x(t) 的点的速度写为

dx(t)dt=Cx(t).\frac{d\boldsymbol x(t)}{dt}=C\boldsymbol x(t).

用 p=(p1,p2,p3)T\boldsymbol p=(p_1,p_2,p_3)^T 的元素 pip_i(i=1,2,3i=1,2,3)和角速度 ω\omega 表示 CC。矩阵 AA 的转置记为 ATA^T。

问3 设 XX 为实反对称矩阵,即元素为实数且满足 XT=−XX^T=-X 的方阵。对于 Y=exp⁡(X)Y=\exp(X),回答以下问题。不仅要给出答案,还要写出推导。其中

exp⁡(X)=∑n=0∞1n!Xn.\exp(X)=\sum_{n=0}^{\infty}\frac1{n!}X^n.

(1) 求 YY 的行列式 det⁡(Y)\det(Y)。可以使用以下事实:实反对称矩阵可对角化,全部非零特征值均为纯虚数,且以 ±λ1,±λ2,…\pm\lambda_1,\pm\lambda_2,\ldots 的形式成对出现。

(2) 求 YYTYY^T 及 YTYY^TY。

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}.

Reference​