牛顿（NEWTON）法求根.pptVIP

§2.4 牛顿迭代法 一、 Newton迭代法的基本思想 二、牛顿法的收敛速度 三、牛顿法的几何意义 四、牛顿迭代法的步骤 下山因子λ的选取：搜索试探 先从λ=1开始，试探|f(xk+1)||f(xk)|是否成立，若不成立，则将λ减半进行试算，也即在集合 中依次挑选下山因子，直到找到某个使单调递减条件成立的下山因子。（书上P27—P28，例27） 注：如果在上述过程中挑选的λ已很小了仍无法使递减条件满足，应考虑重新选择初始值x0，然后进行另一轮的迭代。 例：取x0=0.6，对方程x3-x-1=0求使下山成功的下山因子λ。 解： 迭代公式为： λ=1， =-17.3 =17.6 f(x1)f(x0) */38 计算方法二④⑤ */39 计算方法二④ 学习和了解科学计算的桥梁 上节课回顾 由方程f(x)=0变换为等价方程 x=?(x) （#） 先取一个估计值x0（初始值），若?(x0) =x0，则x*=x0(可能性很小)一般?(x0) ≠x0， 记 x1=?(x0) ， 若x1≠?(x1) ， 记 x2=?(x1) ，再用x2继续试探 记 x3=?(x2) 如此反复计算…… 即形成一迭代公式 xk+1=?(xk) ，(k=0,1,2,…) 简单迭代法的思想 x*是方程f(x)=0的根 x*是方程x=?(x)的根 x*为?(x) 的不动点 即：x*= ?(x*) 当{xk}收敛于x*， x*即是?(x)的不动点，也就是原方程的根。 ?(x)是连续函数，称为迭代函数 1、迭代公式 xk+1=?(xk)收敛 {xk}收敛于x* 2、若 则迭代公式xk+1=?(xk) 收敛 全局收敛性 3、若?(x)在x*邻近有连续的一阶导数，且| ?’(x*)|1,则迭代公式xk+1=?(xk)具有局部收敛性。 注：有局部收敛性的条件| ?’(x*)|1可用| ?’(x0)|1来近似代替。 迭代公式收敛性的判断 误差估计式 当迭代公式xk+1=?(xk) 收敛时，有误差估计式： 即：只要前后两次迭代值的差值足够小，就可使近似值xk达到任意的精度。 | ?’(x*)|≤L1 或| ?’(x0)|≤L1 特别地,当L≤1/2时，有不等式|x*-xk|≤|xk-xk-1|,此时，只要|xk-xk-1| ε，就可以终止迭代，求出满足精度要求的近似根xk。 定理：对于迭代过程xk+1= ?(xk) ，如果?(p)(x) 在所求根x*的邻近连续，并且?’(x*)= ?’’(x*) =...= ?(p-1)(x*) =0，(#) ?(p)(x*)≠0，则该迭代过程在点x* 邻近是P阶收敛的。 该定理告诉我们，迭代过程的收敛速度依赖于迭代函数. 如果选取当x∈[a,b] 时?’(x)≠0，则该迭代过程只能是线性收敛。 因此,简单迭代法的收敛速度一般较慢 . 有的公式未必收敛. 局部收敛时阶的确定 2 2 设 xk 是f(x)=0的一个近似根,把f(x)在xk处作泰勒展开 若取前两项来近似代替f(x)(称为f(x)的线性化)，则得近似的线性方程 设 f ’(xk)≠0 ,令其解为 xk+1 ,得 称其为f(x)=0的牛顿迭代公式。 解出x 一. Newton迭代法的基本思想 4 应用公式（1）来解方程的方法就称为牛顿迭代法。它是解代数方程和超越方程的有效方法之一。 它对应的迭代方程为 显然是f(x)=0的同解方程 故其迭代函数为 迭代公式为: (1) 17 解 因 f ’(x)=3x2+4x+10 ,所以迭代公式为 1.368808189 4 1.368808189 1.369336471 1.411764706 xk 3 2 1 k 例1:用牛顿法求下面方程的根f(x)=x3+2x2+10x-20 从计算结果可以看出,牛顿法的收敛速度是很快的,进行了四次迭代就得到了较满意的结果. 选取x0=1,计算结果列于下表 18 由牛顿迭代公式 例2 计算 的近似值。 ?=10-6