跳到主要内容

埼玉大学 理工学研究科 数理電子情報系専攻 情報システム工学コース 2019年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

  1. 以下の問に答えよ. [Solve the following problems.]

(a) 次の行列 A について考える. m=1,2,3,...m = 1, 2, 3, ... のとき, AmA^m を求めよ. ここで, BBn×nn \times n 正方行列, IIn×nn \times n 単位行列とする. [Let us consider the following matrix A. Obtain AmA^m for m=1,2,3,...m = 1, 2, 3, ... , where BB is an n×nn \times n square matrix, II is the n×nn \times n identity matrix.]

A=(IBBBI+B)A = \begin{pmatrix} I - B & B \\ -B & I + B \end{pmatrix}

(b) 次の行列 C について考える. [Let us consider the following matrix C.]

C=(110112022)C = \begin{pmatrix} 1 & 1 & 0 \\ 1 & -1 & 2 \\ 0 & 2 & -2 \end{pmatrix}

(1) 行列 C の固有値と対応する固有ベクトルをすべて求めよ. [Obtain all eigenvalues and corresponding eigenvectors of the matrix C.]

(2) 集合 {CxxR3}\{Cx | x \in \mathbb{R}^3\} は平面であることを示せ. ここで, R3\mathbb{R}^3 は 3 次元ユークリッド空間を表し, xx は 3 次元列ベクトルを表す. [Show the set {CxxR3}\{Cx | x \in \mathbb{R}^3\} is a plane, where R3\mathbb{R}^3 denotes the three-dimensional Euclidean space and xx denotes the three-dimensional column vector.]

(3) 上の (2) で示した平面と点 p=(1,1,1)Tp = (1, 1, 1)^T の距離を求めよ. ここで, 上付添字 TT はベクトルの転置を表す. [Obtain the distance between the plane shown in (2) and the point p=(1,1,1)Tp = (1, 1, 1)^T , where the superscript TT denotes the transposition of a vector.]

题目描述

  1. 回答下列问题。

(a) 考虑分块矩阵

A=(IBBBI+B),A=\begin{pmatrix} I-B&B\\ -B&I+B \end{pmatrix},

其中 BBn×nn\times n 方阵,IIn×nn\times n 单位矩阵。对

m=1,2,3,,m=1,2,3,\dots,

AmA^m

(b) 考虑矩阵

C=(110112022).C=\begin{pmatrix} 1&1&0\\ 1&-1&2\\ 0&2&-2 \end{pmatrix}.

(1) 求矩阵 CC 的全部特征值及对应的特征向量。

(2) 证明集合

{Cx|xR3}\left\{Cx\,\middle|\,x\in\mathbb{R}^3\right\}

是一个平面,其中 R3\mathbb{R}^3 表示三维欧几里得空间,xx 表示三维列向量。

(3) 求 (2) 中平面与点

p=(1,1,1)Tp=(1,1,1)^T

之间的距离,其中上标 TT 表示向量的转置。

Kai

(a) Let A=(IBBBI+B)A = \begin{pmatrix} I - B & B \\ -B & I + B \end{pmatrix} . Then

A2=(IBBBI+B)(IBBBI+B)=((IB)2B2(IB)B+B(I+B)B(IB)(I+B)BB2+(I+B)2)=(I2B+B2B2BB2+B+B2B+B2BB2B2+I+2B+B2)=(I2B2B2BI+2B)A^2 = \begin{pmatrix} I - B & B \\ -B & I + B \end{pmatrix} \begin{pmatrix} I - B & B \\ -B & I + B \end{pmatrix} = \begin{pmatrix} (I - B)^2 - B^2 & (I - B)B + B(I + B) \\ -B(I - B) - (I + B)B & -B^2 + (I + B)^2 \end{pmatrix} = \begin{pmatrix} I - 2B + B^2 - B^2 & B - B^2 + B + B^2 \\ -B + B^2 - B - B^2 & -B^2 + I + 2B + B^2 \end{pmatrix} = \begin{pmatrix} I - 2B & 2B \\ -2B & I + 2B \end{pmatrix}

In general, we can prove by induction that

Am=(ImBmBmBI+mB)A^m = \begin{pmatrix} I - mB & mB \\ -mB & I + mB \end{pmatrix} .

Base case: m=1m = 1 , it is correct.

Inductive step: Assume Ak=(IkBkBkBI+kB)A^k = \begin{pmatrix} I - kB & kB \\ -kB & I + kB \end{pmatrix} . Then

Ak+1=AkA=(IkBkBkBI+kB)(IBBBI+B)=((IkB)(IB)kB2(IkB)B+kB(I+B)kB(IB)(I+kB)BkB2+(I+kB)(I+B))=(IBkB+kB2kB2BkB2+kB+kB2kB+kB2BkB2kB2+I+B+kB+kB2)=(I(k+1)B(k+1)B(k+1)BI+(k+1)B)A^{k+1} = A^k A = \begin{pmatrix} I - kB & kB \\ -kB & I + kB \end{pmatrix} \begin{pmatrix} I - B & B \\ -B & I + B \end{pmatrix} = \begin{pmatrix} (I - kB)(I - B) - kB^2 & (I - kB)B + kB(I + B) \\ -kB(I - B) - (I + kB)B & -kB^2 + (I + kB)(I + B) \end{pmatrix} = \begin{pmatrix} I - B - kB + kB^2 - kB^2 & B - kB^2 + kB + kB^2 \\ -kB + kB^2 - B - kB^2 & -kB^2 + I + B + kB + kB^2 \end{pmatrix} = \begin{pmatrix} I - (k+1)B & (k+1)B \\ -(k+1)B & I + (k+1)B \end{pmatrix}

(b) (1) To find eigenvalues, we compute the determinant of CλIC - \lambda I :

det(CλI)=det(1λ1011λ2022λ)=(1λ)[(1λ)(2λ)4]1[1(2λ)0]+0=(1λ)(λ2+3λ+24)(2λ)=(1λ)(λ2+3λ2)+2+λ=λ2+3λ2λ33λ2+2λ+2+λ=λ32λ2+6λ=λ(λ2+2λ6)det(C - \lambda I) = det \begin{pmatrix} 1-\lambda & 1 & 0 \\ 1 & -1-\lambda & 2 \\ 0 & 2 & -2-\lambda \end{pmatrix} = (1-\lambda)[(-1-\lambda)(-2-\lambda) - 4] - 1[1(-2-\lambda) - 0] + 0 = (1-\lambda)(\lambda^2 + 3\lambda + 2 - 4) - (-2 - \lambda) = (1-\lambda)(\lambda^2 + 3\lambda - 2) + 2 + \lambda = \lambda^2 + 3\lambda - 2 - \lambda^3 - 3\lambda^2 + 2\lambda + 2 + \lambda = -\lambda^3 - 2\lambda^2 + 6\lambda = -\lambda(\lambda^2 + 2\lambda - 6)

Therefore, the eigenvalues are λ1=0\lambda_1 = 0 , λ2,3=2±44(6)2=1±7\lambda_{2,3} = \frac{-2 \pm \sqrt{4 - 4(-6)}}{2} = -1 \pm \sqrt{7} . For λ1=0\lambda_1 = 0 , we have to solve Cx=0Cx = 0 , i.e.

(110112022)(x1x2x3)=(000)\begin{pmatrix} 1 & 1 & 0 \\ 1 & -1 & 2 \\ 0 & 2 & -2 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

x1+x2=0x_1 + x_2 = 0 , x1x2+2x3=0x_1 - x_2 + 2x_3 = 0 , 2x22x3=02x_2 - 2x_3 = 0 . Thus, x2=x3x_2 = x_3 and x1=x2x_1 = -x_2 . The eigenvector corresponding to λ1=0\lambda_1 = 0 is v1=(111)v_1 = \begin{pmatrix} -1 \\ 1 \\ 1 \end{pmatrix} .

For λ2=1+7\lambda_2=-1+\sqrt7 , a corresponding eigenvector is

v2=(3+71+72).v_2=\begin{pmatrix}3+\sqrt7\\1+\sqrt7\\2\end{pmatrix}.

For λ3=17\lambda_3=-1-\sqrt7 , a corresponding eigenvector is

v3=(37172).v_3=\begin{pmatrix}3-\sqrt7\\1-\sqrt7\\2\end{pmatrix}.

Therefore, the corresponding eigenspaces are span{v1}\operatorname{span}\{v_1\} , span{v2}\operatorname{span}\{v_2\} , and span{v3}\operatorname{span}\{v_3\} , respectively.

(2) Since the eigenvalue 0 exists, CC is not invertible. Then, the column space of CC is of dimension at most 2. We observe that the first two columns of CC are linearly independent. Hence, the column space is of dimension 2. Thus {CxxR3}\{Cx | x \in \mathbb{R}^3\} is a plane.

(3) Let v1=(111)v_1 = \begin{pmatrix} -1 \\ 1 \\ 1 \end{pmatrix} be the normal vector of the plane. The plane is given by v1x=0v_1 \cdot x = 0 . i.e. x1+x2+x3=0-x_1 + x_2 + x_3 = 0 . The distance from the point p=(1,1,1)p = (1, 1, 1) to the plane is given by

1+1+1(1)2+12+12=13=33\frac{|-1 + 1 + 1|}{\sqrt{(-1)^2 + 1^2 + 1^2}} = \frac{1}{\sqrt{3}} = \frac{\sqrt{3}}{3}

.