ラグランジュ誤差界計算機
ラグランジュの剰余定理を使って、テイラー多項式近似の最大誤差を見積もります。
下の4つのパラメータを入力すると、テイラー多項式近似の誤差上限を計算できます。
ラグランジュ誤差界計算機
ラグランジュの剰余定理を使って、テイラー多項式近似の最大誤差を見積もります。
ラグランジュ誤差界の例
次数を上げる、または区間を狭めることで誤差界がどのように縮むかを示す4つの定番例です。
| 関数 / 設定 | 誤差界 | 詳細 |
|---|---|---|
| eˣ, n=3, a=0, x=0.5, M=1.6487 | ≤ 0.004298 | eˣ の 4 階導関数は eˣ のままです。[0,0.5] での最大値は e⁰⋅⁵ ≈ 1.6487。上限 = 1.6487/24 × 0.5⁴。 |
| cos(x), n=2, a=0, x=0.1, M=0.09983 | ≤ 0.00001664 | cos(x) の 3 階導関数は sin(x) です。[0,0.1] での最大値は約 0.09983。上限 = 0.09983/6 × 0.1³。 |
| ln(x), n=3, a=1, x=1.2, M=6 | ≤ 0.0004 | ln(x) の 4 階導関数は 6/x⁴ です。[1,1.2] では x=1 で最大となるため M=6。上限 = 6/24 × 0.2⁴。 |
| √x, n=2, a=4, x=4.1, M=0.01172 | ≤ 0.0000000195 | √x の 3 階導関数は (3/8)x⁻⁵ᴱ² です。x=4 で最大となるため M≈0.01172。上限 = 0.01172/6 × 0.1³。 |
ラグランジュ誤差界計算機について
ラグランジュ誤差界は、テイラー剰余定理またはラグランジュ剰余とも呼ばれ、テイラー多項式が近似する実際の関数からどれだけ外れうるかに厳密な上限を与えます。eˣ、cos(x)、ln(x) のような複雑な関数を次数 n の多項式で置き換えると、打ち切り誤差が生じます。ラグランジュ界は、その誤差が指定区間で最悪どこまで大きくなるかを示すため、精度が重要な場面では欠かせません。
式は |Rₙ(x)| ≤ M × |x − a|ⁿ⁺¹ / (n+1)! です。ここで n はテイラー多項式の次数、a は展開中心（多項式を作る基準点）、x は近似値を評価する具体的な点、M は a と x の間の閉区間における (n+1) 階導関数の絶対値の最大値です。重要なのは、n が大きくなるほど誤差が小さくなることです。分母の階乗は、分子の (x − a) のべきよりはるかに速く増えるからです。
M を見つけることが、この手順で最も考察を要する部分です。まず関数の (n+1) 階導関数を記号的に求め、区間 [a, x]（x < a なら [x, a]）でその絶対値の最大値を見つけます。指数関数や三角関数のような扱いやすい関数では、M は比較的簡単です。eˣ の (n+1) 階導関数は依然として eˣ なので、右端点での eˣ を M に取れば十分です。cos(x) ではすべての導関数が 1 以下に抑えられるため、M = 1 は常に安全です（ただし、もっと厳しい上限を取れることもあります）。その他の関数では、記号微分を行い、区間上の式を簡単に解析すれば足ります。
この上限は、数値解析、科学計算、工学で広く使われます。電卓、計算機代数システム、あるいは多項式で超越関数を評価する組み込みファームウェアでは、結果が所定の小数精度を満たすよう、内部で何らかの形でこの上限が働いています。物理学では、波動関数、ポテンシャル面、確率密度の多項式近似にも同様の精度要件があります。金融では、オプション価格モデルの級数展開も制御された打ち切り誤差に依存します。
よくある誤解は、高次多項式なら必ず誤差が小さくなるというものです。一般には次数を上げるほど上限は締まりますが、|x − a| が大きい場合や導関数の増え方が速い場合には、単純に次数を上げても効果がないことがあります。その場合は、展開中心 a を評価点 x にできるだけ近づけ、必要な許容誤差を下回るまで n を増やすのが最善です。
この計算機はラグランジュ公式の計算を自動化します。M（これは導関数の分析を自分で行う必要があります）、n、a、x を入力すれば、すぐに上限を算出します。結果は保証です。真の絶対誤差 |f(x) − Pₙ(x)| は、表示された値を超えません。
ラグランジュ誤差界計算機の使い方
- 近似したい関数 f(x)、テイラー多項式の次数 n、展開中心 a、評価点 x を確認します。
- f(x) の (n+1) 階導関数を記号的に求め、a と x の間の閉区間でその絶対値の最大値 M を見つけます。
- M、n、a、x を4つの入力欄に入れ、「誤差界を計算」をクリックします。
- 結果を確認します。表示値は |f(x) − Pₙ(x)| の上限です。実際の誤差はそれ以下です。
- 上限がまだ大きすぎる場合は、n を増やすか、x により近い展開中心 a を選んで再計算してください。
よくある質問
ラグランジュ誤差界とは何ですか？
ラグランジュ誤差界は、テイラー多項式近似の誤差が M × |x − a|ⁿ⁺¹ / (n+1)! を超えないことを保証する定理です。ここで M は区間内の (n+1) 階導関数の絶対値の最大値です。打ち切り誤差の厳密で計算可能な最悪値を与えます。
M の値はどうやって求めますか？
関数を n+1 回微分し、その導関数の絶対値を a と x の間のすべての点で評価します。M はその最大値です。eˣ では導関数は常に eˣ なので、M は大きい方の端点での e のべき乗にできます。正弦と余弦では、すべての導関数が 1 以下なので M = 1 は常に有効です（ただし、より厳密にできます）。
次数を上げれば、必ず誤差界は小さくなりますか？
一般にはそうです。多くの標準的な関数や小さな区間では、分母の (n+1)! が分子の |x−a|ⁿ⁺¹ より速く増えるからです。ただし |x−a| が大きい場合や、関数の導関数が急速に増える場合は、次数を上げても必ずしも役立たず、区間分割のような別の方法の方が有効なことがあります。
誤差界と実際の誤差の違いは何ですか？
実際の誤差 |f(x) − Pₙ(x)| は、点 x における関数と多項式の本当の差です。ラグランジュ界は、その誤差に対する保証された上限です。実際の誤差はほとんど常に誤差界より小さく、誤差界は保守的な最悪値の見積もりです。
マクローリン級数にも使えますか？
はい。マクローリン級数は、a = 0 を中心とするテイラー級数にすぎません。「展開中心（a）」欄に 0 を入力し、通常どおり進めてください。式も計算も同じです。
ラグランジュ誤差界の実用例は何ですか？
計算機や数式ライブラリの多項式近似の精度保証、有限要素解析での補間誤差の上限設定、数値積分で求積法が許容誤差を満たすことの確認、制御系で線形化モデルが真の非線形動特性から許容範囲内しかずれないことの検証などに使われます。テイラー展開が厳密な関数の代わりになるあらゆる場面で、ラグランジュ界は実務者や監査担当者が求める厳密な保証を与えます。