線性函數(紅色)會作用估算函數
(藍色)。
梯形公式是數學中数值积分的基础公式之一:
公式由来
由积分中值定理可得
,
但由于ξ其值一般难于确定,故难以准确算出
的值。
如果用两端点
与
的算术平均值估算
,有
,
这就是梯形公式。
类似地,如果用区间中点
其高度
取代
,从而有中矩形公式
。
复合求积公式
每一區間相同
梯形公式的示意圖(長度相同的區間)。
為了計算出更加準確的定積分,可以把積分的區間
分成
份,當中
趨向無限,分割出的每一個區間長度必定要是一樣的,然後就可以應用梯形公式:
![{\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].}](./_assets_/eb734a37dd21ce173a46342d1cc64c92/f7b41bcb81b7c32debfceee32576be9b2de8e098.svg)
亦可以寫成:

當中

其余项为
當區間的長度並不相同時,這一條公式便不能使用。
每一區間並不相同
梯形公式的示意圖(長度不相同的區間)
給予
以及
,定積分就可以估算成
,
當中
.
誤差分析
應用梯形公式的誤差值是真值數字與運用梯形公式結果的差異:
![{\displaystyle {\text{error}}=\int _{a}^{b}f(x)\,dx-{\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]}](./_assets_/eb734a37dd21ce173a46342d1cc64c92/ac04c64166cbfd24aac10c03c629641a2d7c13d4.svg)
如果
中存在一個實數
,那麼

对于中矩形公式,其误差类似的有:
如果被積函數是一個凸函數(亦即有一個正值二階導數),那麼誤差會是一個負數,也代表梯形公式的估算值高估了真實數字。這可以利用一個幾何圖形代去表達:梯形不但覆蓋曲線下的面積更超越其範圍。同樣地,如果被積函數是一個凹函數,梯形公式就會低估其真實數字因為曲線下部份面積沒有被計算在內。如果被積函數中有拐點。它的錯誤是比較難去估計。
一般而言有數種方法可以去分析誤差,例如是:傅利葉級數。
在
的情況下,趨向性的估計誤差是:
![{\displaystyle {\text{error}}=-{\frac {(b-a)^{2}}{12N^{2}}}{\big [}f'(b)-f'(a){\big ]}+O(N^{-3}).}](./_assets_/eb734a37dd21ce173a46342d1cc64c92/be8c3e11b2ab8dfeff4a79a5bea6db9f35a8e2a5.svg)
参考文献
- 《数值分析》,清华大学出版社,李庆扬等编,书号ISBN 978-7-302-18565-9
|
|---|
| 牛頓-寇次公式 | |
|---|
| 高斯求积 |
- Gauss–Hermite quadrature
- Gauss–Jacobi quadrature
- Gauss–Kronrod quadrature formula
- Gauss–Laguerre quadrature
- Gauss–Legendre quadrature
- Chebyshev–Gauss quadrature
|
|---|
| 其他 |
- Barnes–Hut simulation
- Bayesian quadrature
- Clenshaw–Curtis quadrature
- Filon quadrature
- Lebedev quadrature
- Tanh-sinh quadrature
|
|---|