Newton-Cotes求积公式讲课.pptVIP

4.2.3 Newton-Cotes公式的误差分析 因此,Newton-Cotes公式为 称为柯特斯系数.　柯特斯系数不但与被积函数无关，而且与积分区间也无关． Newton-Cotes公式是一类数值积分公式，它是封闭公式（区间端点也是积分节点），它是由拉格朗日插值公式推导来的． 由于是多项式积分，Cotes 系数计算不会遇到实质性困难。 2、常用的Newton-Cotes公式： 这就是抛物线公式，又称辛浦生— Simpson 公式 。几何意义 就是用抛物线下的面积近似曲线f(x)下的面积。 当n=3时，称之为Newton公式。 (4.2.12) 这就是柯特斯公式（Cotes ) 当n 较大时，例如 n=8 时，系数 中出现负数，而且有正有负会使舍入误差增大，数值稳定性较差，因此实际计算并不用 n较大的 公式，而是将区间[a,b] 分割成若干个小区间，对每个或几个小区间应用n 较小的公式去计算。Cotes系数表附后. (4.2.13) 当n = 8时，从表中可以看出出现了负系数，从而影响稳定性(Ak0)和收敛性，因此实用的只是低阶公式。 Cotes 系数仅取决于 n 和 i，可查表得到(给出了n从1～8的柯特斯系数)。 Cotes 系数与 f (x) 及区间[a, b]均无关。 例4.10 分别用梯形公式、辛卜生公式和柯特斯公式计算定积分 的近似值 (计算结果取5位有效数字) (1) 用梯形公式计算 (2) 用辛卜生公式 (3) 用柯特斯公式计算，系数为 积分的准确值为 可见，三个求积公式的精度逐渐提高。 引理：n阶Newton-Cotes公式的代数精度至少是n. 证明：如果f (x)是一个次数不超过n次的多项式，则f (n+1)(x)?0, 其拉格朗日插值公式的插值余项为： 即，Newton-Cotes公式的值精确地等于定积分的值，故n阶Newton-Cotes公式的代数精度至少是n. 故f (x)=Pn(x), 这是对一切x均相等，精确成立．所以， 结论：当n 为奇数时，n阶Newton-Cotes公式的代数精度为n； 　　　当n 为偶数时，n阶Newton-Cotes公式的代数精度为n+1(稍后再证明) 2、Cotes系数特点： 即有 　　由引理知，求积公式至少有n次代数精度,对于1,积分公式总是精确成立, (-1)n+k=(-1)n-k 1.偶阶求积公式的代数精度 作为插值型的求积公式，n 阶的牛顿-柯特斯公式至少具有n 次的插值精度（定理3.2）。从上面几个特殊的公式可以猜想, n为偶数时, 代数精度还可进一步提高， 先看Simpson公式，它是二阶Newton---Cotes公式, 因此至少具有二次代数精度。进一步用 x3 进行检验。 按Simpson公式计算得： 另一方面，直接求积分得： 易验证S=I, 即Simpson 公式对次数不超过三次的多项式均能准确成立。 §4.2 Newton-Cotes求积公式 总结 4.2.3 Newton-Cotes 公式的误差分析 4.2.2 Newton-Cotes求积公式 4.2.1 插值型求积法 数值求积法与代数精度 　　我们的目的就是根据一定原则,　选择求积节点xk和系数Ak，使得求积一般公式(4.2.1)具有较高的精确度, 同时又计算简单。 权Ak仅仅与节点xk的选取有关，而不依赖于被积函数f(x) 的具体形式。 使积分公式具有通用性 右端公式称为左端定积分的某个数值积分公式．其中xk称为积分节点, Ak为求积系数, 也称之为伴随节点xk的权． 　　 (4.2.1) 一、求积公式的代数精度 记 称(4.2.2)为数值求积公式,(4.2.3)为求积公式余项(误差). 构造或确定一个求积公式，要讨论解决的问题有 (i) 确定求积系数Ak和求积节点xk ； (ii)　求积公式的误差估计和收敛性 为了构造形如式(4.2.1)的求积公式,需要提供一种判定求积方法精度高低准则.用什么标准来判定两个节点数相同的求积公式的“好”与“差”呢？通常用“代数精确度”的高低作为求积公式“好”与“差”的一个标准．在后面的讨论中我们将看到，节点相同的求积公式，代数精确度越高，求出的积分近似值精确度一般越好．下面给出代数精确度的定义． 数值求积方法是近似方法，为要保证精度，我们自然希望求积公式能对“尽可能多”的函数准确地成立，这就提出了所谓代数精度的概念．由于闭区间[a,b]上的连续函数可用多项式逼近，所以一个求积公式能对多大次数的多项式 f (x) 成为准确等式，是衡量该公式的精确程度的重要指标，为此给出以下定义。 定义4.1