数学分析(六)-微分中值定理及其应用7：方程的近似解

§ 7 方程的近似解
在实际应用中, 常常求方程
$$f(x)=0$$
的解. 而方程求解的方法主要有两种: 解析法和数值法. 解析法也称为公式法,
得到的解是精确的,
比如一元二次方程的求解公式.然而并不是所有的方程的根都能通过这种方法而求得.法国著名数学家伽罗瓦
(Galois) 在 19 世纪就证明了形如
$$y=a_n{x}^n+a_{n-1}{x}^{n-1}+\cdots +a_0=0\left(a_n\ne 0\right)$$
的代数方程, 当 $n⩾5$ 时,一般不存在求解公式. 因此对于一般方程,
我们必须寻求其他的求解方法.
下面我们介绍一种数值解法—牛顿切线法.
设 $f$ 为 $[a,b]$ 上的二阶可导函数, 满足
$$f^{\prime }(x)\cdot f^{\prime \prime }(x)\ne 0,f(a)\cdot f(b)<0.$$
牛顿切线法的基本思想是构造一收敛点列 { x n } \left\{x_{n}\right\} { xn​}, 使其极限
$\underset{n\to \infty }{\lim }{x}_n=\xi$ 恰好是方程 (1) 的解.
因此当 $n$ 充分大时, $x_n$ 可作为 $\xi$
的近似值.下面分四种情形进行讨论.
(1) 设 $f^{\prime }(x)<0,f^{\prime \prime }(x)>0$. 从而有
$f(a)>0,f(b)<0$, 并设 $f(\xi )=0$. 令
$$x_0=a,x_n=x_{n-1}-\frac{f\left(x_{n-1}\right)}{f^{\prime }\left(x_{n-1}\right)},n=1,2,\cdots .$$
因为 $f^{\prime \prime }(x)>0$, 所以 $f$ 为 $[a,b]$
上的严格凸函数,由定理 6.14
$$f(x)>f(a)+f^{\prime }(a)(x-a),x\in (a,b].$$
设 $x_0=a$, 则 $y=f(x)$ 在点 $a$ 的切线与 $x$ 轴的交点为
$$x_1=a-\frac{f(a)}{f^{\prime }(a)}=x_0-\frac{f\left(x_0\right)}{f^{\prime }\left(x_0\right)}.$$
由 (3) 式可知 $f\left(x_1\right)>0$ (见图 6-17(1)).
以 $\left[x_1,b\right]$ 代替 $[a,b]$ 重复上述步骤可得 $y=f(x)$ 在点
$x_1$ 的切线与 $x$ 轴交点为
$$x_2=x_1-\frac{f\left(x_1\right)}{f^{\prime }\left(x_1\right)},$$
其中 $f\left(x_2\right)>0,a=x_0<x_1<x_2<\xi <b$.
如此继续上述过程可得如 (2) 式确定的点列 { x n } \left\{x_{n}\right\} { xn​}. 显然
{ x n } \left\{x_{n}\right\} { xn​} 严格递增且有上界, 故可设
$\underset{n\to \infty }{\lim }{x}_n=c$. 由于 $f$ 和 $f$