台形則計算機

台形公式計算機は弧長を数値的に見積もります。入力した関数の弧長積分に合成台形公式を適用します。次の手順に従えば、1分以内に最初の見積もりが得られます。

1. 関数を入力します。関数フィールドに f(x) を入力します。例: x^2、sin(x)、ln(x+1)。べき乗には ^ を使い、sin、cos、exp、ln などの標準的な関数名を使います。
2. 開始値を入力します。積分の下限である a を設定します。x = 0 から x = 1 までの曲線区間では、0 を入力します。
3. 終了値を入力します。積分の上限である b を設定します。計算機では b > a が必要です。
4. 分割数を選びます。等しい部分区間の数である n を入力します。滑らかな曲線では、最初は n = 10 で試すのが適切です。
5. 台形公式の方式を選択します。数値積分のインターフェースで台形公式を選び、下に表示されるパネルの公式が適用されるようにします。
6. 「計算」をクリックします。ツールは n + 1 個の標本点で弧長の被積分関数を評価し、端点の重みを適用して、近似値を表示します。
7. 結果を解釈します。表示される数値は、x = a から x = b までの曲線の長さの見積もりであり、x の値と同じ長さの単位です。
8. 精度を上げます。n を2倍にして再計算します。連続する2回の見積もりがほとんど変わらなくなれば、実用上は収束したとみなせます。

各入力が制御する内容は次のとおりです。

f(x) = x²、a = 0、b = 1、n = 10 と入力すると、1.479 に近い弧長の見積もりが得られ、厳密値 1.478943 に近くなります。このガイドの残りでは、計算機が正確に何を計算しているのか、そしてその答えの質をどう判断するかを説明します。

台形公式は数値積分法の一つです。被積分関数の下に台形を描き、曲線下の面積を台形の面積で置き換えることで定積分を近似します。各台形は隣接する2つの標本点から作られ、すべての台形の面積の合計が積分の見積もりとなります。

この名前は幾何学に由来します。2つの標本点 x_i と x_i+1 の間で、被積分関数 g(x) は直線で置き換えられます。その直線の下の領域は台形になります。高さが g(x_i) と g(x_i+1) の2つの平行な垂直辺と、水平方向の幅 h = x_i+1 - x_i を持つ台形です。

1つの台形の面積は h/2 · [g(x_i) + g(x_i+1)] です。n 個の台形すべての面積を足し合わせると、積分の合成台形近似が得られます。「合成（composite）」とは、区間を1つではなく複数のパネルに分けることを意味します。

この公式は数値求積法の一族に属します。透明性が高く評価されています。すべての手順に明確な幾何学的意味があり、唯一のパラメータは分割数 n です。主な考え方は簡単に言えます。より真っすぐで細いパネルほど被積分関数に密接に沿うため、n が大きくなるほど見積もりが改善されます。

弧長の場合、被積分関数は元の関数 f(x) ではありません。弧長の被積分関数 g(x) = √(1 + (f'(x))²) であり、これは曲線に沿って長さがどれだけ速く蓄積されるかを測るものです。この区別を明確にしておくことで、この手法の最も一般的な誤解を防げます。

台形公式には2つの形があります。単一パネル版と、実用で使われる合成版です。どちらも同じ台形の幾何学から導かれます。

単一パネルの公式。区間全体を1つの台形で覆う場合:

幅 (b - a) に2つの端点の高さの平均を掛けます。これはちょうど1つの台形の面積です。

合成公式。[a, b] を幅 h = (b - a)/n の n 個の等しい部分区間に分け、標本点を x₀ = a、x₁ = a + h、...、x_n = b とします:

各記号の意味は決まっています:

• T_n は n 個のパネルによる数値近似です。
• h は部分区間の幅、(b - a)/n です。
• a と b は積分の下限と上限です。
• n は部分区間の数、正の整数です。
• x_i は標本点で、h ずつ離れています。
• g(x_i) はそれらの点での被積分関数の値です。

端点が半分の重みを受けるのは、各内部の標本点が隣接する2つの台形に共有されるのに対し、各端点は1つの台形にしか属さないからです。代数的に、隣接する h/2 · [g(x_i) + g(x_i+1)] の項が結合され、内部の値は完全な重み、外側の2つの値は半分の重みになります。

短い計算例。n = 2 で ∫₀¹ x² dx を近似します。すると h = 0.5、x₀ = 0、x₁ = 0.5、x₂ = 1 で、T₂ = 0.5 · [0/2 + 0.25 + 1/2] = 0.5 · 0.75 = 0.375 となります。厳密な積分は 1/3 ≈ 0.333 なので、粗い2パネルの見積もりは正しい範囲にあり、さらに改善の余地があります。

弧長は積分であるため、台形公式を直接適用できます。曲線 y = f(x) の x = a から x = b までの弧長は次式で与えられます:

計算機は弧長の被積分関数 g(x) = √(1 + (f'(x))²) を作り、合成台形公式に渡します。完全な手順は次のとおりです:

1. y = f(x) から始めます。[a, b] で関数が微分可能であることを確認します。
2. 導関数 f'(x) を求めます。ツールは入力された式を内部で微分します。
3. 導関数を2乗します。(f'(x))² を作ります。傾きの方向に関わらず長さが蓄積されるため、2乗によって符号の情報が消えます。
4. 1を加えます。被積分関数は 1 + (f'(x))² で、1より小さくなることはありません。
5. 平方根を取ります。これにより g(x)、すなわち x の単位あたりの長さの蓄積の局所的な割合が得られます。
6. 区間を分割します。n を選んで h = (b - a)/n を計算します。
7. 被積分関数を評価します。標本点で g(x₀)、...、g(x_n) を計算します。
8. 合成公式を適用します。端点に半分の重みを付けて値を組み合わせます。
9. 結果を解釈します。この合計が曲線区間のおおよその長さです。

導関数を2乗して平方根の中に入れるのには幾何学的な理由があります。幅 dx、高さ dy の微小な曲線部分の長さは、ピタゴラスの定理により √(dx² + dy²) = √(1 + (dy/dx)²) dx です。この局所的な長さを [a, b] で積分すると全長が得られ、台形公式がその積分を近似します。

被積分関数は常に非負です。なぜなら 1 + (f'(x))² が少なくとも 1 だからです。最小値は 1 で、曲線が平らなところで達成されます。傾きが急になるほど g(x) は大きくなり、同じ x 区間でも急な曲線の方が緩やかな曲線より長くなるのはまさにこのためです。

この完全に解いた例では、台形公式が出力するすべての数値を示します。f(x) = x² を [0, 1] で、n = 4 として使います。

1. 導関数。 f'(x) = 2x。

2. 被積分関数。 g(x) = √(1 + (2x)²) = √(1 + 4x²)。

3. 刻み幅。 h = (1 - 0)/4 = 0.25。

4. 標本点。 0、0.25、0.5、0.75、1。

5. 被積分関数の値と重み。

6. 重み付き和。 1.000000/2 + 1.118034 + 1.414214 + 1.802776 + 2.236068/2 = 5.953057。

7. h を掛けます。

8. 厳密値と比較します。厳密な弧長は L = √5/2 + asinh(2)/4 ≈ 1.478943 です。n = 4 の見積もり 1.488264 はそれを 0.009321 上回り、誤差は約 0.63% です。n = 8 で繰り返すと 1.481272 となり、すでに厳密値の 0.002330 以内に入ります。

見積もりが厳密値より大きいことに注目してください。この凸型の被積分関数では、直線の弦が曲線の上に来るため、台形公式は過大評価します。この系統的な振る舞いは正常であり、過大評価は n が大きくなるほど縮小します。

分割数 n と刻み幅 h は h = (b - a)/n で結ばれています。この2つが、被積分関数をどれだけ細かく標本化するか、そして見積もりの精度を制御します。

n が表すもの。n は台形パネルの数であり、正の整数です。標本点の数は常に n + 1 です。なぜなら連続するパネルが境界点を共有するからです。

n = 1 が単一パネル公式である理由。n = 1 のとき、h = b - a で T₁ = h/2 · [g(a) + g(b)] となります。これはこの手法の最も粗い形であり、十分な精度が得られることはまれです。

[0, 1] でのパネル幅。下の表は、n が大きくなるにつれて h がどう縮小するかを示しています:

パネルが小さいほど被積分関数に密接に沿うため、各パネル内の直線近似が真の曲線から外れる度合いが小さくなります。そのため、n を増やすと通常は精度が向上します。代償は計算量です。新しいパネルが1つ増えるごとに関数評価が1回増えます。

分割を増やしても、すべての被積分関数で必ずしも良い結果が得られるわけではありません。滑らかでない被積分関数や振動する被積分関数では不規則な誤差の振る舞いが見られることがあり、浮動小数点の丸めが最終的には非常に大きな n での改善を制限します。ただし、このガイドのような滑らかな曲線では、改善は着実で予測可能です。

実用的な戦略: n = 10 または n = 20 から始め、必要な精度で見積もりが変わらなくなるまで n を2倍にします。この収束チェックは、単一の n の値を信じるよりも信頼できます。

6つの例は、厳密に解ける場合から、弧長が初等的な閉形式を持たない関数までをカバーします。各例で導関数、被積分関数、代入、数値結果を示します。

例1: 直線、f(x) = 2x（[0, 3]）。導関数は f'(x) = 2（定数）なので、被積分関数も g(x) = √(1 + 4) = √5（定数）です。被積分関数が曲がらないため単一パネルで厳密になります: T₁ = 3 · √5 ≈ 6.708204。厳密な長さ 3√5 と一致します。これで次の法則が確認できます。被積分関数が線形または定数の場合、台形公式は厳密です。

例2: 二次関数、f(x) = x²（[0, 1]、n = 4）。f'(x) = 2x、g(x) = √(1 + 4x²) のもとで、前の節の完全な計算により L ≈ 1.488264 が得られ、厳密値 1.478943 と比較できます。小さな差は凸型の被積分関数に由来します。曲線の上に描かれた弦が余分な面積を加えるのです。

例3: 長い二次関数、f(x) = x²（[0, 2]、n = 8）。ここでも h = 0.25 ですが、区間は2倍の長さで、被積分関数はより急です。見積もりは T₈ = 4.656890 です。高解像度の計算では 4.646784 が得られるので、誤差は約 0.010106 です。区間を延ばすと答えが単に比例して伸びるわけではありません。被積分関数が x とともに増大するため、弧長は区間の長さよりも速く伸びます。

例4: 正弦曲線、f(x) = sin x（[0, π]）。導関数は f'(x) = cos x で、g(x) = √(1 + cos²x) となります。n = 8 での見積もりは 3.820198、n = 16 でも 3.820198 で、小数点以下6桁まで変わりません。対称な被積分関数は曲率が穏やかで均等に分布しているため、速く収束します。（角度はラジアンで入力してください。度で入力すると異なる曲線になります。）

例5: 対数、f(x) = ln(x + 1)（[0, 2]、n = 10）。導関数は f'(x) = 1/(x + 1) で、g(x) = √(1 + 1/(x + 1)²) となります。定義域は x > -1 を要求しますが、この区間は満たしています。見積もりは T₁₀ = 2.304215、精密な計算では 2.301988 となり、誤差は約 0.002227 です。x = 0 付近で傾きが最も急なため、長さの大部分は区間の前半で蓄積されます。

例6: 三次関数、f(x) = x³（[0, 1]、n = 10）。ここでは f'(x) = 3x²、g(x) = √(1 + 9x⁴) です。この弧長には初等的な原始関数が存在しません。まさにこういう場合に数値積分が力を発揮します。見積もりは T₁₀ = 1.552609 で、精密値 1.547866 に対する誤差は約 0.004743 です。数値解法がなければ、この長さを得るのは非常に困難でしょう。

上記のすべての結果は合成公式から直接計算し、小数点以下6桁に丸めています。同じ n で計算機で繰り返せば、同じ数字が再現されます。

収束とは、n が大きくなるにつれて近似値が真の値に落ち着くことを意味します。下の表は f(x) = x²（[0, 1]、厳密な長さは 1.478943）を追跡したものです:

n が2倍になるたびに、誤差はおよそ4分の1になります: 0.037181、0.009321、0.002330、0.000582、0.000146、0.000036。これは2次の収束であり、滑らかな被積分関数に対する合成台形公式の特徴です。h を半分にすると、主要な誤差項が約4分の1になります。

「前回からの変化」の列は、計算機のワークフローが推奨する実用的な収束チェックです。n を2倍にしても答えが必要な精度より小さくしか変わらなくなれば、見積もりは安定しています。この曲線では、n = 16 ですでに真の値の 0.000582 以内に長さが確定しています。

合成台形公式の誤差には正確な数学的形があります。[a, b] 上で十分滑らかな被積分関数 g に対して:

ただし ξ は [a, b] 内のある値です。符号が重要です。g'' が正（凸型の被積分関数）のとき誤差は負になり、T_n が積分を過大評価することを意味します。これは x² の例と一致します。すべての見積もりが厳密値 1.478943 を上回っていました。

g''(ξ) を [a, b] 上での最大の大きさ M = max|g''(x)| で置き換えると、保証された上限が得られます:

重要な用語。絶対誤差は |I - T_n|、すなわち真の値と見積もりの距離です。相対誤差はそれを |I| で割ったものです。打ち切り誤差は、積分を有限和で近似することによる誤差であり、有限精度の演算による丸め誤差とは別物です。2階導関数 g'' は被積分関数の曲率を測るもので、これこそが直線パネルでは捉えきれない部分です。

x² の例に対する具体的な上限。[0, 1] では g''(x) = 4/(1 + 4x²)^3/2 で、x = 0 で最大となり M = 4 です。この上限から |I - T₄| ≤ 4/(12 · 16) = 0.020833 が得られます。実際の誤差 0.009321 は余裕を持ってその範囲内です。上限は設計上保守的であり、正確な誤差ではなく上限を約束するものです。

2つの注意点があります。第一に、g は弧長の被積分関数であり、元の f ではありません。平方根を通して2回微分したものを M が表しています。第二に、複雑な関数では誤差の上限の計算が面倒になることがあります。M 自体を最大化する必要があるからです。実際には、n を2倍にする収束チェックの方が有用なツールであり、この計算機は自動的な誤差上限を表示しません。

シンプソンの公式は台形公式の最も近い親戚です。どちらも [a, b] を等しい部分区間に分けますが、台形公式が2点を通る1本の直線の弦を引くのに対し、シンプソンは連続する3点を通る放物線を当てはめます。

滑らかな被積分関数では、シンプソンの4次収束が通常は勝ります。n を2倍にすると誤差がおよそ16分の1になりますが、台形公式では4分の1です。しかし、シンプソンが常に優れているわけではありません。粗い被積分関数や不連続な被積分関数では、どちらの手法も劣化し、余分な精度の次数が現れないことがあります。台形公式は奇数の n も受け付けるため、手早いチェックに便利です。

堅実なワークフローでは両方を使います。透明な基準として台形公式の見積もりを、高い精度のために シンプソンの公式 を使います。両方が必要な桁数で一致すれば信頼度は高くなります。一致しない場合は、どちらかを信じる前に n を増やしてください。

2つの異なる「台形のような」考え方は混同しやすいものです。同じ曲線に対して異なる数値を生むため、その区別は重要です。

方法1: 弧長の被積分関数の台形積分。これがこの計算機の行うことです。被積分関数の曲線の下に台形パネルを描いて、g(x) = √(1 + (f'(x))²) を積分します:

方法2: 元の曲線の折れ線（弦の和）近似。曲線自体の上の隣接する点を直線の線分で結び、その長さを合計します:

弦の和は元の曲線上の点 (x_i, f(x_i)) を結びます。一方、台形パネルは被積分関数 g(x) の下の面積を近似します。どちらの図にも直線が登場しますが、パネルは元の曲線の直線部分を表しているわけではありません。

どちらの手法も n が大きくなるにつれて真の長さに収束し、どちらも正当な数値ツールです。このページの計算機は方法1を実装しています。座標データから弦の和が必要な場合は、座標点からの弧長ツールが代わりに方法2に従います。

単純さと透明性が、精度の最後の桁まで絞り出すことよりも重要な場合、台形公式が適切な選択です。

• 数値積分の学習。すべての手順に目に見える幾何学的意味があり、微積分の授業で標準的な最初の手法となっています。
• 曲線の長さの迅速な見積もり。滑らかな曲線では、少数のパネルで2〜3桁の正しい数字が得られることが多いです。
• 解析解のチェック。手計算で導いた弧長を、独立した数値見積もりと比較します。
• 工学や物理のモデル。ケーブル、アーチ、軌道の長さの見積もり。近似値がより大きな計算に使われます。
• 数値解法の比較。シンプソンの公式や高次の手法が判断される基準となります。
• 滑らかな標本関数。f が滑らかで n を中程度にできる場合、この手法は信頼できます。
• 教育的な検証。予測可能な O(h²) の振る舞いにより、収束チェックが解釈しやすくなります。

被積分関数が急激に変化する場合、折れ点や不連続点がある場合、または低コストで高精度が必要な場合は、別の手法を選んでください。滑らかな曲線では、シンプソンの公式が通常、パネルあたりより多くの桁を提供します。パラメトリック曲線、極座標曲線、空間曲線には、専用のパラメトリック、極座標、または3D弧長計算機を使ってください。同じ数値的アイデアを適切な座標系で適用します。

これらは学生の課題やサポートへの質問で最もよく見られる誤りです。存在を知ってさえいれば、どれも簡単に避けられます。

1. 誤った積分範囲の入力。区間 [0, 2] は [2, 0] と同じ曲線部分ではありません。常に a < b を確認し、範囲が問題と一致していることを確かめてください。
2. 間違った関数の使用。x^2 のつもりで x^3 と打ち間違えると、導関数とその後のすべての値が変わります。
3. 弧長の公式で導関数を忘れる。弧長は √(1 + (f')²) を積分するものであり、f 自体を積分するのではありません。
4. 弧長のために公式を f(x) に直接適用する。パネルは g(x) の下になければなりません。f(x) を積分すると符号付き面積が計算され、長さではありません。
5. 分割数が少なすぎる。曲がった被積分関数では n = 2 で数パーセントずれることがあります。x² の例では n = 2 で 2.5% の誤差がありました。
6. 分割数を0にする。n = 0 は h = (b - a)/0 を意味し、これは定義されません。n は正の整数でなければなりません。
7. 整数でない分割数の入力。n = 2.5 はパネルの数え方として意味がありません。整数に丸めてください。
8. n を標本点の数と混同する。n 個のパネルには n + 1 個の標本点が必要です。n = 4 では4点ではなく5点です。
9. 端点の半分の重みを忘れる。g(x₀) と g(x_n) に完全な重みを与えると、見積もりは h/2 · (g(x₀) + g(x_n)) だけ膨らみます。
10. 間違った刻み幅の使用。h = (b - a)/n であり、b - a でも 1/n でもありません。[0, 2] で n = 8 の場合、h = 0.25 です。
11. 度とラジアンの混同。x を度で測った sin(x) は、ラジアンで測った sin(x) とは異なる曲線であり、弧長もそれに応じて異なります。
12. 関数の定義域の制限の無視。ln(x + 1) は x > -1 を要求します。定義域の外で評価すると無意味な値が生じます。
13. 近似と厳密な積分の混同。T_n は L に近づくものであり、定数のような特殊な被積分関数の場合にのみ L と等しくなります。
14. パネルが多いほど必ず精度が上がると思い込む。滑らかな被積分関数の漸近的な領域では真実ですが、極端な場合には丸め誤差や非滑らかさがこのパターンを崩すことがあります。
15. 収束チェックの無視。1回の見積もりは推測にすぎません。n と 2n での2つの見積もりが一致して初めて根拠となります。
16. シンプソンの公式と台形公式の混同。シンプソンは 1-4-2-4-...-1 の重みを使い、偶数の n が必要です。台形公式は 1/2-1-...-1-1/2 の重みを使い、任意の正の n が使えます。
17. 台形積分と弦の合計の混同。被積分関数の下のパネルは曲線上の弦ではありません。2つの手法は異なる数値を与えます。
18. 過剰な小数精度の報告。誤差が 0.002 程度の場合、小数点以下6桁を書くことは分かっている以上の主張になります。確認された精度に合わせて丸めてください。

弧長の背後にあるのと同じパネルの算術は、滑らかな量を数値的に蓄積しなければならないあらゆる場面に現れます。

• 工学における曲線長の見積もり。ケーブル、懸垂曲線、配管の曲がりの長さは、製造前に設計関数から計算されます。
• 機械設計。カムの輪郭や歯車の歯の曲線は、材料計画のために正確な周長見積もりが必要です。
• 道路や通路の幾何学。道路や鉄道の線形曲線は中心線に沿って測られますが、それらが初等関数であることはまれです。
• コンピュータ支援設計。CAD カーネルは、レンダリング、工具経路の生成、寸法記入のために曲線の長さを見積もります。
• 物理シミュレーション。粒子の軌道や磁力線の長さは、数値モデルから蓄積されます。
• 数値モデリング。確率密度から信号エネルギーまで、閉形式を持たないあらゆる積分は、台形パネルから始めることができます。
• 数学教育。誤差の振る舞いが示しやすいため、この公式は数値解析への標準的な入り口です。
• 曲線の測定と検証。不規則な曲線の物理的測定値は、数値予測と照合されます。

特殊な幾何学にはこの手法が適応します。同じ合成公式がパラメトリック速度 √((x')² + (y')²) や極座標の被積分関数を積分します。ここの計算機は直交座標の場合を扱い、関連ツールが他の座標系をカバーします。CAD ファイルや測量データセットを直接処理するわけではありません。それらはまず関数と範囲に落とし込む必要があります。

これらを手計算で解き、計算機で照合してください。各解答では代入を示しているので、自分の計算ミスを見つけられます。

問題1。n = 1 の台形公式で ∫₀² 3 dx を近似し、結果が厳密かどうかを述べてください。
解答。h = 2、T₁ = 2/2 · [3 + 3] = 6。被積分関数は定数なので、単一の台形は長方形と正確に一致します。厳密な積分も 6 なので、結果は厳密です。

問題2。n = 1 の台形公式で f(x) = x の [0, 1] での弧長を求めてください。
解答。f'(x) = 1 なので g(x) = √2（定数）。T₁ = 1 · √2 ≈ 1.414214。被積分関数が定数のため厳密です。幾何学的には、これは単位正方形の対角線です。

問題3。n = 2 で f(x) = x² の [0, 1] での弧長を近似してください。
解答。h = 0.5。g(0) = 1、g(0.5) = 1.414214、g(1) = 2.236068。T₂ = 0.5 · [1/2 + 1.414214 + 2.236068/2] = 0.5 · 3.032249 = 1.516124。厳密値は 1.478943 なので、誤差は 0.037181 です。

問題4。n = 4 で f(x) = sin x の [0, π] での弧長を近似してください。
解答。g(x) = √(1 + cos²x)、h = π/4。0、π/4、π/2、3π/4、π で評価し、重みを適用すると T₄ ≈ 3.819944。n = 8 では 3.820198 なので、n = 4 の結果はすでに収束値の 0.000254 以内に入っています。

問題5。a = 2、b = 5、n = 6 のときの刻み幅を計算してください。
解答。h = (5 - 2)/6 = 3/6 = 0.5。標本点は 2、2.5、3、3.5、4、4.5、5 で、6つのパネルに対して7点です。

問題6。n = 10 のとき、標本点は何個あり、それぞれの重みはどうなりますか。
解答。標本点は n + 1 = 11 個です。端点 x₀ と x₁₀ は重み 1/2、内部の9点 x₁ から x₉ は重み 1 です。

問題7。f(x) = x² の [0, 1] について、T₂ と T₄ を比較し、変化を述べてください。
解答。T₂ = 1.516124、T₄ = 1.488264 で、変化は 0.027860 です。厳密値は 1.478943 なので、n を2倍にすることで誤差は 0.037181 から 0.009321 に減少しました。約4分の1であり、2次の収束の予測通りです。

問題8。f(x) = x² の [0, 1] で T₈ = 1.481272。厳密値 1.478943 に対する絶対誤差と相対誤差を報告してください。
解答。絶対誤差 = |1.481272 - 1.478943| = 0.002330。相対誤差 = 0.002330 / 1.478943 ≈ 0.001576、約 0.16%。ここでは見積もりと誤差の両方を小数点以下3桁で述べるべきです。1.481272 の余分な桁はまだ信頼できません。