1.
Interpolació
n de Lagrange
Sean 𝑰 = [𝒂, 𝒃], 𝒙𝟎 < 𝒙𝟏 < . . . < 𝒙𝒏 , puntos de I y f una funciónn definida sobre I .
Problemá
tica 1.1. Encontrar un polinomio P de grado menor o igual a n tal que
𝑓 (𝑥𝑖 ) = 𝑃 (𝑥𝑖 ), 𝑖 = 0, . . . , 𝑛
Definición 1.2. Los polinomios de Lagrange asociados a los puntos x i, 𝑖 = 0, . . . , 𝑛, son los
polinomios de grado n definidos por
𝑛
𝑥 − 𝑥𝑗
𝐿𝑖(𝑥) = ∏ ( )
𝑥𝑖 − 𝑥𝑗
𝑗=0
𝑗≠𝑖
Nota
1. Para todo 𝑖 = 0, . . . , 𝑛, Li es un polinomio de grado n.
2. 𝐿𝑖 tiene n raíces reales que son 𝑥0, . . . , 𝑥𝑖 − 1, 𝑥𝑖 + 1, . . . , 𝑥𝑛.
3. Para todo 𝑖, 𝑗 = 0, . . . , 𝑛
𝐿𝑖(𝑥𝑗 ) = 𝛿𝑖𝑗
Donde
donde 𝛿𝑖𝑗 es el símbolo de Kronecker definido por es el símbolo de Kronecker definido por
0, 𝑠𝑖 𝑖 ≠ 𝑗
𝛿𝑖𝑗 = {
1, 𝑠𝑖 𝑖 = 𝑗
1.1. Existencia y unicidad del polinomio de interpolación
Teorema 1.
Dados n + 1 puntos 𝑥𝑖, 𝑖 = 0, . . . , 𝑛, distintos de un intervalo 𝐼 = [𝑎, 𝑏] y 𝑛 + 1 números reales
𝑦𝑖, 𝑖 = 0, . . . , 𝑛, entonces existe un único polinomio Pn de grado menor o igual a n tal que
𝑃𝑛(𝑥𝑖) = 𝑦𝑖 𝑖 = 0, . . . , 𝑛.
Nota 1.3.
En este teorema los números reales 𝑦𝑖, 𝑖 = 0, . . . , 𝑛, representan los valores 𝑓(𝑥𝑖), 𝑖 = 0, . . . , 𝑛,
cuando se trata de la interpolación de la función f y en este caso Pn se llama el polinomio de
interpolación de lagrange de f en los puntos 𝑥𝑖, 𝑖 = 0, . . . , 𝑛.
Demostración. Sea
𝑛
𝑃𝑛 (𝑥) = ∑ 𝑦𝑖 𝐿𝑖 (𝑥)
𝑖=0
𝑃𝑛 (𝑥𝑗 ) = ∑ 𝑦𝑖 𝐿𝑖 (𝑥𝑗 )
𝑖=0
= ∑ 𝑦𝑖 𝛿𝑖𝑗
𝑖=0
= 𝑦𝑗
Es decir Pn interpola f en los puntos 𝑥0 , 𝑥1 , . . . , 𝑥𝑛 .
Mostremos la unicidad. Designamos
◦ d ◦ ◦P el grado de P . Supongamos que Pn y Qn son
dos polinomios tales que 𝑑 𝑃𝑛 ≤ 𝑛 𝑑 𝑄𝑛 ≤ 𝑛 𝑦 𝑃𝑛 (𝑥𝑖 ) = 𝑄𝑛(𝑥𝑖 ), 𝑖 = 0, . . . , 𝑛. Sea el
,
polinomio
𝑊𝑛 (𝑥) = 𝑃𝑛 (𝑥) − 𝑄𝑛 (𝑥)
◦
Notamos que 𝑑 𝑊𝑛 ≤ 𝑛 𝑦 𝑊 (𝑥𝑖 ) = 0 para todo 𝑖 = 0, . . . , 𝑛. Así, Wn es un polinomio
de grado menor o igual a n con (𝑛 + 1) raíces, lo que es imposible a menos que 𝑊𝑛 (𝑥 ) ≡
0 entonces, 𝑄𝑛 (𝑥) = 𝑃𝑛 (𝑥 ) y 𝑃𝑛 es ú
n ico.
Nota
En el teorema anterior el polinomio de interpolación de Lagrange no se presenta en la
forma −
𝑃𝑛 (𝑥) = 𝑎𝑛 𝑥𝑛 + 𝑎𝑛 − 1 𝑥𝑛 1 + . . . + 𝑎0
Eso representa una dificultad numérica mayor por necesitarse muchas operaciones para
calcular Pn (α) para un α real dado. Existen algoritmos que permiten la reducción de los
el cálculos para evaluar Pn (α) .
2. Interpolació
n de Hermite
Vimos que con el conocimiento de las parejas (xi , f (xi )), xi ∈ I = [a, b], i = 0, . . . , n, se
puede hallar el polinomio de interpolació n de grado menor o igual a n. Supongamos ahora
que se quiere encontrar un polinomio “más parecido” a f usando información adicional sobre
f en los puntos xi , i = 0, . . . , n.
Lo má s natural en este caso es pensar en las derivadas de f en los puntos 𝑥𝑖 , 𝑖 =
0, . . . , 𝑛.
El problema se plantea de la siguiente manera. Sean I = [a, b] y x0 < x1 <. . . < xn (n +
1) puntos en I . Sean k0 , k1 , k2 , . . . , kn , n + 1 enteros no negativos y f una función definida
en I tal que los nú m eros
)
𝑓 𝑗 (𝑥𝑖 ), 𝑖 = 0, . . . , 𝑛, = 0, . . . , 𝑘𝑖,
son definidos y conocidos. Denotemos igualmente 𝑚 = 𝑛 + 𝑘0 + . . . + 𝑘𝑛 .
Nota 1.
𝑓𝑗𝑥𝑖 ) denota la j− ésima derivada de f en el punto xi y f (0) = f .
Problemá tica 1.2. Encontrar un polinomio Pm de grado menor o igual a m tal que
Pm (j) (xi ) = f (j) (xi ), i = 0, . . . , n, j = 0, . . . , ki . (1.1)
Definición .Un polinomio como (1.1) se llama polinomio de interpolación de Hermite de la
función f en los puntos xi , i = 0, . . . , n, con respecto a los enteros ki , i = 0, . . . , n.
Teorema. El polinomio de interpolació
n de Hermite existe y es ú
n ico.
Antes de demostrar este teorema introducimos la noción de raíz mú
l tiple y miramos su
relación con las derivadas de polinomios.
Definición . Sea P (x) un polinomio con coeficientes en K (K=R o C) y sean α ∈ K un
escalar y m ∈ N. Se dice que α es una raíz de P de multiplicidad m si (x − α)m divide a P
(x) pero (x − α)m+1 no lo divide.
En otras palabras, α es raíz de multiplicidad m de P si P (x) se escribe en la forma
P (x) = (x − α)m Q(x)
con 𝑄(𝛼) 0.
