台形公式

数学において、台形公式（だいけいこうしき、Trapezoidal rule）もしくは台形則（だいけいそく）は定積分を近似計算するための方法の 1つである。

具体的に言えば、求めたいx-yグラフのy=0を含む面積内に無数の台形を置くと、その台形の面積の集合和は本物の面積に限りなく近い値となる。

概要

台形公式はニュートン・コーツの公式と呼ばれる数値積分の公式のひとつである。被積分関数を一次関数で近似し、台形の面積の公式に帰着させて積分の近似値を求める方法である。一次関数による近似なので精度はそれほど期待できず、二次関数で近似するシンプソンの公式などの方が精度が高い。

シンプソンの公式やその他の類似の手法は、2階連続微分可能な関数に対する台形公式の改良とみなせるが、細かく変動しない荒い関数に対しては台形公式で十分であり、計算も簡単である。

近似誤差

台形公式の利点は、近似誤差が容易に分かることである。

凸関数に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ小さい値になり、 凹関数に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ大きい値になる。

また積分区間が変曲点を含むとき上記の凸部分の誤差と凹部分の誤差が打ち消し合い全体的な誤差は小さくなる。

さらに、台形公式は周期関数をその周期よりも長い区間積分する場合にはきわめて精度が高くなる傾向がある。これはオイラーの和公式との関係をみると良く理解できる。しかしながら非周期関数に対しては一般に、 ガウス求積やクレーンショー・カーチス数値積分則のような非等分点法の方がより精度が高い。

台形による近似

関数f(x) (青線)を線形関数で近似(赤線)。

合成台形公式の図(不等間隔の例).

x ∈(a,b) のとき f(x) > 0 とする。 定積分

$\int _{a}^{b}f(x)\,dx$

の値は、xy座標平面上で y = f(x) と x 軸、 x = a、 x = b とで囲まれた図形の面積になる。一般に積分計算は複雑であることも多いが、近似値が欲しいだけであれば f(x) を一次関数で近似することによりこの図形は台形とみなすことができるので

$\int _{a}^{b}f(x)\,dx\approx (b-a){\frac {f(a)+f(b)}{2}}$

のように、台形の面積公式を用いて簡単に計算できる。

この公式は y = f(x) が直線に近いほど精度がよくなるが、直線とは全く異なる曲線であれば精度が悪くなり欲しい近似値を得ることが難しくなってくる。そこで、この積分をより正確に計算するために、

a = a₀ ≤ a₁ ≤ … ≤ a_n = b

として積分区間 [a, b] をより小さいn個の部分区間 [a₀,a₁], [a₁,a₂], … [a_n−1,a_n] に分け、それぞれの区間で台形公式を用いて近似する。

$\int _{a}^{b}f(x)\,dx$

$=\int _{a_{0}}^{a_{1}}f(x)\,dx+\int _{a_{1}}^{a_{2}}f(x)\,dx+\int _{a_{2}}^{a_{3}}f(x)\,dx$ ・・・ $\int _{a_{n-1}}^{a_{n}}f(x)\,dx$

$=\sum _{k=1}^{n}\int _{a_{k-1}}^{a_{k}}f(x)\,dx\approx \sum _{k=1}^{n}(a_{k}-a_{k-1}){\frac {f(a_{k-1})+f(a_{k})}{2}}$

さらに1区間の幅が　$a_{n}-a_{n-1}={\frac {b-a}{n}}$　と、一定になり

k番目は　$a_{k}=a+k{\frac {b-a}{n}},(k=0\ldots n)$　と置けることで

$\int _{a}^{b}f(x)\,dx\approx \sum _{k=1}^{n}(b-a){\frac {f(a_{k-1})+f(a_{k})}{2n}}$

${\frac {b-a}{2n}}$　は定数より

$={\frac {b-a}{2n}}\sum _{k=1}^{n}{f(a_{k-1})+f(a_{k})}$

$={\frac {b-a}{2n}}\left(f(a_{0})+2f(a_{1})+2f(a_{2})+\cdots +2f(a_{n-1})+f(a_{n})\right)$

$={\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\}$

という公式が得られる。

関連項目

• 数値解析

参考文献

• Burden, Richard L.; J. Douglas Faires (2000). Numerical Analysis (7th Ed. ed.). Brooks/Cole. ISBN 0-534-38216-9 
• Trapezoidal Rule of Integration - Notes, PPT, Mathcad, Maple, Mathematica, Matlab at Holistic Numerical Methods Institute

外部リンク

• Trapezoidal Rule for Numerical Integration （英語）