数值分析牛顿迭代法matlab程序,数值分析4(牛顿迭代法).ppt

To start today, to start now and to iterate, make it a habit. */25 */25 */25 */25 */25 */25 不动点框架: * 收敛性 收敛速度 《数值分析》4 Newton迭代格式 Newton迭代法的收敛性 Newton迭代法收敛速度 弦截法迭代格式 ? ? ? ? * * Nature and Nature' law lay hid in night. God said, "Let Newton be," and all was light. Alexander Pope 给定初值 x0 , 迭代产生数列 x0, x1, x2,···, xn, ··· * 设 x*是方程 f(x)=0 的根, x0是x*的近似值。在 x0 附近对函数做局部线性化 x1比x0更接近于x* x0 x1 x* f(x) = 0 ? ? 化难为易 化繁为简 * 应用——求正数平方根算法 设C > 0, x2 – C = 0 令 f(x) = x2 – C , 则 * 初值: x0=1.5 迭代格式: xn+1=0.5(xn+2/xn) (n = 0,1,2,·····) 例1. 平方根算法求 xn | en | 1.416666666666667 2.45e-003 1.414215686274510 2.12e-006 1.414213562374690 1.59e-012 1.414213562373095 2.22e-016 1.414213562373095 2.22e-016 表1 平方根算法实验 * 收敛性: (1) 符合不动点框架 * (2) 从序列收敛的角度(单调有界序列) 由此可知平方根算法具有 2 阶收敛速度。 思考: 如何求倒数、平方根和立方根？ Newton迭代法的局部收敛性 定理2.7 设 f(x) 在点x*的某邻域内具有二阶连续导数, 且 f(x*)=0和 f′(x*) ≠ 0, 则对充分靠近点x*的初值x0, Newton迭代法至少平方收敛。 所以Newton迭代法至少平方收敛。 ? * 例2.求 x3 +10x – 20 =0 在 x0=1.5 附近的根 解:取 牛顿迭代格式 则有 n xn | en | 0 1.5 1 1.59701492537 0.002452808741981 2 1.59456374876 1.632137654805e-06 3 1.59456211663 7.227551890309e-13 4 1.59456211663 2.220446049250e-16 表2 牛顿迭代法实验 * 注释1: 为了二次收敛有意义我们需要f′(x)相除, 这个假设是关键的。 f(x)=x3 –3x + 2 = 0 在x*=1附近 * x* x0 ? x0 ? x0 Newton方法收敛性依赖于x0 的选取。存在 x0使Newton迭代法陷入死循环。 注释2: * Newton迭代法的变型－弦截法 由于 代入牛顿迭代格式 ? x0 x1 * n xn | en | | en+1 |/| en |1.618 1 -1.5 5.00e-001 2 -2.5 5.00e-001 1.5347 3 -1.83783783783 1.62e-001 0.4978 4 -1.95420890762 4.57e-002 0.8691 5 -2.00552244119 5.52e-003 0.8109 6 -1.99982796307 1.72e-004 0.7742 7 -1.99999936831 6.31e-007 0.7785 8