数值计算方法no非线性方程求根.pptxVIP

数值计算方法No非线性方程求根 2 第二章 非线性方程求根 第1页/共55页 3 本章要点 本章主要介绍非线性方程求根方法,尤其是迭代法 主要方法 二分法、简单迭代法、Newton迭代法、SOR方法和Aitken加速方法 P. 本章作业 第2页/共55页 4 有唯一根 有多根 所有根均为单根 有重根 然后在每个区间上判断是否有根 第3页/共55页 5 若成立 统计根的个数 则所有的有根区间均为单根区间 则继续对分区间,并重新判断 直到找到所有根的所在区间 然后在每个有根区间进行求根 第4页/共55页 6 可得一系列的小区间和中点 第5页/共55页 7 小区间 中点 显然每个小区间都有单根 搜索法—二分法 第6页/共55页 8 例8. 解: 由于 可知方程的解在区间[0,10]内 将区间[0,10]等分成三等份 [0,3.33] [3.33,6.67] [6.67,10] [0,3.33]内至少有一个根 第7页/共55页 9 [5,6.67] [3.33,5] 将[3.33,6.67]再分成两个区间 [5,6.67]内至少有一个根 [3.33,5]内至少有一个根 [0,3.33] [5,6.67] [3.33,5] 因此找到了三个有单根的区间 依此类推结果为 对分 第8页/共55页 10 2.2 迭代法 方程是在科学研究中不可缺少的工具 方程求解是科学计算中一个重要的研究对象 几百年前就已经找到 了代数方程中二次至 五次方程的求解公式 但是,对于更高次数 的代数方程目前仍 无有效的精确解法 对于无规律的非代数方程的求解也无精确解法 因此,研究非线性方程的数值解法成为必然 第9页/共55页 11 设非线性方程 --------(1) 本节主要研究单根区间上的求解方法 第10页/共55页 12 一、简单迭代法(基本迭代法) --------(2) 将非线性方程(1)化为一个同解方程 继续 --------(3) 称(3)式为求解非线性方程(2)的简单迭代法 第11页/共55页 13 则称迭代法(3)收敛,否则称为发散 --------(4) 如果将(2)式表示为 与方程(2)同解 收敛 第12页/共55页 14 例1. 解: (1) 将原方程化为等价方程 发散 第13页/共55页 15 显然迭代法发散 (2) 如果将原方程化为等价方程 第14页/共55页 16 仍取初值 x2 = 0.9644 x3 = 0.9940 x4 = 0.9990 x5 = 0.9998 x6 = 1.0000 x7 = 1.0000 依此类推,得 已经收敛,故原方程的解为 同样的方程 不同的迭代格式 有不同的结果 什么形式的迭代法 能够收敛呢? 迭代函数的构造有关 第15页/共55页 17 定理1. --------(5) --------(6) --------(7) (局部收敛性) 第16页/共55页 18 证: 由条件(1) 由根的存在定理, 由 第17页/共55页 19 由微分中值定理 第18页/共55页 20 证毕. 第19页/共55页 21 定理1指出, 由(6)式,只要 因此,当 迭代就可以终止, 只要构造的迭代函数满足 此时虽收敛但不 一定是唯一根 --------(8) 第20页/共55页 22 例2. 用迭代法求方程的近似解,精确到小数点后6位 解: 本题迭代函数有两种构造形式 因此采用迭代函数 第21页/共55页 23 d1 = 0.1000000 d2 = -0.0105171 d3 = 0.1156e-002 d4 = -0.1265e-003 d5 = 0.1390e-004 d6 = -0.1500e-005 d7 = 0.1000e-006 由于|d7| =0.1000e-0061e-6 因此原方程的解为 x7 = 0.090525 x1 = 0.1000000 x2 = 0.0894829 x3 = 0.0906391 x4 = 0.0905126 x5 = 0.0905265 x6 = 0.0905250 x7 = 0.0905251 第22页/共55页 24 由定理1的(7)式出, 迭代法收敛就越快 定义1. --------(9) 第23页/共55页 25 不可能直接确定 第24页/共55页 26 定理2. 第25页/共55页 27 例3. 为线性收敛