第六章 非线性方程.pptVIP

6.2.5 斯特芬森加速法 §6.3 牛顿法 6.3.1 牛顿迭代公式的构造 6.3.2 牛顿迭代法的收敛性 6.3.3 牛顿迭代法的变形 (2) 若 是 的不动点， 在 的邻域存在且 连续，且 ，则斯特芬森迭代法是局部收 敛的. (3) 若 ，则斯特芬森迭代法是平方收敛的； 若 ，则斯特芬森迭代法是超平方收敛的. 牛顿迭代法用线性方程代替原方程，建立迭代格式求 近似根. 其最大优点是在单根附近的近似根序列具有较高的 收敛速度，还可以用来求代数方程的重根、复根. 设 是非线性方程 的一个近似根，将 在 处作泰勒展开： 得到近似的线性方程 若 ，求此线性方程的根，记为 ，即 将 在 处泰勒展开，得到近似的线性方程 若 ，记此线性方程的根为 ，则 以此类推，一般地，若 ，则 （6.13） 上式称为方程 的牛顿迭代公式，它也是一种 不动点迭代法，其迭代函数为 其中 牛顿法的几何意义是每次 在近似根附近以切线代替 曲线，并以切线与横轴的 交点近似曲线与横轴的交 点，如图: 因此，牛顿法又称为 切线法. 由定理6.4 ，可以得到牛顿迭代法的局部收敛定理. 定理6.6 设 是方程 的根， 在 的邻 域内连续，则牛顿迭代法局部收敛. 且求单根时， 牛顿迭代至少二阶收敛；求重根时，牛顿迭代只 有一阶收敛. 证明： 已知牛顿迭代函数 两边求导得 当 是 的单根时，有 且一般有 .由此可得 由定理6.4知，此时牛顿法至少是二阶局部收敛，且 当 是 的 重根时，有 且 ，此时 故 于是 因为 ，所以 ，故此时牛顿法是一 阶局部收敛的，且 局部收敛性限制了牛顿法的使用范围，为保证牛顿法 的大范围收敛性，须增加一些条件. 定理6.7 给定 ，设 在 上存在二阶连 续导数，且满足条件： (1) (2) (3) 不变号； (4) 初始值 ，且 则可得 (1) 在 内有惟一根 (2) 由初始值产生的牛顿迭代序列收敛于 (3) 即牛顿迭代二阶收敛. 虽然牛顿法的迭代收敛速度快，但每一步迭代都需 计算 和 的值，计算量大，其次对初值选 择要求比较高，偏离精确值很远的初值产生的迭代 序列往往不收敛到所求的根，为了克服这些缺点， 产生了若干牛顿法的变形. 1. 求重根的牛顿变形法 定理6.6表明，当 是 的 重根时， 牛顿迭代一阶收敛，但收敛缓慢，改善如下： 法一：若 已知，构造 （6.14） 此时迭代函数 ，由定理6.6推导过程知， ，故（6.14）式局部收敛且至少二阶收敛， 不过，实际计算时往往很难知道根的重数，此时直接 用上式难以计算，则可考虑法二. 法二： 令 ，若 是 的 重根，则 是 的 重根，从而 是 的单根. 因此，由定理6.6，对 采用牛顿迭代法所得的 近似根序列至少具有2阶收敛性.由于 则方程 的牛顿迭代公式为 （6.15） 这种方法不要求知道具体重根数，但需要求 的二阶导数. 2. 牛顿下山法 牛顿法的收敛性依赖于初值的选取，若初值偏离所求 根较远，则牛顿法可能发散，为了防止迭代过程发散 ，我们对其附加限制，使迭代过程单调，即 （6.16） 构造牛顿下山公式： （6.17） 其中 称为下山因子. 构造牛顿下山公式的重点是选择合适的下山因子. 先取 ，若计算结果满足下山条件式（6.16）， 就继续下一步迭代；若不满足，则减半取下山因子， 取 ，重复上述过程，直至满足下山条件. 一般而言，只要满足下山条件，从任意初值出发都 能得到收敛的近似根序列. 3. 弦截法(割线法) 当 的形式比较复杂难以计算时，可改用差商 形式 带入(6.13)式，得 （6.18） 利用(6.18)式作为迭代格式的方法称为弦截法. 几何意义： 用经过 及 两点的割线 近似曲线 并以割线与横轴的交点近 似曲线与横轴的交点. 弦截法需要两个初值， 避免了牛顿法的求导运算. 第六章 非线性方程求根 引言 基本概念与二分法 不动点迭代法及其收敛性 牛顿法 引言 在科学研究和工程应用中经常遇到求解非线性方程 的问题，对于非线性方程 （6.1） 其中 . 若 为不高于4次的多项式函数 则这样的方程有解析的求根公式可用. 若 为高于4次的多项式函数或者不为多项式函数， 则往往没有求根公式，只能采用适当的数值方法求得 达到精度要求的近似解. §6.1 基本概念与二分法 几个有关方程根的基本概念 定义6.1 对于单变量方程 ，若存在 ，使得 ，则称 为 的零点，或称 为方程 的根. 若 可表示为 其中 或者 为正整数，则称 是 的 重零点，或称 是方程 的 重根. 定理6.1 设函数 则点 是 的 重零点，当且仅当 6.1.2 二分法 用计算机求解非线性方程分两步进行：一是根的隔离，即找出隔离区间；二是根的精确化，即从隔离区间的近似值出