数值分析与2.2-2__方程求根(牛顿法和弦截法) .pptVIP

* 《科学计算方法》 1.Newton迭代格式 2.Newton迭代法的收敛性 3.弦截法 ? ? ? 初值: x1=1.5 迭代格式: xn+1=0.5(xn+2/xn) (n = 1,2,·····) 平方根算法求 xn Error 1.416666666666667 2.45e-003 1.414215686274510 2.12e-006 1.414213562374690 1.59e-012 1.414213562373095 2.22e-016 1.414213562373095 2.22e-016 设 x*是方程 f(x)=0 的根, x0是x*的近似值.在 x0 附近,对函数做局部线性化 x1比x0更接近于x* x0 x1 x* 令f(x) = 0 ? ? (n = 0, 1, 2, ·····) 牛顿迭代格式 给定初值 x0 , 迭代产生数列 x0, x1, x2,·········, xn, ······· 应用——求正数平方根算法 设C 0, x2 – C = 0 令 f(x) = x2 – C , 则 由此可知, 平方根迭代具有 2 阶收敛速度 牛顿迭代法误差分析： Newton迭代法的局部收敛性 定理3 设 f(x) 在点x*的某邻域内具有二阶连续导数,且设 f(x*)=0, f ’(x*) ≠ 0, 则对充分靠近点x*的初值x0, Newton迭代法至少平方收敛. 所以, Newton迭代法至少平方收敛（利用收敛定理） ? 例4.求 f(x)=xex – 1= 0 在 x0=0.5 附近的根 解: 迭代格式为 (n = 0, 1, ·····) f=inline(x*exp(x)-1); f1=inline((x+1)*exp(x)); x0=1.5;er=1;k=0; while er0.00001 x=x0-f(x0)/f1(x0); er=abs(x-x0) x0=x;k=k+1 end x = 0.5671 k=6 er=4.3596e-009 牛顿法的缺陷： 1.被零除错误 2.程序死循环 y = arctan x 例如方程: f(x)=x3 –3x + 2 = 0 在重根x*=1附近,f’(x)近似为零 对 f(x) = arctan x 存在 x0,Newton迭代法陷入死循环 Newton迭代法陷入死循环的另一个例子(p.72例5) 取 x0=0, (n = 0, 1, ·····) f’0, f”0 f ’0, f ”0 f’0, f”0 f’0, f”0 牛顿迭代法收敛的四种情况 定理4: 若函数f(x) 在[a，b] 上满足条件 则方程 f(x) = 0 在[a，b] 上有唯一根 x*，且由初值x0按牛顿迭代公式求得的序列{ xn} 二阶收敛于x*。 （1）f(a) f(b) 0； （2）f ’(x)，f ”(x) 在[a，b]上连续且不变号（恒为正或恒为负）； （3）取x0∈[a，b] 使得 f(x0)f ”(x0) 0。 *