Ch4_Sec_1非线性方程的迭代解法.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *2）若|f (xk+1)| ≥| f (xk) | ，则当?≤??且|f(xk+1)| ε2 ，取x*? xk，计算过程结束；否则若?≤??，而|f(xk+1)| ≥ε2时，则把xk+1加上一个适当选定的小正数，即取xk+1+?作为新的xk值，并转向(3)重复计算；当?＞??；且|f(xk+1)| ≥ε2 ，则将下山因子缩小一半，取?/2代入，并转向（3）重复计算。 下山因子下界 开 始 输入x0, |x - x0|ε1 输出x 结 束 否 是 x0=x λ=λ/2 |f(x)||f(x0)| 是 否 是 否 否 是 例：求方程f (x) = x3 – x – 1 = 0 的根 牛顿法的计算结果 ：kxkf(xk)00.6-1.384117.95716.439211.16923500.45.6147.53.443.0096132462.112.……121.7-1.08×10-9 例：求方程f (x) = x3 – x – 1 = 0 的根 牛顿下山法的计算结果 ：kxk010.611/2513668131132472 |f(x1)||f(x0)| 下山成功 4.1.6 求方程m重根的Newton迭代法 Newton迭代法降为一阶！ 改进办法一 为二阶收敛 改进办法二 此方法不管根的阶数，收敛精度至少二阶！ 例：求方程f (x) = x4-1.4x3-0.48x2+1.408x-0.512=0的根 已知有一个三重根s=0.8. 改进办法的计算结果 ：kxk0110.794 019 93320.799 962 73430.800 019 389 牛顿法的计算结果 ：kxk0110.935 483 871……100.803 649 381 例：求方程f (x) = (x-1)3的根. 牛顿法的计算结果 ：kxk00.810.666720.777830.8519……100.9942 4.1.7 割线法(弦截法) 牛顿法收敛速度快计算f’(xk)(可能复杂), 且其绝对值可能很小 双点弦截法 ： 切线斜率　?　割线斜率 需要2个初值 x0 和 x1. x0 x1 x2 双点弦截法具有局部收敛性且收敛速度的阶至少是1.618 例 求方程x3-x-1=0在x=1.5附近的根. 解 用割线法, 取x0=1.5, x1=1.4, 迭代公式为 由此知, x2=1.33522x3=1.32541…… 用牛顿法 kxk 01.5 1.34783 1.32520 1.32472 1.32472 4.1.8 单点弦截法 ： 牛顿法一步要计算 f 和 f ’，相当于2个函数值。现用f 的值近似 f ’ : x0 x1 切线 割线 切线斜率　?　割线斜率 x2 单点弦截法具有局部收敛性且线性收敛速度 开 始 输入x0, x1,ε x2=x1- f (x1)(x1-x0) /(f(x1)-f(x0)) |x2 – x1|ε 输出x2 结 束 否 是 x1=x2 例 求方程x3-x-1=0在x=1.5附近的根. 解 用割线法, 取x0=1.5, x1=1.4, 迭代公式为 由此知, x2=1.33522x3=1.32532…… 用牛顿法 kxk 01.5 1.34783 1.32520 1.32472 1.32472 小结 一、二分法方法简单，f(x)的性质要求高, 线性收敛主要用于确定初始值x0. 二、迭代法非局部: [a, b],局 部:1. 牛顿法-----非局部收敛, f(x)的性质要求高, 单根附近处平方收敛.2. 割线法-----局部收敛, f(x)的性质要求不高, 线性收敛. 作业 课后习题1、2、4、5、7、9 * * * * * * * * * * * * 证（1）F(x)=x-? (x) ?C[a, b] F(a)=a-?(a)F(b)=b-?(b) (2) (1) 当x?[a, b]时，?(x)?[a, b] 证 (3) 又 迭代法的结束条件 例：求方程x3-3x+1=0在[0, 0.5]内的根，精确到10-5.x1 = ? (0.25) = 0.x2 = ? (x1) = 0.x3 = ? (x2) =0.x4 = ? (x3) =0.x5 = ? (x4) =0.x6 = ? (x5) =0.x7 = ? (x6) =0. 0?(x)≤0.