東京大学 情報理工学研究科 数理情報学 2020年度 第1問
Author
Description
\(n\) を正の整数とし、\(A\) は \(n\) 次の実正方行列、\(B\) は \(n\) 次の正定値実対称行列であるとする。以下の問いに答えよ。
(1) 条件
\[
BC + CB = A
\]
を満たす \(n\) 次の実正方行列 \(C\) がただ一つ存在することを示せ。以降では、この行列 \(C\) を \(C_{A,B}\) で表す。
(2) \(BC_{A,B} = C_{A,B}B\) と \(AB = BA\) が同値であることを示せ。
Kai
(1)
\(B=QDQ^T \ (\mathrm{where} \; D \; \mathrm{is}\; \mathrm{diagonal}, Q^TQ=I)\) とする。
\[
\begin{aligned}
& QDQ^TC+CQDQ^T=A \\
\Leftrightarrow & Q^TQDQ^TCQ+Q^TCQDQ^TQ=Q^TAQ \\
\Leftrightarrow & DQ^TCQ+Q^TCQD=Q^TAQ \\
\Leftrightarrow & DC'+C'D=A' \\
\end{aligned}
\]
よって、対角行列のみを考えれば良い。
\[
\begin{aligned}
& A_{ij}=2b_{ii}c_{ij} \\
\Leftrightarrow & c_{ij}=\frac{A_{ij}}{2b_{ii}}\quad (\because \; b_{ii} \geq 0) \\
\end{aligned}
\]
よって、\(C\) は一意。
(2)
\(BC_{A,B}=C_{A,B}B \Rightarrow AB=BA\) は自明。
\(AB=BA \Rightarrow BC_{A,B}=C_{A,B}B\) は
\[
\begin{aligned}
& AB=BA \\
\Leftrightarrow & CBB=BBC \\
\Leftrightarrow & CQD^2Q^T=QD^2Q^TC \\
\Leftrightarrow & \sum_k (CQ)_{ik}D^2_{kk}Q^T_{kj}=\sum_k Q_{ik}D^2_{kk}(Q^TC)_{kj} \\
\Leftrightarrow & \sum_k (CQ)_{ik}D_{kk}Q^T_{kj}=\sum_k Q_{ik}D_{kk}(Q^TC)_{kj} \\
\Leftrightarrow & CQDQ^T=QDQ^TC \\
\Leftrightarrow & CB=BC
\end{aligned}
\]
より、示される。
Knowledge
正定値性は、基本的には対称行列かエルミート行列 (つまり、自己随伴行列、\(A=A^{\dagger} \Leftrightarrow A^T=\overline{A}\)) のみに限って考えることが多い。
実対称行列の正定値性は以下の4条件が同値。
- (1) \(x^TAx>0 \quad (\forall x\in \mathbb{R}^n\setminus {\boldsymbol{0}})\)
- (2) \(A\) の固有値が全て正 (固有ベクトルに対して、\(\langle Ax,x \rangle = \langle \lambda x,x \rangle =\lambda |x|^2\)、逆も (3) より従う)
- (3) ある直交行列 \(Q(Q^TQ=QQ^T=I)\) と、対角成分が全て正の対角行列 \(D\) を用いて、\(A=QDQ^T\) と表せる
- (4) ある正則行列 \(S\) を用いて、\(A=S^TS\) と表せる
以下の性質が上記の証明で重要。
\(A=A^T\) ならば、\(P^{-1}AP\) と対角化可能。ただし、\(P\) は直交行列で、\(P^TP=PP^T=I\)
- (1) 固有値がすべて異なる場合、固有ベクトル \(\{p_k\}\) は自動的に直交するので、大きさが１になるように選ぶことにより \((r_k=\frac{1}{|p_k|}p_k)\)、\(R=[r_1, r_2, \cdots, r_n]\) は直交行列となり、この \(R\) を用いて、\(R^{-1}AR\) を対角行列にできる。
- (2) 固有値に重複がある場合にも、対称行列では、重複する固有値に属する１次独立な固有ベクトルを重複度分だけ見つけることが常に可能それらをグラム・シュミットの直交化法により正規直交化すれば、他の固有ベクトルも合わせてすべての固有ベクトルが直交することとなり、それらを規格化すればやはり直交行列 \(R\) が得られる。
また、正則であり、逆行列も正定値行列。
\[
\begin{aligned}
A^{-1}=(UDU^{-1})^{-1}=UD^{-1}U^{-1}
\end{aligned}
\]
ちなみに、本問はシルベスター方程式というものらしい。