《第七章前半部分改完》精选课件.ppt

第七章 非线性方程求根 §1 二分法（对分法） §2 迭代法的算法和理论 §3 迭代加速收敛的方法 §4 Newton迭代法 §5 弦割法和抛物线法 * * 一. 研究数值解法的必要性 1 一般方程的根无法用解析表达式给出； 2 三次、四次方程的求根公式较繁。 需要给出求根的近似值的方法。 二. 方程的根 方程f(x)=0 的解 称为方程f(x)=0的根或称为f(x) 的零点。 若f(x) ＝ g(x)，其中m为正整数，g(x)满足 ，显然 为f(x)的零点。这时，称 为f(x)的m重零点，或称 为f(x)=0的m重根。 定理 若f(x)具有m阶连续导数,则 是f(x)的m重零点之充要条件为： 证明 必要性 设 是f(x)的m重零点，则由 当 时 当 k=m时 充分性 设 使得 由Taylor公式得 其中0θ1，令 则有 且 根据定义， 为f(x)的m重零点 三. 根的搜索 求方程根的近似值之前，一般需要首先确定隔（有）根区间[a,b]（在[a,b]上方程仅有一个根）。 方法： 通过函数f(x)的增减性、凹凸性、变号特征等，并结合做草图来确定隔（有）根区间[a,b]。 § 1 二分法（对分法） 基本思想：通过区间逐次分半，将有根区间逐步缩小。 设f(x)在[a,b]上连续，f(a)f(b)0，且在[a,b]内f(x)=0至少有一个实根 。记[a,b]为[ , ] ①计算[ a1 , b1 ]中点 的函数值 若 ＝0，则 若 ，则 令 若 ，则 令 新的有根区间[ , ]的长度 ②再计算[ , ]中点 的函数值 若 ＝0，则 若 ，则 令 若 ，则 令 新的有根区间[ , ]的长度 如此对分下去，则得到一系列有根区间 且 由 得 k=1,2,3 当对分过程无限继续下去，则有根区间 必收缩为一点，即 具体做法 (1)给定ε每步检查 是否成立， 若成立，取 ，否则继续对分。 (2) 令 ，先确定对分次数k，再计算 (3) 误差估计为 优点 对函数性质要求不高（只要函数连续）； 计算简单，且可达到任意精度。 缺点 计算量大； 不能求复根与偶重根。 §2 迭代法的算法和理论 一不动点迭代法 对给定的方程f(x)=0，将其变为等价的方程 构造 k=1,2 称为迭代序列， (x)称为迭代函数。 称为迭代格(公)式或迭代过程。 当 (x)连续时，若 则有 即 故序列 的极限为方程x= (x)(或f(x)=0)的根 若 满足 ， 称 为 (x)的不动点。 即映射关系 将 映射到 本身。 因 求f(x) ) 的零点等价求 的不动点。 也称 k=1,2 为不动点迭代法 （简单迭代法或逐次逼近法）。 迭代序列的收敛性及收敛速度依赖于 迭代函数的选取。 二 不动点迭代法的一般理论 定理（不动点定理） 已知x= (x)，若 且 ①对 ，有 ; ②存在常数0L1,使对 ，有 则 ① (x)在[a,b]上有唯一的不动点 ； ②对任意 k=0,1,2 产生的序列 必定收敛到 (x)的不动点 ； ③ 有误差估 证明 ① 作辅助函数ψ(x) = (x)-x 由于 (x)在[a,b]上连续， 则ψ(x)∈C[a,b]， 且 ψ(a) = (a)-a≥0 ψ(b) = (b)-b≤0 故由连续函数的介值定理， 至少存在 使ψ( )=0 。 即