数值积分 计算方法课件及实验 培训讲解.ppt

2. Simpson公式及其余项(n=2) Cotes系数为 求积公式为 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有 次代数精度。 3 3. Cotes公式及其余项(n=4) Cotes系数为 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度。 注：n?8时，Cotes系数出现负数，会引起误差增大，计算不稳定。 因此，在实际应用中一般不使用高阶Newton-Cotes公式，而是采用低阶复合求积法(下节)。 Cotes系数表： n Ck(n) 1 2 3 4 5 8 1/2 1/2 1/6 4/6 1/6 1/8 3/8 3/8 1/8 7/90 16/45 2/15 16/45 7/90 19/288 25/96 25/144 25/144 25/96 19/288 ………… 989/28350 5888/28350 -928/28350 10496/28350 -4540/28350 … 第三节 复合求积公式 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 ? 复合梯形公式： 在每个 上用梯形公式： = Tn 由介值定理知： 使 即有： 余项： ? 复合 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 例1: 分别利用复合梯形公式和复合Simpson公式计算积分： 积分的相对精确值为 解：设 =0 步长 h=1/8。 =0运算量基本相同 ?复合求积法的余项和收敛阶： 复合梯形（ Trapezoid ）公式的余项： 复合辛甫生（Simpson）公式的余项： 复合柯特斯（Cotes）公式的余项： * 计 算 方 法 第四章 数值积分 孙成立 §第四章 数值积分 §4.1 机械求积公式 §4.2 Newton_Cotes公式 §4.3 变步长求积公式及其加速 收敛技巧 §4.1 机械求积公式 第1节 引言 第2节 数值积分的基本方法 第3节 代数精度法 第4节 插值求积法 §4.1.1 引言 定积分的计算可用著名的牛顿-莱布尼兹公式来计算: 其中 F(x) 是 f (x) 的原函数之一，可用不定积分求得. 被积函数 f (x) 是用函数表格提供; f(x) 极为复杂，求不出原函数; 大量函数的原函数不容易或根本无法求出. 只能运用数值积分, 求积分近似值 . 问题 其中, 称为积分节点， 称为求积系数。 §4.1.2 数值积分的基本方法 就是在区间 [a , b] 内取 n+1 个点 利用被积函数 f (x) 在这 n+1 个点的函数值的某一种线性组合来近似作为待求的定积分. 求积分 此积分的几何意义相当于如图所示的曲边梯形的面积。 用 f (x) 的零次多项式 来近似代替 于是有 简单算例 解 左矩公式 推广: 右矩公式 中矩公式 用 f(x) 的一次多项式 来近似代替 , 于是， 梯形公式 推广: 来近似代替 ,于是， 特别地：当 有 用 f (x) 的二次插值多项式, 推广: Simpson公式 §4.1.3 代数精度法 为了使一个求积公式能对更多的积分具有较好的实际 计算意义,就要求它对尽可能多的被积函数都准确地成立. 因此定义代数精度的概念: 若积分 的数值积分公