非线性方程迭代解法yjs.pptVIP

7.1 根的搜索与二分法 等步长扫描法求有根区间 用计算机求有根区间：等步长扫描法。 设h0是给定的步长，取x0=a, x1=a+h ， 若f(x0)*f(x1)0则扫描成功；否则令x0=x1, x1=x0+h，继续上述方法，直到成功。 如果 x1b 则扫描失败。再将h 缩小， 继续以上步骤。 等步长扫描算法 算法7.1：（求方程 f(x)=0 的有根区间） (1) 输入 a,b,h ；(2) f0 =f(a) ； (3) x=a+h, f1=f(x) , 若xb 输出失败信息，停机。 (4)若 f1=0 , 输出 x ，已算出方程的一个根，停机。 用二分法(将区间对平分)求解。 令a0=a, b0=b, c0=(a0+b0)/2 若 f(a0)f(c0)0 ，则 [a0, c0] 为有根区间，否则 [c0, b0] 为有根区间 记新的有根区间为 [a1, b1] ， 则[a0, b0] 包含[a1, b1] ，且 b1-a1= (a0+b0)/2 设 所求的根为 x*，则 即 取 为 x* 的近似解 求方程f(x)=0的根的二分法算法 迭代法需解决的三个问题 迭代函数的构造 由迭代函数产生的解序列的收敛性 序列的收敛速度和误差估计 迭代法的几何意义 交点的横坐标 例1 试用迭代法求方程 在区间(1，2)内的实根。 解：由 建立迭代关系 k=0,1,2,3……. 但如果由 建立迭代公式 仍取 ，则有 ， 显然结果越来越大， 是发散序列. Steffensen迭代格式几何解释 例2 试用Steffensen算法求解方程 的根。 解法一、取 ，由 解法二、取 ， 对于该迭代函数在一般迭代法中是发散的，而Steffensen格式却是收敛的。 Aitken-Steffensen迭代算法 例 用Newton法计算 。 解： 牛顿迭代法 定理 （收敛的充分条件）设 f ?C2[a, b]，若 (1) f (a) f (b) 0；(2) 在整个[a, b]上 f ”不变号且 f ’(x) ? 0； (3) 选取 x0 ? [a, b] 使得 f (x0) f ”(x0) 0； 则Newton’s Method产生的序列{ xk } 收敛到f (x) 在 [a, b] 的唯一根。 有根 根唯一 产生的序列单调有界，保证收敛。 定理 （局部收敛性）设 f ?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根，且 f ’(x*) ? 0，则存在 x* 的邻域 使得任取初值 ，Newton’s Method产生的序列{ xk } 收敛到x*，且满足 证明：Newton’s Method 事实上是一种特殊的不动点迭代 其中 ，则 收敛 由 Taylor 展开： 只要 f ’(x*) ? 0，则令 可得结论。 在单根 /*simple root */ 附近收敛快 ? Newton’s Method 有 ，只要 就有 p ? 2。至少是二阶收敛。 算法7.5（牛顿法） 取初始点x0，最大迭代次数N和精度要求ε，置k:=0； (2)计算 xk+1= xk- f(xk)/f’(xk) (3)若|xk+1-xk| ε, 则停算； (4)若k=N, 则停算 ；否则，置k:=k+1,转步2 . 程序见P161。 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 牛顿法的优点 牛顿法是目前求解非线性方程 (组) 的