NA-2-1非线性方程.pptVIP

第二章一元非线性方程的求解 一、非线性方程求根的基本问题包括：根的存在性、根的隔离和根的精确化 第二节 二分法 算法 计算结果 误差分析 第三节 不动点迭代 第三节 不动点迭代 第三节 不动点迭代 fi=inline(0.5*sqrt(10-x^3)); x0=1.5;er=1;k=0; while er0.00001 x=fi(x0); er=abs(x-x0); x0=x;k=k+1; end k=16 x0=1.3652 fi=inline(sqrt(10/(4+x))); x0=1.5;er=1;k=0; while er0.00001 x=fi(x0); er=abs(x-x0); x0=x;k=k+1; end k=6 x0=1.3652 迭代序列的收敛性和收敛速度以及误差估计。 例 用收敛的充分条件验证例2.2 的两种迭代格式的收敛性 Newton迭代 Newton迭代 全局收敛性定理 重根情形 提高收敛阶 降低初始点的要求 割线法 五、迭代法和离散动力系统 Newton 法也可以看作一类特殊的加速迭代 取 ?(x) = x-f(x) 定理 设 f ?C2[a, b]，且 f 满足 (1) f (a) f (b) 0 (2) 对? x?[a, b]， f’(x) ? 0 且 f”(x) ? 0 ； 则对任意初始点 x0 ? [a, b]，Newton 法产生的序列收敛到 f 在 [a, b] 的唯一零点 x*。 (3) 数值分析 数值分析 （局部收敛性）设 f ?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根，且 ? 0，则存在 x* 的邻域 使得任取初值 ，Newton’s Method产生的序列{ xk } 收敛到x*，且满足 定理 证明：Newton’s Method 事实上是一种特殊的不动点迭代 其中 ，则 收敛 由 Taylor 展开： 数值分析 数值分析 只要 f ’(x*) ? 0，则令 可得结论。 注：在单根 附近收敛快，是平方收敛的． 由 Taylor 展开： 数值分析 数值分析 y x 0 b a x0 y x 0 b a x0 y x 0 b a x0 y x 0 B a x0 数值分析 数值分析 例用Newton迭代法求 在隔根区间[1.4,1.5] 内的根，要求准确到小数点后第4位。 （1）牛顿迭代公式为 （2） 当 时有， 因 ,故取 ,牛顿迭代法收敛。 例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=1.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=6 er=4.3596e-009 8/21 设 x* 是 f(x) 的 m(m?2) 重根，Newton法是否收敛？ Taylor 展式 Newton 迭代： 线性收敛。 且重数 m 越高，收敛越慢。 提高收敛速度 但 m 通常无法预先知道！ 法一：取 二阶收敛 法二：将求 f(x) 的重根转化为求 另一个函数 的单根。 构造针对 ?(x) 的具有二阶收敛的 Newton 迭代： 令 ，则 x* 是 ?(x) 的单重根。 例：求 sin(x)-x/6=0 的正根。 Newton 下山法： ?k ?k 为数列 中满足 的最大数。 算法 (Newton下山法) 给定初始点 x0，精度要求 ? 1. 如果 | f (xk)| ? , 停机，输出 xk 2. 计算 ，? =1 3. 如果| f (xk+?dk)| | f(xk)|，令 xk+1= xk+?dk ，返回第1步； 否则? 折半，重新计算第3步 Newton 法的收敛依赖于初始点的选取。 数值分析 数值分析 1、优点：牛顿迭代法