在数值分析 上,梯形法则 和辛卜生法则 均是数值积分 的方法。它们都是计算定积分 的。
这两种方法都属于牛顿-寇次公式 。它们以函数于等距
n
+
1
{\displaystyle n+1}
点的值,取得一个
n
{\displaystyle n}
次的多项式 来近似原来的函数,再行求积。
原函数(蓝色)近似为红色的线性函数
多重梯形法则 梯形法则 是:
∫
a
b
f
(
x
)
d
x
≈
(
b
−
a
)
f
(
a
)
+
f
(
b
)
2
.
{\displaystyle \int _{a}^{b}f(x)\,dx\approx (b-a){\frac {f(a)+f(b)}{2}}.}
这等同将被积函数近似为直线 函数,被积的部分近似为梯形 。
要求得较准确的数值,可以将要求积的区间分成多个小区间,再个别估计,即:
∫
a
b
f
(
x
)
d
x
≈
b
−
a
n
(
f
(
a
)
+
f
(
b
)
2
+
∑
k
=
1
n
−
1
f
(
a
+
k
b
−
a
n
)
)
.
{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\frac {b-a}{n}}\left({f(a)+f(b) \over 2}+\sum _{k=1}^{n-1}f\left(a+k{\frac {b-a}{n}}\right)\right).}
可改写成
∫
a
b
f
(
x
)
d
x
≈
b
−
a
2
n
(
f
(
x
0
)
+
2
f
(
x
1
)
+
2
f
(
x
2
)
+
⋯
+
2
f
(
x
n
−
1
)
+
f
(
x
n
)
)
{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\frac {b-a}{2n}}\left(f(x_{0})+2f(x_{1})+2f(x_{2})+\cdots +2f(x_{n-1})+f(x_{n})\right)}
其中
对
k
=
0
,
1
,
…
,
n
{\displaystyle k=0,1,\dots ,n}
,
x
k
=
a
+
k
b
−
a
n
,
{\displaystyle x_{k}=a+k{\frac {b-a}{n}},}
。
辛普森法则 (Simpson's rule,又称森逊法则 )是:
∫
a
b
f
(
x
)
d
x
≈
b
−
a
6
[
f
(
a
)
+
4
f
(
a
+
b
2
)
+
f
(
b
)
]
.
{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\frac {b-a}{6}}\left[f(a)+4f\left({\frac {a+b}{2}}\right)+f(b)\right].}
同样地,辛普森法则也有多重的版本:
∫
a
b
f
(
x
)
d
x
≈
h
3
⋅
[
f
(
x
0
)
+
2
∑
k
=
1
n
−
1
f
(
x
k
)
+
4
∑
k
=
1
n
f
(
x
k
−
1
+
x
k
2
)
+
f
(
x
n
)
]
{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\frac {h}{3}}\cdot \left[f(x_{0})+2\sum _{k=1}^{n-1}f(x_{k})+4\sum _{k=1}^{n}f\left({\frac {x_{k-1}+x_{k}}{2}}\right)+f(x_{n})\right]}
h
=
b
−
a
n
,
x
k
=
a
+
k
⋅
h
.
{\displaystyle h={\frac {b-a}{n}},\ x_{k}=a+k\cdot h.}
或写成
∫
a
b
f
(
x
)
d
x
≈
h
3
[
f
(
x
0
)
+
4
f
(
x
1
)
+
2
f
(
x
2
)
+
4
f
(
x
3
)
+
2
f
(
x
4
)
+
⋯
+
4
f
(
x
n
−
1
)
+
f
(
x
n
)
]
{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\frac {h}{3}}{\bigg [}f(x_{0})+4f(x_{1})+2f(x_{2})+4f(x_{3})+2f(x_{4})+\cdots +4f(x_{n-1})+f(x_{n}){\bigg ]}}
牛顿-寇次公式(Newton-Cotes rule / Newton-Cotes formula)以Roger Cotes和艾萨克·牛顿 命名。其内容是:
∫
a
b
f
(
x
)
d
x
≈
∑
i
=
0
n
w
i
f
(
x
i
)
{\displaystyle \int _{a}^{b}f(x)\,dx\approx \sum _{i=0}^{n}w_{i}\,f(x_{i})}
其中对
k
=
0
,
1
,
…
,
n
{\displaystyle k=0,1,\dots ,n}
,
w
i
{\displaystyle w_{i}}
是常数(由
n
{\displaystyle n}
的值决定),
x
k
=
a
+
k
b
−
a
n
{\displaystyle x_{k}=a+k{\frac {b-a}{n}}}
。
梯形法则和辛卜生法则便是
n
=
1
,
2
{\displaystyle n=1,2}
的情况。
亦有不采用在边界点来估计的版本,即取
x
k
=
a
+
k
b
−
a
n
+
1
{\displaystyle x_{k}=a+k{\frac {b-a}{n+1}}}
。
对于次数较高的多项式而有很大误差(龙格现象 ),不如高斯积分法 。
M. Abramowitz and I. A. Stegun, eds. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables . New York: Dover, 1972. (See Section 25.4.)
George E. Forsythe, Michael A. Malcolm, and Cleve B. Moler. Computer Methods for Mathematical Computations . Englewood Cliffs, NJ: Prentice-Hall, 1977. (See Section 5.1.)
William H. Press, Brian P. Flannery, Saul A. Teukolsky, William T. Vetterling. Numerical Recipes in C . Cambridge, UK: Cambridge University Press, 1988. (See Section 4.1.)
Josef Stoer and Roland Bulirsch. Introduction to Numerical Analysis . New York: Springer-Verlag, 1980. (See Section 3.1.)