梯形公式
∫
a
b
f
(
x
)
d
x
≈
b
−
a
2
[
f
(
a
)
+
f
(
b
)
]
=
T
\int_a^bf(x)dx\approx \frac{b-a}{2} [f(a)+f(b)]=T
∫abf(x)dx≈2b−a[f(a)+f(b)]=T
余项:
E
T
(
f
)
=
−
(
b
−
a
)
3
12
f
′
′
(
η
)
,
η
∈
(
a
,
b
)
E_T(f)=-\frac{(b-a)^3}{12}f''(\eta),\eta \in (a,b)
ET(f)=−12(b−a)3f′′(η),η∈(a,b)
Simpson(抛物线)公式
∫
a
b
f
(
x
)
d
x
≈
b
−
a
6
[
f
(
a
)
+
4
f
(
a
+
b
2
)
+
f
(
b
)
]
=
S
\int_a^bf(x)dx\approx \frac{b-a}{6} [f(a)+4f(\frac{a+b}{2})+f(b)]=S
∫abf(x)dx≈6b−a[f(a)+4f(2a+b)+f(b)]=S
余项:
E
S
(
f
)
=
−
(
b
−
a
)
5
2880
f
(
4
)
(
η
)
,
η
∈
(
a
,
b
)
E_S(f)=-\frac{(b-a)^5}{2880}f^{(4)}(\eta),\eta \in (a,b)
ES(f)=−2880(b−a)5f(4)(η),η∈(a,b)
Cotes公式
∫ a b f ( x ) d x ≈ b − a 90 [ 7 f ( x 0 ) + 32 f ( x 1 ) + 12 f ( x 2 ) + 32 f ( x 3 ) + 7 f ( x 4 ) ] \int_a^bf(x)dx\approx \frac{b-a}{90}[7f(x_0)+32f(x_1)+12f(x_2)+32f(x_3)+7f(x_4)] ∫abf(x)dx≈90b−a[7f(x0)+32f(x1)+12f(x2)+32f(x3)+7f(x4)]
Newton-Cotes代数精度
当n为奇数时,Newton-Cotes公式具有n次代数精度
当n为偶数时,Newton-Cotes公式具有n+1次代数精度
复化梯形公式
将区间[a,b]作n等分
T
n
=
h
2
[
f
(
a
)
+
f
(
b
)
+
2
∑
k
=
1
n
−
1
f
(
a
+
k
h
)
]
T_n=\frac{h}{2}[f(a)+f(b)+2\sum_{k=1}^{n-1}f(a+kh)]
Tn=2h[f(a)+f(b)+2k=1∑n−1f(a+kh)]
h
=
b
−
a
n
h=\frac{b-a}{n}
h=nb−a
另外:
T
2
n
=
1
2
(
T
n
+
U
n
)
T_{2n}=\frac{1}{2}(T_n+U_n)
T2n=21(Tn+Un)
其中:
U
n
=
h
∑
k
=
0
n
−
1
f
(
x
k
+
1
2
)
U_n=h\sum_{k=0}^{n-1}f(x_{k+\frac{1}{2}})
Un=hk=0∑n−1f(xk+21)
复化抛物线公式
将区间[a,b]作n等分
S
n
=
h
6
[
f
(
a
)
+
4
∑
k
=
0
n
−
1
f
(
x
k
+
1
2
)
+
2
∑
k
=
1
n
−
1
f
(
x
k
)
+
f
(
b
)
]
S_n=\frac{h}{6}[f(a)+4\sum_{k=0}^{n-1}f(x_{k+\frac{1}{2}})+2\sum_{k=1}^{n-1}f(x_k)+f(b)]
Sn=6h[f(a)+4k=0∑n−1f(xk+21)+2k=1∑n−1f(xk)+f(b)]
h
=
b
−
a
n
h=\frac{b-a}{n}
h=nb−a
另外:
S
n
=
1
3
T
n
+
2
3
U
n
S_n=\frac{1}{3}T_n+\frac{2}{3}U_n
Sn=31Tn+32Un
Romberg积分法
其中,
T
0
(
k
)
T_0^{(k)}
T0(k)的计算方式与复化梯形公式一致。之后,
T
m
(
k
)
T_m^{(k)}
Tm(k)的求解遵循下述公式:
T
m
k
=
T
m
−
1
k
+
1
−
T
m
−
1
k
4
m
−
1
T_m^k=\frac{T_{m-1}^{k+1}-T_{m-1}^k}{4^m-1}
Tmk=4m−1Tm−1k+1−Tm−1k