next up previous contents
Next: 4.3 三重対角化の手法 Up: 4.2 解法についての概観 Previous: 4.2.5.0.1 計算上の注意:

4.2.6 シフト法

既に得られている近似固有値の精度を改良したい、あるいは、ある指定した 値に最も近い固有値を求める方法である。

行列 $A$ の固有値 $\lambda_i$ に対する近似固有値 $\lambda_i'$ が分かっ ているとしよう。この時、

\begin{displaymath}
A' = A - \lambda_i' I
\end{displaymath}

の固有値は

\begin{displaymath}
\lambda_1-\lambda_i', \lambda_2-\lambda_i', \cdots, \lambda_n-\lambda_i'
\end{displaymath}

となる。$\lambda_i'$$\lambda_i$ の近似値ということで、絶対値が最小 なのは $\lambda_i-\lambda_i'$ であると期待できる。よって、$A'$ に対し て逆反復法を適用すれば、この値(それを $\Delta\lambda$ とおこう)が高精 度に計算できる。こうして

\begin{displaymath}
\lambda_i := \lambda_i' + \Delta\lambda
\end{displaymath}

により $\lambda_i$ が求まる。


\begin{displaymath}
A'=
\left(
\begin{array}{ccccc}
4& 1& & & \bigzerou \\
...
...\\
& & 1& 4& 1 \\
\bigzerol& & & 1& 4
\end{array}\right).
\end{displaymath}


next up previous contents
Next: 4.3 三重対角化の手法 Up: 4.2 解法についての概観 Previous: 4.2.5.0.1 計算上の注意:
Masashi Katsurada
平成17年6月2日