next up previous
: 重根近くでの収束 : ニュートン・ラフソン法(その2) : ニュートン・ラフソン法(その2)

2次収束

ニュートン法は根 ($x_s; f(x_s)=0$) に十分近く、また根での微係数が有限の時 ( $f'(x_s) \neq 0$)には極めて収束が早い。実際根の近くで $x_s$ から微小距離 $\epsilon$ 離れた点で Taylor 展開すると、関数と導関数は近似 的に

$\displaystyle f(x_s+\epsilon)$ $\textstyle =$ $\displaystyle f(x_s)+\epsilon f'(x_s) +\frac{\epsilon^2}{2}f''(x_s) + \cdots$ (1)
$\displaystyle f'(x_s+\epsilon)$ $\textstyle =$ $\displaystyle f'(x_s)+\epsilon f''(x_s) + \cdots$ (2)

で表される。これらを $\epsilon_{i} \equiv x_{i} - x_{s}$ として Newton-Raphson の公式
\begin{displaymath}
x_{i+1}=x_i - \frac{f(x_i)}{f'(x_{i})}
\end{displaymath} (3)

に代入し、$f(x_s)=0$ を使うと、
\begin{displaymath}
\epsilon_{i+1} \approx - \epsilon_i^2 \frac{f''(x_s)}{2f'(x_s)}
\end{displaymath} (4)

となるので、2次収束で極めて早く収束する。



Kiyohide Nomura 平成15年6月25日