数值计算方法(第2章)教学材料.ppt

Steffensen迭代格式几何解释 Steffensen迭代算法 Steffensen迭代算法 2.2.2 Steffensen加速收敛法 迭代法收敛的阶 定义2.2.1 设序列 收敛到 ，若有实数 和非零常数C，使得 其中， ，则称该序列是p 阶收敛的，C 称为渐进误差常数。 迭代法收敛的阶 当p=1时，称为线性收敛(0c1)；(达朗贝尔判别法) 当p1时，称为超线性收敛； 当p=2时，称为平方收敛或二次收敛。 2.3 Newton迭代法 设x * 是方程f (x ) = 0的根，又x0 为x * 附近的一个值 ，将f (x ) 在x0附近做泰勒展式 令 ，则 Newton迭代法 去掉 的二次项，有： 即 以x1代替x0重复以上的过程，继续下去得： Newton迭代法 以此产生的序列{Xn}得到f(x)=0的近似 解，称为Newton法，又叫切线法。 Newton迭代法几何解释 几何意义 Newton迭代法收敛性 定理2.3.1 设函数 ，且满足 若初值 时，由Newton法产生的序列收敛到 在Δ上的唯一根,并且收敛阶至少为2, 。 证明迭代序列的收敛性 例题 例2.3.1 用Newton法求 的近似解。 解：由零点定理。 例2.3.2 用Newton法计算 。 解： 压缩映像原理 则 ， 即 由条件2） 是 上的连续函数 是 上的连续函数。 故由零点定理 在 上至少有一根 压缩映像原理 再证根的唯一性 设有 均为方程的根 则 因为 0L1 ，所以只可能 ， 即根是唯一的。 压缩映像原理 最后证迭代序列的收敛性 与n 无关，而0L1 即 压缩映像原理 误差估计 若 满足定理2.2.1条件，则 这是事后估计，也就是停机标准。L越小，收敛速度越快。 这是事前估计。选取n，预先估计迭代次数。 例题 例2.2.2 证明函数 在区间[1，2]上满足迭代收敛条件。 证明： 例题 例题 若取迭代函数 ， 不满足压缩映像原理，故不能肯定 收敛到方程的根。 简单迭代收敛情况的几何解释 压缩映像原理推论 推论 设迭代函数 在闭区间 上满足 定理2.2.2 设迭代函数 在 上连续可微 2.2.2 Steffensen加速收敛法 迭代法收敛的阶 定义2.2.1 设序列 收敛到 ，若有实数 和非零常数C，使得 其中， ，则称该序列是p 阶收敛的，C 称为渐进误差常数。 迭代法收敛的阶 当p=1时，称为线性收敛(0c1)； 当p1时，称为超线性收敛； 当p=2时，称为平方收敛或二次收敛。 迭代法收敛的阶 定理2.2.2’ 设 是方程 的不动点， 若为足够小的正数 。如果 且 ，则从任意 出发，由 产生的序列 收敛到 ，当 时敛速是线性的。 迭代法收敛的阶 证明： 满足压缩映像原理 迭代法收敛的阶 敛速是线性的 线性收敛到 。 Steffensen迭代格式 由线性收敛知 当n充分大时有 即