4.5.2.0.1.1 ステップ1

($A^{k+1}$ の QR 分解) まず

[

l]主張1 $P_k\DefEq Q_0Q_1\cdots Q_{k-1}Q_k$ とおくと

$$A_k=P_{k-1}^\ast A P_{k-1}.$$ (4.10)

実際、

$$\begin{eqnarray*} A_1&=&R_0 Q_0=(Q_0^\ast A)Q_0=Q_0^\ast A Q_0,\\ A_2&=&R_1 Q... ...\cdots Q_1^\ast Q_0^\ast A Q_0 Q_1\cdots Q_{k-1})Q_k,\\ \cdots \end{eqnarray*}$$

となる (帰納法で証明すべきかもしれないが)。

[

l]主張2 $U_k\DefEq R_k R_{k-1}\cdots R_1 R_0$ とおくと $A^{k+1}=P_k U_k$ であり、 これは $A^{k+1}$ の QR 分解である。

まず

$$\begin{eqnarray*} P_k U_k&=&(Q_0 Q_1\cdots Q_{k-1}Q_k)(R_k R_{k-1}\cdots R_1 R_... ...{k-1})(Q_kR_k) (R_{k-1}\cdots R_1 R_0)\\ &=&P_{k-1}A_k U_{k-1} \end{eqnarray*}$$

で、主張1より $P_{k-1}A_k=A P_{k-1}$ だから、

$$P_k U_k=A P_{k-1} U_{k-1}.$$

この式を再帰的に用いて、

$$P_k U_k= A P_{k-1} U_{k-1} = A (A P_{k-1} U_{k-2}) \cdots = A^k P_0 U_0 = A^k Q_0 R_0 = A^k A =A^{k+1}.$$

一方、$P_k$ は unitary 行列の積であるから unitary 行列、$U_k$ は 対角線分が正である上三角行列の積であるから対角成分が正である上三角行列で ある。ゆえに $A^{k+1}=P_k U_k$ は $A^{k+1}$ の QR 分解である。