「台形公式」の版間の差分
Shunta suzuki (会話 | 投稿記録) 編集の要約なし |
編集の要約なし |
||
1行目: | 1行目: | ||
[[数学]]において、'''台形公式'''(だいけいこうしき、 |
[[数学]]において、'''台形公式'''(だいけいこうしき、{{lang-en|Trapezoidal rule}})もしくは'''台形則'''(だいけいそく)は[[積分|定積分]]を[[近似]]計算するための方法、すなわち[[数値積分]]のひとつである。これはは[[ニュートン・コーツの公式]]の1次の場合である。被積分関数を[[区分線形関数]]で近似し、[[台形]]の面積の公式に帰着させて積分の近似値を求める。 |
||
具体的に言えば、求めたいx-yグラフのy=0を含む面積内に無数の[[台形]]を置くと、その台形の面積の集合和は本物の面積に限りなく近い値となる。 |
具体的に言えば、求めたい''x'' -''y'' グラフの''y'' = 0を含む面積内に無数の[[台形]]を置くと、その台形の面積の集合和は本物の面積に限りなく近い値となる。 |
||
⚫ | |||
== 概要 == |
|||
台形公式は[[ニュートン・コーツの公式]]と呼ばれる[[数値積分]]の公式のひとつである。被積分関数を[[一次関数]]で近似し、[[台形]]の面積の公式に帰着させて積分の近似値を求める方法である。一次関数による近似なので精度はそれほど期待できず、[[二次関数]]で近似する[[シンプソンの公式]]などの方が精度が高い。 |
|||
⚫ | |||
== 近似誤差 == |
== 近似誤差 == |
||
台形公式の利点は、近似誤差が容易に分かることである。 |
台形公式の利点は、近似誤差が容易に分かることである。 |
||
[[凸関数]]に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ小さい値になり、 |
[[凸関数]]に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ小さい値になり、[[凹関数]]に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ大きい値になる。また積分区間が[[変曲点]]を含むとき上記の凸部分の誤差と凹部分の誤差が打ち消し合い''全体的な誤差''は小さくなる。 |
||
[[凹関数]]に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ大きい値になる。 |
|||
また積分区間が[[変曲点]]を含むとき上記の凸部分の誤差と凹部分の誤差が打ち消し合い'''全体的な誤差'''は小さくなる。 |
|||
さらに、台形公式は[[周期関数]]をその周期よりも長い区間積分する場合にはきわめて精度が高くなる傾向がある。これは[[オイラーの和公式]]との関係をみると良く理解できる。しかしながら非周期関数に対しては一般に、 [[ガウス求積]]や[[クレーンショー・カーチス数値積分則]]のような非等分点法の方がより精度が高い。 |
さらに、台形公式は[[周期関数]]をその周期よりも長い区間積分する場合にはきわめて精度が高くなる傾向がある。これは[[オイラーの和公式]]との関係をみると良く理解できる。しかしながら非周期関数に対しては一般に、 [[ガウス求積]]や[[クレーンショー・カーチス数値積分則]]のような非等分点法の方がより精度が高い。 |
||
{{節stub}}<!--誤差の大きさは区間幅に対して何次のオーダーか、など--> |
|||
== 台形による近似 == |
== 台形による近似 == |
||
[[Image:trapezoidal_rule_illustration.png|right|thumb|関数''f''(''x'') |
[[Image:trapezoidal_rule_illustration.png|right|thumb|関数''f'' (''x'')(青線)を線形関数(赤線)で近似]] |
||
[[Image:Composite trapezoidal rule illustration.png|right|thumb|合成台形公式の図 |
[[Image:Composite trapezoidal rule illustration.png|right|thumb|合成台形公式の図(不等間隔の例)]] |
||
⚫ | |||
⚫ | |||
⚫ | |||
⚫ | |||
定積分 |
|||
⚫ | |||
⚫ | |||
:<math> \int_{a}^{b} f(x)\, dx \approx (b-a)\frac{f(a) + f(b)}{2}</math> |
:<math> \int_{a}^{b} f(x)\, dx \approx (b-a)\frac{f(a) + f(b)}{2}</math> |
||
のように、台形の面積公式を用いて簡単に計算できる。 |
のように、台形の面積公式を用いて簡単に計算できる。 |
||
この公式は ''y'' = ''f''(''x'') が直線に近いほど精度がよくなるが、直線とは全く異なる曲線であれば精度が悪くなり欲しい近似値を得ることが難しくなってくる。そこで、この積分をより正確に計算するために、 |
この公式は曲線 ''y'' = ''f'' (''x'') が直線に近いほど精度がよくなるが、直線とは全く異なる曲線であれば精度が悪くなり欲しい近似値を得ることが難しくなってくる。そこで、この積分をより正確に計算するために、 |
||
:''a'' = ''a''<sub>0</sub> |
:''a'' = ''a''<sub>0</sub> < ''a''<sub>1</sub> < … < ''a''<sub>''n''</sub> = ''b'' |
||
として積分区間 [''a'', ''b''] をより小さい''n''個の部分区間 [''a''<sub>0</sub>,''a''<sub>1</sub>], [''a''<sub>1</sub>,''a''<sub>2</sub>], … [''a''<sub>''n''−1</sub>,''a''<sub>''n''</sub>] に分け、それぞれの区間で台形 |
として積分区間 [''a'', ''b''] をより小さい''n'' 個の部分区間 [''a''<sub>0</sub>,''a''<sub>1</sub>], [''a''<sub>1</sub>,''a''<sub>2</sub>], … , [''a''<sub>''n''−1</sub>,''a''<sub>''n''</sub>] に分け、それぞれの区間で台形とみなして面積を近似計算する。 |
||
⚫ | |||
⚫ | |||
&= \sum_{k=1}^{n} \int_{a_{k-1}}^{a_k} f(x)\, \mathrm{d}x \\ |
|||
⚫ | |||
=== 区間幅が等間隔の場合 === |
|||
さらに一区間の幅が |
|||
⚫ | |||
と一定になるように''k'' 番目の区切りを |
|||
⚫ | |||
ととることで |
|||
:<math>\begin{align} |
|||
⚫ | |||
&= \frac{b-a}{2n} \sum_{k=1}^{n}(f(a_{k-1}) + f(a_k))\\ |
|||
⚫ | |||
⚫ | |||
\end{align}</math> |
|||
⚫ | |||
区間分割が十分多い、すなわち''n'' が十分大きい場合、上式右辺の第1項は |
|||
:<math>\frac{b-a}{n}\cdot\frac{f(a) + f(b)}{2} \rightarrow 0 \quad (n\rightarrow\infty)</math> |
|||
となるため、端点''x'' = ''a'', ''b'' における関数''f'' の値は積分の結果に大きな影響を及ぼさず、長方形近似した場合とほぼ同じとなる(誤差のオーダーが同じ)。 |
|||
⚫ | |||
⚫ | |||
⚫ | |||
さらに1[[区間 (数学)|区間]]の幅が <math> a_n - a_{n-1} = \frac{b-a}{n}</math> と、一定になり |
|||
⚫ | |||
⚫ | |||
:<math>\frac{b-a}{2n}</math> は定数より |
|||
⚫ | |||
⚫ | |||
⚫ | |||
⚫ | |||
==関連項目== |
==関連項目== |
2012年9月20日 (木) 00:43時点における版
数学において、台形公式(だいけいこうしき、英語: Trapezoidal rule)もしくは台形則(だいけいそく)は定積分を近似計算するための方法、すなわち数値積分のひとつである。これははニュートン・コーツの公式の1次の場合である。被積分関数を区分線形関数で近似し、台形の面積の公式に帰着させて積分の近似値を求める。
具体的に言えば、求めたいx -y グラフのy = 0を含む面積内に無数の台形を置くと、その台形の面積の集合和は本物の面積に限りなく近い値となる。
一次関数による近似なので精度はそれほど期待できず、二次関数で近似するシンプソンの公式などの方が精度が高い。シンプソンの公式やその他の類似の手法は、2階連続微分可能な関数に対する台形公式の改良とみなせるが、細かく変動しない関数に対しては台形公式で十分であり、計算も簡単である。
近似誤差
台形公式の利点は、近似誤差が容易に分かることである。
凸関数に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ小さい値になり、凹関数に対してこの公式で積分を求めると、結果は実際の値よりも台形と実際の関数曲線の差分の分だけ大きい値になる。また積分区間が変曲点を含むとき上記の凸部分の誤差と凹部分の誤差が打ち消し合い全体的な誤差は小さくなる。
さらに、台形公式は周期関数をその周期よりも長い区間積分する場合にはきわめて精度が高くなる傾向がある。これはオイラーの和公式との関係をみると良く理解できる。しかしながら非周期関数に対しては一般に、 ガウス求積やクレーンショー・カーチス数値積分則のような非等分点法の方がより精度が高い。
この節の加筆が望まれています。 |
台形による近似
x ∈(a,b) のとき f (x) > 0 とする。定積分
の値は、xy 座標平面上で曲線 y = f (x) と x 軸、 x = a、 x = b とで囲まれた図形の面積になる。一般に積分計算は複雑であることも多いが、近似値が欲しいだけであれば f (x) を一次関数で近似することによりこの図形は台形とみなすことができるので
のように、台形の面積公式を用いて簡単に計算できる。
この公式は曲線 y = f (x) が直線に近いほど精度がよくなるが、直線とは全く異なる曲線であれば精度が悪くなり欲しい近似値を得ることが難しくなってくる。そこで、この積分をより正確に計算するために、
- a = a0 < a1 < … < an = b
として積分区間 [a, b] をより小さいn 個の部分区間 [a0,a1], [a1,a2], … , [an−1,an] に分け、それぞれの区間で台形とみなして面積を近似計算する。
区間幅が等間隔の場合
さらに一区間の幅が
と一定になるようにk 番目の区切りを
ととることで
という公式が得られる。
区間分割が十分多い、すなわちn が十分大きい場合、上式右辺の第1項は
となるため、端点x = a, b における関数f の値は積分の結果に大きな影響を及ぼさず、長方形近似した場合とほぼ同じとなる(誤差のオーダーが同じ)。
関連項目
参考文献
- 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