复化科特斯公式matlab_基于牛顿—科特斯积分的误差分析

[摘 要]本文首先简单介绍了牛顿科特斯公式的基本思想及算法步骤，然后讨论了牛顿科特斯公式的阶数及复化的子区间个数对误差的影响，得出布尔求积公式(即阶数取4)，复化子区间个数取90左右较为理想。

[关键词]积分 复化牛顿-科特斯积分 误差

[中图分类号] G642 [文献标识码] A [文章编号] 2095-3437(2013)18-0152-02

一、引言

我们知道利用插值多项式来构造数值求积公式是最常用的一种方法，为了便于计算与应用，常将积分区间n等分，其中的每个节点作为求积节点，这样构造出来的插值型求积公式就称为牛顿-科特斯(Newton-Cotes)求积公式，这里的n称为牛顿-科特斯公式的阶数。当n=1时，该公式即为梯形求积公式；当n=2时，为辛普森求积公式；当n=3时，为3/8辛普森求积公式；当n=4时，为布尔求积公式。

由文[1]我们知道，当 n≤7 时，牛顿-布尔公式是稳定的。而当 n≥8 时，出现负数，稳定性得不到保证。而且当n较大时，由于Runge现象，收敛性也无法保证。[2]故一般不采用高阶的牛顿-科特斯求积公式。为了提高精度我们通常把积分区间分成若干子区间，然后在每个子区间上应用低阶牛顿-科特斯求积公式求积分，即为复化求积法。[3]

本文借助Matlab[4][5]符号计算系统，首先讨论不同的方法(即阶数的不同)对积分的精度与速度的影响，其次，讨论复化的子区间段数对积分误差的影响。

二、复化的牛顿科特斯求积算法实现

在积分区间[a，b]上取n+1个等距节点xk=a+kn(k=0，1