《计算方法》5 数值积分（国家级精品课程）.pptVIP

上式称为梯形求积公式,也称两点公式，记为 梯形公式的余项为 积分第二中值定理 梯形(trapezia)公式具有 次代数精度。 故 均差性质 1 设在区间[a, b]上函数f (x)连续，而函数?(x)可积且不变号，则在开区间(a, b)内至少存在一点ξ，使 积分第二中值定理： 均差性质： 2. Simpson公式及其余项 Cotes系数为 求积公式为 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有 次代数精度。 3 3. Cotes公式及其余项 Cotes系数为 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度。 第3节 复化求积公式 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复化求积公式。 ? 复化梯形公式： 在每个 上用梯形公式： = Tn 由介值定理知： 使 即有： 余项： ? 复化 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 例1: 分别利用复化梯形公式和复化Simpson公式计算积分： 积分的相对精确值为 解：设 =0 步长 h=1/8。 =0运算量基本相同 ?复化求积法的余项和收敛阶： 复化梯形（ Trapezoid ）公式的余项： 复化辛甫生（Simpson）公式的余项： 复化柯特斯（Cotes）公式的余项： * 计 算 方 法 数学与统计学院 第五章 数值积分 计算方法课程组 §5 数值积分 §5.1 机械求积公式 §5.2 Newton_Cotes公式 §5.3 变步长求积公式及其加速 收敛技巧 §5.4 Gauss公式 §5.1 机械求积公式 第1节 引言 第2节 数值积分的基本方法 第3节 代数精度法 第4节 插值求积法 §5.1.1 引言 定积分的计算可用著名的牛顿-莱布尼兹公式来计算: 其中 F(x) 是 f (x) 的原函数之一，可用不定积分求得. 被积函数 f (x) 是用函数表格提供; f(x) 极为复杂，求不出原函数; 大量函数的原函数不容易或根本无法求出. 只能运用数值积分, 求积分近似值 . 问题 其中, 称为积分节点， 称为求积系数。 §5.1.2 数值积分的基本方法 就是在区间 [a , b] 内取 n+1 个点 利用被积函数 f (x) 在这 n+1 个点的函数值的某一种线性组合来近似作为待求的定积分. 其中, 称为积分节点， 称为求积系数。 §5.1.2 数值积分的基本方法 因此，数值积分公式关键在于积分节点 的选取 和积分系数 的决定，其中 与被积函数 f (x) 无关。 称为机械求积公式。 求积分 简单算例 求积分 此积分的几何意义相当于如图所示的曲边梯形的面积。 用 f (x) 的零次多项式 来近似代替 于是有 简单算例 解 左矩公式 推广: 右矩公式 中矩公式 用 f(x) 的一次多项式 来近似代替 , 于是， 梯形公式 推广: 来近似代替 ,于是， 特别地：当 有 用 f (x) 的二次插值多项式, 推广: Simpson公式 §5.1.3 代数精度法 为了使一个求积公式能对更多的积分具有较好的实际 计算意义,就要求它对尽可能多的被积函数都准确地成立. 因此定义代数精度的概念: 若积分 的数值积分公式 对于任意 多项式都精确成立， 但对 不精确成立， 则称该数值积分公式具m 次代数精确度。 对于[a, b]上线性插值，如图所示有 考察其代数精度。 梯形公式 算例: f(x)