数值分析—— 方程迭代求根

最新推荐文章于 2024-07-19 13:59:11 发布

小刀来啦

最新推荐文章于 2024-07-19 13:59:11 发布

阅读量746

点赞数 1

分类专栏：数学原理数值分析笔记文章标签：算法

本文链接：https://blog.csdn.net/weixin_40519529/article/details/124573258

版权

6 篇文章 2 订阅

订阅专栏

3 篇文章 0 订阅

订阅专栏

前言

有时候我们需要知道非线性方程在某个区间内是否存在根，存在的话有几个等问题，直接求解往往是困难，可以通过迭代的思想去逼近根，比如经典的二分法和牛顿迭代法。

具体内容详见 《数值分析》李庆扬、王能超、易大义（华中科大出版社）第六章

设有待求根方程 $f (x)$ ，将其写成隐式方程格式： $x=\phi(x)$ ，给出根的某个猜想值 $x_0$ 代入得到 $x_1=\phi(x_0)$ ，重复这个过程， $x_{k+1}=\phi(x_k),(k=0,1,...,n)$ ，若 ${ x_k \}$ 存在极限 $x^*=\lim_{k \rightarrow \infty} x_k$ ，则称此迭代过程收敛， $x^*$ 即为方程的根。
为了避免迭代过程发散，对迭代方程 $x=\phi(x)$ 作出以下约束：若存在定数 $0 < L < 1$ ，使对于任意 $\in [a,b]$ 均有 $|\phi(x)^{'}| \leq L$ 则该迭代过程收敛，且存在唯一根。

证明存在零点：

由假设可知，对于 $\in [a,b]$ ，总有 $\phi(x) \in [a,b]$ ，记 $f(x)=\phi(x)-x$ ，有 $f(a)=\phi(a)-a \geq0, f(b)=\phi(b)-b \leq 0$ ，即 $f (x)$ 在 $[a, b]$ 上有根。
证明零点唯一：

设存在两个零点， $x^*=\phi(x^*), \bar{x}=\phi(\bar{x})$ ， $x^*-\bar{x}=\phi(x^*)-\phi(\bar{x})=\phi^{'}(\xi)(x^*-\bar{x})$ ，由于 $\phi^{'}(\xi) \leq 1$ ，则 $(x^*-\bar{x})$ 只能为0，即只有一个零点。
证明收敛：

$|x^*-x_k|=|\phi(x^*)-\phi(x_k)|=|\phi^{'}(\xi)||x^*-x_k| \leq L|x^*-x_k| \leq ... \leq L^k |x^*-x_0|$ ，又 $\lim_{k \rightarrow \infty} L^k=0$ ，当k足够大时，迭代过程收敛。

给出误差估计式：
$|x_k-x^*| \leq \frac{L^k}{1-L}|x_1-x_0|\\|x_k-x^*| \leq \frac{1}{1-L}|x_{k+1}-x_k|$

即若两次迭代的结果相差足够小，就可以保证结果的精度。

记 $x_0$ 是迭代初值，用迭代公式迭代一次得到： $x_1=\phi(x_0)$ ，此时有微分中值定理： $x_1-x^*=\phi(x_0)-\phi(x^*)=\phi^{'}(\xi)(x_0-x^*)$ 若 $\phi^{'}(\xi)$ 是一个小量，用L去近似代替，则有 $x_1-x^* \approx L(x_0-x^*)$ ，得 $x^*=\frac{1}{1-L}x_1-\frac{L}{1-L}x_0$ 即 $x_2=\frac{1}{1-L}x_1-\frac{L}{1-L}x_0=x_1+\frac{L}{1-L}(x_1-x_0)$ 是比 $x_1$ 更好的迭代值。

上面的改进格式中存在近似值L不好估计，若再通过迭代公式迭代一次即可联系消去L，设 $x_1=\phi(x_0),x_2=\phi(x_1)$ ，有 $x_1-x^* \approx L(x_0-x^*)\\ x_2-x^* \approx L(x_1-x^*)$
联立可得， $\frac{x_1-x^*}{x_2-x^*} \approx \frac{x_0-x^*}{x_1-x^*}$
解出 $x^*$ ， $x^* \approx \frac{x_0x_2-x_1^2}{x_0-2x_1+x_2}=x_2-\frac{(x_2-x_1)^2}{x_0+2x_1+x_2}$ ，同理， $x_3=x_2-\frac{(x_2-x_1)^2}{x_0+2x_1+x_2}$ 是更好的迭代值。