2-2数值分析.pptVIP

6 弦截法 弦截法几何含义 提高收敛速度 但m通常无法预先知道！ 法一：取 二阶收敛 法二：将求f(x)的重根转化为求另一个函数 的单根. 构造针对?(x) 的具有二阶收敛的Newton迭代： 令 ，则x*是?(x)的单重根. 例 取初始点x0=1.5, 分别用Newton法和求重根的两种方法计算 f(x)=(x+1)(x-1)2=0 解 (1) Newton法迭代公式 1.0001 …… 1.0012 1.0024 1.0048 1.0096 xk 13 …… 9 8 7 6 k 1.0191 1.0378 1.0744 1.1441 1.2727 1.5 xk 5 4 3 2 1 0 k (2) 方法一：取 , m=2 迭代公式为 方法二: 迭代公式为 取 1.0000 0.9996 0.9608 1.5 xk 方法二 1.0000 1.0005 1.0455 1.5 xk 方法一 3 2 1 0 k 两种改进方法的结果比较 用Newton迭代公式求解, 只能是线性收敛, 而改进的两种方法都具有二阶收敛, 所以计算速度要快得多. 但是很多问题实现不知道根的重数，所以方法一有时并不实用. Newton 法的收敛依赖于初始点的选取. 5 Newton下山法 如果x0偏离所求根x*较远,则Newton法可能发散,为防止迭代发散,我们对迭代过程再附加一项要求,即具有单调性: 满足这项要求的算法称为下山法 ?k为数列 中满足 的最大数. 弦截法迭代格式： k = 1, 2, 3, . . . 注：弦截法需要提供两个迭代初始值 目的：避免计算 Newton 法中的导数，且具有较高的收敛性（超线性收敛） * 数值分析 第二章 解非线性方程的数值方法 一、 二分法 二、 迭代法 三、 Newton法 对给定方程f(x)=0,可以用各种方法转化成等价方程 二、 迭代法 1 迭代法的基本思想 若x*是f(x)的根,即若 ,则有 称x*为函数 的一个不动点. 设x0是一个近似解，可以构造序列 迭代法(2.2)称为简单迭代法或单点迭代法. 称函数 为迭代函数. 简单迭代法(2.2), 若迭代序列{xk}保持有界,全在 定义域内称为适定的;若进一步有 称为是收敛的. 若 收敛，即存在x*使得 则由?的连续性和 可得 x* = ? (x*)，即x*是?的不动点，也就是f (x)的零点。 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 例 求x3-x-1=0在1.5附近的根x* 解 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 ? x2 ? 2 收敛定理和误差估计 定理1 设 在[a, b]上有连续的一阶导数,且 (1) 有 (2) 则有 (1) 函数 在[a ,b]上存在唯一的不动点x* 由迭代公式(2.2)得到的序列 (2) 都收敛到方程的根x* (3) (4) 证明 (1) 先证明存在性,作函数 由条件(1)可知 由连续函数根的存在定理可知,必有 使得h(x*)=0，即 再证明唯一性,设 也是一个解,即 那么 因为L1, 所以有 (2) 由条件(2)和Lagrange中值定理得 因为L1, 所以当 时, 有 (3)和(4) 利用(2)的方法 令 即分别得 定理1的几点说明 (i) 通常将条件(1)称为映内性； (ii) 条件(2)也可以改为：存在常数L且0L1 使得 结论仍然成立, 这个条件通常称为压缩性； (iii) 结论(3)说明 的误差大概是常数 乘上Lk, 但是一般L未知; (iv) 结论(4)说明 与 有关 因此得到迭代法的终止条件 3 一般迭代法的算法 算法： (1) 取初始点x0,最大迭代次数N和精度要求 并置k=0； (2) 计算 (3) 若