数值分析__第7章非线性方程和方程组的数值解法.ppt

第7章 非线性方程与方程组的数值解法;7.1 方程求根与二分法 ;若 可分解为 ;它在整个 轴上有无穷多个解，若 取值范围不同，解也 不同，因此讨论非线性方程(1.1)的求解必须强调 的定 义域，即 的求解区间; 迭代法要求先给出根 的一个近似，若 且 ，根据连续函数性质可知 在 内至少有一个实根，这时称 为方程(1.1)的有根区间.; 例1 求方程 的有根区间.; 检查 与 是否同号， 如果同号，说明所求的根 在 的右侧,这时令 否则 必在 的左侧， 这时令 见图7-1. ; 对压缩了的有根区间 又可施行同样的手续，即用 中点 将区间 再分为两半，然后通过根 的搜索判定所求的根在 的哪一侧，从而又确定一个新的 有根区间 ，其长度是 的一半. ; 就是说，如果二分过程无限地继续下去，这些区间最终必收缩于一点 ，该点显然就是所求的根.; 由于 ; 例2 求方程 ; 计算结果如表7-2. ; 二分法是计算机上的一种常用算法，计算步骤为： ;此时中点 即为所求近似根. ;7.2 不动点迭代法及其收敛性 ;如此反复迭代计算 ; 方程 的求根问题在 平面上就是要确定曲 线 与直线 的交点 ;则点 的横坐标为 ，; 例3 求方程 ;各步迭代的结果见表7-3. ;但若采用方程（2.3）的另一种等价形式; 7.2.2 不动点的存在性与迭代法的收敛性 ; 因 ，; 再证唯一性. ;（2.5）; 再证明估计式（2.5），;在上式令 ，注意到 即得式（2.5）.; 由此可见，只要相邻两次计算结果的偏差 ; 例3中，当 时， , 在区间 中， ，故（2.7）成立.; 7.2.3 局部收敛性与收敛阶 ; 证明 由连续函数的性质，存在 的某个邻域 使对于任意 成立; 讨论迭代序列的收敛速度. ;取 ，对上述4种迭代法，计算三步所得的结果如下表. ; 从计算结果看到迭代法（1）及（2）均不收敛，且它们均不满足定理3中的局部收敛条件.; 定义2 设迭代过程 收敛于方程 的根 ，如果迭代误差 当 时成立下列 渐近关系式; 定理4 对于迭代过程 ，如果 在所 求根 的邻近连续，并且 ;注意到 ，; 在例4中，迭代法（3）的 ，故它只是线性 收敛.;7.3 迭代收敛的加速方法 ;由于 ; 一般情形是由 计算 ，; 7.3.2 斯蒂芬森迭代法; 它的理解为，要求 的根 ，; 实际上（3.3）是将不动点迭代法（2.2）计算两步合 并成一步得到的，可将它写成另一种不动点迭代 ; 定理5 若 为（3.5）定义的迭代函数 的不动点， 则 为 的不动点.反之，若 为 的不动点， 设 存在， ，则 是 的不动点， 且斯蒂芬森迭代法（3.3）是2阶收敛的.; 至于原来已收敛的迭代法（2.2），由定理5可知它可达到2阶收敛. ; 更进一步还可知若（2.2）为 阶收敛，则（3.3）为; 若取 迭代16次得 ，有六位有效数 字.;7.4 牛顿法 ;这是个线性方程，记其根为 ，; 设 是根 的某个近似值，过曲线 上横坐标为 的点 引切线，并将该切线与 轴的交点的横坐标 作为 的新的近似值. ;由于 ;故由（2.9）可得 ; 牛顿法的计算步骤： ; 此处 是允许误差，而 ;或者 ，则方法失败；; 7.4.2 牛顿法应用举例 ;以上两式相除得 ; 解 取初值 ，对 按(4.5)式迭代3次便得到精度为 的结