数值计算方法--第18-1讲--牛顿法及其改进解析.pptVIP

牛顿法的计算步骤： 插 值 法 主讲教师：刘春凤 8/eol/jpk/course/welcome.jsp?courseId=1220 非线性方程（组）数值解法 主讲教师：杨爱民 第 七章 方程求根与二分法 牛顿法及改进的牛顿法 迭代法及其收敛性 牛顿下山法与弦截法 非线性方程组的牛顿法 1 3 4 2 5 牛顿迭代法原理 牛顿迭代法的几何意义 简化的牛顿法与牛顿重根法 取x0作为初始近似值，将 f(x)在x0做Taylor展开: 重复上述过程 ? 作为第一次近似值 Newton 迭代公式 基本思路：将非线性方程 f (x)=0 线性化。 牛顿迭代法（切线法）原理 牛顿迭代法原理验证 简单验证： 牛顿迭代法几何解释 牛顿法的几何解释 x y x* x0 x 1 x 2 故：牛顿法也称为切线法 牛顿迭代法（切线法）例题 解： 则由牛顿迭代公式可得： 例1. 牛顿迭代法流程图 (1) 选定初值x0 ,计算f (x0) , f ?(x0) (2) 按公式 迭代 得新的近似值xk+1 (3) 对于给定的允许精度?，如果 则终止迭代，取 ;否则k=k+1，再转 步骤(2)计算 允许精度 最大迭代次数 迭代信息 牛顿迭代法典型例题 解 例2. 牛顿迭代法典型例题 程序 Clear[f,x,xx] f[x_]:=x*Exp[x]-1 g[x_]:=x-f[x]/f[x] x=0.5; xx=0.51; Do[x=g[x];Print[k, ,x],{k,1,8}] Abs[g[xx]]-1//N Clear[x] FindRoot[f[x]==0,{x,3}] 迭代解 1. 0.57102 2. 0.567156 3. 0.567143 4. 0.567143 5. 0.567143 6. 0.567143 7. 0.567143 8. 0.567143 精确解 {x - 0.567143} 牛顿迭代法如何收敛？ 结论：Newton法的收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 ? 疑问：Newton法是否永远收敛呢？如果不是，那么与谁有直接关系呢？ 牛顿法的特点 优点: 收敛快! 缺点: 牛顿法的特点 牛顿迭代法不足 牛顿法的不足： 1）牛顿迭代法存在从一个根跳到另一个根的情况。 2）如果f(x)=0没有实根，则牛顿迭代序列不收敛。 3）每次迭代都要求导，程序易中断。 简化的牛顿迭代法 x y x* x0 x 1 x 2 斜率恒定不变 思考题 如何改进牛顿法？ 简化的牛顿法 修改为 简化的牛顿迭代法 牛顿重根法 牛顿重根法 且若三种方法都收敛，则一般重根法收敛最快， 牛顿法次之，简化迭代法最慢。 牛顿迭代法典型例题 例3 解 三种迭代公式分别为 牛顿迭代法典型例题 三种迭代公式迭代结果分别为: 牛顿法迭代解 1. 0.9 2. 0.93235 3. 0.95446 4. 0.97773 5. 1.00000 6. 1.00000 7. 1.00000 8. 1.00000 精确解 {x - 1.00000} 重根法迭代解 1. 0.9 2. 0.99706 3. 0.99997 4. 1.00000 5. 1.00000 6. 1.00000 7. 1.00000 8. 1.00000 精确解 {x - 1.00000} 简化法迭代解 1. 0.9 2. 0.92583 3. 0.94806 4. 0.96752 5. 0.98001 6. 0.99952 7. 1.00000 8. 1.00000 精确解 {x - 1.00000} 思考题 牛顿迭代法及其改进都无法回避计算机中求导的问题，那有没有其它算法可以回避呢？