跳到主要内容

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

Author

KardeniaPoyu

Description

以下の問いに答えよ。

(1) 行列 A=(211121112)A = \begin{pmatrix} 2 & -1 & -1 \\ -1 & 2 & -1 \\ -1 & -1 & 2 \end{pmatrix} の固有値を求めよ。

(2) 行列 AA の固有空間の基底を求めよ。

(3) PTAPP^T AP が対角行列になるような直交行列 PP を求めよ。

(4) x=(x1,x2,x3)T\boldsymbol{x} = (x_1, x_2, x_3)^T を三次元実ベクトルとする。f(x)=xTAxf(\boldsymbol{x}) = \boldsymbol{x}^T A \boldsymbol{x} の最小値を与える最適解的集合は、三次元空間における直線になっている。この直線を表す方程式を示せ。

(5) 実数の成分からなる対称正方行列に対して、すべての固有値が異なるならば、すべての固有ベクトルは相互に直交することを証明せよ。

Kai

(1)

行列 AA の固有方程式 λIA=0|\lambda I - A| = 0 を解く。

λIA=λ2111λ2111λ2|\lambda I - A| = \begin{vmatrix} \lambda - 2 & 1 & 1 \\ 1 & \lambda - 2 & 1 \\ 1 & 1 & \lambda - 2 \end{vmatrix}

各列を第1列に加えると、

λIA=λ11λλ21λ1λ2=λ1111λ2111λ2\begin{aligned} |\lambda I - A| &= \begin{vmatrix} \lambda & 1 & 1 \\ \lambda & \lambda - 2 & 1 \\ \lambda & 1 & \lambda - 2 \end{vmatrix} \\ &= \lambda \begin{vmatrix} 1 & 1 & 1 \\ 1 & \lambda - 2 & 1 \\ 1 & 1 & \lambda - 2 \end{vmatrix} \end{aligned}

第2行および第3行から第1行を引くと、

=λ1110λ3000λ3=λ(λ3)2=0\begin{aligned} &= \lambda \begin{vmatrix} 1 & 1 & 1 \\ 0 & \lambda - 3 & 0 \\ 0 & 0 & \lambda - 3 \end{vmatrix} \\ &= \lambda (\lambda - 3)^2 = 0 \end{aligned}

ゆえに、行列 AA の固有値は λ=0,3\lambda = 0, 333 は2重根) である。


(別解:サラスの公式による展開)

λIA=λ2111λ2111λ2=(λ2)3+(111)+(111)(λ2)1111(λ2)1(λ2)1=(λ36λ2+12λ8)+23(λ2)=λ36λ2+12λ63λ+6=λ36λ2+9λ=λ(λ26λ+9)=λ(λ3)2\begin{aligned} |\lambda I - A| &= \begin{vmatrix} \lambda - 2 & 1 & 1 \\ 1 & \lambda - 2 & 1 \\ 1 & 1 & \lambda - 2 \end{vmatrix} \\ &= (\lambda - 2)^3 + (1 \cdot 1 \cdot 1) + (1 \cdot 1 \cdot 1) - (\lambda - 2) \cdot 1 \cdot 1 - 1 \cdot 1 \cdot (\lambda - 2) - 1 \cdot (\lambda - 2) \cdot 1 \\ &= (\lambda^3 - 6\lambda^2 + 12\lambda - 8) + 2 - 3(\lambda - 2) \\ &= \lambda^3 - 6\lambda^2 + 12\lambda - 6 - 3\lambda + 6 \\ &= \lambda^3 - 6\lambda^2 + 9\lambda \\ &= \lambda(\lambda^2 - 6\lambda + 9) \\ &= \lambda(\lambda - 3)^2 \end{aligned}

(2)

固有値 λ\lambda に属する固有空間 W(λ)W(\lambda) は、方程式 (λIA)x=0(\lambda I - A)\boldsymbol{x} = \boldsymbol{0} の解空間である。

λ=0\lambda = 0 のとき

(λIA)x=(211121112)(x1x2x3)=(000)\begin{aligned} (\lambda I - A)\boldsymbol{x} = \begin{pmatrix} -2 & 1 & 1 \\ 1 & -2 & 1 \\ 1 & 1 & -2 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \end{aligned}

行基本変形により x1x2=0x_1 - x_2 = 0 かつ x2x3=0x_2 - x_3 = 0 を得る。すなわち x1=x2=x3x_1 = x_2 = x_3。 したがって、W(0)W(0) の基底は以下の通りとなる。

{(111)}\begin{aligned} \left\{ \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix} \right\} \end{aligned}

λ=3\lambda = 3 のとき

(3IA)x=(111111111)(x1x2x3)=(000)\begin{aligned} (3I - A)\boldsymbol{x} = \begin{pmatrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \end{aligned}

これは x1+x2+x3=0x_1 + x_2 + x_3 = 0 と同値である。x2=s,x3=tx_2 = s, x_3 = t とおくと x1=stx_1 = -s - t となる。 したがって、W(3)W(3) の基底は以下の通りとなる。

{(110),(101)}\begin{aligned} \left\{ \begin{pmatrix} -1 \\ 1 \\ 0 \end{pmatrix}, \begin{pmatrix} -1 \\ 0 \\ 1 \end{pmatrix} \right\} \end{aligned}

(3)

実対称行列は直交行列によって対角化可能である。各固有空間における正規直交基底を構成する。

  • W(0)W(0) の正規直交基底: v1=(1,1,1)T\boldsymbol{v}_1 = (1, 1, 1)^T を正規化すると、

    u1=13(111)\boldsymbol{u}_1 = \frac{1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}

    となる。

  • W(3)W(3) の正規直交基底: v2=(1,1,0)T,v3=(1,0,1)T\boldsymbol{v}_2 = (-1, 1, 0)^T, \boldsymbol{v}_3 = (-1, 0, 1)^T とし、グラム・シュミットの直交化法を用いる。 まず v2\boldsymbol{v}_2 を正規化すると、

    u2=12(110)\boldsymbol{u}_2 = \frac{1}{\sqrt{2}} \begin{pmatrix} -1 \\ 1 \\ 0 \end{pmatrix}

    となる。次に v3\boldsymbol{v}_3u2\boldsymbol{u}_2 に対する直交成分 v3\boldsymbol{v}_3' を求めると、

    v3=v3(v3u2)u2=(101)12(110)=(1/21/21)\begin{aligned} \boldsymbol{v}_3' &= \boldsymbol{v}_3 - (\boldsymbol{v}_3 \cdot \boldsymbol{u}_2)\boldsymbol{u}_2 \\ &= \begin{pmatrix} -1 \\ 0 \\ 1 \end{pmatrix} - \frac{1}{2}\begin{pmatrix} -1 \\ 1 \\ 0 \end{pmatrix} = \begin{pmatrix} -1/2 \\ -1/2 \\ 1 \end{pmatrix} \end{aligned}

    これを正規化すると、

    u3=16(112)\boldsymbol{u}_3 = \frac{1}{\sqrt{6}} \begin{pmatrix} -1 \\ -1 \\ 2 \end{pmatrix}

    を得る。

以上より、求める直交行列 PP は以下の通りである。

P=(u1u2u3)=(13121613121613026)P = (\boldsymbol{u}_1 \quad \boldsymbol{u}_2 \quad \boldsymbol{u}_3) = \begin{pmatrix} \frac{1}{\sqrt{3}} & -\frac{1}{\sqrt{2}} & -\frac{1}{\sqrt{6}} \\ \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{2}} & -\frac{1}{\sqrt{6}} \\ \frac{1}{\sqrt{3}} & 0 & \frac{2}{\sqrt{6}} \end{pmatrix}

(4)

直交行列 PP を用いて x=Py\boldsymbol{x} = P\boldsymbol{y}y=(y1,y2,y3)T\boldsymbol{y} = (y_1, y_2, y_3)^T)と変数変換を行うと、二次形式 f(x)f(\boldsymbol{x}) は以下のように対角化される。

f(x)=yT(PTAP)y=0y12+3y22+3y32\begin{aligned} f(\boldsymbol{x}) &= \boldsymbol{y}^T (P^T A P) \boldsymbol{y} \\ &= 0y_1^2 + 3y_2^2 + 3y_3^2 \end{aligned}

y1,y2,y3y_1, y_2, y_3 は実数であるため、f(x)=3y22+3y320f(\boldsymbol{x}) = 3y_2^2 + 3y_3^2 \geq 0 が成り立つ。 したがって、f(x)f(\boldsymbol{x}) の最小値は 00 であり、条件は y2=y3=0y_2 = y_3 = 0 である。このとき、

x=P(y100)=y1u1=y13(111)\boldsymbol{x} = P \begin{pmatrix} y_1 \\ 0 \\ 0 \end{pmatrix} = y_1 \boldsymbol{u}_1 = \frac{y_1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}

これは固有空間 W(0)W(0) に属する任意のベクトルを表す。よって、最適解の集合がなす直線の方程式は x1=x2=x3x_1 = x_2 = x_3 である。

English:

By applying the orthogonal transformation x=Py\boldsymbol{x} = P\boldsymbol{y} using the matrix PP from (3), the quadratic form f(x)f(\boldsymbol{x}) is diagonalized as:

f(x)=yT(PTAP)y=0y12+3y22+3y32\begin{aligned} f(\boldsymbol{x}) &= \boldsymbol{y}^T (P^T A P) \boldsymbol{y} \\ &= 0y_1^2 + 3y_2^2 + 3y_3^2 \end{aligned}

Since yiRy_i \in \mathbb{R}, f(x)0f(\boldsymbol{x}) \geq 0 holds. The minimum value 00 is achieved when y2=y3=0y_2 = y_3 = 0, which implies x\boldsymbol{x} is any multiple of the eigenvector u1\boldsymbol{u}_1:

x=P(y100)=y1u1=y13(111)\boldsymbol{x} = P \begin{pmatrix} y_1 \\ 0 \\ 0 \end{pmatrix} = y_1 \boldsymbol{u}_1 = \frac{y_1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}

Thus, the set of optimal solutions is the line x1=x2=x3x_1 = x_2 = x_3.

(5)

実対称行列 AA (AT=AA^T = A) の相異なる固有値を λ1,λ2\lambda_1, \lambda_2、対応する固有ベクトルを x1,x2\boldsymbol{x}_1, \boldsymbol{x}_2 とする。

Ax1=λ1x1(i)Ax2=λ2x2(ii)\begin{aligned} A \boldsymbol{x}_1 &= \lambda_1 \boldsymbol{x}_1 \quad \dots (i) \\ A \boldsymbol{x}_2 &= \lambda_2 \boldsymbol{x}_2 \quad \dots (ii) \end{aligned}

(i)(i) の両辺を転置すると、x1TAT=λ1x1T\boldsymbol{x}_1^T A^T = \lambda_1 \boldsymbol{x}_1^T となる。AT=AA^T = A より、

x1TA=λ1x1T\begin{aligned} \boldsymbol{x}_1^T A = \lambda_1 \boldsymbol{x}_1^T \end{aligned}

この両辺に右から x2\boldsymbol{x}_2 を掛けると、

x1TAx2=λ1x1Tx2(iii)\begin{aligned} \boldsymbol{x}_1^T A \boldsymbol{x}_2 = \lambda_1 \boldsymbol{x}_1^T \boldsymbol{x}_2 \quad \dots (iii) \end{aligned}

一方、式 (iii)(iii) の左辺に式 (ii)(ii) を代入すると、

x1T(Ax2)=x1T(λ2x2)=λ2x1Tx2(iv)\begin{aligned} \boldsymbol{x}_1^T (A \boldsymbol{x}_2) = \boldsymbol{x}_1^T (\lambda_2 \boldsymbol{x}_2) = \lambda_2 \boldsymbol{x}_1^T \boldsymbol{x}_2 \quad \dots (iv) \end{aligned}

(iii)(iii)(iv)(iv) より、

(λ1λ2)x1Tx2=0\begin{aligned} (\lambda_1 - \lambda_2) \boldsymbol{x}_1^T \boldsymbol{x}_2 = 0 \end{aligned}

λ1λ2\lambda_1 \neq \lambda_2 より λ1λ20\lambda_1 - \lambda_2 \neq 0 であるから、

x1Tx2=0\begin{aligned} \boldsymbol{x}_1^T \boldsymbol{x}_2 = 0 \end{aligned}

が成り立つ。これは x1\boldsymbol{x}_1x2\boldsymbol{x}_2 が直交することを意味する。(証明終)

English:

Let λ1,λ2\lambda_1, \lambda_2 be distinct eigenvalues of a real symmetric matrix AA, with eigenvectors x1,x2\boldsymbol{x}_1, \boldsymbol{x}_2.

Ax1=λ1x1(i)Ax2=λ2x2(ii)\begin{aligned} A\boldsymbol{x}_1 &= \lambda_1 \boldsymbol{x}_1 \quad \dots (i) \\ A\boldsymbol{x}_2 &= \lambda_2 \boldsymbol{x}_2 \quad \dots (ii) \end{aligned}

Taking the transpose of (i) and using AT=AA^T = A, we have x1TA=λ1x1T\boldsymbol{x}_1^T A = \lambda_1 \boldsymbol{x}_1^T. Multiplying by x2\boldsymbol{x}_2 from the right yields:

x1TAx2=λ1x1Tx2(iii)\begin{aligned} \boldsymbol{x}_1^T A \boldsymbol{x}_2 = \lambda_1 \boldsymbol{x}_1^T \boldsymbol{x}_2 \quad \dots (iii) \end{aligned}

Simultaneously, from (ii), we have:

x1TAx2=λ2x1Tx2(iv)\begin{aligned} \boldsymbol{x}_1^T A \boldsymbol{x}_2 = \lambda_2 \boldsymbol{x}_1^T \boldsymbol{x}_2 \quad \dots (iv) \end{aligned}

From (iii) and (iv), it follows that:

(λ1λ2)x1Tx2=0\begin{aligned} (\lambda_1 - \lambda_2) \boldsymbol{x}_1^T \boldsymbol{x}_2 = 0 \end{aligned}

Since λ1λ2\lambda_1 \neq \lambda_2, we must have:

x1Tx2=0\begin{aligned} \boldsymbol{x}_1^T \boldsymbol{x}_2 = 0 \end{aligned}

proving orthogonality. (Q.E.D.)