跳到主要内容

東京大学 新領域創成科学研究科 メディカル情報生命専攻 2025年8月実施 問題7

Author

KardeniaPoyu

Description

行列 A=(1110)A = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} について、以下の問いに答えよ。

(1) AA の固有値、固有ベクトルを求めよ。
(2) limnAn(xy)=(00)\lim_{n \to \infty} A^n \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} を満たす、非ゼロの (xy)\begin{pmatrix} x \\ y \end{pmatrix} を求めよ。

m×nm \times n 実行列 AAl×nl \times n 実行列 BB に対して、Ax=0A \boldsymbol{x} = \boldsymbol{0} を満たすベクトル x\boldsymbol{x} は必ず Bx=0B \boldsymbol{x} = \boldsymbol{0} を満たすとき、以下の二つを示せ。

(3) rank Arank B\text{rank } A \ge \text{rank } B
(4) B=CAB = CA を満たす l×ml \times m 実行列 CC が存在する。

Kai

(1)

行列 AA の固有値を λ\lambda とすると、特性方程式は次のように表される。

det(AλI)=1λ11λ=λ(1λ)1=λ2λ1=0\begin{aligned} \det(A - \lambda I) &= \begin{vmatrix} 1 - \lambda & 1 \\ 1 & -\lambda \end{vmatrix} \\ &= -\lambda(1 - \lambda) - 1 \\ &= \lambda^2 - \lambda - 1 = 0 \end{aligned}

これを解くと、固有値は λ=1±52\lambda = \frac{1 \pm \sqrt{5}}{2} となる。

次に、それぞれの固有値に対する固有ベクトルを求める。 λ1=1+52,λ2=152\lambda_1 = \frac{1 + \sqrt{5}}{2}, \lambda_2 = \frac{1 - \sqrt{5}}{2} とおく。

λ1=1+52\lambda_1 = \frac{1 + \sqrt{5}}{2} のとき (Aλ1I)v=0(A - \lambda_1 I)\boldsymbol{v} = \boldsymbol{0} より、

(11+52111+52)(v1v2)=(00)\begin{pmatrix} 1 - \frac{1 + \sqrt{5}}{2} & 1 \\ 1 & -\frac{1 + \sqrt{5}}{2} \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}
(152111+52)(v1v2)=(00)\begin{pmatrix} \frac{1 - \sqrt{5}}{2} & 1 \\ 1 & -\frac{1 + \sqrt{5}}{2} \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}

よって、v11+52v2=0v_1 - \frac{1 + \sqrt{5}}{2} v_2 = 0 が得られる。 したがって、対応する固有ベクトルは c1(1+521)(c10)c_1 \begin{pmatrix} \frac{1 + \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_1 \neq 0) となる。

λ2=152\lambda_2 = \frac{1 - \sqrt{5}}{2} のとき (Aλ2I)v=0(A - \lambda_2 I)\boldsymbol{v} = \boldsymbol{0} より、

(115211152)(v1v2)=(00)\begin{pmatrix} 1 - \frac{1 - \sqrt{5}}{2} & 1 \\ 1 & -\frac{1 - \sqrt{5}}{2} \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}
(1+5211152)(v1v2)=(00)\begin{pmatrix} \frac{1 + \sqrt{5}}{2} & 1 \\ 1 & -\frac{1 - \sqrt{5}}{2} \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}

よって、v1152v2=0v_1 - \frac{1 - \sqrt{5}}{2} v_2 = 0 が得られる。 したがって、対応する固有ベクトルは c2(1521)(c20)c_2 \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_2 \neq 0) となる。

以上より、固有値と固有ベクトルは以下の通りである。

  • 固有値 1+52\frac{1 + \sqrt{5}}{2} に対する固有ベクトル: c1(1+521)(c10)c_1 \begin{pmatrix} \frac{1 + \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_1 \neq 0)
  • 固有値 152\frac{1 - \sqrt{5}}{2} に対する固有ベクトル: c2(1521)(c20)c_2 \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_2 \neq 0)

(2)

(1) で求めた固有ベクトル v1=(1+521),v2=(1521)\boldsymbol{v}_1 = \begin{pmatrix} \frac{1 + \sqrt{5}}{2} \\ 1 \end{pmatrix}, \boldsymbol{v}_2 = \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} は一次独立であるため、これらを基底として (xy)\begin{pmatrix} x \\ y \end{pmatrix} を表すことができる。 実数 α,β\alpha, \beta を用いて、

(xy)=αv1+βv2\begin{pmatrix} x \\ y \end{pmatrix} = \alpha \boldsymbol{v}_1 + \beta \boldsymbol{v}_2

と表すとする。このとき、An(xy)A^n \begin{pmatrix} x \\ y \end{pmatrix} は次のように計算できる。

An(xy)=An(αv1+βv2)=αAnv1+βAnv2=αλ1nv1+βλ2nv2\begin{aligned} A^n \begin{pmatrix} x \\ y \end{pmatrix} &= A^n (\alpha \boldsymbol{v}_1 + \beta \boldsymbol{v}_2) \\ &= \alpha A^n \boldsymbol{v}_1 + \beta A^n \boldsymbol{v}_2 \\ &= \alpha \lambda_1^n \boldsymbol{v}_1 + \beta \lambda_2^n \boldsymbol{v}_2 \end{aligned}

ここで、nn \to \infty の極限を考える。 λ1=1+521.618>1\lambda_1 = \frac{1 + \sqrt{5}}{2} \approx 1.618 > 1 であるため、λ1n\lambda_1^n は発散する。 一方、λ2=1520.618\lambda_2 = \frac{1 - \sqrt{5}}{2} \approx -0.618 であり、λ2<1|\lambda_2| < 1 であるため、λ2n0(n)\lambda_2^n \to 0 \quad (n \to \infty) となる。

したがって、limnAn(xy)=(00)\lim_{n \to \infty} A^n \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} が成り立つためには、発散する項の係数 α\alpha00 でなければならない。

α=0\alpha = 0 のとき、

(xy)=βv2=β(1521)\begin{pmatrix} x \\ y \end{pmatrix} = \beta \boldsymbol{v}_2 = \beta \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix}

となる。(xy)\begin{pmatrix} x \\ y \end{pmatrix} は非ゼロであるため、β0\beta \neq 0 である。 よって、求める非ゼロのベクトル (xy)\begin{pmatrix} x \\ y \end{pmatrix} は、

c(1521)(c0)c \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c \neq 0)

である。

(3)

AAm×nm \times n 行列、BBl×nl \times n 行列とする。 仮定より、Ax=0A\boldsymbol{x} = \boldsymbol{0} を満たすベクトル x\boldsymbol{x} は必ず Bx=0B\boldsymbol{x} = \boldsymbol{0} を満たす。 すなわち、行列 AA の核 (Null space) を ker(A)\ker(A)、行列 BB の核を ker(B)\ker(B) とすると、

ker(A)ker(B)\ker(A) \subseteq \ker(B)

が成り立つ。 これより、それぞれの核の次元 (Nullity) について以下の不等式が成り立つ。

dim(ker(A))dim(ker(B))\dim(\ker(A)) \le \dim(\ker(B))

ここで、次元定理 (Rank-Nullity Theorem) より、任意の k×nk \times n 行列 MM について、

rank(M)+dim(ker(M))=n\text{rank}(M) + \dim(\ker(M)) = n

が成り立つ。行列 AABB はともに nn 列の行列であるため、

rank(A)+dim(ker(A))=n\text{rank}(A) + \dim(\ker(A)) = n
rank(B)+dim(ker(B))=n\text{rank}(B) + \dim(\ker(B)) = n

となる。これらを変形して不等式に代入すると、

nrank(A)nrank(B)n - \text{rank}(A) \le n - \text{rank}(B)
rank(A)rank(B)\text{rank}(A) \ge \text{rank}(B)

が示された。

English:

Let AA be an m×nm \times n matrix and BB be an l×nl \times n matrix. By assumption, any vector x\boldsymbol{x} satisfying Ax=0A\boldsymbol{x} = \boldsymbol{0} also satisfies Bx=0B\boldsymbol{x} = \boldsymbol{0}. That is, if we denote the null space (kernel) of AA and BB as ker(A)\ker(A) and ker(B)\ker(B) respectively, then:

ker(A)ker(B)\ker(A) \subseteq \ker(B)

This implies the following inequality for the dimensions of the kernels (nullities):

dim(ker(A))dim(ker(B))\dim(\ker(A)) \le \dim(\ker(B))

From the Rank-Nullity Theorem, for any k×nk \times n matrix MM:

rank(M)+dim(ker(M))=n\text{rank}(M) + \dim(\ker(M)) = n

Since matrices AA and BB both have nn columns, we have:

rank(A)+dim(ker(A))=nrank(B)+dim(ker(B))=n\begin{aligned} \text{rank}(A) + \dim(\ker(A)) &= n \\ \text{rank}(B) + \dim(\ker(B)) &= n \end{aligned}

Substituting these into the inequality:

nrank(A)nrank(B)    rank(A)rank(B)n - \text{rank}(A) \le n - \text{rank}(B) \implies \text{rank}(A) \ge \text{rank}(B)

This completes the proof.

(4)

仮定より ker(A)ker(B)\ker(A) \subseteq \ker(B) である。 任意の部分空間 VV について、その直交補空間を VV^\perp と表す。包含関係の直交補空間をとると、包含関係が逆転するため、

ker(B)ker(A)\ker(B)^\perp \subseteq \ker(A)^\perp

が成り立つ。 任意の行列について、その核の直交補空間は行空間 (Row space) と一致する。すなわち、行列 MM の行空間を R(M)R(M) とすると、ker(M)=R(M)\ker(M)^\perp = R(M) である。 これを用いると、上の包含関係は次のように表せる。

R(B)R(A)R(B) \subseteq R(A)

これは、行列 BB のすべての行ベクトルが、行列 AA の行ベクトルの線形結合として表せることを意味する。 行列 AAjj 番目の行ベクトルを ajT(j=1,,m)\boldsymbol{a}_j^T \quad (j=1, \dots, m) とし、行列 BBii 番目の行ベクトルを biT(i=1,,l)\boldsymbol{b}_i^T \quad (i=1, \dots, l) とすると、 任意の ii について、あるスカラー cijc_{ij} が存在して、

biT=j=1mcijajT\boldsymbol{b}_i^T = \sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T

と表せる。 ここで、cijc_{ij}(i,j)(i, j) 成分に持つ l×ml \times m 行列を CC と定義する。 このとき、行列の積 CACA の第 ii 行は j=1mcijajT\sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T となり、これは biT\boldsymbol{b}_i^T に等しい。 よって、B=CAB = CA が成り立つ。 以上より、B=CAB = CA を満たす l×ml \times m 実行列 CC が存在することが示された。

English:

By assumption, ker(A)ker(B)\ker(A) \subseteq \ker(B). For any subspace VV, let VV^\perp denote its orthogonal complement. Taking the orthogonal complement of both sides reverses the inclusion relation:

ker(B)ker(A)\ker(B)^\perp \subseteq \ker(A)^\perp

For any matrix, the orthogonal complement of its null space is its row space. That is, if R(M)R(M) is the row space of matrix MM, then ker(M)=R(M)\ker(M)^\perp = R(M). Thus, the inclusion can be rewritten as:

R(B)R(A)R(B) \subseteq R(A)

This means every row vector of matrix BB can be expressed as a linear combination of the row vectors of matrix AA. Let ajT(j=1,,m)\boldsymbol{a}_j^T \quad (j=1, \dots, m) be the jj-th row vector of AA, and biT(i=1,,l)\boldsymbol{b}_i^T \quad (i=1, \dots, l) be the ii-th row vector of BB. For each ii, there exist scalars cijc_{ij} such that:

biT=j=1mcijajT\boldsymbol{b}_i^T = \sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T

Define an l×ml \times m matrix CC where the (i,j)(i, j)-th entry is cijc_{ij}. Then, the ii-th row of the product CACA is j=1mcijajT\sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T, which equals biT\boldsymbol{b}_i^T. Therefore, B=CAB = CA holds.

This proves that there exists an l×ml \times m real matrix CC satisfying B=CAB = CA.