数值分析第八章.pptVIP

例 用牛顿法求方程 的根, 取 x0=1.5. 解： 迭代公式为 代入初值得： 可见牛顿法收敛很快. 定理 （收敛的充分条件）设 f ?C2[a, b]，若 (1) f (a) f (b) 0；(2) 在整个[a, b]上 f ??不变号且 f ?(x) ? 0；(3) 选取 x0 ? [a, b] 使得 f (x0) f ??(x0) 0； 则Newton’s Method产生的序列{ xk } 收敛到 f (x) 在 [a, b] 的唯一根. 有根 根唯一 产生的序列单调有界，保证收敛. （局部收敛性）设 f ?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根，且 f ?(x*) ? 0，则存在 x* 的邻域 使得任取初值 Newton’s Method产生的序列{ xk } 收敛到x*，且满足 定理 证明：Newton’s Method 事实上是一种特殊的不动点迭代其中 则 收敛 ? 由 Taylor 展开： 只要 f ?(x*) ? 0, 则令 可得结论. 在单根 附近收敛快 再证： 注：Newton’s Method 收敛性依赖于x0 的选取. x* x0 ? x0 ? x0 ? 当x*是方程的m重根时, 牛顿迭代法是线性收敛的, 且 则称 x*为方程 f (x)=0的m重根. ? 设 ? 求重根的改进Newton法 此时迭代函数 则 ? 迭代法具有2阶收敛，但要知道x*的重数m. 其中m为 f (x)=0的根的重数. ? 构造求重根的另一迭代法 则 故x*是? (x)=0的单根. 对它用牛顿法，其迭代函数为 从而可构造迭代法 它是二阶收敛的. 例 方程 的根 是二重根，用上述三种方法求根. 解 先求出三种方法的迭代公式 (3) 方法3 取初值 x0=1.5, 计算结果见下表 (1) 牛顿法 (2) 方法2 三种方法数值结果 1 2 3 方法(1) 方法(2) 方法(3) 1.458333333 1.436607143 1.425497619 1.416666667 1.414215686 1.414213562 1.411764706 1.414211438 1.414213562 ? 计算三步, 方法(2)及(3)均达到10位有效数字, 而用牛顿法只有线性收敛, 要达到同样精度需迭代30次. ? 弦截法(割线法, 线性插值法) ? Newton’s Method 一步要计算 f 和 f ?，相当于2个函数值，比较费时. 现用 f 的值近似 f ?，可少算一个函数值. ? 在Newton迭代法中用差商 代替导数 得到迭代公式 弦截法 ? 须两个初始点x0, x1, 通常可取有根区间的两个端点. Secant Method 过曲线上两点 的直线为 它与x轴的交点为 以割线近似曲线，故弦截法又称为割线法. 设 f (x)在 x*邻近二阶连续可导，且 f (x*)=0, f ?(x*)?0, 则存在? 0，当x0 , x1?[ x*? ?，x*+ ? ]时，由弦截法产生的数列{xn}收敛于x*, 且收敛阶为1.618. (方程 的正根) 定理 例 用弦截法求下列方程的根. 解 因为 f (1.5)0, f (2)0, [1.5, 2]为有根区间. 取 按弦截法计算公式求, 结果如下 ? 此计算结果表明, 迭代5次所得近似解精确到8位有效数字. 它的收敛速度比Newton法慢. §4 抛物线法(Müller法） 基本思想: 过 作一抛物线 P2(x), 该抛物线与x轴的一个交点 xn+1作为根x*的下一个近似. Secant and Muller’s Method Muller’s Method y(x) x Secant line x1 Fit a parabola (quadratic) to exact curve x2 x3 Parabola 过 或 其中 的抛物线 P2(x) 该抛物线与x轴的交点为， 在xn , xn?1 , xn?2三个近似根中, 自然假定xn更接近所求的根x*, 这时为了保证精度，选上式中较接近xn的一个值