11 非线性方程求根.pptVIP

引理2.1 如果 ,满足条件: ; (2) 则 在 [a, b] 有唯一的不动点 x*. 不动点迭代序列的收敛速度 定理2.4 如果 ,满足条件: ; (2) 则对任意的 x0∈ [a, b] , 迭代格式 产生的序列 { xn }收敛到不动点 x*,且有 不动点迭代序列的收敛速度 数列的 r 阶收敛概念 设 , 若存在 a0 , r0 使得 则称数列{xn} r 阶收敛. 特别: (1) 收敛阶r=1时,称为线性收敛; (2) 收敛阶r1时,称为超收敛; (3) 收敛阶r=2 时,称为平方收敛。 序列的收敛阶数越高,收敛速度越快。 f’0, f”0 f’0, f”0 f’0, f”0 f’0, f”0 牛顿迭代法收敛的四种情况 * * * * * * 非线性方程求根 何国良 hegl@uestc.edu.cn 数学科学学院 参考文献 [1]李庆扬 关治 白峰杉, 数值计算原理(清华) [2]蔡大用 白峰杉, 现代科学计算 [3] 李庆扬 等, 数值分析 [4]Numerical Analysis (Seventh Edition) 数值分析 （第七版 影印版） [5]David Kincaid,数值分析(第三版) [6] John H. Mathews,数值方法(MATLAB版) 模型误差: 建立数学模型时所引起的误差； 误差分类： 舍入误差:计算机表示的数的位数有限，通常用四舍五入的办法取近似值，由此引起的误差. 截断误差:求解数学模型时，用简单代替复杂,或者用有限过程代替无限过程所引起的误差 观测误差:测量工具的限制或在数据的获取时随机因素所引起的物理量的误差； 数值计算中的一些基本概念 假设某一数据的准确值为 x*,其近似值 为 x，则称 而称 为 x 的相对误差 误差的有关概念 e(x)= x - x* 为 x 的绝对误差。 避免绝对值小的数作除数 这一原则主要指尽量避免除数绝对值远远小于被除数绝对值的除法。 设 （ x≠0），如果 x 的绝对值远小于 y 的绝对值，由于 二、数值计算中的一些基本原则 避免两个相近的数据相减 如果 y ≈ x，现分析两个数的近似数作减法所得结果的误差. 设 z= y – x，则利用误差估计 有相对误差估计 当 y ≈ x 时，有 z≈0，计算结果的相对误差限可能很大，导致数值计算结果的有效数字位数减少。 要防止大数“吃掉”小数 一个绝对值很大的数和一个绝对值很小的数直接相加时，很可能发生所谓“大数吃小数”的现象。 例如，a= 1013，b= 4，设想这两个数在具有12位浮点数计算机系统（12位有效位数系）中相加 a + b= 1013 + 4=1. 0000000000000 ×1013 +0. 0000000000004 ×1013 实际加法操作如下 尽量减少计算工作量 在考虑算法时应注意简化计算步骤，减少运算次数 。 计算工作量小的算法不仅节约运行时间，而且使误差积累小。 算法一： S 0= a0 ， Sk= Sk－1 + ak xk ，（ k= 1，2，…，n ） Pn（x）= Sn 例2 设计算法用于计算多项式 计算一个 n 次多项值需要用 2n 次乘法。 尽量减少计算工作量 算法二： Sn= an， Sk－1= a k －1+ xSk，（k= n，n－1，…，1）， Pn（x）= S0 另一种典型算法是秦九韶算法 计算一个 n 次多项值需要用 n 次乘法。 选用数值稳定性好的算法 不同的算法在执行过程中对数据误差的影响是不一样的。舍入误差对计算结果影响不大的算法被称为数值稳定的算法. 例3 利用递推式计算定积分 ( n= 0， 1， 2， …， 20 )的值。 得递推关系式 算法一： 其中 利用递推式可得20个数据如下表： S 1 0. 36787944117144 S 11 0. 07735222935878 S 2 0. 26424111765712 S 12 0. 07177324769464 S 3 0. 20727664702865 S 13 0