Polinomul de interpolare Lagrange

Polinomul de interpolare Lagrange pentru 4 puncte: ((−9, 5), (−4, 2), (−1, −2), (7, 9)), L(x), care este suma polinoamelor de bază scalate y₀l₀(x), y₁l₁(x), y₂l₂(x) and y₃l₃(x). Polinomul de interpolare trece prin toate cele 4 puncte, iar fiecare polinom scalat de bază trece prin punctul său de control respectiv, în cazul în care este 0, unde x corespunde celorlalte 3 puncte. Codul Matlab pentru acest exemplu este disponibil în .

Definiție

Fie un set de k + 1 puncte de date, diferite între ele:

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

Polinomul de interpolare Lagrange este combinația liniară

$L(x):=\sum _{j=0}^k{y}_j{\ell }_j(x)$

de polinoame Lagrange de bază

${\ell }_j(x):=\prod _{\begin{array}{c}0\le m\le k\\ m\ne j\end{array}}\frac{x-x_m}{x_j-x_m}=\frac{(x-x_0)}{(x_j-x_0)}\cdots \frac{(x-x_{j-1})}{(x_j-x_{j-1})}\frac{(x-x_{j+1})}{(x_j-x_{j+1})}\cdots \frac{(x-x_k)}{(x_j-x_k)}.$

Deși numit după Joseph Louis Lagrange în 1795, a fost descoperit pentru prima data în 1779 de către Edward Waring și a fost publicat în 1783 de Leonhard Euler.

Având în vedere ipoteza inițială că $x_i$ sunt diferite între ele, această expresie este întotdeauna bine definită.

Se verifică imediat că polinomul interpolează corect funcția, adică: $P(x_i)$=$y_{i}i=0,\dots ,k$, pentru orice i=1..n.

Exemple

1

Polinomul de interpolare al funcției tangentă

Să găsim o formulă de interpolare pentru funcția f(x) = tan(x) dată de următoarele seturi de valori:

$\begin{array}{cccccccc}{x}_0& =-1.5& & & & & f(x_0)& =-14.1014\\ x_1& =-0.75& & & & & f(x_1)& =-0.931596\\ x_2& =0& & & & & f(x_2)& =0\\ x_3& =0.75& & & & & f(x_3)& =0.931596\\ x_4& =1.5& & & & & f(x_4)& =14.1014.\end{array}$

Polinoamele de bază sunt:

${\ell }_0(x)=\frac{x-x_1}{x_0-x_1}\cdot \frac{x-x_2}{x_0-x_2}\cdot \frac{x-x_3}{x_0-x_3}\cdot \frac{x-x_4}{x_0-x_4}=\frac{1}{243}x(2x-3)(4x-3)(4x+3)$

${\ell }_1(x)=\frac{x-x_0}{x_1-x_0}\cdot \frac{x-x_2}{x_1-x_2}\cdot \frac{x-x_3}{x_1-x_3}\cdot \frac{x-x_4}{x_1-x_4}=-\frac{8}{243}x(2x-3)(2x+3)(4x-3)$

${\ell }_2(x)=\frac{x-x_0}{x_2-x_0}\cdot \frac{x-x_1}{x_2-x_1}\cdot \frac{x-x_3}{x_2-x_3}\cdot \frac{x-x_4}{x_2-x_4}=\frac{3}{243}(2x+3)(4x+3)(4x-3)(2x-3)$

${\ell }_3(x)=\frac{x-x_0}{x_3-x_0}\cdot \frac{x-x_1}{x_3-x_1}\cdot \frac{x-x_2}{x_3-x_2}\cdot \frac{x-x_4}{x_3-x_4}=-\frac{8}{243}x(2x-3)(2x+3)(4x+3)$

${\ell }_4(x)=\frac{x-x_0}{x_4-x_0}\cdot \frac{x-x_1}{x_4-x_1}\cdot \frac{x-x_2}{x_4-x_2}\cdot \frac{x-x_3}{x_4-x_3}=\frac{1}{243}x(2x+3)(4x-3)(4x+3).$

Deci polinomul de interpolare este:

$\begin{array}{cc}L(x)& =\frac{1}{243}(f(x_0)x(2x-3)(4x-3)(4x+3)\\ & -8f(x_1)x(2x-3)(2x+3)(4x-3)\\ & +3f(x_2)(2x+3)(4x+3)(4x-3)(2x-3)\\ & -8f(x_3)x(2x-3)(2x+3)(4x+3)\\ & +f(x_4)x(2x+3)(4x-3)(4x+3))\\ & =4.834848x^3-1.477474x.\end{array}$

2

Să interpolăm funcția f(x) = x² pe domeniul 1 ? x ? 3, prin următoarele 3 puncte:

$\begin{array}{cccccc}{x}_0& =1& & & f(x_0)& =1\\ x_1& =2& & & f(x_1)& =4\\ x_2& =3& & & f(x_2)& =9.\end{array}$

Polinomul este:

$\begin{array}{cc}L(x)& =1\cdot \frac{x-2}{1-2}\cdot \frac{x-3}{1-3}+4\cdot \frac{x-1}{2-1}\cdot \frac{x-3}{2-3}+9\cdot \frac{x-1}{3-1}\cdot \frac{x-2}{3-2}\\ & =x^2.\end{array}$

3

Să interpolăm funcția f(x) = x³ pe domeniul 1 < x < 3, prin punctele:

$x_0=1$	$f(x_0)=1$
$x_1=2$	$f(x_1)=8$
$x_2=3$	$f(x_2)=27$

Polinomul este:

$\begin{array}{cc}L(x)& =1\cdot \frac{x-2}{1-2}\cdot \frac{x-3}{1-3}+8\cdot \frac{x-1}{2-1}\cdot \frac{x-3}{2-3}+27\cdot \frac{x-1}{3-1}\cdot \frac{x-2}{3-2}\\ & =6x^2-11x+6.\end{array}$

Interpolarea baricentrică

Exemplu de divergență al polinomului de interpolare Lagrange

Forma Lagrange de interpolare polinomului arată caracterul liniar al polinomului de interpolare și unicitatea acestui polinom. De aceea, este de preferat în probe și argumente teoretice. Dar, după cum se poate observa din construcții, de fiecare dată când un nod x_k se modifică, toate polinoame Lagrange de bază trebuie să fie recalculate. O formă mai bună a polinomului de interpolare în practică este forma baricentrică de interpolare Lagrange formula Newton a polinomului.

Utilizând

$\ell (x)=(x-x_0)(x-x_1)\cdots (x-x_k)$

putem rescrie polinoamele de bază Lagrange ca

${\ell }_j(x)=\frac{\ell (x)}{x-x_j}\frac{1}{\prod _{i=0,i\ne j}^k(x_j-x_i)}$

sau, prin definirea ponderilor baricentrice

$w_j=\frac{1}{\prod _{i=0,i\ne j}^k(x_j-x_i)}$

putem scrie pur și simplu

${\ell }_j(x)=\ell (x)\frac{w_j}{x-x_j}$

care este denumit în mod obișnuit ca prima formă a formulei de interpolare baricentrică.

Avantajul este că această reprezentare polinomul de interpolare poate fi acum evaluat ca

$L(x)=\ell (x)\sum _{j=0}^k\frac{w_j}{x-x_j}{y}_j$

care, în cazul în care ponderile $w_j$ au fost pre-calculate, are nevoie doar de $mathcalO(n)$ (operații de evaluare $\ell (x)$ și ponderile $w_j/(x-x_j)$), spre deosebire de $𝒪(n^2)$ pentru evaluarea polinoamelor Lagrange de bază $el{l}_j(x)$ individual.

Formula de interpolare baricentrică poate fi, de asemenea, ușor de actualizat pentru a include un nod nou $x_{k+1}$ prin împărțirea nodurilor $w_j$, $j=0\dots k$ la$(x_j-x_{k+1})$ și construirea noului $w_{k+1}$ ca mai sus.

Putem simplifica și mai mult prima formă prin luarea în considerare prima interpolare baricentrică a funcției constante $g(x)\equiv 1$:

$g(x)=\ell (x)\sum _{j=0}^k\frac{w_j}{x-x_j}.$

Împărțirea $L(x)$ la $g(x)$ nu modifică interpolarea, dar conduce la rezultatul

$L(x)=\frac{\sum _{j=0}^k\frac{w_j}{x-x_j}{y}_j}{\sum _{j=0}^k\frac{w_j}{x-x_j}}$

care este menționat ca forma a doua sau adevarata forma a formulei de interpolare baricentrică. Această formă are avantajul că $\ell (x)$ nu trebuie să fie evaluate pentru fiecare evaluare a $L(x)$.

Referințe

Bibliografie

• Constantin Ilioi, Probleme de optimizare și algoritmi de aproximare a soluțiilor, Editura Academiei Republicii Socialiste România, București, 1980.

Legături externe

• www.utgjiu.ro/math/mbuneci/book/mn2009.pdf
• http://cs.upm.ro/~bela.finta/.files/carti/Numerika.pdf Format:Webarchive