logo

ラグランジュの公式の導出 📂数値解析

ラグランジュの公式の導出

公式 1

異なるx0,,xnx_{0} , \cdots , x_{n}のデータ(x0,y0),,(xn,yn)(x_{0}, y_{0}) , \cdots , (x_{n} , y_{n})についてli(x):=ij(xxjxixj)\displaystyle l_{i} (x) := \prod_{i \ne j} \left( {{ x - x_{j} } \over { x_{i} - x_{j} }} \right)とすると、 pn(x)=i=0nyili(X) p_{n} (x) = \sum_{i=0}^{n} y_{i} l_{i} (X)

説明

ラグランジュの公式は、多項式補間を見つける方法の中で最もシンプルな公式だ。

導出

戦略: lil_{i}がインデックスに対するクロネッカーのデルタ関数であることを示す。


li(xi)=ij(xixjxixj)=1 l_{i} (x_{i}) = \prod_{i \ne j} \left( {{ x_{i} - x_{j} } \over { x_{i} - x_{j} }} \right) = 1

li(xj)=ij(xjxjxixj)=0 l_{i} (x_{j}) = \prod_{i \ne j} \left( {{ x_{j} - x_{j} } \over { x_{i} - x_{j} }} \right) = 0 整理するとli(xj)=δijl_{i}(x_{j}) = \delta_{ij}である。 pn(x)=y0l0(x)+y1l1(x)+ynln(X) p_{n}(x) = y_{0} l_{0}(x) + y_{1} l_{1}(x) + \cdots y_{n} l_{n}(X) これを設定すると、全てのi=0,1,,ni=0,1, \cdots , nについて、 pn(xi)=0++yi1++0=yi p_{n}(x_{i}) =0 + \cdots + y_{i} \cdot 1 + \cdots + 0 = y_{i} が成立する。


  1. Atkinson. (1989). An Introduction to Numerical Analysis(2nd Edition): p134. ↩︎