数值分析3-1研.ppt

* 例3-3： 用迭代法求方程 在区间[0,1] 内的一个实根，精确到4位有效数字。 解： 当 满足定理3.1的条件① 3.2.2 不动点的存在性与迭代法的收敛性 * 定理3.1的条件②未得到满足。 缩小有根区间。用二分法。 ，条件②不满足 ，有根区间为[0,0.5]. 3.2.2 不动点的存在性与迭代法的收敛性 * 条件②不满足 条件②满足 当 条件①不满足 3.2.2 不动点的存在性与迭代法的收敛性 * 条件②满足 条件①满足 根据定理3.1，对任意 x0∈ [0.4444,0.4839]，此迭代格式收敛。 对于定理3.1，迭代初值必须同时满足两个条件，迭代格式才收敛。收敛条件要求较高。 3.2.2 不动点的存在性与迭代法的收敛性 * fi=inline(1/(x+1)^2 ); x0=0.46;er=1;k=0; while er0.00005 x=fi(x0) er=abs(x-x0); x0=x;k=k+1; end 若取x0＝0.46，k=13, x=0.46558659888486 若取x0＝0.45，k=15, x=0.46558867795970 3.2.2 不动点的存在性与迭代法的收敛性 不动点迭代法-迭代格式的构造 迭代的敛散性条件 迭代法的局部收敛性 迭代法的收敛速度 牛顿法 牛顿法的局部收敛性 3.2 迭代法 3.2.3 迭代法的局部收敛性 定义 对于方程 x=φ(x) ，若在 x*的某邻域 * 内，对任意 x0∈ R，迭代格式 都收敛，则该迭代格式在 x*的附近 局部收敛。 * 若x0靠近根x*，可用 来判断迭代格式局部收敛或发散。 3.2.3 迭代法的局部收敛性 定理3.2 设方程 x=φ(x) 有根x*，且在其某邻域 内, φ(x) 存在一阶连续导数，则 （1）当 时，此迭代格式局部收敛； （2）当 时，此迭代格式发散。 回忆：与定理3.1的区别 * 3.2.3 迭代法的局部收敛性 证明：(1) 设 时， 要证： 即： 定理3.1的条件得到满足 迭代格式局部收敛 * 由定理 3.2知，迭代格式 局部收敛 迭代格式局部收敛 例: 用迭代法求方程 在x0=0.4附近的一个实根。 解： * fi=inline(1/(x+1)^2 ); x0=0.4;er=1;k=0; while er0.00005 x=fi(x0) er=abs(x-x0); x0=x;k=k+1; end 若取x0＝0.4，k=19, x=0.465583582752150 3.2.3 迭代法的局部收敛性 不动点迭代法-迭代格式的构造 迭代的敛散性条件 迭代法的局部收敛性 迭代法的收敛速度 牛顿法 牛顿法的局部收敛性 3.2 迭代法 3.2.4 迭代法的收敛速度 * 定义： 设序列{xk} 收敛于x*， ek= xk - x*， (k=0,1,2,…) 。若有非零常数 c 和正整数 p，使 则称序列{xk} p阶收敛。 线性收敛：p=1, 且 超线性收敛：p1。其中，p=2，平方收敛。 序列的收敛阶数越高,收敛速度越快。 迭代格式产生的序列 p 阶收敛，则称此迭代格式 p阶收敛。 * 例： 2个迭代格式分别是线性收敛和平方收敛： 其中， ，若要求精度 估计迭代格式收敛所需的迭代次数。 3.2.4 迭代法的收敛速度 * 解： （1）由 应迭代32次 （2）由 应迭代4次 * 定理3.3：设φ(x) 在x* 附近的某邻域内有p ( p ≥1 )阶连续导数，且 则对一个任意靠近x* 的初值，迭代公式 p 阶收敛，且 如果p＝1，则 要求 3.2.4 迭代法的收敛速度 * 3.2.4 迭代法的收敛速度 证明：(1)设 根据定理3.2，可以断定迭代过程 具有局部收敛性。 将 在根 x*处做泰勒展开，利用条件 则有： 迭代过程为p阶收敛 * (2) 若 线性收敛。 定理3.1已证明 迭代法的收敛速度和迭代函数的选取有关 如何构造高阶迭代函数？ 3.2.4 迭代法的收敛速度 不动点迭代法-迭代格式的构造 迭代的敛散性条件 迭代法的局部收敛性 迭代法的收敛速度 牛顿法 牛顿法的局部收敛性 3.3 牛顿法 * 3.3 牛顿法 简单迭代法（不动点迭代法） 不易得到收敛速度快的迭代函数，且可能发散 牛顿法：用线性方程近似代替非线性方程 f(x)=0 可以方便地得到快速收敛的迭代公式 * 设方程 f (x) = 0