跳到主要内容

明治大学 先端数理科学研究科 現象数理学専攻 2019年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

(1) 実数を成分とする行列 AAtA=A{}^t A = -A を満たすとする。ただし, tA{}^t AAA の転置行列を表す。

(a) AA の固有値は純虚数であることを示せ。

(b) AA の相異なる固有値 λ,μ\lambda, \mu に対し, λ\lambda に対する固有ベクトル xxμ\mu に対する固有ベクトル yy は直交することを示せ。

(2) A=(011101110)A = \begin{pmatrix} 0 & 1 & 1 \\ -1 & 0 & 1 \\ -1 & -1 & 0 \end{pmatrix} とする。

(a) AA の固有値をすべて求めよ。

(b) AA の各固有値に対する固有ベクトルを1つずつ求めよ。

(c) AA を対角化するユニタリ行列 UU を1つ求めよ。

题目描述

(1) 设实矩阵 AA 满足

tA=A,{}^tA=-A,

其中 tA{}^tA 表示 AA 的转置矩阵。

(a) 证明 AA 的特征值均为纯虚数。

(b) 对于 AA 的两个不同特征值 λ,μ\lambda,\mu,设 xx 是对应于 λ\lambda 的特征向量,yy 是对应于 μ\mu 的特征向量,证明 xxyy 正交。

(2) 设

A=(011101110).A=\begin{pmatrix} 0&1&1\\ -1&0&1\\ -1&-1&0 \end{pmatrix}.

(a) 求 AA 的全部特征值。

(b) 对 AA 的每个特征值,各求一个对应的特征向量。

(c) 求一个将 AA 对角化的酉矩阵 UU

Kai

(1) の解答

(a) AA の固有値は純虚数であることの証明

AA の固有値を λ\lambda 、対応する固有ベクトルを x(x0)x (x \neq \boldsymbol{0}) とする。定義より、 Ax=λxAx = \lambda x が成り立つ。 この式の両辺のエルミート共役(随伴)をとると、 (Ax)H=(λx)H(Ax)^H = (\lambda x)^H となる。 これは、 xHAH=λˉxHx^H A^H = \bar{\lambda} x^H となる。

行列 AA は実数成分を持つため、 AH=(tAˉ)=tAA^H = ({}^t \bar{A}) = {}^t A である。 また、問題の条件より tA=A{}^t A = -A なので、 AH=AA^H = -A となる。 これを代入すると、 xHA=λˉxH-x^H A = \bar{\lambda} x^H を得る。

元の式 Ax=λxAx = \lambda x の両辺に左から xHx^H をかけると、

xHAx=xH(λx)=λ(xHx)(1)x^H A x = x^H (\lambda x) = \lambda (x^H x) \quad (1)

随伴をとった式 xHA=λˉxH-x^H A = \bar{\lambda} x^H の両辺に右から xx をかけると、

xHAx=(λˉxH)x=λˉ(xHx)(2)-x^H A x = (\bar{\lambda} x^H) x = \bar{\lambda} (x^H x) \quad (2)

(1) と (2) より、 λ(xHx)=λˉ(xHx)\lambda (x^H x) = -\bar{\lambda} (x^H x) となる。 これを整理すると、 (λ+λˉ)(xHx)=0(\lambda + \bar{\lambda})(x^H x) = 0 となる。

xx は固有ベクトルなので x0x \neq \boldsymbol{0} であり、 xHx=x2>0x^H x = ||x||^2 > 0 である。 したがって、 λ+λˉ=0\lambda + \bar{\lambda} = 0 でなければならない。 λ=a+bi\lambda = a + bi ( a,ba, b は実数) とおくと、 λˉ=abi\bar{\lambda} = a - bi なので、 λ+λˉ=2a=0\lambda + \bar{\lambda} = 2a = 0 となり、 a=0a=0 が導かれる。 よって、 λ=bi\lambda = bi となり、 AA の固有値は純虚数である( 00 を含む)。

(b) 固有ベクトルの直交性の証明

AA の相異なる固有値 λ,μ\lambda, \mu に対応する固有ベクトルをそれぞれ x,yx, y とする。 Ax=λx(3)Ax = \lambda x \quad (3) Ay=μy(4)Ay = \mu y \quad (4)

(3) のエルミート共役をとると、 xHAH=λˉxHx^H A^H = \bar{\lambda} x^H となる。(a) と同様に AH=AA^H = -A であるから、 xHA=λˉxH-x^H A = \bar{\lambda} x^H となる。 この式の両辺に右から yy をかけると、

xHAy=λˉ(xHy)(5)-x^H A y = \bar{\lambda} (x^H y) \quad (5)

(4) を (5) の左辺に代入すると、

xH(μy)=λˉ(xHy)-x^H (\mu y) = \bar{\lambda} (x^H y)
μ(xHy)=λˉ(xHy)-\mu (x^H y) = \bar{\lambda} (x^H y)
(λˉ+μ)(xHy)=0(\bar{\lambda} + \mu) (x^H y) = 0

(a) の結果から、固有値は純虚数である。 λ=ik1,μ=ik2\lambda = ik_1, \mu = ik_2 ( k1,k2Rk_1, k_2 \in \mathbb{R} ) とおける。 すると、 λˉ=ik1=λ\bar{\lambda} = -ik_1 = -\lambda となる。 したがって、 (λ+μ)(xHy)=0(-\lambda + \mu)(x^H y) = 0 となる。

問題の条件より λμ\lambda \neq \mu なので、 λ+μ0-\lambda + \mu \neq 0 である。 よって、 xHy=0x^H y = 0 でなければならない。これは、ベクトル xxyy が直交することを示している。

(2) の解答

A=(011101110)A = \begin{pmatrix} 0 & 1 & 1 \\ -1 & 0 & 1 \\ -1 & -1 & 0 \end{pmatrix} とする。

(a) 固有値の計算

固有方程式 det(AλI)=0\det(A - \lambda I) = 0 を解く。

det(AλI)=λ111λ111λ=λ(λ2+1)1(λ+1)+1(1λ)=λ3λλ1+1λ=λ33λ\det(A - \lambda I) = \begin{vmatrix} -\lambda & 1 & 1 \\ -1 & -\lambda & 1 \\ -1 & -1 & -\lambda \end{vmatrix} = -\lambda(\lambda^2+1) - 1(\lambda+1) + 1(1-\lambda) = -\lambda^3 - \lambda - \lambda - 1 + 1 - \lambda = -\lambda^3 - 3\lambda

よって、固有方程式は λ(λ2+3)=0-\lambda(\lambda^2 + 3) = 0 である。 これを解くと、固有値は λ1=0,λ2=i3,λ3=i3\lambda_1 = 0, \lambda_2 = i\sqrt{3}, \lambda_3 = -i\sqrt{3} となる。

(b) 固有ベクトルの計算

  • ** λ1=0\lambda_1 = 0 の場合**: (A0I)x1=0(A-0I)x_1 = \boldsymbol{0} を解く。
(011101110)(xyz)=(000)    {y+z=0x+z=0xy=0\begin{pmatrix} 0 & 1 & 1 \\ -1 & 0 & 1 \\ -1 & -1 & 0 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \implies \begin{cases} y+z=0 \\ -x+z=0 \\ -x-y=0 \end{cases}

z=tz=t とおくと、 x=t,y=tx=t, y=-t となる。 t=1t=1 とすると、固有ベクトル x1=(111)x_1 = \begin{pmatrix} 1 \\ -1 \\ 1 \end{pmatrix} を得る。

  • ** λ2=i3\lambda_2 = i\sqrt{3} の場合**: (Ai3I)x2=0(A-i\sqrt{3}I)x_2 = \boldsymbol{0} を解く。
(i3111i3111i3)(xyz)=(000)\begin{pmatrix} -i\sqrt{3} & 1 & 1 \\ -1 & -i\sqrt{3} & 1 \\ -1 & -1 & -i\sqrt{3} \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

1行目から2行目を引くと (1i3)x+(1+i3)y=0(1-i\sqrt{3})x + (1+i\sqrt{3})y = 0x=1+i3x=1+i\sqrt{3} とすると y=(1i3)=1+i3y=-(1-i\sqrt{3})=-1+i\sqrt{3} 。これを3行目の式 xyi3z=0-x-y-i\sqrt{3}z=0 に代入すると、 (1+i3)(1+i3)i3z=0    2i3i3z=0    z=2-(1+i\sqrt{3}) - (-1+i\sqrt{3}) - i\sqrt{3}z = 0 \implies -2i\sqrt{3} - i\sqrt{3}z = 0 \implies z=-2 。よって固有ベクトル x2=(1+i31+i32)x_2 = \begin{pmatrix} 1+i\sqrt{3} \\ -1+i\sqrt{3} \\ -2 \end{pmatrix} を得る。

  • ** λ3=i3\lambda_3 = -i\sqrt{3} の場合**: λ3=λˉ2\lambda_3 = \bar{\lambda}_2 であり、行列 AA が実数行列なので、対応する固有ベクトルは x2x_2 の複素共役となる。 よって、固有ベクトル x3=xˉ2=(1i31i32)x_3 = \bar{x}_2 = \begin{pmatrix} 1-i\sqrt{3} \\ -1-i\sqrt{3} \\ -2 \end{pmatrix} を得る。

(c) ユニタリ行列の計算

(b)で求めた固有ベクトルを正規化し、それらを列ベクトルとするユニタリ行列 UU を構成する。

  • x1x_1 のノルム: x1=12+(1)2+12=3||x_1|| = \sqrt{1^2+(-1)^2+1^2} = \sqrt{3}
  • x2x_2 のノルム: x2=1+i32+1+i32+22=(1+3)+(1+3)+4=12=23||x_2|| = \sqrt{|1+i\sqrt{3}|^2 + |-1+i\sqrt{3}|^2 + |-2|^2} = \sqrt{(1+3)+(1+3)+4} = \sqrt{12} = 2\sqrt{3}
  • x3x_3 のノルム: x3=xˉ2=x2=23||x_3|| = ||\bar{x}_2|| = ||x_2|| = 2\sqrt{3}

正規化された固有ベクトル u1,u2,u3u_1, u_2, u_3 は、

u1=13(111),u2=123(1+i31+i32),u3=123(1i31i32)u_1 = \frac{1}{\sqrt{3}} \begin{pmatrix} 1 \\ -1 \\ 1 \end{pmatrix}, \quad u_2 = \frac{1}{2\sqrt{3}} \begin{pmatrix} 1+i\sqrt{3} \\ -1+i\sqrt{3} \\ -2 \end{pmatrix}, \quad u_3 = \frac{1}{2\sqrt{3}} \begin{pmatrix} 1-i\sqrt{3} \\ -1-i\sqrt{3} \\ -2 \end{pmatrix}

これらを列ベクトルとして並べ、ユニタリ行列 UU を得る。

U=[u1,u2,u3]=(131+i3231i323131+i3231i32313223223)=(131+i3231i323131+i3231i323131313)U = [u_1, u_2, u_3] = \begin{pmatrix} \frac{1}{\sqrt{3}} & \frac{1+i\sqrt{3}}{2\sqrt{3}} & \frac{1-i\sqrt{3}}{2\sqrt{3}} \\ -\frac{1}{\sqrt{3}} & \frac{-1+i\sqrt{3}}{2\sqrt{3}} & \frac{-1-i\sqrt{3}}{2\sqrt{3}} \\ \frac{1}{\sqrt{3}} & -\frac{2}{2\sqrt{3}} & -\frac{2}{2\sqrt{3}} \end{pmatrix} = \begin{pmatrix} \frac{1}{\sqrt{3}} & \frac{1+i\sqrt{3}}{2\sqrt{3}} & \frac{1-i\sqrt{3}}{2\sqrt{3}} \\ -\frac{1}{\sqrt{3}} & \frac{-1+i\sqrt{3}}{2\sqrt{3}} & \frac{-1-i\sqrt{3}}{2\sqrt{3}} \\ \frac{1}{\sqrt{3}} & -\frac{1}{\sqrt{3}} & -\frac{1}{\sqrt{3}} \end{pmatrix}

この行列 UU によって AA は対角化され、 UHAU=(0000i3000i3)U^H A U = \begin{pmatrix} 0 & 0 & 0 \\ 0 & i\sqrt{3} & 0 \\ 0 & 0 & -i\sqrt{3} \end{pmatrix} となる。