Polinomul de interpolare Newton

În analiză numerică, polinomul de interpolare Newton, numit după inventatorul său Isaac Newton, este polinomul de interpolare, exprimat sub forma Newton, folosind diferențe divizate.

Definiție

Având un set de k + 1 puncte de date diferite între ele:

$(x_0,y_0),\dots ,(x_k,y_k)$

polinomul de interpolare în forma Newton este o combinație liniară din polinoamele Newton polinoame de bază

$N(x):=\sum _{j=0}^k{a}_j{n}_j(x)$

unde polinoamele Newton de bază sunt definite astfel:

$n_j(x):=\prod _{i=0}^{j-1}(x-x_i)$

pentru $j>0$ și $n_0(x)\equiv 1$. Coeficienții sunt definiți ca:

$a_j:=[y_0,\dots ,y_j]$

unde

$[y_0,\dots ,y_j]$

sunt diferențele divizate.

Astfel, polinomul Newton poate fi scris ca:

$N(x)=[y_0]+[y_0,y_1](x-x_0)+\cdots +[y_0,\dots ,y_k](x-x_0)(x-x_1)\cdots (x-x_{k-1}).$

"Polinomul Newton 'de mai sus poate fi exprimat într-o formă simplificată atunci când $x_0,x_1,\dots ,x_k$ sunt aranjate consecutiv, la distanțe egale. Introducând notația $h=x_{i+1}-x_i$ pentru fiecare $i=0,1,\dots ,k-1$ și $x=x_0+sh$, diferența $x-x_i$ poate fi scrisă ca $(s-i)h$. Deci, "polinomul Newton" de mai sus devine:

$N(x)=[y_0]+[y_0,y_1]sh+\cdots +[y_0,\dots ,y_k]s(s-1)\cdots (s-k+1)h^k$

$=\sum _{i=0}^{k}s(s-1)\cdots (s-i+1)h^i[y_0,\dots ,y_i]=\sum _{i=0}^k(\genfrac{}{}{0ex}{}{s}{i})i!h^i[y_0,\dots ,y_i]$

$N(x)=\sum _{i=0}^k(\genfrac{}{}{0ex}{}{s}{i})i!h^i[y_0,\dots ,y_i]$

se numește formula diferențelor divizate ale lui Newton".

În cazul în care nodurile sunt ordonate ca $x_k,x_{k-1},\dots ,x_0$, polinomul Newton devine:

$N(x)=[y_k]+[y_k,y_{k-1}](x-x_k)+\cdots +[y_k,\dots ,y_0](x-x_k)(x-x_{k-1})\cdots (x-x_1)$

Dacă$x_k,x_{k-1},\dots ,x_0$ sunt la fel de distanțate cu x=$x_k+sh$ and $x_i=x_k-(k-i)h$ for $i=0,1,\dots ,k$, atunci

$N(x)=[y_k]+[y_k,y_{k-1}]sh+\cdots +[y_k,\dots ,y_0]s(s+1)\cdots (s+k-1)h^k=\sum _{i=0}^k{(-1)}^i(\genfrac{}{}{0ex}{}{-s}{i})i!h^i[y_k,\dots ,y_{k-i}]$

$N(x)=\sum _{i=0}^k{(-1)}^i(\genfrac{}{}{0ex}{}{-s}{i})i!h^i[y_k,\dots ,y_{k-i}]$

se numește formula diferențelor divizate inversate ale lui Newton".

Exemplu

Diferențele divizate pot fi scrise în forma unui tabel. De exemplu, pentru o funcție $f$ este de a fi interpolate pe puncte $x_0,ldots,x_n$. Scriem

$\begin{array}{cccc}{x}_0& f(x_0)& & \\ & & \frac{f(x_1)-f(x_0)}{x_1-x_0}& \\ x_1& f(x_1)& & \frac{\frac{f(x_2)-f(x_1)}{x_2-x_1}-\frac{f(x_1)-f(x_0)}{x_1-x_0}}{x_2-x_0}\\ & & \frac{f(x_2)-f(x_1)}{x_2-x_1}& \\ x_2& f(x_2)& & ⋮\\ & & ⋮& \\ ⋮& & & ⋮\\ & & ⋮& \\ x_n& f(x_n)& & \end{array}$

Atunci polinomul de interpolare este format ca mai sus folosind mențiunile cel mai de sus din fiecare coloană ca coeficienți.

De exemplu, să presupunem că trebuie să construim polinomul de interpolare pentru $f(x)=\tan x$ folosind diferențele divizate, la punctele

$x_0=-1.5$	$x_1=-0.75$	$x_2=0$	$x_3=0.75$	$x_4=1.5$
$f(x_0)=-14.1014$	$f(x_1)=-0.931596$	$f(x_2)=0$	$f(x_3)=0.931596$	$f(x_4)=14.1014$

Pentru utilizarea unei precizii de 6 zecimale, vom construi tabelul

$\begin{array}{cccccc}-1.5& -14.1014& & & & \\ & & 17.5597& & & \\ -0.75& -0.931596& & -10.8784& & \\ & & 1.24213& & 4.83484& \\ 0& 0& & 0& & 0\\ & & 1.24213& & 4.83484& \\ 0.75& 0.931596& & 10.8784& & \\ & & 17.5597& & & \\ 1.5& 14.1014& & & & \end{array}$

Astfel, polinomul de interpolare este:$-14.1014+17.5597(x+1.5)-10.8784(x+1.5)(x+0.75)$

$+4.83484(x+1.5)(x+0.75)(x)+0(x+1.5)(x+0.75)(x)(x-0.75)$

$=-0.00005-1.4775x-0.00001x^2+4.83484x^3$

Având în vedere o acuratețe mai mare în tabel, coeficienții primul și al treilea vor fi egali cu zero.

Bibliografie

• Constantin Ilioi, Probleme de optimizare și algoritmi de aproximare a soluțiilor, Editura Academiei Republicii Socialiste România, București, 1980.
• www.utgjiu.ro/math/mbuneci/book/mn2009.pdf/ Metode numerice - Aspecte teoretice și practice, Mădălina Roxana Buneci, Editura Academică Brâncuși, Târgu Jiu, 2009
• http://cs.upm.ro/~bela.finta/.files/carti/Numerika.pdf Format:Webarchive
• www.vpetrehus.home.ro/Lectii_AN.pdf/ Lecții de analiză numerică, Viorel Petrehus, Universitatea Tehnică de Construcții București, 2010