Interpolació polinòmica de Newton

L'interpolació polinòmica de Newton és un mètode interpolació polinòmica. Encara que només existeix un únic polinomi que interpola una sèrie de punts, hi ha diferents formes de calcular-lo. Aquest mètode és útil per a situacions que requereixen un nombre baix de punts per a interpolar, ja que a mesura que creix el nombre de punts, també ho fa el grau del polinomi.

Hi ha certs avantatges en l'ús d'aquest polinomi respecte al polinomi interpolador de Lagrange. Per exemple, si fos necessari afegir algun nou punt o nus a la funció, només caldria calcular aquest últim punt, donada la relació de recurrència existent i demostrada anteriorment.

Definició de pendent

El primer pas per a trobar la fórmula de la interpolació és definir el pendent d'ordre $n$ de manera recursiva:

• $f_0(x_i)$: terme i-èsim de la seqüència
• $f_1(x_0,x_1)=\frac{f_0(x_1)-f_0(x_0)}{x_1-x_0}$
• $f_2(x_0,x_1,x_2)=\frac{f_1(x_1,x_2)-f_1(x_0,x_1)}{x_2-x_0}$

En general:

$f_i(x_0,x_1,\dots ,x_{i-1},x_i)=\frac{f_{i-1}(x_1,\dots ,x_{i-1},x_i)-f_{i-1}(x_0,x_1,\dots ,x_{i-1})}{x_i-x_0}$,

on $x_i-x_j$ representa la distància entre dos elements (per exemple, es pot tenir l'element en $x=3$ i $x=5$ però desconéixer el valor de la seqüència en $x=4$).

Es pot apreciar com en la definició general s'utilitza el pendent del pas anterior, $f_{i-1}(x_1,\dots ,x_{i-1},x_i)$, a la qual se li resta el pendent previ del mateix ordre, és a dir, el subíndex dels termes es decrementa en $1$, com si es desplaçara, per a obtenir $f_{i-1}(x_0,x_1,\dots ,x_{i-1})$.

Cal notar també que encara que el terme inicial sempre siga $x_0$, aquest pot ser en realitat qualsevol altre, per exemple, es pot definir $f_1(x_{i-1},x_i)$ de manera anàloga al cas mostrat anteriorment.

Definició del polinomi

Una vegada coneixem el pendent, ja és possible definir el polinomi de grau $n$ de manera també recursiva:

• $p_0(x)=f_0(x_0)=y_0$. Es defineix així ja que aquest valor és l'únic que s'ajusta a la seqüència original per al primer terme.
• $p_1(x)=p_0(x)+f_1(x_0,x_1)*(x-x_0)$.
• $p_2(x)=p_1(x)+f_2(x_0,x_1,x_2)*(x-x_0)*(x-x_1)$.

En general:

$p_i(x)=p_{i-1}(x)+f_i(x_0,x_1,\dots ,x_{i-1},x_i)\prod _{j=0}^{i-1}(x-x_j)$

Exemples

Posem com a exemple la seqüència $f_0$ tal que $f_0(1)=6,f_0(2)=9,f_0(3)=2$ i $f_0(4)=5$, és a dir, són els termes $6,9,2,5$ per a $x_0=1$ fins a $x_3=4$.

S'obté les pendents d'ordre $1$ de la següent manera:

• $f_1(x_0,x_1)=\frac{f_0(x_1)-f_0(x_0)}{x_1-x_0}=\frac{9-6}{2-1}=3$
• $f_1(x_1,x_2)=\frac{f_0(x_2)-f_0(x_1)}{x_2-x_1}=\frac{2-9}{3-2}=-7$
• $f_1(x_2,x_3)=\frac{f_0(x_3)-f_0(x_2)}{x_3-x_2}=\frac{5-2}{4-3}=3$

Una vegada tenim les pendents d'ordre $1$, és possible obtenir les de següent ordre:

• $f_2(x_0,x_1,x_2)=\frac{f_1(x_1,x_2)-f_1(x_0,x_1)}{x_2-x_0}=\frac{-7-3}{3-1}=-5$
• $f_2(x_1,x_2,x_3)=\frac{f_1(x_2,x_3)-f_1(x_1,x_2)}{x_3-x_1}=\frac{3-(-7)}{4-2}=5$

Finalment, definim el pendent d'ordre $3$:

• $f_3(x_0,x_1,x_2,x_3)=\frac{f_2(x_1,x_2,x_3)-f_2(x_0,x_1,x_2)}{x_3-x_0}=\frac{5-(-5)}{4-1}=\frac{10}{3}$

Una vegada tenim el pendent, podem definir els successius polinomis:

• $p_0(x)=6$.
• $p_1(x)=6+3(x-1)$
• $p_2(x)=6+3(x-1)-5(x-1)(x-2)$.
• $p_3(x)=6+3(x-1)-5(x-1)(x-2)+\frac{10}{3}(x-1)(x-2)(x-3)$

Podem provar, per exemple, la interpolació lineal pel valor $1.5$, que resulta ser $p_1(1.5)=6+3(1.5-1)=7.5$. Efectivament, al ser una recta, podem veure que aquest valor és igual a $3+\frac{3+6}{2}=7.5$, el punt mitjà entre tots dos (més el punt inicial, $3$).