数值分析积分上.ppt

第七章 数值积分与微 分（上） 第七章目录 第七章目录 序(1) 序(2) §1 数值积分的基本概念 构造数值求积公式的基本思想（续） 构造数值求积公式的基本思想（续1） 1.2 代数精度 代数精度(续1） 代数精度（续2） 代数精度举例 代数精度举例（续1） 代数精度举例（续2） 待定系数法注释 1.3 插值型求积公式 插值型求积公式（续） 插值型求积公式代数精度定理 定理2（续） 插值型求积公式举例 §2 牛顿一柯特斯(Newton-Cotes)公式 牛顿一柯特斯(Newton-Cotes)公式（续） 牛顿一柯特斯(Newton-Cotes)公式（续1） 牛顿一柯特斯(Newton-Cotes)公式（续2） 柯特斯系数的性质 柯特斯系数的性质（续） N为偶时的牛—柯公式的代数精度证明 N-C公式应用举例 N-C公式应用举例（续） 2.2 几种低价N-C求积公式的余项 辛卜生公式误差估计式的 推导 辛卜生公式误差估计式的 推导（续） 2.3 牛顿一柯特斯公式的稳定性和收敛性 牛顿一柯公式的稳定性和收敛性（续） §3 复化求积公式 复化梯形公式 复化梯形公式的截断误差 复化梯形公式的数值稳定性讨论 3.2 复化Simpson公式和复化Cotes公式 复化Simpson公式的截断误差 复化Cotes公式 复化求积公式举例 例6(续） 复化求积公式举例（续） 若用复化Simpson公式 例6、例7说明 例6、例7说明（续） §4 逐次分半算法(变步长方法) 4.1 梯形法的递推公式 梯形公式的逐次分半算法（续1） 梯形公式的逐次分半算法（续1） 梯形公式的逐次分半算法（续2） 梯形公式的逐次分半算法（续3） 复化梯形公式的停止计算控制 复化simpson的停止计算控制 4.2 Simpson公式的逐次分半法 Simpson公式的逐次分半法（续） 梯形公式的逐次分半法举例 例8（续） 例8说明 例8说明（续1） 例8说明（续2） 例8说明（续3） 例8说明（续4） §5 龙贝格（Romberg）求积公式 李查逊（Richardson）外推法 李查逊（Richardson）外推法（续） 5.2 龙贝格（Romberg）求积公式 龙贝格(Romberg)求积公式（续1） 龙贝格(Romberg)求积公式（续2） 外推公式 外推公式（续） Romberg方法举例 Romberg方法小结 由于 为复化Cotes序列， 即由Simpson序列可构造出收敛更快的Cotes序列 。 并且我们的具体做法都是利用控制结束的误差式，构成新的，收敛更快的序列，而由前面的推导可知，下面这些公式具有如下规律性： 类似地，也可以推导出： 5.1 外推法 从上面例，我们看到复化梯形序列{T2m}收敛较慢，而利用梯形序列这些较粗略的近似值，重新进行线性组合得到的结果收敛更快，更准确。这种利用若干精略近似值推算更精确的近似值的方法，称为外推法。 下面再举例说明： 从一个基本公式出发，利用加速收敛的技巧，可以构造 出收敛速度更快的近似序列，下面的李查逊（Richardson） 外推法就是这样一种方法。 设用某种数值方法求积分I的近似值，一般设近似值是 步长的函数，记为I1(h)，若有相应的误差关系式： 回到这里的积分问题我们可由变步长公式计算复化梯形 T2m，得到梯形序列{T2m}，它收敛较慢，能否由它构造出 新的，收敛更快的序列？ 将上述Richardson外推法应用到{T2m}，就可构造出计 算简便，收敛却很快的数值积分公式。 T2k与积分之间有下面的误差关系式: 这个基本关系式中系数a1,a2,…,ak,…均与h无关，称为 Euler-Maclaurin公式，利用它可构造出收敛更快的序列 （同前一样）即首先构成这里的具体的误差关系式： 而辛卜生公式至少具有三次代数精度，因此对上述三次多项式H(x) 应准确成立，即有： 其插值 余项为： 因此，辛卜生公式的误差就是对上述误差公式的积分： 3. 柯特斯公式（6-10）的余项为： 根据定理2，牛顿一柯特斯公式（6-7）对f (x)=1精确成立，即： 由此可得： 下面来分析f (xk) 的误差对数值求积结果的影响。 设f (xk)有误差?k，并设 : 则由此引起的计算误差为： 关于收敛性可