2012数值分析-lec10非线性方程迭代解法.pptVIP

Newton 迭代法 Newton 法的几何解释 割线法的几何表示 抛物线法 抛物线法计算公式 Y x0 X x* x1 x2 x3 Y P0 P2 P1 割线法的收敛性 引理：设f(s)=0，在 s 的某领域 [s- δ, s+δ] 内 f″(x) 连续，f′(x)≠0，又设 xk-1，xk∈[s-δ, s+δ] 且 xk-1, xk, s 互异，记 ek=s-xk, 则有 其中 xk＋1 是由割线法产生，ηk, ξk 在 min(xk-1, xk, s) 与 max(xk-1, xk, s) 之间。 定理：设 f(s)=0,在 s 的某领域 Iδ= [s- δ, s+δ] 内 f″(x) 连续，f′(x)≠0，则存在ε0，当 x-1, x0∈Iδ时，则由割线法产生的序列 {xk} 收敛于 s ，且收敛速度的阶至少为 1.618。 单点割线法 定理4.9：单点割线法有一阶的收敛速度 Newton 迭代法的变体：单点割线法（一阶） Newton迭代法需要计算f(x)的一阶导数，对复杂的函数，特别是多元隐函数，求函数值是计算量比较大或者比较繁琐，尽量减少函数值的计算！ 在Newton迭代法中用 来近似f(x)在xk处的一阶导数 由此得到的算法叫单点割线法。 定理4.9：单点割线法有一阶的收敛速度 牛顿下山法： 目的是解决初值的选取范围太小这一困难。 构造迭代格式为： 其中的参数满足： 这个方法称为牛顿下山法。其中的参数称为下山因子： 通常取 ，然后逐步减半。 牛顿下山法当 时，只有线性收敛速度，但对初值的选取却放的相当宽。 作业 课后习题1、2、4、5、7、9 * 数值分析 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang2012 答疑时间：星期三下午2:00－5:00 答疑地点：主216 第十讲非线性方程的迭代解法 第三章非线性方程与非线性方程组的迭代解法 这一部分的主要任务是解 其中f（x）是一个一元非线性函数。 非线性科学是当今科学发展的一个重要研究方向，很多实际工程物理问题都归结为非线性方程（组）的求解。 例如： 无穷组解 无解 一个解 两个解 四个解 非线性方程的求根非常复杂。 求根问题包括下面三个问题： 根的存在性：即f(x)=0有没有根？若有， 有几个根？ 哪儿有根？确定有根区间 根的近似求解 常用的求非线性方程根的方法 二分法（对分法、搜索法） 不动点法 （简单迭代法、压缩映象法）及其加速算法 Newton方法及其变体 根，重根的定义： 根：如果存在常数s，使得 f(s)=0则称s是f(x)=0的根（零点）； 重根：如果 ， 称s为m重根。特别地对 f(x) 是多项式，则有 其中 。 二分法（对分法） 理论依据：设函数 f(x)∈C[a,b], f(a)f(b)0，则 f(x)在区间[a,b]上有一实根s使 f(s)=0。 二分法工作原理：每次把半分根所在的区间[ak, bk]分成两个等长的子区间[ak, 0.5(ak+ bk) ] 和 [0.5(ak+ bk) ，bk]，并剔除根所不在的区间，直至根所在的区间变得充分小！即|bk-ak|ε。 计算框图为： 对分法（二分法）：利用连续函数的性质进行对分 算法： 1.令a0=a, b0=b, 给定绝对误差限ε和最大迭代次数N，计算f(a)和f(b)的值。如果f(a)f(b)0，运算停止，输出计算失败；如果|f(a)| ε，计算结束，输出a即是所求解；如果|f(b)| ε，计算结束，输出b即是所求解;如果都不是转入下一步； 2.开始循环迭代，k=0,1,…; 计算x0=0.5(a0+b0)和f(x0)。如果 |f(x0)| ε，计算结束，输出x0即是所求解 ；否则如果f(a0)f(x0)0， 令b0=x0； 否则，令a0=x0 ; 3.令k=k+1，如果b0-a0ε 而且kN, 则转入步骤2; 4.如果kN, 输出计算失败，停止计算。 简单迭代法及其收敛速度 迭代法的构想 从一个初值x0出发，计算 如果 {xk}收敛，即存在x*, 使得 则由 得 即 x*是 φ(x)的不动点，也就是 f(x) 的根。 问题： {xk} 收敛吗？怎样实现φ(x) 简单迭代法的收敛速度 Steffensen加速收敛方法 非线性问题的最简单解