Mètode Lagrange

Prendrem:

$P_{n} (x) = \sum_{k = 0}^{n} f_{k} \cdot l_{k} (x)$

on $l_{k} (x) = \frac{\prod_{i \neq k} (x - x_{i})}{\prod_{i \neq k} (x_{k} - x_{i})} =$ $= \frac{(x - x_{0}) \cdot (x - x_{1}) \dots \hat{(x - x_{k})} \dots (x - x_{n})}{(x_{k} - x_{0}) \cdot (x_{k} - x_{1}) \dots \hat{(x_{k} - x_{k})} \dots (x_{k} - x_{n})}$

on el barret indica que no apareix aquest terme en el producte.

És a dir, en el numerador multipliquem per tots els polinomis de grau $1$ excepte el $k$ -èsim i després dividim pel producte de totes les diferències excepte la $k$ -èsima.Efectivament, si avaluem el polinomi en algun dels $x_{k}$ , obtenim el valor $f_{k}$ .

Al polinomi calculat d'aquesta manera se l'anomena polinomi de Lagrange.

Exemple

Anem a veure un exemple per tenir més clar. Considerem la taula de valors:

$x_{k}$	$1$	$2$	$4$	$5$
$f_{k}$	$0$	$2$	$12$	$21$

i volem calcular el valor $f (3)$ a partir d'interpolació. Tenim doncs $n + 1 = 4$ punts, per tant $n = 3$ , la interpolació serà cúbica. El polinomi s'escriu:

$P_{3} (x) = f_{0} \cdot l_{0} (x) + f_{1} \cdot l_{1} (x) + f_{2} \cdot l_{2} (x) + f_{3} \cdot l_{3} (x) =$

$= 2 \cdot l_{1} (x) + 12 \cdot l_{2} (x) + 21 \cdot l_{3} (x)$

Hem de calcular ara els polinomis de Lagrange associats,

$\begin{aligned} l_{0} (x) = & \frac{(x - x_{1}) \cdot (x - x_{2}) \cdot (x - x_{3})}{(x_{0} - x_{1}) \cdot (x_{0} - x_{2}) \cdot (x_{0} - x_{3})} \\ = & \frac{(x - 2) \cdot (x - 4) \cdot (x - 5)}{(1 - 2) \cdot (1 - 4) \cdot (1 - 5)} = \frac{- 1}{12} (- 40 + 38 x - 11 x^{2} + x^{3}) \\ l_{1} (x) = & \frac{(x - x_{0}) \cdot (x - x_{2}) \cdot (x - x_{3})}{(x_{1} - x_{0}) \cdot (x_{1} - x_{2}) \cdot (x_{1} - x_{3})} \\ = & \frac{(x - 1) \cdot (x - 4) \cdot (x - 5)}{(2 - 1) \cdot (2 - 4) \cdot (2 - 5)} = \frac{1}{6} (- 20 + 29 x - 10 x^{2} + x^{3}) \\ l_{2} (x) = & \frac{(x - x_{0}) \cdot (x - x_{1}) \cdot (x - x_{3})}{(x_{2} - x_{0}) \cdot (x_{2} - x_{1}) \cdot (x_{2} - x_{3})} \\ = & \frac{(x - 1) \cdot (x - 2) \cdot (x - 5)}{(4 - 1) \cdot (4 - 2) \cdot (4 - 5)} = \frac{- 1}{6} (- 10 + 17 x - 8 x^{2} + x^{3}) \\ l_{3} (x) = & \frac{(x - x_{0}) \cdot (x - x_{1}) \cdot (x - x_{2})}{(x_{3} - x_{0}) \cdot (x_{3} - x_{1}) \cdot (x_{3} - x_{2})} \\ = & \frac{(x - 1) \cdot (x - 2) \cdot (x - 4)}{(5 - 1) \cdot (5 - 2) \cdot (5 - 4)} = \frac{1}{12} (- 8 + 14 x - 7 x^{2} + x^{3}) \end{aligned}$

Així doncs el polinomi interpolador és:

$P_{3} (x) = 2 \cdot l_{1} (x) + 12 \cdot l_{2} (x) + 21 \cdot l_{3} (x) = \frac{1}{12} (- 8 + 2 x + 5 x^{2} + x^{3})$

D'aquesta manera $f (3) = \frac{35}{6} = 5.8333 \dots$