第04章非线性方程和方程组的数值解法要素.ppt

4.1 引言 * 第4章 非线性方程的数值解法 则可用搜索法求有根区间. ? ? + + f(x)的符号 ?1 0 1 2 x 求根问题的三个方面：存在性，分布，精确化。 4.3 不动点迭代法 4.3.1不动点和不动点迭代法 1.5 1.35721 1.33086 1.32588 1.32494 1.32476 1.32473 1.32472 0 1 2 3 4 5 6 7 xk k Fixed-Point Iteration x y y = x x y y = x x y y = x x y y = x x* x* x* x* y=?(x) y= ?(x) y= ?(x) y= ?(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? Fixed-Point Iteration 定理 考虑方程 x = ?(x), ?(x)?C[a, b], 若 ( I ) 当 x?[a, b] 时， ?(x)?[a, b]； ( II ) ? 0 ? L 1 使 | ? ′(x) | ? L 1 对 ? x?[a, b] 成立。 则任取 x0?[a, b]，由 xk+1 = ?(xk) 得到的序列 收敛于?(x) 在[a, b]上的唯一不动点。并且有误差估计式： ? ? ( k = 1, 2, … ) 且存在极限 k 所替代 由微分学中值定理，条件(2)常常用条件： 此时 试在1.5附近讨论根的存在惟一性，并构造两种 不同的收敛迭代格式，再用其中一种收敛迭代格式 计算该方程在1.5附近的一个根( ) 3、已知一个三次方程为 解：设 当 时， 是严格单调增加的，而 所以 在（1,2）内有唯一的根。对此建立两种迭代格式： 和 （2） （1） 对于格式（1） 当 时， 并且 ，故 单调增加，有 故当 时 从而对任意初值 ，迭代格式（1）收敛。 对于格式（2） 当 时， 又由于 ，故 单调下降，有 从而对任意初值 ，迭代格式（2）收敛。 当 时， 用格式（1）计算得到： x0=1.5, x1=1.35721, x2=1.33086, x3=1.32588, x4=1.32494, x5=1.32476, x6=1.32473 4.3.3、局部收敛性 2 1.75 1.73475 1.732631 ? 迭代法(3) 2 1.5 2 1.5 ? 迭代法(2) 2 3 9 87 ? 迭代法(1) x0 x1 x2 x3 ? xk 2 1.75 1.732143 1.732051 ? 0 1 2 3 ? 迭代法(4) k 4.5 Newton迭代法和割线法 原理：将非线性方程线性化 —— Taylor 展开 /* Taylor’s expansion */ 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间。 将 (x* ? x0)2 看成高阶小量，则有： 线性 /* linear */ x y x* x0 一般地 Newton - Raphson Method 定理 （局部收敛性）设 f ?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根，且 f ’(x*) ? 0，则存在 x* 的邻域 使得任取初值 ，Newton’s Method产生的序列{ xk } 收敛到x*，且满足 二、牛顿法应用举例 0.5 0.57102 0.56716 0.56714 0 1 2 3 xk k 10 10.750000 10.723837 10.723805 10.723805 0 1 2 3 4 xk k 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 *