東京大学 新領域創成科学研究科 メディカル情報生命専攻 2025年8月実施 問題7
Author
KardeniaPoyu
Description
行列 A = ( 1 1 1 0 ) A = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} A = ( 1 1 1 0 ) について、以下の問いに答えよ。
(1) A A A の固有値、固有ベクトルを求めよ。
(2) lim n → ∞ A n ( x y ) = ( 0 0 ) \lim_{n \to \infty} A^n \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} lim n → ∞ A n ( x y ) = ( 0 0 ) を満たす、非ゼロの ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) を求めよ。
m × n m \times n m × n 実行列 A A A 、l × n l \times n l × n 実行列 B B B に対して、A x = 0 A \boldsymbol{x} = \boldsymbol{0} A x = 0 を満たすベクトル x \boldsymbol{x} x は必ず B x = 0 B \boldsymbol{x} = \boldsymbol{0} B x = 0 を満たすとき、以下の二つを示せ。
(3) rank A ≥ rank B \text{rank } A \ge \text{rank } B rank A ≥ rank B
(4) B = C A B = CA B = C A を満たす l × m l \times m l × m 実行列 C C C が存在する。
Kai
(1)
行列 A A A の固有値を λ \lambda λ とすると、特性方程式は次のように表される。
det ( A − λ I ) = ∣ 1 − λ 1 1 − λ ∣ = − λ ( 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} det ( A − λ I ) = 1 − λ 1 1 − λ = − λ ( 1 − λ ) − 1 = λ 2 − λ − 1 = 0
これを解くと、固有値は λ = 1 ± 5 2 \lambda = \frac{1 \pm \sqrt{5}}{2} λ = 2 1 ± 5 となる。
次に、それぞれの固有値に対する固有ベクトルを求める。
λ 1 = 1 + 5 2 , λ 2 = 1 − 5 2 \lambda_1 = \frac{1 + \sqrt{5}}{2}, \lambda_2 = \frac{1 - \sqrt{5}}{2} λ 1 = 2 1 + 5 , λ 2 = 2 1 − 5 とおく。
λ 1 = 1 + 5 2 \lambda_1 = \frac{1 + \sqrt{5}}{2} λ 1 = 2 1 + 5 のとき
( A − λ 1 I ) v = 0 (A - \lambda_1 I)\boldsymbol{v} = \boldsymbol{0} ( A − λ 1 I ) v = 0 より、
( 1 − 1 + 5 2 1 1 − 1 + 5 2 ) ( v 1 v 2 ) = ( 0 0 ) \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 − 2 1 + 5 1 1 − 2 1 + 5 ) ( v 1 v 2 ) = ( 0 0 )
( 1 − 5 2 1 1 − 1 + 5 2 ) ( v 1 v 2 ) = ( 0 0 ) \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} ( 2 1 − 5 1 1 − 2 1 + 5 ) ( v 1 v 2 ) = ( 0 0 )
よって、v 1 − 1 + 5 2 v 2 = 0 v_1 - \frac{1 + \sqrt{5}}{2} v_2 = 0 v 1 − 2 1 + 5 v 2 = 0 が得られる。
したがって、対応する固有ベクトルは c 1 ( 1 + 5 2 1 ) ( c 1 ≠ 0 ) c_1 \begin{pmatrix} \frac{1 + \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_1 \neq 0) c 1 ( 2 1 + 5 1 ) ( c 1 = 0 ) となる。
λ 2 = 1 − 5 2 \lambda_2 = \frac{1 - \sqrt{5}}{2} λ 2 = 2 1 − 5 のとき
( A − λ 2 I ) v = 0 (A - \lambda_2 I)\boldsymbol{v} = \boldsymbol{0} ( A − λ 2 I ) v = 0 より、
( 1 − 1 − 5 2 1 1 − 1 − 5 2 ) ( v 1 v 2 ) = ( 0 0 ) \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 − 2 1 − 5 1 1 − 2 1 − 5 ) ( v 1 v 2 ) = ( 0 0 )
( 1 + 5 2 1 1 − 1 − 5 2 ) ( v 1 v 2 ) = ( 0 0 ) \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} ( 2 1 + 5 1 1 − 2 1 − 5 ) ( v 1 v 2 ) = ( 0 0 )
よって、v 1 − 1 − 5 2 v 2 = 0 v_1 - \frac{1 - \sqrt{5}}{2} v_2 = 0 v 1 − 2 1 − 5 v 2 = 0 が得られる。
したがって、対応する固有ベクトルは c 2 ( 1 − 5 2 1 ) ( c 2 ≠ 0 ) c_2 \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_2 \neq 0) c 2 ( 2 1 − 5 1 ) ( c 2 = 0 ) となる。
以上より、固有値と固有ベクトルは以下の通りである。
固有値 1 + 5 2 \frac{1 + \sqrt{5}}{2} 2 1 + 5 に対する固有ベクトル: c 1 ( 1 + 5 2 1 ) ( c 1 ≠ 0 ) c_1 \begin{pmatrix} \frac{1 + \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_1 \neq 0) c 1 ( 2 1 + 5 1 ) ( c 1 = 0 )
固有値 1 − 5 2 \frac{1 - \sqrt{5}}{2} 2 1 − 5 に対する固有ベクトル: c 2 ( 1 − 5 2 1 ) ( c 2 ≠ 0 ) c_2 \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c_2 \neq 0) c 2 ( 2 1 − 5 1 ) ( c 2 = 0 )
(2)
(1) で求めた固有ベクトル v 1 = ( 1 + 5 2 1 ) , v 2 = ( 1 − 5 2 1 ) \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} v 1 = ( 2 1 + 5 1 ) , v 2 = ( 2 1 − 5 1 ) は一次独立であるため、これらを基底として ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) を表すことができる。
実数 α , β \alpha, \beta α , β を用いて、
( x y ) = α v 1 + β v 2 \begin{pmatrix} x \\ y \end{pmatrix} = \alpha \boldsymbol{v}_1 + \beta \boldsymbol{v}_2 ( x y ) = α v 1 + β v 2
と表すとする。このとき、A n ( x y ) A^n \begin{pmatrix} x \\ y \end{pmatrix} A n ( x y ) は次のように計算できる。
A n ( x y ) = A n ( α v 1 + β v 2 ) = α A n v 1 + β A n v 2 = α λ 1 n v 1 + β λ 2 n v 2 \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} A n ( x y ) = A n ( α v 1 + β v 2 ) = α A n v 1 + β A n v 2 = α λ 1 n v 1 + β λ 2 n v 2
ここで、n → ∞ n \to \infty n → ∞ の極限を考える。
λ 1 = 1 + 5 2 ≈ 1.618 > 1 \lambda_1 = \frac{1 + \sqrt{5}}{2} \approx 1.618 > 1 λ 1 = 2 1 + 5 ≈ 1.618 > 1 であるため、λ 1 n \lambda_1^n λ 1 n は発散する。
一方、λ 2 = 1 − 5 2 ≈ − 0.618 \lambda_2 = \frac{1 - \sqrt{5}}{2} \approx -0.618 λ 2 = 2 1 − 5 ≈ − 0.618 であり、∣ λ 2 ∣ < 1 |\lambda_2| < 1 ∣ λ 2 ∣ < 1 であるため、λ 2 n → 0 ( n → ∞ ) \lambda_2^n \to 0 \quad (n \to \infty) λ 2 n → 0 ( n → ∞ ) となる。
したがって、lim n → ∞ A n ( x y ) = ( 0 0 ) \lim_{n \to \infty} A^n \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} lim n → ∞ A n ( x y ) = ( 0 0 ) が成り立つためには、発散する項の係数 α \alpha α が 0 0 0 でなければならない。
α = 0 \alpha = 0 α = 0 のとき、
( x y ) = β v 2 = β ( 1 − 5 2 1 ) \begin{pmatrix} x \\ y \end{pmatrix} = \beta \boldsymbol{v}_2 = \beta \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} ( x y ) = β v 2 = β ( 2 1 − 5 1 )
となる。( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) は非ゼロであるため、β ≠ 0 \beta \neq 0 β = 0 である。
よって、求める非ゼロのベクトル ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) は、
c ( 1 − 5 2 1 ) ( c ≠ 0 ) c \begin{pmatrix} \frac{1 - \sqrt{5}}{2} \\ 1 \end{pmatrix} \quad (c \neq 0) c ( 2 1 − 5 1 ) ( c = 0 )
である。
(3)
A A A を m × n m \times n m × n 行列、B B B を l × n l \times n l × n 行列とする。
仮定より、A x = 0 A\boldsymbol{x} = \boldsymbol{0} A x = 0 を満たすベクトル x \boldsymbol{x} x は必ず B x = 0 B\boldsymbol{x} = \boldsymbol{0} B x = 0 を満たす。
すなわち、行列 A A A の核 (Null space) を ker ( A ) \ker(A) ker ( A ) 、行列 B B B の核を ker ( B ) \ker(B) ker ( B ) とすると、
ker ( A ) ⊆ ker ( B ) \ker(A) \subseteq \ker(B) ker ( A ) ⊆ ker ( B )
が成り立つ。
これより、それぞれの核の次元 (Nullity) について以下の不等式が成り立つ。
dim ( ker ( A ) ) ≤ dim ( ker ( B ) ) \dim(\ker(A)) \le \dim(\ker(B)) dim ( ker ( A )) ≤ dim ( ker ( B ))
ここで、次元定理 (Rank-Nullity Theorem) より、任意の k × n k \times n k × n 行列 M M M について、
rank ( M ) + dim ( ker ( M ) ) = n \text{rank}(M) + \dim(\ker(M)) = n rank ( M ) + dim ( ker ( M )) = n
が成り立つ。行列 A A A と B B B はともに n n n 列の行列であるため、
rank ( A ) + dim ( ker ( A ) ) = n \text{rank}(A) + \dim(\ker(A)) = n rank ( A ) + dim ( ker ( A )) = n
rank ( B ) + dim ( ker ( B ) ) = n \text{rank}(B) + \dim(\ker(B)) = n rank ( B ) + dim ( ker ( B )) = n
となる。これらを変形して不等式に代入すると、
n − rank ( A ) ≤ n − rank ( B ) n - \text{rank}(A) \le n - \text{rank}(B) n − rank ( A ) ≤ n − rank ( B )
rank ( A ) ≥ rank ( B ) \text{rank}(A) \ge \text{rank}(B) rank ( A ) ≥ rank ( B )
が示された。
English:
Let A A A be an m × n m \times n m × n matrix and B B B be an l × n l \times n l × n matrix.
By assumption, any vector x \boldsymbol{x} x satisfying A x = 0 A\boldsymbol{x} = \boldsymbol{0} A x = 0 also satisfies B x = 0 B\boldsymbol{x} = \boldsymbol{0} B x = 0 .
That is, if we denote the null space (kernel) of A A A and B B B as ker ( A ) \ker(A) ker ( A ) and ker ( B ) \ker(B) ker ( B ) respectively, then:
ker ( A ) ⊆ ker ( B ) \ker(A) \subseteq \ker(B) ker ( A ) ⊆ 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)) dim ( ker ( A )) ≤ dim ( ker ( B ))
From the Rank-Nullity Theorem , for any k × n k \times n k × n matrix M M M :
rank ( M ) + dim ( ker ( M ) ) = n \text{rank}(M) + \dim(\ker(M)) = n rank ( M ) + dim ( ker ( M )) = n
Since matrices A A A and B B B both have n n n columns, we have:
rank ( A ) + dim ( ker ( A ) ) = n rank ( B ) + dim ( ker ( B ) ) = n \begin{aligned}
\text{rank}(A) + \dim(\ker(A)) &= n \\
\text{rank}(B) + \dim(\ker(B)) &= n
\end{aligned} rank ( A ) + dim ( ker ( A )) rank ( B ) + dim ( ker ( B )) = n = n
Substituting these into the inequality:
n − rank ( A ) ≤ n − rank ( B ) ⟹ rank ( A ) ≥ rank ( B ) n - \text{rank}(A) \le n - \text{rank}(B) \implies \text{rank}(A) \ge \text{rank}(B) n − rank ( A ) ≤ n − rank ( B ) ⟹ rank ( A ) ≥ rank ( B )
This completes the proof.
(4)
仮定より ker ( A ) ⊆ ker ( B ) \ker(A) \subseteq \ker(B) ker ( A ) ⊆ ker ( B ) である。
任意の部分空間 V V V について、その直交補空間を V ⊥ V^\perp V ⊥ と表す。包含関係の直交補空間をとると、包含関係が逆転するため、
ker ( B ) ⊥ ⊆ ker ( A ) ⊥ \ker(B)^\perp \subseteq \ker(A)^\perp ker ( B ) ⊥ ⊆ ker ( A ) ⊥
が成り立つ。
任意の行列について、その核の直交補空間は行空間 (Row space) と一致する。すなわち、行列 M M M の行空間を R ( M ) R(M) R ( M ) とすると、ker ( M ) ⊥ = R ( M ) \ker(M)^\perp = R(M) ker ( M ) ⊥ = R ( M ) である。
これを用いると、上の包含関係は次のように表せる。
R ( B ) ⊆ R ( A ) R(B) \subseteq R(A) R ( B ) ⊆ R ( A )
これは、行列 B B B のすべての行ベクトルが、行列 A A A の行ベクトルの線形結合として表せることを意味する。
行列 A A A の j j j 番目の行ベクトルを a j T ( j = 1 , … , m ) \boldsymbol{a}_j^T \quad (j=1, \dots, m) a j T ( j = 1 , … , m ) とし、行列 B B B の i i i 番目の行ベクトルを b i T ( i = 1 , … , l ) \boldsymbol{b}_i^T \quad (i=1, \dots, l) b i T ( i = 1 , … , l ) とすると、
任意の i i i について、あるスカラー c i j c_{ij} c ij が存在して、
b i T = ∑ j = 1 m c i j a j T \boldsymbol{b}_i^T = \sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T b i T = j = 1 ∑ m c ij a j T
と表せる。
ここで、c i j c_{ij} c ij を ( i , j ) (i, j) ( i , j ) 成分に持つ l × m l \times m l × m 行列を C C C と定義する。
このとき、行列の積 C A CA C A の第 i i i 行は ∑ j = 1 m c i j a j T \sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T ∑ j = 1 m c ij a j T となり、これは b i T \boldsymbol{b}_i^T b i T に等しい。
よって、B = C A B = CA B = C A が成り立つ。
以上より、B = C A B = CA B = C A を満たす l × m l \times m l × m 実行列 C C C が存在することが示された。
English:
By assumption, ker ( A ) ⊆ ker ( B ) \ker(A) \subseteq \ker(B) ker ( A ) ⊆ ker ( B ) .
For any subspace V V V , let V ⊥ V^\perp V ⊥ 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 ker ( B ) ⊥ ⊆ ker ( A ) ⊥
For any matrix, the orthogonal complement of its null space is its row space. That is, if R ( M ) R(M) R ( M ) is the row space of matrix M M M , then ker ( M ) ⊥ = R ( M ) \ker(M)^\perp = R(M) ker ( M ) ⊥ = R ( M ) .
Thus, the inclusion can be rewritten as:
R ( B ) ⊆ R ( A ) R(B) \subseteq R(A) R ( B ) ⊆ R ( A )
This means every row vector of matrix B B B can be expressed as a linear combination of the row vectors of matrix A A A .
Let a j T ( j = 1 , … , m ) \boldsymbol{a}_j^T \quad (j=1, \dots, m) a j T ( j = 1 , … , m ) be the j j j -th row vector of A A A , and b i T ( i = 1 , … , l ) \boldsymbol{b}_i^T \quad (i=1, \dots, l) b i T ( i = 1 , … , l ) be the i i i -th row vector of B B B .
For each i i i , there exist scalars c i j c_{ij} c ij such that:
b i T = ∑ j = 1 m c i j a j T \boldsymbol{b}_i^T = \sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T b i T = j = 1 ∑ m c ij a j T
Define an l × m l \times m l × m matrix C C C where the ( i , j ) (i, j) ( i , j ) -th entry is c i j c_{ij} c ij .
Then, the i i i -th row of the product C A CA C A is ∑ j = 1 m c i j a j T \sum_{j=1}^{m} c_{ij} \boldsymbol{a}_j^T ∑ j = 1 m c ij a j T , which equals b i T \boldsymbol{b}_i^T b i T .
Therefore, B = C A B = CA B = C A holds.
This proves that there exists an l × m l \times m l × m real matrix C C C satisfying B = C A B = CA B = C A .