Polinomul de interpolare Lagrange pentru 4 puncte: ((−9, 5) , (−4, 2) , (−1, −2) , (7, 9) ), L (x ), care este suma polinoamelor de bază scalate y0 l 0 (x ) , y1 l 1 (x ) , y2 l 2 (x ) and y3 l 3 (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 [1] Fie un set de k + 1 puncte de date, diferite între ele:
(
x
0
,
y
0
)
,
…
,
(
x
j
,
y
j
)
,
…
,
(
x
k
,
y
k
)
{\displaystyle (x_{0},y_{0}),\ldots ,(x_{j},y_{j}),\ldots ,(x_{k},y_{k})}
Polinomul de interpolare Lagrange este combinația liniară
L
(
x
)
:=
∑
j
=
0
k
y
j
ℓ
j
(
x
)
{\displaystyle L(x):=\sum _{j=0}^{k}y_{j}\ell _{j}(x)}
de polinoame Lagrange de bază
ℓ
j
(
x
)
:=
∏
0
≤
m
≤
k
m
≠
j
x
−
x
m
x
j
−
x
m
=
(
x
−
x
0
)
(
x
j
−
x
0
)
⋯
(
x
−
x
j
−
1
)
(
x
j
−
x
j
−
1
)
(
x
−
x
j
+
1
)
(
x
j
−
x
j
+
1
)
⋯
(
x
−
x
k
)
(
x
j
−
x
k
)
.
{\displaystyle \ell _{j}(x):=\prod _{\begin{smallmatrix}0\leq m\leq k\\m\neq j\end{smallmatrix}}{\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
{\displaystyle x_{i}}
bine definită .
Se verifică imediat că polinomul interpolează corect funcția, adică:
P
(
x
i
)
{\displaystyle P(x_{i})}
y
i
i
=
0
,
…
,
k
{\displaystyle y_{i}\qquad i=0,\ldots ,k}
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:
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.
{\displaystyle {\begin{aligned}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{aligned}}}
Polinoamele de bază sunt:
ℓ
0
(
x
)
=
x
−
x
1
x
0
−
x
1
⋅
x
−
x
2
x
0
−
x
2
⋅
x
−
x
3
x
0
−
x
3
⋅
x
−
x
4
x
0
−
x
4
=
1
243
x
(
2
x
−
3
)
(
4
x
−
3
)
(
4
x
+
3
)
{\displaystyle \ell _{0}(x)={x-x_{1} \over x_{0}-x_{1}}\cdot {x-x_{2} \over x_{0}-x_{2}}\cdot {x-x_{3} \over x_{0}-x_{3}}\cdot {x-x_{4} \over x_{0}-x_{4}}={1 \over 243}x(2x-3)(4x-3)(4x+3)}
ℓ
1
(
x
)
=
x
−
x
0
x
1
−
x
0
⋅
x
−
x
2
x
1
−
x
2
⋅
x
−
x
3
x
1
−
x
3
⋅
x
−
x
4
x
1
−
x
4
=
−
8
243
x
(
2
x
−
3
)
(
2
x
+
3
)
(
4
x
−
3
)
{\displaystyle \ell _{1}(x)={x-x_{0} \over x_{1}-x_{0}}\cdot {x-x_{2} \over x_{1}-x_{2}}\cdot {x-x_{3} \over x_{1}-x_{3}}\cdot {x-x_{4} \over x_{1}-x_{4}}={}-{8 \over 243}x(2x-3)(2x+3)(4x-3)}
ℓ
2
(
x
)
=
x
−
x
0
x
2
−
x
0
⋅
x
−
x
1
x
2
−
x
1
⋅
x
−
x
3
x
2
−
x
3
⋅
x
−
x
4
x
2
−
x
4
=
3
243
(
2
x
+
3
)
(
4
x
+
3
)
(
4
x
−
3
)
(
2
x
−
3
)
{\displaystyle \ell _{2}(x)={x-x_{0} \over x_{2}-x_{0}}\cdot {x-x_{1} \over x_{2}-x_{1}}\cdot {x-x_{3} \over x_{2}-x_{3}}\cdot {x-x_{4} \over x_{2}-x_{4}}={3 \over 243}(2x+3)(4x+3)(4x-3)(2x-3)}
ℓ
3
(
x
)
=
x
−
x
0
x
3
−
x
0
⋅
x
−
x
1
x
3
−
x
1
⋅
x
−
x
2
x
3
−
x
2
⋅
x
−
x
4
x
3
−
x
4
=
−
8
243
x
(
2
x
−
3
)
(
2
x
+
3
)
(
4
x
+
3
)
{\displaystyle \ell _{3}(x)={x-x_{0} \over x_{3}-x_{0}}\cdot {x-x_{1} \over x_{3}-x_{1}}\cdot {x-x_{2} \over x_{3}-x_{2}}\cdot {x-x_{4} \over x_{3}-x_{4}}=-{8 \over 243}x(2x-3)(2x+3)(4x+3)}
ℓ
4
(
x
)
=
x
−
x
0
x
4
−
x
0
⋅
x
−
x
1
x
4
−
x
1
⋅
x
−
x
2
x
4
−
x
2
⋅
x
−
x
3
x
4
−
x
3
=
1
243
x
(
2
x
+
3
)
(
4
x
−
3
)
(
4
x
+
3
)
.
{\displaystyle \ell _{4}(x)={x-x_{0} \over x_{4}-x_{0}}\cdot {x-x_{1} \over x_{4}-x_{1}}\cdot {x-x_{2} \over x_{4}-x_{2}}\cdot {x-x_{3} \over x_{4}-x_{3}}={1 \over 243}x(2x+3)(4x-3)(4x+3).}
Deci polinomul de interpolare este:
L
(
x
)
=
1
243
(
f
(
x
0
)
x
(
2
x
−
3
)
(
4
x
−
3
)
(
4
x
+
3
)
−
8
f
(
x
1
)
x
(
2
x
−
3
)
(
2
x
+
3
)
(
4
x
−
3
)
+
3
f
(
x
2
)
(
2
x
+
3
)
(
4
x
+
3
)
(
4
x
−
3
)
(
2
x
−
3
)
−
8
f
(
x
3
)
x
(
2
x
−
3
)
(
2
x
+
3
)
(
4
x
+
3
)
+
f
(
x
4
)
x
(
2
x
+
3
)
(
4
x
−
3
)
(
4
x
+
3
)
)
=
4.834848
x
3
−
1.477474
x
.
{\displaystyle {\begin{aligned}L(x)&={1 \over 243}{\Big (}f(x_{0})x(2x-3)(4x-3)(4x+3)\\&{}\qquad {}-8f(x_{1})x(2x-3)(2x+3)(4x-3)\\&{}\qquad {}+3f(x_{2})(2x+3)(4x+3)(4x-3)(2x-3)\\&{}\qquad {}-8f(x_{3})x(2x-3)(2x+3)(4x+3)\\&{}\qquad {}+f(x_{4})x(2x+3)(4x-3)(4x+3){\Big )}\\&=4.834848x^{3}-1.477474x.\end{aligned}}}
Să interpolăm funcția f (x ) = x 2 pe domeniul 1 ? x ? 3, prin următoarele 3 puncte:
x
0
=
1
f
(
x
0
)
=
1
x
1
=
2
f
(
x
1
)
=
4
x
2
=
3
f
(
x
2
)
=
9.
{\displaystyle {\begin{aligned}x_{0}&=1&&&f(x_{0})&=1\\x_{1}&=2&&&f(x_{1})&=4\\x_{2}&=3&&&f(x_{2})&=9.\end{aligned}}}
Polinomul este:
L
(
x
)
=
1
⋅
x
−
2
1
−
2
⋅
x
−
3
1
−
3
+
4
⋅
x
−
1
2
−
1
⋅
x
−
3
2
−
3
+
9
⋅
x
−
1
3
−
1
⋅
x
−
2
3
−
2
=
x
2
.
{\displaystyle {\begin{aligned}L(x)&={1}\cdot {x-2 \over 1-2}\cdot {x-3 \over 1-3}+{4}\cdot {x-1 \over 2-1}\cdot {x-3 \over 2-3}+{9}\cdot {x-1 \over 3-1}\cdot {x-2 \over 3-2}\\[10pt]&=x^{2}.\end{aligned}}}
Să interpolăm funcția f (x ) = x 3 pe domeniul 1 < x < 3, prin punctele:
x
0
=
1
{\displaystyle x_{0}=1\,}
f
(
x
0
)
=
1
{\displaystyle f(x_{0})=1\,}
x
1
=
2
{\displaystyle x_{1}=2\,}
f
(
x
1
)
=
8
{\displaystyle f(x_{1})=8\,}
x
2
=
3
{\displaystyle x_{2}=3\,}
f
(
x
2
)
=
27
{\displaystyle f(x_{2})=27\,}
Polinomul este:
L
(
x
)
=
1
⋅
x
−
2
1
−
2
⋅
x
−
3
1
−
3
+
8
⋅
x
−
1
2
−
1
⋅
x
−
3
2
−
3
+
27
⋅
x
−
1
3
−
1
⋅
x
−
2
3
−
2
=
6
x
2
−
11
x
+
6.
{\displaystyle {\begin{aligned}L(x)&={1}\cdot {x-2 \over 1-2}\cdot {x-3 \over 1-3}+{8}\cdot {x-1 \over 2-1}\cdot {x-3 \over 2-3}+{27}\cdot {x-1 \over 3-1}\cdot {x-2 \over 3-2}\\[8pt]&=6x^{2}-11x+6.\end{aligned}}}
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
Utilizând
ℓ
(
x
)
=
(
x
−
x
0
)
(
x
−
x
1
)
⋯
(
x
−
x
k
)
{\displaystyle \ell (x)=(x-x_{0})(x-x_{1})\cdots (x-x_{k})}
putem rescrie polinoamele de bază Lagrange ca
ℓ
j
(
x
)
=
ℓ
(
x
)
x
−
x
j
1
∏
i
=
0
,
i
≠
j
k
(
x
j
−
x
i
)
{\displaystyle \ell _{j}(x)={\frac {\ell (x)}{x-x_{j}}}{\frac {1}{\prod _{i=0,i\neq j}^{k}(x_{j}-x_{i})}}}
sau, prin definirea ponderilor baricentrice [2]
w
j
=
1
∏
i
=
0
,
i
≠
j
k
(
x
j
−
x
i
)
{\displaystyle w_{j}={\frac {1}{\prod _{i=0,i\neq j}^{k}(x_{j}-x_{i})}}}
putem scrie pur și simplu
ℓ
j
(
x
)
=
ℓ
(
x
)
w
j
x
−
x
j
{\displaystyle \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
)
=
ℓ
(
x
)
∑
j
=
0
k
w
j
x
−
x
j
y
j
{\displaystyle L(x)=\ell (x)\sum _{j=0}^{k}{\frac {w_{j}}{x-x_{j}}}y_{j}}
care, în cazul în care ponderile
w
j
{\displaystyle w_{j}}
m
a
t
h
c
a
l
O
(
n
)
{\displaystyle \ mathcalO(n)}
ℓ
(
x
)
{\displaystyle \ell (x)}
w
j
/
(
x
−
x
j
)
{\displaystyle w_{j}/(x-x_{j})}
O
(
n
2
)
{\displaystyle {\mathcal {O}}(n^{2})}
e
l
l
j
(
x
)
{\displaystyle \ ell_{j}(x)}
Formula de interpolare baricentrică poate fi, de asemenea, ușor de actualizat pentru a include un nod nou
x
k
+
1
{\displaystyle x_{k+1}}
w
j
{\displaystyle w_{j}}
j
=
0
…
k
{\displaystyle j=0\dots k}
(
x
j
−
x
k
+
1
)
{\displaystyle (x_{j}-x_{k+1})}
w
k
+
1
{\displaystyle w_{k+1}}
Putem simplifica și mai mult prima formă prin luarea în considerare prima interpolare baricentrică a funcției constante
g
(
x
)
≡
1
{\displaystyle g(x)\equiv 1}
g
(
x
)
=
ℓ
(
x
)
∑
j
=
0
k
w
j
x
−
x
j
.
{\displaystyle g(x)=\ell (x)\sum _{j=0}^{k}{\frac {w_{j}}{x-x_{j}}}.}
Împărțirea
L
(
x
)
{\displaystyle L(x)}
g
(
x
)
{\displaystyle g(x)}
L
(
x
)
=
∑
j
=
0
k
w
j
x
−
x
j
y
j
∑
j
=
0
k
w
j
x
−
x
j
{\displaystyle 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ă
ℓ
(
x
)
{\displaystyle \ell (x)}
L
(
x
)
{\displaystyle L(x)}
Constantin Ilioi, Probleme de optimizare și algoritmi de aproximare a soluțiilor , Editura Academiei Republicii Socialiste România, București, 1980. 
![{\displaystyle {\begin{aligned}L(x)&={1}\cdot {x-2 \over 1-2}\cdot {x-3 \over 1-3}+{4}\cdot {x-1 \over 2-1}\cdot {x-3 \over 2-3}+{9}\cdot {x-1 \over 3-1}\cdot {x-2 \over 3-2}\\[10pt]&=x^{2}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/274ccc9dab43ec4880e953e43ebcd5a8470cf12a)
![{\displaystyle {\begin{aligned}L(x)&={1}\cdot {x-2 \over 1-2}\cdot {x-3 \over 1-3}+{8}\cdot {x-1 \over 2-1}\cdot {x-3 \over 2-3}+{27}\cdot {x-1 \over 3-1}\cdot {x-2 \over 3-2}\\[8pt]&=6x^{2}-11x+6.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/14d94ef9be0953a87421f693085e703fb618de34)
