《计算方法》第五章 数值积分.ppt

上式称为梯形求积公式,也称两点公式，记为 梯形公式的余项为 积分第二中值定理 梯形 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 a b f b f a 的代数精度。 解：逐次检查公式是否精确成立 代入 L0 1： 代入 L1 x ： 代入 L2 x2 ： ? 得代数精度 1 梯形公式 f x a b f b f a 试确定下面积分公式中的参数使其代数精确度尽量高. 算例: 试确定下面积分公式中的参数使其代数精确度尽量高. 解: 算例: 令 因此 由此得该积分公式具有 3 次代数精确度. 类似地, 可以证明矩形公式具 0 次代数精度 可以证明 Simpson公式具 3 次代数精度. 利用插值多项式来构造数值求积公式,具体步骤如下: 近似计算 §5.1.2 插值求积法 利用插值多项式 , 则定积分容易计算。 在 [a, b] 上取 a ? x0 x1 … xn ? b，做 f 的 n 次插值多 项式 ， 即得到 利用插值多项式来构造数值求积公式