《数值计算方法》第二版习题2第15题

《数值计算方法》第二版习题2第$15$题

题目：

设函数$f(x)$在$[a,b]$上至少三阶连续可微，$p\in (a,b)$为$f(x)$的一个$m$重零点，求一个$\lambda$值使改进的$Newton$下山法
$$x_{k+1}=x_k-\lambda \frac{f(x_k)}{f^{\prime }(x_k)}$$
至少是二阶收敛的。

求解：

首先，对于一个迭代过程$x_{k+1}=\phi (x_k)$，如果${\phi }^{(p)}(x)$在所求根$x^{\ast }$的附近连续，并且
$${\phi }^{\prime }(x^{\ast })={\phi }^{\prime \prime }(x^{\ast })=\cdots ={\phi }^{(p-1)}(x^{\ast })=0,{\phi }^{(p)}(x^{\ast })\ne 0$$
则该迭代过程在点$x^{\ast }$附近是$p$阶收敛的。

此定理的证明可由$\phi (x)$在根$x^{\ast }$附近做$Taylor$展开来证明，此处不详述。

因此，要证明一个迭代过程至少是二阶收敛的，我们只需证明$\phi (x)$在其根$x^{\ast }$处的一阶导数为$0$即可。

对于此题，我们首先将迭代过程标准化，即表示为$x_{k+1}=\phi (x)$的形式。因此我们得到:
$$\phi (x)=x-\lambda \frac{f(x)}{f^{\prime }(x)}\text{\hspace{1em}(1)}$$
而$p$为$f(x)=0$在$(a,b)$的一个根，因此若要证明改进的$Newton$下山法至少是二阶收敛的，只需要证明${\phi }^{\prime }(p)=0$即可。

由$p$是$f(x)$的$m$重零点，因此我们可设$f(x)=(x-p{)}^{m}g(x)$，其中$g(x)$满足$g(p)\ne 0$。

对$\phi (x)$求一阶导有：
$$\begin{array}{rl}{\phi }^{\prime }(x)& =1-\lambda (\frac{[f^{\prime }(x){]}^2-f(x)f^{\prime \prime }(x)}{[f^{\prime }(x){]}^2})\\ & =1-\lambda +\lambda \frac{f(x)f^{\prime \prime }(x)}{[f^{\prime }(x){]}^2}\end{array}\text{\hspace{1em}(2)}$$
由于$f^{(j)}(p)=0$，其中$j\in [0,m-1]$。因此，${\phi }^{\prime }(p)$中最后一项分子分母均为$0$，故必须表示为如下极限形式：
$${\phi }^{\prime }(p)=1-\lambda +\lambda \underset{x\to p}{\lim }\frac{f(x)f^{\prime \prime }(x)}{[f^{\prime }(x){]}^2}\text{\hspace{1em}(3)}$$
对$f(x)$分别求一阶和二阶导有：
$$\begin{gathered} \begin{array}{rl}{f}^{\prime }(x)& =m(x-p{)}^{m-1}g(x)+(x-p{)}^m{g}^{\prime }(x)\\ & =(x-p{)}^{m-1}[mg(x)+(x-p)g^{\prime }(x)]\end{array} \\ \begin{array}{rl}{f}^{\prime \prime }(x)& =m[(m-1)(x-p{)}^{m-2}g(x)+(x-p{)}^{m-1}{g}^{\prime }(x)]+m(x-p{)}^{m-1}{g}^{\prime }(x)+(x-p{)}^m{g}^{\prime \prime }(x)\\ & =m(m-1)(x-p{)}^{m-2}g(x)+2m(x-p{)}^{m-1}{g}^{\prime }(x)+(x-p{)}^m{g}^{\prime \prime }(x)\\ & =(x-p{)}^{m-2}[m(m-1)g(x)+2m(x-p)g^{\prime }(x)+(x-p{)}^2{g}^{\prime \prime }(x)]\end{array}\text{\hspace{1em}(4)} \end{gathered}$$
将式(4)代入式(3)，并将分子分母约分和代入，有：
$$\begin{array}{rl}{\phi }^{\prime }(p)& =1-\lambda +\lambda \underset{x\to p}{\lim }\frac{g(x)[m(m-1)g(x)+2m(x-p)g^{\prime }(x)+(x-p{)}^2{g}^{\prime \prime }(x)]}{[mg(x)+(x-p)g^{\prime }(x){]}^2}\\ & =1-\lambda +\lambda \underset{x\to p}{\lim }\frac{m(m-1)g^2(x)}{m^2{g}^2(x)}\\ & =1-\lambda +\lambda \frac{m-1}{m}\\ & =1-\frac{\lambda }{m}\end{array}\text{\hspace{1em}(5)}$$
我们令${\phi }^{\prime }(p)=0$，即可得到$\lambda =m$。

因此，当取$\lambda =m$时，改进的$Newton$下山法至少是二阶收敛的。