第7单元 非线性方程迭代解法yjs.pptVIP

7.1 根的搜索与二分法 首先确定有限区间：依据零点定理。 设 ，且 ，则方程 在区间 上至少有一个根。如果 在 上恒正或恒负，则此根唯一。 等步长扫描法求有根区间 用计算机求有根区间：等步长扫描法。 设h0是给定的步长，取x0=a,x1=a+h ， 若f(x0)*f(x1)0则扫描成功；否则令x0=x1, x1=x0+h，继续上述方法，直到成功。 如果 x1b 则扫描失败。再将h 缩小， 继续以上步骤。 等步长扫描算法 算法：（求方程 f(x)=0 的有根区间） (1) 输入 a,b,h ；(2) f0 =f(a) ； (3) x=a+h, f1=f(x) , 若xb 输出失败信息，停机。 (4)若 f1=0 , 输出 x ，已算出方程的一个根，停机。 用二分法(将区间对平分)求解。 令a1=a, b1=b, c1=(a1+b1)/2 若 f(a1)f(c1)0 ，则 [a1, c1] 为有根区间，否则 [c1, b1] 为有根区间 记新的有根区间为 [a2, b2] ， 则[a1, b1] 包含[a2, b2] ，且 b2-a2= (a1+b1)/2 设 所求的根为 x*，则 即 取 为 x* 的近似解 求方程f(x)=0的根的二分法算法 迭代法需解决的三个问题 迭代函数的构造 由迭代函数产生的解序列的收敛性 序列的收敛速度和误差估计 迭代法的几何意义 交点的横坐标 例1 试用迭代法求方程 在区间(1，2)内的实根。 解：由 建立迭代关系 k=0,1,2,3……. 计算结果如下: 但如果由 建立迭代公式 仍取 ，则有 ， 显然结果越来越大， 是发散序列. 例2 试用Steffensen算法求解方程 解法一、取 ，由 解法二、取 ， 对于该迭代函数在一般迭代法中是发散的，而Steffensen格式却是收敛的。 Steffensen迭代格式几何解释 Steffensen迭代算法 例 用Newton法计算 。 解： 牛顿法迭代公式 证明：Newton’s Method 事实上是一种特殊的不动点迭代 其中 ，则 收敛 由 Taylor 展开： 只要 f ’(x*) ? 0，则令 可得结论。 在单根 /*simple root */ 附近收敛快 ? Newton’s Method 有 ，只要 就有 p ? 2。至少是二阶收敛。 算法7.5（牛顿法） 取初始点x0，最大迭代次数N和精度要求eps，置k:=0； (2)计算 xk+1= xk- f(xk)/f’(xk) (3)若|xk+1-xk|eps, 则停算； (4)若k=N,则停算 ；否则，置k:=k+1,转步2 . 程序见P161。 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 科大研究生学位课程 数值分析 第7章 非线性方程迭代解法 1．根的存在性。方程有没有根？如果有根，有几个根？ 2．这些根大致在哪里？如何把根隔离开来？ 3．根的精确化 本章介绍求解非线性方程 f(x)=0 的几种常见和有效的数值方法. 其中?(x)是高次多项式函数或超越函数.如 ?(x)=3x5-2x4+8x2-7x+1 ?(x)=e2x+1-xln(sinx)-2 等等. 求非线性方程 确定方程的有根区间 计算根的近似值 f(x)=0的根的方法 分为两步： (5) 若 f0 f1 0 .输出 a,x, [a,