next up previous contents
Next: 4.5.3 QR 法の歴史 Up: 4.5.2.0.1 証明 Previous: 4.5.2.0.1.3 ステップ3

4.5.2.0.1.4 ステップ4
(4.10) と (4.15) から

\begin{displaymath}
A_{k+1}=P_k^\ast A P_k=
D_1^k D_2 S_k^\ast (Q^\ast A Q)S_k D_2^\ast D_1^{k\ast}.
\end{displaymath}


\begin{displaymath}
A=M \Lambda M^{-1}=(Q R)\Lambda (R^{-1}Q^\ast)
\end{displaymath}

より

\begin{displaymath}
Q^\ast A Q=R\Lambda R^{-1}.
\end{displaymath}

この行列は (右辺に注目することで) 上三角行列で、かつ対角線上に $\lambda_1$, $\cdots$, $\lambda_n$ がこの順にならんでいる。 $k\to\infty$ のとき、$S_k\to I$ だから

\begin{displaymath}
A_{k+1}\to D_1^k (D_2 R\Lambda R^{-1}D_2^\ast) D_1^{\ast k}
\end{displaymath}

この右辺は上三角行列で、対角成分は $\lambda_1$, $\lambda_2$, $\cdots$, $\lambda_n$ である。

この定理から $k\to\infty$ とするときの収束の速さは

\begin{displaymath}
\left(\frac{\lambda_i}{\lambda_j}\right)^k
\end{displaymath}

という因子に関わることが分かる。


next up previous contents
Next: 4.5.3 QR 法の歴史 Up: 4.5.2.0.1 証明 Previous: 4.5.2.0.1.3 ステップ3
Masashi Katsurada
平成17年6月2日