72第二节 牛顿—柯特斯公式.pptVIP

数学学院 信息与计算科学系 第二节 牛顿—柯特斯公式 一、牛顿—柯特斯公式 设将积分区间[a, b]划分成 n 等分,步长 h = 求积节点取为xk = a + kh (k = 0,1,…,n),由此构造插值型求积公式, 则其求积系数为 (k=0,1, …, n) 引入变换 x=a+th, 则有 dx=hdt 及系数 (k=0,1,…,n) 记 (k=0,1,…,n) 则有 于是得求积公式 称为n 阶牛顿－柯特斯公式, 其中 称为柯特斯系数. 显然, 柯特斯系数与被积函数 f (x) 和积分区间[a, b]无关, 且为多项式积分. (k=0,1,…,n) 柯特斯系数表(书p176有n=8的表) n=2时公式称为辛卜生(Simpson)公式, 即 n=1时的柯特斯系数为 又称为抛物线公式. 相应的牛顿-柯特斯公式公式为梯形公式, 即 也称为柯特斯公式. 其中，xk=a+kh, (k=0,1,2,3,4), n=4 时的牛顿-柯特斯公式为 定理7.3 n阶牛顿-柯特斯公式的代数精度至少为 当n为偶数 当n为奇数 证 只要验证当n为偶数时，牛顿-柯特斯公式对f(x)=xn+1的余项为0即可 由于 ，从而有余项 这里用到xj=a+jh和变换 x=a+th及其 dx=hdt. 因为被积函数 令n=2k，k为正整数，再令t=u+k，则有 是奇函数，故在对称区间上的积分为0，即Rn(f)=0. 这就证明了n阶牛顿-柯特斯公式在n为偶数的时候代数精度至少为n +1，从而定理得证. 所以抛物线公式的代数精度是3. 抛物线公式(Simpson 公式)是n=2 时的牛顿-柯特斯公式，故其代数精度至少为3，但由于 (2) 若 f(x)?C4[a,b], 则Simpson 公式余项为 　　定理7.4 (1) 若 f(x)?C2[a, b], 梯形公式余项为 (3) 若 f(x)?C6[a,b], 则柯特斯公式余项为 二、误差估计 这里余项就是截断误差. 被积函数中的因子(x-a)(x-b)在区间[a, b]上不变号, 由积分中值定理, 在[a, b]内至少存在一点?, 使 证 (1) (2) 进一步还可以证明辛卜生公式的余项为 (3) 柯特斯公式的余项为 解 由梯形公式得 由辛卜生公式得 例1 分别用梯形公式、辛卜生公式和柯特斯公式计算积分 由柯特斯公式得 积分的精确值为 解 ①由梯形公式得 故误差估计为 例2 分别用梯形公式和抛物线公式计算积分，并估计误差. 由于 　　②由抛物线(辛卜生)公式得 故误差估计为 由于 积分的精确值为 因为牛顿-柯特斯公式对 f (x) = 1精确成立, 即 由此可得 设 f (xk) 有误差 ?k , 则计算误差为 三、牛顿-柯特斯公式的稳定性和收敛性 只要f (xk) 取得足够精确, 初始数据的误差对计算结果影响不大, 方法是稳定的. 当 全为正时 , 从而 当 有正有负时 , 因为 而 可能会很大, f (xk) 可以取得足够精确,但初始数据的误差对计算结果影响会很大, 方法可能是不稳定的. 数学学院 信息与计算科学系