Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

対称行列の直交行列による対角化

証明

(1)    (2)(1) \implies (2)について:

AAは複素数の範囲では重複をふくめてnn個の固有値をもつという定理があり、また実対称行列の固有値はすべて実数であるという定理も存在するため、AAの固有値はすべて実数である。

三角化に関する定理より、直交行列PPが存在し、

P1AP=(λ1Oλn)P^{-1} A P=\left(\begin{array}{ccc} \lambda_1 & & * \\ & \ddots & \\ O & & \lambda_n \end{array}\right)

と三角化できる。

仮定よりPT=P1P^T = P^{-1}であり、AT=AA^T=Aであるから、

(P1AP)T=(PTAP)T=PTAT(PT)T=P1AP(P^{-1} A P)^T = (P^T A P)^T\\ = P^T A^T (P^T)^T\\ = P^{-1} A P\\

すなわち、三角行列P1APP^{-1} A Pも対称行列であるため、非対角要素はすべて0であることがわかる

(2)    (1)(2) \implies (1)について

AAは直交行列PPによって

P1AP=(λ1OOλn)P^{-1} A P=\left(\begin{array}{lll} \lambda_1 & & O \\ & \ddots & \\ O & & \lambda_n \end{array}\right)

と対角化できるとする。

対角行列であるから、(P1AP)T=P1AP(P^{-1} A P)^T = P^{-1} A Pを満たす。

P1=PTP^{-1} = P^Tであるから、

(P1AP)T=(PTAP)T=PTATP=P1ATP(P^{-1} A P)^T = (P^T A P)^T = P^T A^T P = P^{-1} A^T P

であり、(P1AP)T=P1ATP(P^{-1} A P)^T = P^{-1} A^T Pに左からPP、右からP1P^{-1}を掛けることでA=ATA=A^Tが得られ、AAは対称行列であることがわかる

解答
AtE=t0101t010t=(t1)(t+1)2|A-t E|=\left|\begin{array}{ccc} -t & 0 & 1 \\ 0 & -1-t & 0 \\ 1 & 0 & -t \end{array}\right|=-(t-1)(t+1)^2

ゆえに固有値は 1,11,-1 (重複度2) である. 次に固有値 1,11,-1 に属する固有空間 V(1),V(1)V(1), V(-1) を求める。固有値1のとき,

(1010110101)(xyz)=(x+z2yxz)=(000)\left(\begin{array}{rcr} -1 & 0 & 1 \\ 0 & -1-1 & 0 \\ 1 & 0 & -1 \end{array}\right)\left(\begin{array}{l} x \\ y \\ z \end{array}\right)=\left(\begin{array}{c} -x+z \\ -2 y \\ x-z \end{array}\right)=\left(\begin{array}{l} 0 \\ 0 \\ 0 \end{array}\right)

すなわち, x=z,y=0x=z, y=0 であるから

V(1)={c(101)c は任意 }V(1)=\left\{\left.c\left(\begin{array}{l} 1 \\ 0 \\ 1 \end{array}\right) \right\rvert\, c \text { は任意 }\right\}

となる. 固有値 -1 のとき,

((1)0101(1)010(1))(xyz)=(x+z0x+z)=(000)\left(\begin{array}{ccc} -(-1) & 0 & 1 \\ 0 & -1-(-1) & 0 \\ 1 & 0 & -(-1) \end{array}\right)\left(\begin{array}{l} x \\ y \\ z \end{array}\right)=\left(\begin{array}{c} x+z \\ 0 \\ x+z \end{array}\right)=\left(\begin{array}{l} 0 \\ 0 \\ 0 \end{array}\right)

すなわち, x+z=0,yx+z=0, y は任意, であるから

V(1)={c1(010)+c2(101)c1,c2 は任意 }V(-1)=\left\{\left.c_1\left(\begin{array}{l} 0 \\ 1 \\ 0 \end{array}\right)+c_2\left(\begin{array}{r} 1 \\ 0 \\ -1 \end{array}\right) \right\rvert\, c_1, c_2 \text { は任意 }\right\}

となる. 以上の固有ベクトルの集合のうちから正規直交基底を選んで並べると対角化を与える直交行列になる。 そこで

a1=(101),a2=(010),a3=(101)\boldsymbol{a}_1=\left(\begin{array}{l}1 \\ 0 \\ 1\end{array}\right), \quad \boldsymbol{a}_2=\left(\begin{array}{l}0 \\ 1 \\ 0\end{array}\right), \quad \boldsymbol{a}_3=\left(\begin{array}{r}1 \\ 0 \\ -1\end{array}\right)

と おいてシュミットの方法で正規直交化を行う。

定理より V(1)V(1)V(1)V(-1) の元は互いに直交しているから、それぞれで正規直交化を行えばよい。今回の場合、 a2a_2a3a_3 はたまたま直交しているので長さのみ調節すればよい。こうして次の正規直交基底を得る。

v1=12(101),v2=(010),v3=12(101)\boldsymbol{v}_1=\frac{1}{\sqrt{2}}\left(\begin{array}{l} 1 \\ 0 \\ 1 \end{array}\right), \quad \boldsymbol{v}_2=\left(\begin{array}{l} 0 \\ 1 \\ 0 \end{array}\right), \quad \boldsymbol{v}_3=\frac{1}{\sqrt{2}}\left(\begin{array}{r} 1 \\ 0 \\ -1 \end{array}\right)

これらを列ベクトルにもつ行列

P=(1/201/20101/201/2)P=\left(\begin{array}{ccc} 1 / \sqrt{2} & 0 & 1 / \sqrt{2} \\ 0 & 1 & 0 \\ 1 / \sqrt{2} & 0 & -1 / \sqrt{2} \end{array}\right)

は定理より直交行列であり、 行列 AA の対角化

P1AP=(100010001)P^{-1} A P=\left(\begin{array}{rrr} 1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & -1 \end{array}\right)

が得られる。