科学与计算工程办法第5章.pdfVIP

第五章 数值积分 §5.0 引言 §5.1 机械求积公式 §5.2 Newton-Cotes公式 §5.3 变步长求积公式及其加速收敛技巧 §5.4 Gauss公式 §5.5 小结 1 §5.0 引 言 1. 定积分的计算可用著名的牛顿-莱布尼兹公式来计算: b ∫a f (x )dx F (b )=−F (a ) 其中 F(x)是f (x) 的原函数之一，可用不定积分求得。 然而在实际问题中，往往碰到以下问题： (a) 被积函数f (x)是用函数表格提供的； (b) 被积函数表达式极为复杂，求不出原函数，或求出原函数后，由于 形式复杂不利于计算； (c) 大量函数的原函数不容易或根本无法求出，例如 1 2 ∫0 e−x dx ，概率积分 1 sinx ∫0 dx , 正弦型积分 x π 2 4Ir x 2 ⎛ ⎞ 2 H (x ) 2 2 ∫0 1 =−⎜ ⎟ sin θ dθ 回路磁场强度公式 r −x r ⎝ ⎠ 等根本无法用初等函数来表示其原函数，因而也就无法精确计算其 定积分，只能运用数值积分。 2 所谓数值积分就是求积分近似值的方法。 而数值积分只需计算 f (x) 在节点 xi (i 1,2, ,n) 上的值，计算方便 且适合于在计算机上机械地实现。 2 §5.1 机械求积公式 1 数值积分的基本思想 b 区间[a,b] 上的定积分 ∫a f (x )dx ，就是在区间[a,b] 内取 n+ 1 个点 x ,x ,,xn ，利用被积函数f (x)在这 n+ 1 个点的函数值的某一种线性组合 0 1 来近似作为待求定积分的值，即 b n f (x )dx ≈ A f (x ) ∫a ∑ k k k 0 右端公式称为左边定积分的某个数值积分公式。 其中,x 称为积分节点，A 称为求积系数。 k k 因此，一个数值积分公式关键在于积分节点x 的选取和积分系数A 的决定， k k 其中Ak 与被积函数f(x) 无关。称为机械求积公式。 1.1 简单算例说明 x 1 例 1 求积分 ∫x f (x )dx 0 此积分的几何意义相当于如下图所示的曲边梯形的面积。 解：（1）用f (x) 的零次多项式 y L 0 (x ) f (x 0 ) 来近似代替 f (x ) ,于是， x x 1