非线性方程求根

非线性方程求根

隔根区间

即位只包含一个根的区间

描图法

画出 $y=f(x)$ 的函数图像大致判断根所在的区间

可将函数变形，更好绘图

比如将 $f(x)=x^3-x-1$ 变形为 $x^3=x+1$ 作图，交点所在位置即为根所在区间

逐步搜索法

1. 确定实根所在区间 $[a,b]$
2. 确定步长 $h=(b-a)/n$
3. 用零点定理判断每个小区间是否存在根

• 对于多项式存在以下结论

例题

二分法

通过计算隔根区间的中点，逐步将隔根区间缩小，从而可得方程的近似根数列 { x n } \{x_n\} {xn​}

步骤

令 $f(a)<0,f(b)>0$

1. 若 $f(\frac{a+b}{2})=0$ ，则 $x^{\ast }=(a+b)/2$
2. 若 $f(\frac{a+b}{2})<0$ ，隔根区间变为 $[(a+b)/2,b]$
3. 若 $f(\frac{a+b}{2})>0$ ，则 $[a,(a+b)/2]$\

讨论

• $x\ast \approx (a_k+b_k)/2,$ $k$ 足够大
• $b_k-a_K=(b-a)/2^k$
• 误差为 $|x\ast -x_k|\le (b_k-a_k)/2=(b-a)/2^{k+1}$
• 若事先给定精度要求 $eps$ ，只需要 $(b-a)/2^{k+1}<eps$
• 不能求复根及偶数重根

例题

迭代法

简单迭代法

将 $f(x)=0$ 转化为 $x=g(x)$

假设 $g(x)$ 为连续函数

任取一个初值 $x_0$ ，代入 $g(x)$

$x_1=g(x_0)$

$x_2=g(x_1)$

$\cdots \cdots$

$x_{k+1}=g(x_k)$

该过程被称为简单迭代法

收敛性

$g(x)$ 迭代函数，$x_0$ 为初始近似值，$x_k$ 为第 $k$ 步迭代值

如果存在一点 $x^{\ast }$ 使得 { x k } 0 ∞ \{x_k\}^{\infty}_{0} {xk​}0∞​

$$\underset{k\to +\infty }{\lim }{x}_k=x^{\ast }$$

则称迭代法收敛，否则称为发散

例题

迭代法收敛的充分条件

将 $x=g(x)$ 建图模拟

定理1

设迭代函数 $g(x)$ 在 $[a,b]$ 上连续且满足

1. 当 $x\in [a,b]$ 时， $a\le g(x)\le b$
2. 存在一个正数 $L$ , 满足 $0<L<1$ , 且 $\forall x\in [a,b]$ ，有 $|g^{，}(x)|\le L$

则

1. 方程 $x=g(x)$ 在 $[a,b]$ 内有唯一解 $x^* $
2. 对于任意 $x_0\in [a,b]$，$x=g(x)$ 均收敛于 $x^{\ast }$
3. $|x_k-x^{\ast }|\le \frac{1}{1-L}|x_{k+1}-x_k|$
4. $|x_k-x^{\ast }|\le \frac{L^K}{1-L}|x_1-x_0|$

例题

牛顿迭代法

$f(x)$ 在区间 $[a,b] $ 可导，且对于 $x\in [a,b]$， 有 $f^{\prime }(x)\ne 0$

牛顿迭代法：$x_{k+1}=x_k-\frac{f(x_k)}{f^{\prime }(x_k)}$

其实就是泰勒展开的前两项

收敛性

用定理一可证明

用定理四可证明

弦割法

将牛顿迭代式中

$f^{\prime }(x_k)\approx \frac{x_k-x_{k-1}}{f(x_k)-f(x_{k-1})}$