跳到主要内容

東京工業大学 工学院 情報通信系 2019年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

H2.

1) 2 つの正方行列 A,BA, B に対して,

B=P1AP(H2.1)B = P^{-1}AP \tag{H2.1}

が成り立つ正則な行列 PP が存在するとき,行列 A,BA, B は相似であるとよぶ.このとき,以下の問に答えよ.

a) 行列 A,BA, B が相似であるとき,これらが同じ固有値を持つことを示せ. b) 式 (H2.1) において A=(4321)A = \begin{pmatrix} 4 & -3 \\ 2 & -1 \end{pmatrix} , P=(2513)P = \begin{pmatrix} 2 & 5 \\ 1 & 3 \end{pmatrix} であるとき,行列 A,BA, B のそれぞれについて,すべての固有値とこれらに対応する長さ 1 の固有ベクトルを求めよ.

2) 実変数 x,y,zx, y, z の二次形式 f(x,y,z)=2x2+2y2+4z2+2xy+2yz+2zxf(x, y, z) = 2x^2 + 2y^2 + 4z^2 + 2xy + 2yz + 2zx が与えられるとき,以下の問に答えよ.

a) f(x,y,z)f(x, y, z) を,実対称行列 TT を用いて (xyz)T(xyz)\begin{pmatrix} x & y & z \end{pmatrix} T \begin{pmatrix} x \\ y \\ z \end{pmatrix} と表すとき, TT を求めよ. b) f(x,y,z)f(x, y, z) は,直交行列 QQ を用いた変換 (xyz)=Q(stu)\begin{pmatrix} x \\ y \\ z \end{pmatrix} = Q \begin{pmatrix} s \\ t \\ u \end{pmatrix} を用いて, 標準形 a1s2+a2t2+a3u2a_1s^2 + a_2t^2 + a_3u^2 に変換できる.ただし, a1,a2,a3a_1, a_2, a_3 は実数で, a1a2a3a_1 \le a_2 \le a_3 とする.このとき, QQ 及び, a1,a2,a3a_1, a_2, a_3 を求めよ. c) 二次曲面 f(x,y,z)23x23y+23z8=0f(x, y, z) - 2\sqrt{3}x - 2\sqrt{3}y + 2\sqrt{3}z - 8 = 0 が楕円面であることを示せ.

题目描述

  1. 对两个方阵 A,BA,B,若存在可逆矩阵 PP 使
B=P1AP,(H2.1)B=P^{-1}AP, \tag{H2.1}

则称 A,BA,B 相似。

  1. 证明相似矩阵 A,BA,B 具有相同的特征值。
  2. 在 (H2.1) 中取
A=(4321),P=(2513).A=\begin{pmatrix}4&-3\\2&-1\end{pmatrix}, \qquad P=\begin{pmatrix}2&5\\1&3\end{pmatrix}.

分别求 AABB 的全部特征值,并对每个特征值求一个对应的长度为 11 的特征向量。

  1. 给定实变量二次形式
f(x,y,z)=2x2+2y2+4z2+2xy+2yz+2zx.f(x,y,z) =2x^2+2y^2+4z^2+2xy+2yz+2zx.
  1. 求实对称矩阵 TT,使
f(x,y,z)=(xyz)T(xyz).f(x,y,z) = \begin{pmatrix}x&y&z\end{pmatrix} T \begin{pmatrix}x\\y\\z\end{pmatrix}.
  1. 通过正交变换
(xyz)=Q(stu)\begin{pmatrix}x\\y\\z\end{pmatrix} = Q\begin{pmatrix}s\\t\\u\end{pmatrix}

ff 化为标准形

a1s2+a2t2+a3u2,a1a2a3.a_1s^2+a_2t^2+a_3u^2, \qquad a_1\leq a_2\leq a_3.

求正交矩阵 QQ 以及 a1,a2,a3a_1,a_2,a_3。 3. 证明二次曲面

f(x,y,z)23x23y+23z8=0f(x,y,z) -2\sqrt3x-2\sqrt3y+2\sqrt3z-8=0

是椭球面。

Kai

問題の解答と詳細な解説

(1) 行列の相似と固有値・固有ベクトル

(a) 相似な行列が同じ固有値を持つことの証明

行列 AABB が相似であるとき、正則行列 PP を用いて B=P1APB = P^{-1}AP と表されます。 ある行列の固有値は、その固有方程式(特性方程式) λIX=0|\lambda I - X| = 0 の解として求められます。ここで、 II は単位行列です。

行列 BB の特性多項式 λIB|\lambda I - B| を変形して、行列 AA の特性多項式と一致することを示します。

λIB=λIP1AP=P1(λI)PP1AP(λI はスカラー行列なので可換)=P1(λIA)P\begin{aligned} |\lambda I - B| &= |\lambda I - P^{-1}AP| \\ &= |P^{-1}(\lambda I)P - P^{-1}AP| \quad (\lambda I \text{ はスカラー行列なので可換}) \\ &= |P^{-1}(\lambda I - A)P| \end{aligned}

行列式の性質 XYZ=XYZ|XYZ| = |X||Y||Z| を用いると、

P1(λIA)P=P1λIAP=P1λIAP=λIA\begin{aligned} |P^{-1}(\lambda I - A)P| &= |P^{-1}| |\lambda I - A| |P| \\ &= |P|^{-1} |\lambda I - A| |P| \\ &= |\lambda I - A| \end{aligned}

となります。

このように、行列 AA と行列 BB の特性多項式は完全に一致する( λIA=λIB|\lambda I - A| = |\lambda I - B| )ため、方程式の解である「固有値」もすべて等しくなります。(証明終)

(b) 行列 A,BA, B の固有値と長さ 1 の固有ベクトルの導出

与えられた行列は A=(4321)A = \begin{pmatrix} 4 & -3 \\ 2 & -1 \end{pmatrix} , P=(2513)P = \begin{pmatrix} 2 & 5 \\ 1 & 3 \end{pmatrix} です。

【行列 AA について】

  1. 固有値の計算 AA の特性方程式 AλI=0|A - \lambda I| = 0 を解きます。
4λ321λ=(4λ)(1λ)(3×2)=λ23λ4+6=λ23λ+2=0\begin{vmatrix} 4 - \lambda & -3 \\ 2 & -1 - \lambda \end{vmatrix} = (4 - \lambda)(-1 - \lambda) - (-3 \times 2) = \lambda^2 - 3\lambda - 4 + 6 = \lambda^2 - 3\lambda + 2 = 0
(λ1)(λ2)=0(\lambda - 1)(\lambda - 2) = 0

よって、 AA の固有値は λ=1,2\lambda = 1, 2 です。

  1. 固有ベクトルの計算 ・ λ=1\lambda = 1 のとき: (AI)x=(3322)(xy)=(00)(A - I)\mathbf{x} = \begin{pmatrix} 3 & -3 \\ 2 & -2 \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} より、 x=yx = y 。 固有ベクトルを定数 cc を用いて c(11)c\begin{pmatrix} 1 \\ 1 \end{pmatrix} と表せます。長さ(ノルム)を 1 にするため、 12+12=2\sqrt{1^2 + 1^2} = \sqrt{2} で割ります。 λ=1\lambda=1 に対応する長さ1の固有ベクトル: ±12(11)\pm\frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 1 \end{pmatrix}

    λ=2\lambda = 2 のとき: (A2I)x=(2323)(xy)=(00)(A - 2I)\mathbf{x} = \begin{pmatrix} 2 & -3 \\ 2 & -3 \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} より、 2x=3y2x = 3y 。 固有ベクトルは c(32)c\begin{pmatrix} 3 \\ 2 \end{pmatrix} と表せます。ノルム 32+22=13\sqrt{3^2 + 2^2} = \sqrt{13} で割ります。 λ=2\lambda=2 に対応する長さ1の固有ベクトル: ±113(32)\pm\frac{1}{\sqrt{13}}\begin{pmatrix} 3 \\ 2 \end{pmatrix}

【行列 BB について】

(a) の結果より、 BB の固有値も AA と同じく λ=1,2\lambda = 1, 2 です。 By=λyB\mathbf{y} = \lambda\mathbf{y}(P1AP)y=λy(P^{-1}AP)\mathbf{y} = \lambda\mathbf{y} であり、左から PP を掛けると A(Py)=λ(Py)A(P\mathbf{y}) = \lambda(P\mathbf{y}) となります。 すなわち、 BB の固有ベクトル y\mathbf{y} は、 y=P1x\mathbf{y} = P^{-1}\mathbf{x}x\mathbf{x}AA の固有ベクトル)として求められます。

まず、 P1P^{-1} を求めます。 P=2×35×1=1|P| = 2\times3 - 5\times1 = 1 より、

P1=11(3512)=(3512)P^{-1} = \frac{1}{1}\begin{pmatrix} 3 & -5 \\ -1 & 2 \end{pmatrix} = \begin{pmatrix} 3 & -5 \\ -1 & 2 \end{pmatrix}

λ=1\lambda = 1 のとき: ベースとなるベクトルとして x=(11)\mathbf{x} = \begin{pmatrix} 1 \\ 1 \end{pmatrix} を用います。 y=P1x=(3512)(11)=(21)\mathbf{y} = P^{-1}\mathbf{x} = \begin{pmatrix} 3 & -5 \\ -1 & 2 \end{pmatrix}\begin{pmatrix} 1 \\ 1 \end{pmatrix} = \begin{pmatrix} -2 \\ 1 \end{pmatrix} これを正規化(ノルム (2)2+12=5\sqrt{(-2)^2 + 1^2} = \sqrt{5} )します。 λ=1\lambda=1 に対応する長さ1の固有ベクトル: ±15(21)\pm\frac{1}{\sqrt{5}}\begin{pmatrix} -2 \\ 1 \end{pmatrix}

λ=2\lambda = 2 のとき: ベースとなるベクトルとして x=(32)\mathbf{x} = \begin{pmatrix} 3 \\ 2 \end{pmatrix} を用います。 y=P1x=(3512)(32)=(9103+4)=(11)\mathbf{y} = P^{-1}\mathbf{x} = \begin{pmatrix} 3 & -5 \\ -1 & 2 \end{pmatrix}\begin{pmatrix} 3 \\ 2 \end{pmatrix} = \begin{pmatrix} 9 - 10 \\ -3 + 4 \end{pmatrix} = \begin{pmatrix} -1 \\ 1 \end{pmatrix} これを正規化(ノルム (1)2+12=2\sqrt{(-1)^2 + 1^2} = \sqrt{2} )します。 λ=2\lambda=2 に対応する長さ1の固有ベクトル: ±12(11)\pm\frac{1}{\sqrt{2}}\begin{pmatrix} -1 \\ 1 \end{pmatrix}

(※ 固有ベクトルは定数倍の自由度があるため、符号が反転したマイナスのものも正解となります。)

(2) 二次形式と直交変換

(a) 実対称行列 TT の導出

二次形式 f(x,y,z)=2x2+2y2+4z2+2xy+2yz+2zxf(x, y, z) = 2x^2 + 2y^2 + 4z^2 + 2xy + 2yz + 2zx を行列で表します。 各変数の2乗の係数が対角成分に、交差項( xy,yz,zxxy, yz, zx )の係数の半分が非対角成分に対応します。

T=(x2の係数xyの係数/2zxの係数/2xyの係数/2y2の係数yzの係数/2zxの係数/2yzの係数/2z2の係数)=(211121114)T = \begin{pmatrix} x^2の係数 & xyの係数/2 & zxの係数/2 \\ xyの係数/2 & y^2の係数 & yzの係数/2 \\ zxの係数/2 & yzの係数/2 & z^2の係数 \end{pmatrix} = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 4 \end{pmatrix}

よって、 T=(211121114)T = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 4 \end{pmatrix}

(b) 直交行列 QQ と標準形の係数 a1,a2,a3a_1, a_2, a_3 の導出

二次形式の標準形の係数 a1,a2,a3a_1, a_2, a_3 は、対称行列 TT の固有値です。

  1. TT の固有値の計算
TλI=2λ1112λ1114λ=0|T - \lambda I| = \begin{vmatrix} 2-\lambda & 1 & 1 \\ 1 & 2-\lambda & 1 \\ 1 & 1 & 4-\lambda \end{vmatrix} = 0

計算を簡略化するため、1行目から2行目を引くと、

1λ1+λ012λ1114λ=(1λ)11012λ1114λ\begin{vmatrix} 1-\lambda & -1+\lambda & 0 \\ 1 & 2-\lambda & 1 \\ 1 & 1 & 4-\lambda \end{vmatrix} = (1-\lambda)\begin{vmatrix} 1 & -1 & 0 \\ 1 & 2-\lambda & 1 \\ 1 & 1 & 4-\lambda \end{vmatrix}

1列目に2列目を足すと、

(1λ)0103λ2λ1214λ=(1λ)13λ124λ=(1λ){(3λ)(4λ)2}(1-\lambda)\begin{vmatrix} 0 & -1 & 0 \\ 3-\lambda & 2-\lambda & 1 \\ 2 & 1 & 4-\lambda \end{vmatrix} = (1-\lambda) \cdot 1 \cdot \begin{vmatrix} 3-\lambda & 1 \\ 2 & 4-\lambda \end{vmatrix} = (1-\lambda)\{ (3-\lambda)(4-\lambda) - 2 \}
=(1λ)(λ27λ+10)=(1λ)(λ2)(λ5)=0= (1-\lambda)(\lambda^2 - 7\lambda + 10) = (1-\lambda)(\lambda - 2)(\lambda - 5) = 0

よって、固有値は λ=1,2,5\lambda = 1, 2, 5 です。 a1a2a3a_1 \le a_2 \le a_3 の条件より、 a1=1,a2=2,a3=5a_1 = 1, a_2 = 2, a_3 = 5 です。

  1. 直交行列 QQ の計算(固有ベクトルの導出)

対称行列の異なる固有値に属する固有ベクトルは互いに直交します。それぞれの単位固有ベクトルを求めます。

λ1=1\lambda_1 = 1 のとき: TI=(111111113){x+y+z=0x+y+3z=0z=0,y=xT - I = \begin{pmatrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 3 \end{pmatrix} \rightarrow \begin{cases} x + y + z = 0 \\ x + y + 3z = 0 \end{cases} \rightarrow z = 0, y = -x 単位固有ベクトル u1=12(110)\mathbf{u}_1 = \frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ -1 \\ 0 \end{pmatrix}

λ2=2\lambda_2 = 2 のとき: T2I=(011101112){y+z=0x+z=0x=y=zT - 2I = \begin{pmatrix} 0 & 1 & 1 \\ 1 & 0 & 1 \\ 1 & 1 & 2 \end{pmatrix} \rightarrow \begin{cases} y + z = 0 \\ x + z = 0 \end{cases} \rightarrow x = y = -z 単位固有ベクトル u2=13(111)\mathbf{u}_2 = \frac{1}{\sqrt{3}}\begin{pmatrix} 1 \\ 1 \\ -1 \end{pmatrix}

λ3=5\lambda_3 = 5 のとき: T5I=(311131111){3x+y+z=0x3y+z=0T - 5I = \begin{pmatrix} -3 & 1 & 1 \\ 1 & -3 & 1 \\ 1 & 1 & -1 \end{pmatrix} \rightarrow \begin{cases} -3x + y + z = 0 \\ x - 3y + z = 0 \end{cases} 第1式から第2式を引くと 4x+4y=0x=y-4x + 4y = 0 \rightarrow x = y 。代入すると z=2xz = 2x 。 単位固有ベクトル u3=16(112)\mathbf{u}_3 = \frac{1}{\sqrt{6}}\begin{pmatrix} 1 \\ 1 \\ 2 \end{pmatrix}

これらを列ベクトルとして並べて直交行列 QQ を作ります。 Q=(12131612131601326)Q = \begin{pmatrix} \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{6}} \\ -\frac{1}{\sqrt{2}} & \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{6}} \\ 0 & -\frac{1}{\sqrt{3}} & \frac{2}{\sqrt{6}} \end{pmatrix} (※各列ベクトルは符号が反転していても直交行列として成立します。)

(c) 二次曲面が楕円面であることの証明

与えられた方程式は f(x,y,z)23x23y+23z8=0f(x, y, z) - 2\sqrt{3}x - 2\sqrt{3}y + 2\sqrt{3}z - 8 = 0 です。 (b) で求めた変数変換 (xyz)=Q(stu)\begin{pmatrix} x \\ y \\ z \end{pmatrix} = Q \begin{pmatrix} s \\ t \\ u \end{pmatrix} を適用します。

まず、二次形式の部分は f(x,y,z)=s2+2t2+5u2f(x,y,z) = s^2 + 2t^2 + 5u^2 に標準化されます。

次に、一次の項(線形項)を変換します。一次項を行列表記すると、

(232323)(xyz)=(232323)Q(stu)\begin{pmatrix} -2\sqrt{3} & -2\sqrt{3} & 2\sqrt{3} \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} -2\sqrt{3} & -2\sqrt{3} & 2\sqrt{3} \end{pmatrix} Q \begin{pmatrix} s \\ t \\ u \end{pmatrix}

係数ベクトルと行列 QQ の積を計算します。

(232323)(12131612131601326)\begin{pmatrix} -2\sqrt{3} & -2\sqrt{3} & 2\sqrt{3} \end{pmatrix} \begin{pmatrix} \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{6}} \\ -\frac{1}{\sqrt{2}} & \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{6}} \\ 0 & -\frac{1}{\sqrt{3}} & \frac{2}{\sqrt{6}} \end{pmatrix}
  • ss の係数: 232+232+0=0\frac{-2\sqrt{3}}{\sqrt{2}} + \frac{2\sqrt{3}}{\sqrt{2}} + 0 = 0
  • tt の係数: 233+233+233=222=6\frac{-2\sqrt{3}}{\sqrt{3}} + \frac{-2\sqrt{3}}{\sqrt{3}} + \frac{-2\sqrt{3}}{\sqrt{3}} = -2 - 2 - 2 = -6
  • uu の係数: 236+236+436=43+436=0\frac{-2\sqrt{3}}{\sqrt{6}} + \frac{-2\sqrt{3}}{\sqrt{6}} + \frac{4\sqrt{3}}{\sqrt{6}} = \frac{-4\sqrt{3} + 4\sqrt{3}}{\sqrt{6}} = 0

したがって、線形項の変換後は 6t-6t となります。 方程式全体を s,t,us, t, u で書き直すと次のようになります。

s2+2t2+5u26t8=0s^2 + 2t^2 + 5u^2 - 6t - 8 = 0

この式について、 tt に関して平方完成を行います。

s2+2(t23t)+5u2=8s^2 + 2(t^2 - 3t) + 5u^2 = 8
s2+2(t32)22(32)2+5u2=8s^2 + 2\left(t - \frac{3}{2}\right)^2 - 2 \left(\frac{3}{2}\right)^2 + 5u^2 = 8
s2+2(t32)292+5u2=8s^2 + 2\left(t - \frac{3}{2}\right)^2 - \frac{9}{2} + 5u^2 = 8
s2+2(t32)2+5u2=252s^2 + 2\left(t - \frac{3}{2}\right)^2 + 5u^2 = \frac{25}{2}

両辺を 252\frac{25}{2} で割って標準形にします。

s225/2+(t32)225/4+u25/2=1\frac{s^2}{25/2} + \frac{\left(t - \frac{3}{2}\right)^2}{25/4} + \frac{u^2}{5/2} = 1

この方程式は、定数 A2,B2,C2>0A^2, B^2, C^2 > 0 を用いて X2A2+Y2B2+Z2C2=1\frac{X^2}{A^2} + \frac{Y^2}{B^2} + \frac{Z^2}{C^2} = 1 の形(各変数の2乗の係数がすべて正)で表されるため、楕円面(ellipsoid)の標準形であることがわかります。

よって、与えられた二次曲面は楕円面であることが示されました。(証明終)

图像边界:原始 JSON 引用了 question_images/东京工业大学_情报通信_2019_20260208_093713_461467.png,但该文件没有随 JSON 一并提供。本文只保留能够由文字和公式确认的内容。