4非线性方程求解4非线性方求解程求解.ppt

华长生制作 第四章 非线性方程数值求解 3)构 造 原 理 定理1.(根的存在定理) ??????这个原理指出了根的存在区间可由两端点处的函数值是否反号确定，将含根区间分为两个长度相等的子区间后，在这两个子区间上也可利用零点原理确定根在那个子区间上，如此继续下去就达到将含根区间逐步缩小的目的，此时，在这一些相互包含的子区间中构造收敛的数列 将它收敛于根 ，见下图?? 　　如精度要求提高,则上式的关键项ln(1/ε), 由于210=1024,如要求误差缩小0. 001,则要求多计算10次。 第四章 非线性方程数值求解 迭代法是求解非线性方程近似根的一种方法，这种方法的关键是确定迭代函数?(x). 简单迭代法:用直接的方法从原方程中隐含的求出x，从而确定迭代函数?(x)，这种迭代法收敛速度较慢，迭代次数多. Newton迭代法采用另一种迭代格式, 具有较快的收敛速度，由牛顿迭代法可以得到很多其他迭代格式。 需要讨论的问题 首先期望每个xk都在? (x)的定义域上,保持有界而且收敛到精确解; 改进条件： 改进条件： 第四章 非线性方程数值求解 3、牛顿迭代法的步骤 步一、准备。选定初始近似值 ，计算 步二、迭代。按公式 迭代一次，得到新的近似值 ，计算 步三、控制。如果 满足 或 .则终止迭代，以 作为所求的根；否则转步四。此处 是允许误差, 5. Newton迭代法的应用----------开方公式 对于给定正数 应用牛顿迭代法解二次方程 可导出求开方值 的计算公式 设 是 的某个近似值，则 自然也是一个近似值，上式表明，它们两者的算术平均值将是更好的近似值。 定理 开方公式对于任意给定的初值 均为平方收敛。 4、弦截法 6. Newton法的改进(III) : 牛顿下山法 一般地说，牛顿法的收敛性依赖于初值 的选取，如果 偏离 较远，则牛顿法可能发散。 为了防止发散，通常对迭代过程再附加一项要求，即保证函数值单调下降： 满足这项要求的算法称为下山法。 牛顿下山法采用以下迭代公式： 其中 称为下山因子。 11 定理4.4：对于迭代过程 ，如果 在所求根 的邻近连续，并且 : (*) , 则该迭代过程在点 邻近是P阶收敛的。 计算流程: 选取初值 对 做 如果 ,跳出循环 否则,置 ,继续循环 输出 将方程变形为 ,这儿 , 在 内 ,所以迭代是收敛的。 例1:用迭代法求方程 的根,精确到0.001。 解:设, , ,在 区间内有根。 ,则 取 ， ,迭代结束。 。 § 4.3 Newton迭代法 §4.3　Newton 迭代法 将f(x)在点xn作Taylor展开: ——Taylor展开线性化 f(x)=0 近似于 f(xn)+ f′(xn)(x-xn)=0 (1) 从(1)解出x, 记为xn+1 ,则 1. Newton迭代公式建立 4 它对应的迭代方程为 显然是f(x)=0的同解方程，故其迭代函数为 在 f(x)=0的根 x* 的某个邻域 内, 在 x* 的邻域R 内，对任意初值 ，应用公式（2）来解方程的方法就称为牛顿迭代法。它是解代数方程和超越方程的有效方法之一. 2. Newton迭代法的几何意义 与x轴(y=0)的交点x，作为下一个迭代点xn+1 ,即 用f(x)在 xn 处的切线 Newton迭代法又称切线法. 15 16 步四、修改。如果迭代次数达到预定指定的次数N，或者 ,则方法失败；否则以 代替 转步二继续迭代。