数值积分.docVIP

数值积分 §5.0 引言 §5.1 机械求积公式 §5.2 Newton-Cotes公式 §5.3 变步长求积公式及其加速收敛技巧 §5.4 Gauss公式 0 引 言 定积分的计算可用著名的牛顿-莱布尼兹公式来计算: 其中F(x)是f(x)的原函数之一，可用不定积分求得。 然而在实际问题中，往往碰到以下问题： 被积函数f(x)是用函数表格提供的； 被积函数表达式极为复杂，求不出原函数，或求出原函数后，由于形式复杂不利于计算； 大量函数的原函数不容易或根本无法求出，例如 概率积分，正弦型积分 (回路磁场强度公式) 等根本无法用初等函数来表示其原函数，因而也就无法精确计算其定积分，只能运用数值积分。 2 所谓数值积分就是求积分近似值的方法。 §1 机械求积公式 1 数值积分的基本思想 区间[a,b]上的定积分，就是在区间[a,b]内取n+1个点，利用被积函数f(x)在这n+1个点的函数值的某一种线性组合来近似作为待求定积分的值，即 右端公式称为左边定积分的某个数值积分公式。 其中,xk称为积分节点，Ak称为求积系数。 因此，一个数值积分公式关键在于积分节点xk的选取和积分系数Ak的决定，其中Ak与被积函数f(x)无关。称为机械求积公式。 1.1 简单算例说明 例1 求积分 此积分的几何意义相当于如下图所示的曲边梯形的面积。 解：（1） 用f(x)的零次多项式 来近似代替,于是， （为左矩公式） 推广: （为右矩公式） （为中矩公式） （2） 用f(x)的一次多项式 来近似代替,于是， （为梯形公式） （3） 用f(x)的二次插值多项式,其中 来近似代替,于是， 特别地：当时，有 （为Simpson公式） 2 代数精确度 定义：若积分的数值积分公式对于任意一个次数不高于m次的多项式都精确成立，且存在一个m+1次多项式使之不精确成立，则称该数值积分公式的代数精确度为m。 对于代数精确度为m的求积公式，若f(x)是不超过m次的代数多项式，求积公式是精确成立的。 2.1 算例 例1: 有积分公式: 求该积分公式的代数精确度。 这个求积公式的几何意义是曲边梯形的面积近似地用两个梯形面积来代替。 解：(1)取f(x)=1，定积分, 而数值积分，两端相等； (2)取f(x)=x，定积分, 而数值积分，两端相等； (3)取，定积分, 而数值积分，两端不相等； 只要取f(x)=1，f(x)=x验证了上述求积公式精确成立，就意味着对于任意一个一次多项式，求积公式都是精确成立的；而取时求积公式不精确成立，也就是存在一个二次多项式使求积公式不精确成立；故该求积公式的代数精确度为1。 例2：在如下求积公式中，求积分节点和相应的求积系数使其代数精确度尽可能高。 解：(1) f(x)=1, ,而数值积分为； 得到方程； (2) f(x)=x，，而数值积分为； 得到方程； (3) ，，而数值积分为；得到方程； (4) ，，而数值积分为；得到方程； 综合上述方程： 解得: 。 于是我们得到积分公式 。 再取，有， 而数值积分为，两式不相等，求积公式不精确成立了。所以，该积分公式的代数精确度为3。 §2 Newton-Cotes公式 1 公式的推导 Newton-Cotes公式是由拉格朗日插值公式推导出来的数值积分公式。 将区间[a,b]等分n等份,记,分点为, k=0,1,...,n,这n+1个节点上的函数值为, 从而区间[a,b]上的拉格朗日插值多项式为 其中, 为插值基本多项式,与函数f(x)无关,k=0,1,...,n。 由于插值结点是等距节点，故插值多项式可以进一步化简: 因为 ， 故 ， 因 ， 作积分变量代换，， 当时，t=0；当x=b时，t=n； 故 记，我们称为柯特斯（Cotes）系数，它不仅与函数f(x)无关，而且与积 分区间[a,b]无关。 例如： 当n=1时 (梯形积分公式中的系数) , ; 当 n=2时 (抛物线积分公式中的系数) , , ; 当n=3时 (3/8积分公式中的系数) ， ， ， ； 于是，由柯特斯（Cotes）系数公式出发，我们得到n阶Newton—Cotes公式： 。 2 低阶公式及其余项 常用的Newton—Cotes公式 a) 梯形公式 n=1时，积分节点为，，则数值积分公式为：