Newton-Cotes求积公式（一） | 一般形式 +公式稳定性 +截断误差与代数精度

1. Newton-Cotes公式的一般形式

对于$I={\int }_a^{b}f(x)dx$，将区间$[a,b]$分成n等份，步长为：
$$h=\frac{(b-a)}{n}$$
故求积节点$x_k=a+kh(k=0,1,2,\cdots ,n)$，通过这$n+1$个节点构造一个n次代数多项式$p_n(x)\cong f(x)$。令
$$x=a+th(t=0,1,2,\cdots ,k,\cdots ,n)$$
作积分变换
$$dx=hdt$$

$$x=a,t=0;x=b,t=n$$

由$A_k={\int }_a^b{l}_k(x)dx$计算求积系数$A_k(k=0,1,2,\cdots ,n)$，即：
$$\begin{gathered} A_k={\int }_a^b{l}_k(x)dx={\int }_a^b\prod _{j=0,j\ne k}^n\frac{x-x_j}{x_k-x_j}dx={\int }_0^n\prod _{j=0,j\ne k}^n\frac{(a+th)-(a+jh)}{(a+kh)-(a+jh)}hdt \\ ={\int }_0^n\prod _{j=0,j\ne k}^n\frac{t-j}{k-j}hdt=h{\int }_0^n\prod _{j=0,j\ne k}^n\frac{t-j}{k-j}dt \end{gathered}$$
因为
$$\begin{gathered} \prod _{j=0,j\ne k}^n(k-j)=k(k-1)\cdots (k-(k-1))\cdot (k-(k+1))\cdots (k-n) \\ =k(k-1)\cdots 1\cdot (-1)(-2)\cdots (-(n-k))=k!(-1{)}^{n-k}(n-k)! \end{gathered}$$
所以
$$A_k=\frac{h(-1{)}^{n-k}}{k!(n-k)!}{\int }_0^n\prod _{j=0,j\ne k}^n(t-j)dt$$
令
$$C_k=\frac{A_k}{b-a}=\frac{A_k}{nh}=\frac{(-1{)}^{n-k}}{nk!(n-k)!}{\int }_0^n\prod _{j=0,j\ne k}^n(t-j)dt\text{\hspace{1em}(1)}$$
称$C_k$为Cotes（柯特斯）系数。故
$${\int }_a^{b}f(x)dx\cong (b-a)\sum _{k=0}^n{C}_{k}f(x_k)\text{\hspace{1em}(2)}$$
（1）式即为n阶的Newton-Cotes公式的一般形式。

如果给定n，则$C_k$可由（1）式计算出来。下标给出了$n=1\sim 10$的Cotes系数。应该注意

1）表中的系数$c_k$要乘倍率K，才是公式（1）和（2）中的$C_k$。例如，当n=3时，$C_0=\frac{1}{8},C_1=\frac{3}{8},C_2=\frac{3}{8},C_3=\frac{1}{8}$。

2）当n大于等于6时，表中未列出的系数可根据对称性得到。例如，当n=6时，$C_6=\frac{41}{840}$。

根据上述规则，可以写出$n=4$的Newton-Cotes公式如下：
$${\int }_a^{b}f(x)dx\cong \frac{1}{90}(b-a)[7f(x_0)+32f(x_1)+12f(x_2)+32f(x_3)+7f(x_4)]$$
其中，$x_0=a,x_4=b$.

2. Newton-Cotes公式的稳定性

在Newton-Cotes公式（2）中，设计算函数值$f(x_k)(k=1,2,\cdots ,n)$所产生的舍入误差为${ϵ}_k(k=1,2,\cdots ,n)$，则用Newton-Cotes公式计算积分的误差为：
$$\eta =(b-a)\sum _{k=0}^n{C}_{k}f(x_k)-(b-a)\sum _{k=0}^n{C}_k[f(x_k)+{ϵ}_k]=(b-a)\sum _{j=0}^n{C}_k{ϵ}_k$$
令$ϵ=ma{x}_{0\le k\le n}|{ϵ}_k|$，当Cotes系数全为正时$(n\le 7,n=9)$，则
$$|\eta |\le (b-a)\sum _{k=0}^n|C_k{ϵ}_k|\le (b-a)ϵ\sum _{k=0}^n|C_k|=(b-a)ϵ$$
所以，Newton-Cotes公式是数值稳定的。如果Cotes系数有正有负，则${\sum }_{k=0}^n|C_k|>1$

随着n的增大，该值变得越来越大，并且对$|\eta |$的影响越来越大，容易出现数值不稳定现象。因此，高阶$(n\ge 8)$的Newton-Cotes公式在实际中很少应用。

3. 截断误差与代数精度

1. 截断误差

Newton-Cotes公式是插值型求积公式，故其余项R是Lagrange插值多项式余项的积分：
$$R={\int }_a^{b}f(x)dx-\sum _{i=0}^n{A}_{i}f(x_i)={\int }_a^b\frac{f^{(n+1)}(\xi )}{(n+1)!}\prod _{i=0}^n(x-x_i)dx,\xi \in [a,b]$$
即
$$R={\int }_a^b\frac{f^{(n+1)}(\xi )}{(n+1)!}\prod _{i=0}^n(x-x_i)dx,\xi \in [a,b]$$
做积分变换$x=a+th,dx=hdt,x=a,t=0;x=b,t=n$，有：
$$R={\int }_0^a\frac{f^{(n+1)}(\xi )}{(n+1)!}(a+th-a-0h)(a+th-a-1h)\cdots (a+th-a-nh)hdt$$
故
$$R=\frac{h^{n+2}}{(n+1)!}{\int }_0^n{f}^{(n+1)}(\xi )\prod _{k=0}^n(t-k)dt,\xi \in [a,b]$$
在推导几个低阶Newton-Cotes公式的截断误差时，会用到广义积分中值定理。即：

定理1：设函数$g(x)$即$f(x)g(x)$在[a,b]上可积，函数$f(x)$在$[a,b]$上连续，且$m\le f(x)\le M$（m，M为常数），函数$g(x)$在$[a,b]$上不变号，即始终有$g(x)\le 0$或$g(x)\ge 0$，则
$${\int }_a^{b}f(x)g(x)dx=f(\xi ){\int }_a^{b}g(x)dx,\xi \in [a,b]$$

2. 代数精度

定理2：具有n+1个求积节点的Newton-Cotes公式在n为偶数时，其代数精度为n+1次。