新第六讲方程求根.ppt

简介（Introduction） 我们知道在实际应用中有许多非线性方程的例子，例如 （1）在光的衍射理论（the theory of diffraction of light)中,我们需要求x-tanx=0的根 （2）在行星轨道（ planetary orbits）的计算中,对任意的a和b，我们需要求x-asinx=b的根 （3） 在数学中，需要求n次多项式xn + a1 xn-1+...+an-1 x + an ＝0的根 举例 1.用一般迭代法求x3-x-1=0的正实根x* 例2 用一般迭代法求方程x-lnx＝2在区间（2，?）内的根，要求|xk-xk-1|/|xk|=10-8 将方程化为等价方程：x＝2＋lnx k xi 0 3.000000000 1 3.098612289 2 3.130954362 3 3.141337866 4 3.144648781 5 3.145702209 6 3.146037143 另一种迭代格式： 由此可见，对同一个非线性方程的迭代格式，在收敛的情形下，有的收敛快，有的收敛慢。 例题1 用Newton法求方程 的根,要求 例题2 求函数 的正实根 精度要求： 初值x0＝8.0 时,计算的是单根, The iterative number is 28,The numerical solution is 7.600001481 初值x0＝1.0 ,计算的是重根, The iterative number is 1356,The numerical solution is 1.198631981 小结 (1) 当f (x)充分光滑且 x* 是f (x) =0的单根时，牛顿法在x*的附近至少是平方收敛的。 (2) 当f (x)充分光滑且 x* 是f (x) =0的重根时，牛顿法在x*的附近是线性收敛的。 (3) Newton法在区间[a,b]上的收敛性依赖于初值x0 的选取。 结论：Newton法的收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 ? 关于初值的问题： 一般来说采用试探法，但对于一些实际问题初值的选择并不困难，它是明确的。 关于重根的问题： 重根时的改进： 关于导数计算的处理： §6.4 割线法 迭代公式： 例6.9(P154) 31 20 36 5 5 26 迭代次数 离散2 离散1 离散 改进2 改进1 牛顿 方法 说明：epsl=1e-8,x0=1.5,m= [n,x]=newtonfa(ff14,1.5i,0.0001); n = 5 x = -0.0000 + 1.0000i 说明牛顿法也可求复数根 迭代格式一： 迭代格式二： 取初值x0＝0.0,计算如下： 对迭代格式一: the iterative number is 27, the numerical solution is 0.442852706 对迭代格式二: the iterative number is 3, the numerical solution is 0.442854401 从图形中我们可以看出： 在x=7和x=8 之间有一单根； 在x=1和x=2 之间有一重根。 取初值x0＝8.0,用牛顿迭代公式计算如下： 取初值x0＝1.0,用牛顿迭代公式计算如下： 一般 3、迭代法 第六章：方程求根 第六章：方程求根 第六章：方程求根 第六章：方程求根 第六章：方程求根 第六章：方程求根 第六章：方程求根 4、迭代法的收敛性 第六章：方程求根 第六章：方程求根 第六章：方程求根 第六章：方程求根 即 第六章：方程求根 由 第六章：方程求根 第六章：方程求根 容易得到：g′(x)在包含x*的某邻域U(x*) 内 连续，且|g′(x*)|1 解：令f(x)=x-lnx-2 f(2)0,f(4)0,故方程在（2,4）内至少有一个根 因此，? x0?（2，?），xk+1＝2＋lnxk产生的序列? xk ?收敛于X* 取初值x0＝3.0,计算结果如下： 7 3.146143611 8 3.146177452 9 3.146188209 10 3.146191628 11 3.146192714 12 3.146193060 13 3.146193169 14 3.146193204 0 3.000000000 1 3.147918433 2 3.146193441 3 3.146193221 §6.3 Newton迭代法 做迭代函数 Newton迭代公式： 第六章：方程求根 取x0