[2018年最新整理]4-牛顿迭代法.pptVIP

*/26 * */24 */24 */24 */24 返回 */24 牛顿迭代法 Newton迭代格式 Newton迭代法的收敛性 弦截法迭代格式 数值实验题介绍 初值: x1=1.5 迭代格式: xn+1=0.5(xn+2/xn) (n = 1,2,·····) 平方根算法求 xn Error 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 牛顿(Newton)迭代法 基本思想： 将方程 f(x)=0中函数 f(x)线性化,以线性方程的解逼近非线性方程的解. 设函数 f(x)在有根区间[a,b]二次连续可微,则 f(x)在x0处的泰勒展开式为: 只取关于x线性项，有 设 ，其解记为 (*) 牛顿迭代法的几何意义 y x O x0 x1 x2 图2.3 牛顿迭代法在单变量情况下又称为切线法. (*) 迭代格式(*)称为牛顿迭代法. 构造迭代格式： 设C 0, x2 – C = 0 令 f(x) = x2 – C , 则 应用——求正数平方根算法 例1 设C0，证明由迭代格式 ( n= 0，1，……)产生的迭代序列 {xn}，对任意的x00，均收敛于 ;且具有 2 阶收敛速度。 分析：由迭代格式，有 证明：由迭代格式，有 等式两端同减 ，配方得 同理有 将上面两式相除有 反复递推，得 令 则有 化简得 由此可知, 平方根迭代具有 2 阶收敛速度 Newton迭代法的局部收敛性 定理 2.7 设 f(x) 在点x*的某邻域内具有二阶连续导数,且设 f(x*)=0, f ’(x*) ≠ 0, 则对充分靠近点x*的初值x0, Newton迭代法至少平方收敛. 所以, Newton迭代法至少平方收敛。 ? 例2. 求 f(x)=xex – 1= 0 在 x0=0.5 附近的根 解: 迭代格式为 (n = 0, 1, ·····) f=inline(x*exp(x)-1); f1=inline((x+1)*exp(x)); x0=0.5;er=1;k=0; while er0.00001 x=x0-f(x0)/f1(x0); er=abs(x-x0) x0=x;k=k+1 end x = 0.5671 k=4 er=1.2347e-010 缺陷 1.被零除错误 2.程序死循环 y = arctan x 方程: f(x)=x3 –3x + 2 = 0 在重根x*=1附近,f’(x)近似为零 对 f(x) = arctan x 存在 x0,Newton迭代法陷入死循环 例3 用牛顿迭代法解方程 f（x）= x3 – x – 3= 0 . x1= – 3， x2= – 1.9615， x3= – 1.1472 ， x4= – 0.0066， …… 数列中的项以四项为一个周期重复（死循环）. 取 x0=0, (n = 0, 1, ·····) 例4 用牛顿迭代法解方程 f（x）= x e – x=0. 初值取x0= 2 ， x1= 4， x2= 5.33333， ……， x15= 19 …… f（x15）=0.0000000536 f’0, f”0 f’0, f”0 f’0, f”0 f’0, f”0 牛顿迭代法收敛的四种情况 定理: 若函数f(x) 在[a，b] 上满足条件 则方程 f(x) = 0 在[a，b] 上有唯一根 x*，且由初值x0按牛顿迭代公式求得的序列{ xn} 二阶收敛于x*。 （1）f(a) f(b) 0； （2）f ’(x)，f ”(x) 在[a，b]上连续且不变号（恒为正或恒为负）； （3）取x0∈[a，b] 使得 f(x0)f ”(x0) 0。 设x*是方程 f(x)=0 的根, x0和x1是x*附近的两个点. Newton迭代法的变形－弦截法 曲线 y=f(x) 在点(x0, f(x0))和点 (x1,f(x1)