ch02b 牛顿-柯特斯公式.pptVIP

第二章 数值积分 牛顿-柯特斯(Newton-Cotes)公式 及其复化求积公式 牛顿-柯特斯公式 等距节点的插值型求积公式称为牛顿-柯特斯公式： 令 x = a + t h 得： 插值型求积公式 其中 牛顿-柯特斯公式（续） 注：Cotes 系数仅取决于 n 和 i，可通过查表得到。与被积函数 f (x) 及积分区间 [a, b] 均无关。 牛顿-柯特斯公式： 几个常见公式 代数精度 = 1 梯形求积公式 代数精度 = 3 抛物线求积公式 辛甫生(Simpson)求积公式 n = 4: 柯特斯(Cotes)求积公式 柯特斯系数表 系数特点和稳定性 柯特斯系数具有以下特点： (3) 当 n ? 8 时，出现负数，稳定性得不到保证。而且当 n 较大时，由于Runge现象，收敛性也无法保证。 故一般不采用高阶的牛顿-柯特斯求积公式。 当 n ? 7 时，牛顿-柯特斯公式是稳定的。 牛顿-柯特斯公式的代数精度 定理 当等距节点的个数为奇数n+1(n 为偶数)时， 牛顿－柯特斯公式至少有 n+1 次代数精度。 证：只要证明当 n 为偶数时，公式对f (x)＝xn+1精确成立。 数值积分的余项 积分中值定理 若f(x), g(x)均在[a,b]上连续, 且g(x)在[a,b]不变号, 则存在??[a,b] 使 插值型求积公式的余项： 左矩形公式 而 x-a 在[a,b]上不变号（非负），由积分中值定理知 证明：设f (x)在[a,b]上连续， 梯形公式 梯形公式的余项证明： 中值定理 Simpson公式的余项证明： 三次Hermite插值 Simpson公式余项(证明: 用积分中值定理+Hermite插值) Simpson公式 Newton-Cotes余项的一般形式 定理 (1) 若 n 为偶数， f (x) ?Cn+2[a, b] ，则存在 ? ?(a, b) 使得 (2) 若 n 为奇数， f (x) ?Cn+1[a, b] ，则存在 ? ?(a, b) 使得 举例（一） 与精确值 0.6321...... 相比得误差分别为 0.0518 和 0.0002。 复化求积公式 提高积分计算精度的常用两种方法 用 复化公式 用 非等距节点 复化求积公式：将积分区间分割成多个小区间，然后在每个小区间上使用低次牛顿－柯特斯求积公式。 将[a, b] 分成 n 等分 [xi , xi+1] ，其中节点 复化梯形公式 复化梯形公式: 余项： 复化辛甫生(simpson)公式 复化simpson公式: 余项： 复化柯特斯公式 复化cotes公式: 余项： 举例（二） 解： xi 0 1/8 2/8 3/8 4/8 5/8 6/8 7/8 1.0 f (xi ) 1 0.997 0.990 0.977 0.954 0.936 0.909 0.877 0.841 比较精确值0.9460831... 定义：若一个复化积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 收敛速度与误差估计： ?i ?(xi, xi+1 ) h 很小时的误差 同理 问题: 给定精度 ?，如何取 n ? 即：取 n = 409