ビジネス

ラグランジュ補間の計算方法は?具体例や手順も解説!(2点補間・3点補間:例題・分点など)

ラグランジュ補間の意味と計算結果
当サイトでは記事内に広告を含みます

ラグランジュ補間の計算方法は?具体例や手順も解説!(2点補間・3点補間・例題・分点など)

ラグランジュ補間は、いくつかの既知の点を通る多項式を作り、未知の位置にある値を推定したいときに役立つ計算方法です。

2点から直線を作る方法を発展させた考え方であり、3点以上では曲線的な変化も扱えます。

数式だけを見ると難しく感じやすい一方、各点に対応する基底多項式を順番に作れば、計算の流れは整理できます。

この記事では、ラグランジュ補間の基本式、2点補間と3点補間の例題、分点の考え方、計算時の注意点をわかりやすく解説します。

ラグランジュ補間の意味と計算結果

ラグランジュ補間の意味と計算結果

それではまずラグランジュ補間の意味と、どのような結果を得られるのかについて解説していきます。

既知の点から未知の値を求める考え方

ラグランジュ補間とは、座標平面上でわかっている複数の点をすべて通る多項式を利用し、途中の値や近い範囲の未知の値を求める手法です。

たとえば、横軸を時間、縦軸を温度としたとき、複数時点の温度が記録されていれば、その間の時刻における温度を推定できます。

点が2つなら一次式、3つなら二次式、4つなら三次式を作るのが基本です。

つまり、n個の点があれば、通常はn-1次以下の補間多項式を1つ作れます。

ここで大切なのは、補間で求める対象が原則として既知データの間にある値だという点です。

既知の範囲を大きく外れた位置を求める場合は外挿となり、実際の変化と離れる可能性が高くなります。

ラグランジュ補間は、与えられた点を必ず通る多項式を直接組み立てる方法です。

未知の値を推定するときは、点の並び、求めたい位置、データが対象範囲内かどうかを最初に確認しましょう。

補間多項式の基本式

n個のデータ点を、xの値とyの値の組として考えます。

i番目の点を、x_i と y_i と表すと、ラグランジュ補間多項式は各点のy_iに基底多項式を掛けて足し合わせる形になります。

P(x) = y_0L_0(x) + y_1L_1(x) + … + y_nL_n(x)

L_i(x) = 他のすべてのx座標との差を使って作る多項式

基底多項式L_i(x)には、対応する点では1になり、それ以外の既知点では0になる特徴があります。

この性質があるため、xがx_0ならy_0だけが残り、xがx_1ならy_1だけが残る仕組みです。

複雑に見える式でも、各点の寄与を必要な場所だけ残すための重み付けと考えると理解しやすくなります。

直線補間との違い

2点だけを用いるラグランジュ補間は、一般に知られる線形補間と同じ結果になります。

2点を結ぶ直線上で値を求めるため、変化がほぼ一定と見なせる場面に向いています。

一方で3点補間では二次式を用いるため、データに山や谷のような曲がりがある場合にも対応可能です。

ただし、点の数を増やせば必ず精度が上がるとは限りません。

高次の多項式は端の付近で大きく振れやすく、観測誤差も増幅されることがあります。

そのため、実務では必要以上に高い次数を使わず、求めたい位置の近くにある点を選ぶことが重要です。

比較項目 2点補間 3点補間
使用する点の数 2点 3点
補間式の次数 一次式 二次式
グラフの形 直線 放物線状の曲線
向くデータ 変化がほぼ一定のデータ 緩やかな曲がりがあるデータ
計算の難しさ 比較的簡単 分数計算に注意が必要

2点補間の計算手順

続いては2点補間の計算手順を確認していきます。

2点補間で使う式

2点補間では、2つの点をAとBとします。

Aを x_0 = a、y_0 = p、Bを x_1 = b、y_1 = q と置くと、補間式は次のように表せます。

P(x) = p × {x-b}÷{a-b} + q × {x-a}÷{b-a}

この式では、最初の項が点Aのy値に対応し、次の項が点Bのy値に対応しています。

xにaを代入すると最初の係数が1、次の係数が0になります。

xにbを代入した場合は反対になり、必ず元の2点を通ることが確認できます。

式を丸暗記するより、片方の点では1、もう片方の点では0になる係数を作ると理解すると、計算ミスを減らせます。

2点補間の具体例

点Aを x = 1、y = 4、点Bを x = 5、y = 12 とします。

このとき、x = 3におけるyの値をラグランジュ補間で求めます。

最初の項は、4 × {3-5}÷{1-5}です。

これは4 × 2÷4となるため、計算結果は2です。

次の項は、12 × {3-1}÷{5-1}です。

こちらは12 × 2÷4となり、計算結果は6になります。

よって、P(3)は2+6で8です。

元の点を結ぶ直線では、xが1から5まで4増える間にyが4から12まで8増えています。

xが1から3まで2増える場合は、yが4増えて8となるため、ラグランジュ補間の結果と直線的な考え方が一致します。

2点補間は、基準点からの距離の比率を使って値を配分する計算です。

求めたいxが2点のちょうど中間なら、yの値も両端の平均になるとは限りませんが、直線上では平均値になります。

分点として考える方法

2点補間は、内分点や外分点の考え方と深く関係しています。

求めたいxが2つのx座標の間にあれば、x方向の位置に応じてyを分ける計算になります。

たとえばxが1と5の間で3なら、両端からの距離はどちらも2です。

そのため、yの値は4と12を等しい割合で混ぜた8になります。

xが2なら、1からの距離は1、5からの距離は3です。

この場合は点Aの影響が大きく、P(2)は6になります。

距離が近い点ほど影響が大きくなるため、分点の位置と重みの関係を意識すると式の意味が見えやすくなるでしょう。

求めるx 点Aからの距離 点Bからの距離 補間値
1 0 4 4
2 1 3 6
3 2 2 8
4 3 1 10
5 4 0 12

3点補間の計算方法

続いては3点補間の計算方法を確認していきます。

3つの基底多項式

3点補間では、点を x_0、x_1、x_2 の3つ用意します。

それぞれのy座標を y_0、y_1、y_2 とすると、二次式の補間多項式が作れます。

最初の点に対応する項では、xから残り2点のx座標を引き、それぞれの差で割ります。

2番目と3番目の項も同様に作り、最後に足し合わせます。

P(x) = y_0 × {x-x_1}{x-x_2}÷{x_0-x_1}{x_0-x_2}

+ y_1 × {x-x_0}{x-x_2}÷{x_1-x_0}{x_1-x_2}

+ y_2 × {x-x_0}{x-x_1}÷{x_2-x_0}{x_2-x_1}

式の構造を見ると、各項は対応する点では1になり、ほかの2点では0になるよう設計されています。

そのため、既知点の座標を代入すれば、必ず与えられたy値に戻ります。

分母の符号を間違えやすいので、分子と分母に含まれる差の順番をそろえることが計算の基本です。

3点補間の例題

点Aを x = 0、y = 1、点Bを x = 1、y = 2、点Cを x = 2、y = 5 とします。

x = 1.5における値を求めてみましょう。

点Aに対応する項は、1 × {1.5-1}{1.5-2}÷{0-1}{0-2}です。

計算すると、0.5×-0.5を2で割る形になり、-0.125です。

点Bに対応する項は、2 × {1.5-0}{1.5-2}÷{1-0}{1-2}です。

結果は2 × 0.75となり、1.5です。

点Cに対応する項は、5 × {1.5-0}{1.5-1}÷{2-0}{2-1}です。

こちらは5 × 0.375となり、1.875です。

3つを合計すると、P(1.5)は3.25になります。

一次補間では直線的にしか変化を表せませんが、この例ではyが1、2、5と増え方を変えています。

その曲がりを反映できる点が、3点補間で二次式を使う大きな特徴です。

二次式へ展開する場合

ラグランジュ形式のまま計算すれば補間値は求められますが、式を展開して二次式にする場面もあります。

たとえば、複数のxについて繰り返し値を求めたい場合は、P(x)を ax²+bx+c の形に整理しておくと便利です。

先ほどの3点では、補間多項式を整理すると x²+1 になります。

実際にxへ0、1、2を代入すると、それぞれ1、2、5となり、元の点と一致します。

xへ1.5を代入すれば、1.5²+1で3.25です。

ただし、展開は途中の符号や係数で間違えやすいため、最初はラグランジュ形式のまま代入する方法がおすすめです。

展開後には既知点を代入し、作成した式がすべてのデータ点を通るかを必ず検算しましょう。

分点と補間値の関係

続いては分点と補間値の関係を確認していきます。

内分する位置での重み

内分とは、2つの点の間を一定の比で分ける位置を考えることです。

2点補間において、求めたいxがAとBの間にあるなら、その位置は内分点として扱えます。

Aから近い位置ではAのy値の影響が強く、Bに近づくほどBの影響が増えます。

たとえばAからBへ全体の4分の1だけ進んだ位置では、Aの重みは4分の3、Bの重みは4分の1です。

このように、補間値は距離に応じた重み付き平均として理解できます。

ただし、重みは自分に近い点ではなく、反対側までの距離で表されることに注意が必要です。

Aに近いほどAの係数が大きくなるよう、式全体は自然に組まれています。

外分と外挿の注意点

求めたいxが2点の外側にある場合、式そのものは同じように計算できます。

しかし、これは補間ではなく外挿に近い考え方です。

たとえば、測定値がx = 1からx = 5までしかないのに、x = 20の値を求める場合、直線や曲線の傾向がその先も続く保証はありません。

3点以上を使う高次の多項式では、外側で急激に大きな値や小さな値になることもあります。

計算結果が出たとしても、その数値が現実的かを別途判断する必要があります。

補間は観測された範囲の中で使うほど信頼しやすくなります。

範囲外の予測を行うなら、補間多項式だけに頼らず、物理的な制約や追加データも確認してください。

等間隔データと不等間隔データ

x座標が0、1、2のように等間隔で並んでいるデータは、途中計算を比較的整理しやすい特徴があります。

一方で、x座標が0、0.7、2.4のように不等間隔でも、ラグランジュ補間は利用できます。

等間隔であることは必須条件ではありません。

ただし、不等間隔では差の計算が複雑になり、近いx座標同士では分母が小さくなります。

分母が非常に小さいと、わずかな測定誤差が結果に大きく影響することがあります。

データの配置を見る際は、点が極端に密集していないか、求めたい位置の近くに適切な点があるかも確認しましょう。

データの状態 補間時の特徴 確認したい点
等間隔 計算過程を整理しやすい 変化が急な区間がないか
不等間隔 任意の位置に対応可能 分母が小さくなりすぎないか
点が密集 局所的な推定に使いやすい 丸め誤差の影響
点が広く離れる 大まかな傾向を見やすい 途中の変化を見落とさないか

計算ミスを防ぐ確認項目

続いてはラグランジュ補間で計算ミスを防ぐ確認項目を確認していきます。

x座標の重複確認

ラグランジュ補間では、使用する各点のx座標がすべて異なっていなければなりません。

x座標が同じ点を2つ入れると、分母に0が現れて計算できなくなります。

同じxに異なるyがある場合は、そもそも1つの関数として扱えないデータです。

同じxに同じyが重複しているだけなら、同一データを重ねて使う意味はありません。

計算前に表へ整理し、xの重複がないかを確認するだけでもトラブルを防げます。

分母が0になるケースは入力データの見直しが必要な合図と覚えておくとよいでしょう。

符号と分数の整理

3点補間以上では、負の数を含む差の計算が増えます。

たとえば x_0-x_1 と x_1-x_0 は絶対値が同じでも符号が逆です。

分子だけ順番を変えたり、分母だけ順番を変えたりすると、項全体の符号が変わってしまいます。

それぞれの基底多項式では、分子と分母に書く順番を対応させることが大切です。

途中で小数に直すと丸め誤差が入りやすいため、可能なら分数のまま進め、最後に小数へ直すと安全です。

電卓や表計算ソフトを使う場合も、手計算で概算を出しておくと、異常な結果に気付きやすくなります。

補間値が既知のy値より大きすぎる、または小さすぎるときは、すぐに誤りと決めつける必要はありません。

ただし、求める位置が既知点の間にあり、データも穏やかに変化しているなら、符号や入力順を再確認する価値があります。

既知点を使った検算

ラグランジュ補間の最も確実な検算は、作成した式に既知のx座標を代入することです。

x_0を代入したときにy_0、x_1を代入したときにy_1となれば、少なくとも通過条件は満たしています。

2点補間なら2点、3点補間なら3点すべてで確認します。

さらに、求めるxが2点の中間にある場合は、結果が大まかにどの程度になるかを先に予想しておくとよいでしょう。

単調に増えているデータなら、内分点の補間値が両端の値から大きく外れることは通常ありません。

このような感覚的な確認と代入による検算を組み合わせると、数式の正しさと結果の自然さを両方確認できます。

ラグランジュ補間の活用場面

続いてはラグランジュ補間の活用場面を確認していきます。

測定データの中間値推定

ラグランジュ補間は、実験、計測、品質管理などで得た離散的なデータの中間値を推定する場面で使われます。

温度と体積、電圧と電流、濃度と吸光度など、連続的に変化する量を扱う場合に考え方を応用できます。

実際のデータには測定誤差が含まれるため、すべての点を厳密に通す補間が常に最適とは限りません。

それでも、少数の基準点から近い範囲の値を知りたい場合には、有力な選択肢になります。

点を増やすほどよいと考えず、求めたい場所の周囲にある信頼できる少数のデータを選ぶことが実用上のポイントです。

数値計算とプログラミング

プログラミングでは、ラグランジュ補間は数値計算の学習題材としてよく使われます。

配列にx座標とy座標を保存し、各点について分子と分母を掛け合わせれば、任意のxの値を計算できます。

少ない点で一度だけ計算するなら、基本式をそのまま実装しても十分です。

一方、多数の点を使って何度も計算する場合は、計算量や数値の安定性を考える必要があります。

そのような場合には、ニュートン補間やスプライン補間が選ばれることもあります。

ラグランジュ補間を学ぶことで、データ点と多項式の関係、数値誤差、補間法の選び方への理解が深まります。

他の補間方法との使い分け

線形補間は、少ない計算で安定した結果を得やすく、多くの実務で利用されています。

ラグランジュ補間は、複数点をまとめて通る多項式を明示的に作りたい場合に便利です。

滑らかな曲線を長い区間で作りたい場合は、区間ごとに低次多項式を使うスプライン補間が適することがあります。

ノイズを含む測定データの傾向を知りたい場合は、すべての点を通す補間ではなく、回帰分析を検討する方法もあります。

目的が中間値の再現なのか、将来値の予測なのか、滑らかな曲線の作成なのかによって、適した手法は変わります。

ラグランジュ補間は、少数の既知点を正確に通過する式を作る目的に特に向いています。

ラグランジュ補間の計算方法まとめ

ラグランジュ補間は、複数の既知点を通る多項式を作り、途中にある未知の値を求める方法です。

2点補間では一次式となり、直線補間と同じ考え方で計算できます。

3点補間では二次式を使えるため、データの曲がりを反映した推定が可能です。

計算では各点に対応する基底多項式を作り、対応するy値を掛けて合計します。

分子と分母の順番をそろえ、x座標の重複がないかを確認し、既知点を代入して検算することが重要です。

求める位置が既知データの範囲内にあるほど、補間の結果は使いやすくなります。

まずは2点補間で距離と重みの関係を理解し、次に3点補間の式へ進むと、ラグランジュ補間を無理なく使いこなせるでしょう。