跳到主要内容

東京工業大学 情報理工学院 数理・計算科学系 2020年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

A\boldsymbol{A}4×44 \times 4 行列, b\boldsymbol{b} を 4 次元列ベクトル, I\boldsymbol{I}4×44 \times 4 単位行列とする.これらを左から順に並べた 4×94 \times 9 行列 (A b I)(\boldsymbol{A} \ \boldsymbol{b} \ \boldsymbol{I}) の左側から 4×44 \times 4 行列 P\boldsymbol{P} を乗じたとき

(111101000011412100001532210000111101)\begin{pmatrix} 1 & -1 & 1 & -1 & 0 & -1 & 0 & 0 & 0 \\ 0 & 1 & 1 & 4 & 1 & 2 & 1 & 0 & 0 \\ 0 & 0 & 1 & 5 & 3 & 2 & 2 & 1 & 0 \\ 0 & 0 & 0 & 1 & 1 & 1 & 1 & 0 & 1 \end{pmatrix}

を得る.以下の問に答えよ.

(1) 4×44 \times 4 行列 P\boldsymbol{P} を記せ.

(2) A\boldsymbol{A} の行列式を計算せよ.

(3) 線形方程式 Ax=b\boldsymbol{Ax} = \boldsymbol{b} の解 x\boldsymbol{x} を計算せよ.

(4) A\boldsymbol{A} の逆行列を計算せよ.

题目描述

A\boldsymbol A4×44\times4 矩阵,b\boldsymbol b 是四维列向量,I\boldsymbol I4×44\times4 单位矩阵。把三者从左到右拼成 4×94\times9 矩阵 (A b I)(\boldsymbol A\ \boldsymbol b\ \boldsymbol I),再从左侧乘以某个 4×44\times4 矩阵 P\boldsymbol P,所得矩阵为

(111101000011412100001532210000111101).\begin{pmatrix} 1&-1&1&-1&0&-1&0&0&0\\ 0&1&1&4&1&2&1&0&0\\ 0&0&1&5&3&2&2&1&0\\ 0&0&0&1&1&1&1&0&1 \end{pmatrix}.
  1. 写出矩阵 P\boldsymbol P
  2. 计算 detA\det\boldsymbol A
  3. 求线性方程组 Ax=b\boldsymbol A\boldsymbol x=\boldsymbol b 的解 x\boldsymbol x
  4. 计算 A1\boldsymbol A^{-1}

Kai

解答

与えられた 4×94\times9 行列を、最初の4列、次の1列、最後の4列に分けて

(U  c  P)(U\ \ c\ \ P)

と書く。仮定 P(A  b  I)=(PA  Pb  P)P(A\ \ b\ \ I)=(PA\ \ Pb\ \ P) と比較すると

PA=U,Pb=cPA=U,\qquad Pb=c

である。

(1)

最後の4列を読み取れば

P=[1000210022101101].\boxed{P= \begin{bmatrix} -1&0&0&0\\ 2&1&0&0\\ 2&2&1&0\\ 1&1&0&1 \end{bmatrix}}.

(2)

最初の4列は

U=[1111011400150001]U=\begin{bmatrix} 1&-1&1&-1\\ 0&1&1&4\\ 0&0&1&5\\ 0&0&0&1 \end{bmatrix}

である。 PA=UPA=U から

det(P)det(A)=det(U).\det(P)\det(A)=\det(U).

両行列は三角行列で、 det(P)=1\det(P)=-1det(U)=1\det(U)=1 だから

det(A)=1.\boxed{\det(A)=-1}.

(3)

Ax=bAx=b に左から PP を掛けると

Ux=c,c=[0131]Ux=c,\qquad c=\begin{bmatrix}0\\1\\3\\1\end{bmatrix}

となる。後退代入により順に

x4=1,x3=2,x2=1,x1=2x_4=1,\qquad x_3=-2,\qquad x_2=-1,\qquad x_1=2

を得る。したがって

x=[2121].\boxed{x=\begin{bmatrix}2\\-1\\-2\\1\end{bmatrix}}.

(4)

PA=UPA=U より A1=U1PA^{-1}=U^{-1}P である。ここで

U1=[1127011100150001].U^{-1}= \begin{bmatrix} 1&1&-2&7\\ 0&1&-1&1\\ 0&0&1&-5\\ 0&0&0&1 \end{bmatrix}.

したがって積を計算すると

A1=[4427101133151101].\boxed{A^{-1}= \begin{bmatrix} 4&4&-2&7\\ 1&0&-1&1\\ -3&-3&1&-5\\ 1&1&0&1 \end{bmatrix}}.

実際、 A=P1UA=P^{-1}U なので (P1U)(U1P)=I(P^{-1}U)(U^{-1}P)=I となり、この結果を検算できる。