第2章 方程求根.ppt

本章讨论非线性方程 的求根问题， function [xvect,xdif,fx,nit]=bisect(a,b,toll,nmax,fun) err=toll+1; nit=0; xvect=[]; fx=[]; xdif=[]; while (nit nmax err toll) nit=nit+1; c=(a+b)/2; x=c; fc=feval(fun,x); xvect=[xvect;x]; fx=[fx;fc]; x=a; if (fc*feval(fun,x) 0) a=c; else b=c; end; err=abs(b-a); xdif=[xdif;err]; end return 由迭代公式可得 据此反复递推，得到 于是对任意正整数，有 在上式令 ，注意到 即得到结果。证毕 根据定理2的结论，对于给定的计算精度，迭代次数是可以 预先确定的，但由于公式中含有常数 ，使得计算迭代次数 较为复杂，根据估计式 我们得到： 令 ，得到 由此可知，只要相邻两次计算结果的偏差 足够小即可 保证近似值 有足够的精度。 对于定理中的条件2，在实际使用时，如 果 且对任意的 有 则由中值定理可知 有 它表明定理中条件2可由 替代。 定义1：设 有不动点 ，如果存在 的某个领 域 ，对任意 ，迭代公式 产生的序列 收敛到 ，则称该迭代法局部收 敛。 全局收敛与局部收敛 例1： 利用牛顿迭代法求解 f(x)=ex-1.5-tan-1x 的零点。初始点x0=-7.0  例1： 利用牛顿迭代法求解 f(x)=ex-1.5-tan-1x 的零点。初始点x0=-7.0  p 阶收敛 迭代法 就收敛, 即要求 定理指出, 只要 因此,当 迭代就可以终止, 只要构造的迭代函数满足 此时虽收敛 但不一定是唯一根 对于预先给定的误差? 可以作为方程的近似解 由于 因此 为有根区间 由于 则 用迭代法求方程的近似解,精确到小数点后6位 本题迭代函数有两种构造形式 因此采用迭代函数 例如: 解: 时 §2.2 迭代法—算例分析3 取初值 d1 = 0.1000000 d2 = -0.0105171 d3 = 0.1156e-002 d4 = -0.1265e-003 d5 = 0.1390e-004 d6 = -0.1500e-005 d7 = 0.1000e-006 由于|d7| =0.1000e-0061e-6 因此原方程的解为 x7 = 0.090525 x1 = 0.1000000 x2 = 0.0894829 x3 = 0.0906391 x4 = 0.0905126 x5 = 0.0905265 x6 = 0.0905250 x7 = 0.0905251 定义 设迭代过程 收敛于方程 的根 ，如果迭代误差 ,当 时 成立下列渐进关系式 则称该迭代过程是 阶收敛的. 特别地， 时称线性收敛， 时称超线性收敛， 时称平方收敛。 §2.2 迭代法的收敛速度 局部收敛性 设 为 的不动点， 在 的某个领域连续，且 ，则迭代法 收敛。 定理: 设 为 的不动点， 在 的某个领域连续，且 ，则迭代法 收敛。 证明：由连续函数的性质，存在 的某个领域 ，使对任意 成立 ， 此外，对于任意 ，总有 ， 这是因为 于是可断定迭代过程对于任意的初值收敛. 定理: 定理的条件保证了不动点迭代的全局收敛性。 即迭代的收敛性与初始点的选取无关。 这种在 x* 的邻域内具有的收敛性称为局部收敛性。 定理中的条件 | ?’(x) | ? L 1 可以适当放宽 (2’) ?’(x) 在 x* 的