Unidad 6: Interpolación

Carlos Alliera (calliera@[Link])

17 de mayo de 2021

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 1 / 15

 Introducción

Problema
Tenemos n + 1 puntos en el plano (x0 , y0 ), ..., (xn , yn ) tales que x0 < x1 < ... < xn
queremos construir un polinomio P ∈ Kn [x] tal que

P (xk ) = yk , 0≤k≤n

Al polinomio que verifica lo pedido se lo llama Polinomio interpolador.

Visto de otra forma, buscamos los coeficientes ak ∈ K de

P (x) = a0 + a1 x + a2 x2 + · · · + an−1 xn−1 + an xn

que pueden verse como el resultado del siguiente sistema lineal:

a0 + a1 x0 + · · · + an xn

 0 = y0
 a + a x + · · · + a xn = y

 0 1 1 n 1 1
..


 .
a0 + a1 xn + · · · + an xn

n = yn

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 2 / 15

 Introducción

Problema
Tenemos n + 1 puntos en el plano (x0 , y0 ), ..., (xn , yn ) tales que x0 < x1 < ... < xn
queremos construir un polinomio P ∈ Kn [x] tal que

P (xk ) = yk , 0≤k≤n

Al polinomio que verifica lo pedido se lo llama Polinomio interpolador.

Visto de otra forma, buscamos los coeficientes ak ∈ K de

P (x) = a0 + a1 x + a2 x2 + · · · + an−1 xn−1 + an xn

que pueden verse como el resultado del siguiente sistema lineal:

a0 + a1 x0 + · · · + an xn

 0 = y0
 a + a x + · · · + a xn = y

 0 1 1 n 1 1
..


 .
a0 + a1 xn + · · · + an xn

n = yn

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 2 / 15

 Introducción

Problema
Tenemos n + 1 puntos en el plano (x0 , y0 ), ..., (xn , yn ) tales que x0 < x1 < ... < xn
queremos construir un polinomio P ∈ Kn [x] tal que

P (xk ) = yk , 0≤k≤n

Al polinomio que verifica lo pedido se lo llama Polinomio interpolador.

Visto de otra forma, buscamos los coeficientes ak ∈ K de

P (x) = a0 + a1 x + a2 x2 + · · · + an−1 xn−1 + an xn

que pueden verse como el resultado del siguiente sistema lineal:

a0 + a1 x0 + · · · + an xn

 0 = y0
 a + a x + · · · + a xn = y

 0 1 1 n 1 1
..


 .
a0 + a1 xn + · · · + an xn

n = yn

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 2 / 15

 Introducción

Problema
Tenemos n + 1 puntos en el plano (x0 , y0 ), ..., (xn , yn ) tales que x0 < x1 < ... < xn
queremos construir un polinomio P ∈ Kn [x] tal que

P (xk ) = yk , 0≤k≤n

Al polinomio que verifica lo pedido se lo llama Polinomio interpolador.

Visto de otra forma, buscamos los coeficientes ak ∈ K de

P (x) = a0 + a1 x + a2 x2 + · · · + an−1 xn−1 + an xn

que pueden verse como el resultado del siguiente sistema lineal:

a0 + a1 x0 + · · · + an xn

 0 = y0
 a + a x + · · · + a xn = y

 0 1 1 n 1 1
..


 .
a0 + a1 xn + · · · + an xn

n = yn

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 2 / 15

 Introducción

Problema
Tenemos n + 1 puntos en el plano (x0 , y0 ), ..., (xn , yn ) tales que x0 < x1 < ... < xn
queremos construir un polinomio P ∈ Kn [x] tal que

P (xk ) = yk , 0≤k≤n

Al polinomio que verifica lo pedido se lo llama Polinomio interpolador.

Visto de otra forma, buscamos los coeficientes ak ∈ K de

P (x) = a0 + a1 x + a2 x2 + · · · + an−1 xn−1 + an xn

que pueden verse como el resultado del siguiente sistema lineal:

a0 + a1 x0 + · · · + an xn

 0 = y0
 a + a x + · · · + a xn = y

 0 1 1 n 1 1
..


 .
a0 + a1 xn + · · · + an xn

n = yn

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 2 / 15

 Método de Coeficientes Indeterminados

La Matriz de Vandermonde

Para resolver el problema antes planteado, consideramos un sistema lineal de la

forma V · −
→
a = y donde V ∈ K(n+1)×(n+1) , a, y ∈ Kn+1 que tiene esta forma:

1 x0 · · · xn
á ë á ë á ë
0 a0 y0
1 x1 · · · xn 1 a1 y1
.. .. .. .. · .. = ..
. . . . . .
1 xn · · · xn n an yn

donde la matriz V se conoce como Matriz de Vandermonde1

Propiedad
El sistema V · −
→a =− →y tiene solución unica cuando los xk son todos distintos. Es
decir, V es inversible si xk 6= xj si k 6= j.

1
Músico y químico frances 1735-1796
Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 3 / 15
 Método de Coeficientes Indeterminados

La Matriz de Vandermonde

Para resolver el problema antes planteado, consideramos un sistema lineal de la

forma V · −
→
a = y donde V ∈ K(n+1)×(n+1) , a, y ∈ Kn+1 que tiene esta forma:

1 x0 · · · xn
á ë á ë á ë
0 a0 y0
1 x1 · · · xn 1 a1 y1
.. .. .. .. · .. = ..
. . . . . .
1 xn · · · xn n an yn

donde la matriz V se conoce como Matriz de Vandermonde1

Propiedad
El sistema V · −
→a =− →y tiene solución unica cuando los xk son todos distintos. Es
decir, V es inversible si xk 6= xj si k 6= j.

1
Músico y químico frances 1735-1796
Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 3 / 15
 Método de Coeficientes Indeterminados

La Matriz de Vandermonde

Para resolver el problema antes planteado, consideramos un sistema lineal de la

forma V · −
→
a = y donde V ∈ K(n+1)×(n+1) , a, y ∈ Kn+1 que tiene esta forma:

1 x0 · · · xn
á ë á ë á ë
0 a0 y0
1 x1 · · · xn 1 a1 y1
.. .. .. .. · .. = ..
. . . . . .
1 xn · · · xn n an yn

donde la matriz V se conoce como Matriz de Vandermonde1

Propiedad
El sistema V · −
→a =− →y tiene solución unica cuando los xk son todos distintos. Es
decir, V es inversible si xk 6= xj si k 6= j.

1
Músico y químico frances 1735-1796
Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 3 / 15
 Método de Coeficientes Indeterminados

La Matriz de Vandermonde

Para resolver el problema antes planteado, consideramos un sistema lineal de la

forma V · −
→
a = y donde V ∈ K(n+1)×(n+1) , a, y ∈ Kn+1 que tiene esta forma:

1 x0 · · · xn
á ë á ë á ë
0 a0 y0
1 x1 · · · xn 1 a1 y1
.. .. .. .. · .. = ..
. . . . . .
1 xn · · · xn n an yn

donde la matriz V se conoce como Matriz de Vandermonde1

Propiedad
El sistema V · −
→a =− →y tiene solución unica cuando los xk son todos distintos. Es
decir, V es inversible si xk 6= xj si k 6= j.

1
Músico y químico frances 1735-1796
Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 3 / 15
 Método de Coeficientes Indeterminados

Ejemplo 1

Calcule el polinomio de menor grado que interpole los puntos:

(−1; −7), (1; 3), (2; 5)

Al tratarse de 3 puntos, entonces buscamos un polinomio de grado a lo sumo 2.

Directamente planteamos el modelo matricial V · −
→
a =−→y:
Ñ é Ñ é Ñ é
1 −1 1 a0 −7
1 1 1 · a1 = 3
1 2 4 a2 5

Al resolver, se obtiene:

a0 = −1, a1 = 5, a2 = −1

∴ P (x) = −1 + 5x − x2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 4 / 15

 Método de Coeficientes Indeterminados

Ejemplo 1

Calcule el polinomio de menor grado que interpole los puntos:

(−1; −7), (1; 3), (2; 5)

Al tratarse de 3 puntos, entonces buscamos un polinomio de grado a lo sumo 2.

Directamente planteamos el modelo matricial V · −
→
a =−→y:
Ñ é Ñ é Ñ é
1 −1 1 a0 −7
1 1 1 · a1 = 3
1 2 4 a2 5

Al resolver, se obtiene:

a0 = −1, a1 = 5, a2 = −1

∴ P (x) = −1 + 5x − x2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 4 / 15

 Método de Coeficientes Indeterminados

Ejemplo 1

Calcule el polinomio de menor grado que interpole los puntos:

(−1; −7), (1; 3), (2; 5)

Al tratarse de 3 puntos, entonces buscamos un polinomio de grado a lo sumo 2.

Directamente planteamos el modelo matricial V · −
→
a =−→y:
Ñ é Ñ é Ñ é
1 −1 1 a0 −7
1 1 1 · a1 = 3
1 2 4 a2 5

Al resolver, se obtiene:

a0 = −1, a1 = 5, a2 = −1

∴ P (x) = −1 + 5x − x2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 4 / 15

 Método de Coeficientes Indeterminados

Ejemplo 1

Calcule el polinomio de menor grado que interpole los puntos:

(−1; −7), (1; 3), (2; 5)

Al tratarse de 3 puntos, entonces buscamos un polinomio de grado a lo sumo 2.

Directamente planteamos el modelo matricial V · −
→
a =−→y:
Ñ é Ñ é Ñ é
1 −1 1 a0 −7
1 1 1 · a1 = 3
1 2 4 a2 5

Al resolver, se obtiene:

a0 = −1, a1 = 5, a2 = −1

∴ P (x) = −1 + 5x − x2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 4 / 15

 Método de Coeficientes Indeterminados

Ejemplo 1

Calcule el polinomio de menor grado que interpole los puntos:

(−1; −7), (1; 3), (2; 5)

Al tratarse de 3 puntos, entonces buscamos un polinomio de grado a lo sumo 2.

Directamente planteamos el modelo matricial V · −
→
a =−→y:
Ñ é Ñ é Ñ é
1 −1 1 a0 −7
1 1 1 · a1 = 3
1 2 4 a2 5

Al resolver, se obtiene:

a0 = −1, a1 = 5, a2 = −1

∴ P (x) = −1 + 5x − x2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 4 / 15

 Método de Coeficientes Indeterminados

Ejemplo 2: Condiciones sobre la derivada

Hallar el polinomio de menor grado tal que

P (−1) = 13, P (2) = 11, P (3) = 51, P 0 (1) = 8

En este caso buscamos un polinomio de grado a lo sumo 3 que cumpla lo pedido. Si

se plantea:
P (X) = a0 + a1 x + a2 x2 + a3 x3
nos queda el siguiente sistema de ecuaciones lineales:


 a0 − a1 + a2 − a3 = 13
a0 + 2a1 + 4a2 + 8a3 = 11


 a0 + 3a1 + 9a2 + 27a3 = 51
a1 + 2a2 + 3a3 = 8


Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 5 / 15

 Método de Coeficientes Indeterminados

Ejemplo 2: Condiciones sobre la derivada

Hallar el polinomio de menor grado tal que

P (−1) = 13, P (2) = 11, P (3) = 51, P 0 (1) = 8

En este caso buscamos un polinomio de grado a lo sumo 3 que cumpla lo pedido. Si

se plantea:
P (X) = a0 + a1 x + a2 x2 + a3 x3
nos queda el siguiente sistema de ecuaciones lineales:

 a0 − a1 + a2 − a3 = 13

a0 + 2a1 + 4a2 + 8a3 = 11


 a0 + 3a1 + 9a2 + 27a3 = 51
a1 + 2a2 + 3a3 = 8


Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 5 / 15

 Método de Coeficientes Indeterminados

Ejemplo 2

En este caso nos queda el siguiente planteo matricial:

Ü ê Ü ê Ü ê
1 −1 1 −1 a0 13
1 2 4 8 a1 11
· =
1 3 9 27 a2 51
0 1 2 3 a3 8

tras aplicar un método de resolución de nuestra preferencia, se tiene:

Ü ê Ü ê
a0 −9
a1 2
= ⇒ P (x) = 2x3 + 2x − 9
a2 0
a3 2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 6 / 15

 Método de Coeficientes Indeterminados

Ejemplo 2

En este caso nos queda el siguiente planteo matricial:

Ü ê Ü ê Ü ê
1 −1 1 −1 a0 13
1 2 4 8 a1 11
· =
1 3 9 27 a2 51
0 1 2 3 a3 8

tras aplicar un método de resolución de nuestra preferencia, se tiene:

Ü ê Ü ê
a0 −9
a1 2
= ⇒ P (x) = 2x3 + 2x − 9
a2 0
a3 2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 6 / 15

 Interpolación de Lagrange

Teorema
Dados x0 , x1 , ..., xn ∈ K distintos y y0 , y1 , ..., yn ∈ K. Existe un único
polinomio de grado a lo sumo n tal que

P (xj ) = yj ∀ j = 0, ..., n

Base de Lagrange
Dados n + 1 pares (xj , yj ), j = 0, 1, ..., n se definen los polinomios `j
tales que
`j (xj ) = 1 & `j (xk ) = 0 si k 6= j

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 7 / 15

 Interpolación de Lagrange

Teorema
Dados x0 , x1 , ..., xn ∈ K distintos y y0 , y1 , ..., yn ∈ K. Existe un único
polinomio de grado a lo sumo n tal que

P (xj ) = yj ∀ j = 0, ..., n

Base de Lagrange
Dados n + 1 pares (xj , yj ), j = 0, 1, ..., n se definen los polinomios `j
tales que
`j (xj ) = 1 & `j (xk ) = 0 si k 6= j

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 7 / 15

 Interpolación de Lagrange

Teorema
Dados x0 , x1 , ..., xn ∈ K distintos y y0 , y1 , ..., yn ∈ K. Existe un único
polinomio de grado a lo sumo n tal que

P (xj ) = yj ∀ j = 0, ..., n

Base de Lagrange
Dados n + 1 pares (xj , yj ), j = 0, 1, ..., n se definen los polinomios `j
tales que
`j (xj ) = 1 & `j (xk ) = 0 si k 6= j

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 7 / 15

 Interpolación de Lagrange

Base de Lagrange

Así nos quedan n + 1 polinomios de grado n definidos de esta manera:

n
Q
(x − xk )
k6=j
`j (x) = n
Q
(xj − xk )
k6=j

Estos polinomios forman la Base de Lagrange:

B := {`0 , `1 , ..., `n }

y dependen de las xj de los puntos a interpolar.

Polinomio Interpolador de Lagrange

n
X
P (x) = yk `k (x)
k=0

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 8 / 15

 Interpolación de Lagrange

Base de Lagrange

Así nos quedan n + 1 polinomios de grado n definidos de esta manera:

n
Q
(x − xk )
k6=j
`j (x) = n
Q
(xj − xk )
k6=j

Estos polinomios forman la Base de Lagrange:

B := {`0 , `1 , ..., `n }

y dependen de las xj de los puntos a interpolar.

Polinomio Interpolador de Lagrange

n
X
P (x) = yk `k (x)
k=0

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 8 / 15

 Interpolación de Lagrange

Base de Lagrange

Así nos quedan n + 1 polinomios de grado n definidos de esta manera:

n
Q
(x − xk )
k6=j
`j (x) = n
Q
(xj − xk )
k6=j

Estos polinomios forman la Base de Lagrange:

B := {`0 , `1 , ..., `n }

y dependen de las xj de los puntos a interpolar.

Polinomio Interpolador de Lagrange

n
X
P (x) = yk `k (x)
k=0

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 8 / 15

 Interpolación de Lagrange

Ejemplo 3

Hallar el polinomio de grado mínimo tal que P (xk ) = yk :

Valores xk 1 -1 3 5
Valores yk 4 20 -4 -52

Calculemos la base de Lagrange:

(x + 1)(x − 3)(x − 5) x3 − 7x2 + 7x + 15

`0 (x) = =
(1 − (−1))(1 − 3)(1 − 5) 16

(x − 1)(x − 3)(x − 5) x3 − 9x2 + 23x − 15

`1 (x) = =
(−1 − 1)(−1 − 3)(−1 − 5) −48

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 9 / 15

 Interpolación de Lagrange

Ejemplo 3

Hallar el polinomio de grado mínimo tal que P (xk ) = yk :

Valores xk 1 -1 3 5
Valores yk 4 20 -4 -52

Calculemos la base de Lagrange:

(x + 1)(x − 3)(x − 5) x3 − 7x2 + 7x + 15

`0 (x) = =
(1 − (−1))(1 − 3)(1 − 5) 16

(x − 1)(x − 3)(x − 5) x3 − 9x2 + 23x − 15

`1 (x) = =
(−1 − 1)(−1 − 3)(−1 − 5) −48

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 9 / 15

 Interpolación de Lagrange

Ejemplo 3

Hallar el polinomio de grado mínimo tal que P (xk ) = yk :

Valores xk 1 -1 3 5
Valores yk 4 20 -4 -52

Calculemos la base de Lagrange:

(x + 1)(x − 3)(x − 5) x3 − 7x2 + 7x + 15

`0 (x) = =
(1 − (−1))(1 − 3)(1 − 5) 16

(x − 1)(x − 3)(x − 5) x3 − 9x2 + 23x − 15

`1 (x) = =
(−1 − 1)(−1 − 3)(−1 − 5) −48

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 9 / 15

 Interpolación de Lagrange

Ejemplo 3

(x − 1)(x + 1)(x − 5) x3 − 5x2 − x + 5

`2 (x) = =
(3 − 1)(3 − (−1))(3 − 5) −16
(x − 1)(x − (−1))(x − 3) x3 − 3x2 − x + 3
`3 (x) = =
(5 − 1)(5 − (−1))(5 − 3) 48
entonces,

P (x) = 4`0 (x) + 20`1 (x) − 4`2 (x) − 52`3 (x) = −x3 + 4x2 − 7x + 8

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 10 / 15

 Interpolación de Lagrange

Ejemplo 3

(x − 1)(x + 1)(x − 5) x3 − 5x2 − x + 5

`2 (x) = =
(3 − 1)(3 − (−1))(3 − 5) −16
(x − 1)(x − (−1))(x − 3) x3 − 3x2 − x + 3
`3 (x) = =
(5 − 1)(5 − (−1))(5 − 3) 48
entonces,

P (x) = 4`0 (x) + 20`1 (x) − 4`2 (x) − 52`3 (x) = −x3 + 4x2 − 7x + 8

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 10 / 15

 Interpolación de Lagrange

Ejemplo 3

(x − 1)(x + 1)(x − 5) x3 − 5x2 − x + 5

`2 (x) = =
(3 − 1)(3 − (−1))(3 − 5) −16
(x − 1)(x − (−1))(x − 3) x3 − 3x2 − x + 3
`3 (x) = =
(5 − 1)(5 − (−1))(5 − 3) 48
entonces,

P (x) = 4`0 (x) + 20`1 (x) − 4`2 (x) − 52`3 (x) = −x3 + 4x2 − 7x + 8

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 10 / 15

 Interpolación de Lagrange

Ejemplo 3

(x − 1)(x + 1)(x − 5) x3 − 5x2 − x + 5

`2 (x) = =
(3 − 1)(3 − (−1))(3 − 5) −16
(x − 1)(x − (−1))(x − 3) x3 − 3x2 − x + 3
`3 (x) = =
(5 − 1)(5 − (−1))(5 − 3) 48
entonces,

P (x) = 4`0 (x) + 20`1 (x) − 4`2 (x) − 52`3 (x) = −x3 + 4x2 − 7x + 8

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 10 / 15

 Interpolación de Lagrange

Interpolación de funciones

Muchas veces se busca interpolar una función f definida en un intervalo

[a, b] ⊂ Dom(f ), es decir, dados x0 , x1 , ..., xn ∈ [a, b] n + 1 valores distintos,
buscamos un polinomio de grado a lo sumo n que interpole a f en esos valores:

P (xk ) = f (xk )

Observación
Si f es un polinomio de grado a lo sumo n, y el polinomio P interpola a f en n + 1
puntos distintos de un intervalo, entonces

P (x) = f (x) ∀ x

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 11 / 15

 Interpolación de Lagrange

Interpolación de funciones

Muchas veces se busca interpolar una función f definida en un intervalo

[a, b] ⊂ Dom(f ), es decir, dados x0 , x1 , ..., xn ∈ [a, b] n + 1 valores distintos,
buscamos un polinomio de grado a lo sumo n que interpole a f en esos valores:

P (xk ) = f (xk )

Observación
Si f es un polinomio de grado a lo sumo n, y el polinomio P interpola a f en n + 1
puntos distintos de un intervalo, entonces

P (x) = f (x) ∀ x

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 11 / 15

 Interpolación de Lagrange

Interpolación de funciones

Muchas veces se busca interpolar una función f definida en un intervalo

[a, b] ⊂ Dom(f ), es decir, dados x0 , x1 , ..., xn ∈ [a, b] n + 1 valores distintos,
buscamos un polinomio de grado a lo sumo n que interpole a f en esos valores:

P (xk ) = f (xk )

Observación
Si f es un polinomio de grado a lo sumo n, y el polinomio P interpola a f en n + 1
puntos distintos de un intervalo, entonces

P (x) = f (x) ∀ x

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 11 / 15

 Interpolación de Lagrange

Interpolación de funciones

Muchas veces se busca interpolar una función f definida en un intervalo

[a, b] ⊂ Dom(f ), es decir, dados x0 , x1 , ..., xn ∈ [a, b] n + 1 valores distintos,
buscamos un polinomio de grado a lo sumo n que interpole a f en esos valores:

P (xk ) = f (xk )

Observación
Si f es un polinomio de grado a lo sumo n, y el polinomio P interpola a f en n + 1
puntos distintos de un intervalo, entonces

P (x) = f (x) ∀ x

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 11 / 15

 Interpolación de Newton

Consideremos x0 , x1 , ..., xn ∈ [a, b] y una función f : [a, b] → R.

Diferencias Divididas
Primera Diferencia Dividida
f (x1 ) − f (x0 )
f [x0 , x1 ] =
x1 − x0

Segunda Diferencia Dividida

f [x1 , x2 ] − f [x0 , x1 ]
f [x0 , x1 , x2 ] =
x2 − x0

...
k−ésima Diferencia Dividida
f [x1 , x2 , ..., xk ] − f [x0 , x1 , ..., xk−1 ]
f [x0 , ..., xk ] =
xk − x0

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 12 / 15

 Interpolación de Newton

Consideremos x0 , x1 , ..., xn ∈ [a, b] y una función f : [a, b] → R.

Diferencias Divididas
Primera Diferencia Dividida
f (x1 ) − f (x0 )
f [x0 , x1 ] =
x1 − x0

Segunda Diferencia Dividida

f [x1 , x2 ] − f [x0 , x1 ]
f [x0 , x1 , x2 ] =
x2 − x0

...
k−ésima Diferencia Dividida
f [x1 , x2 , ..., xk ] − f [x0 , x1 , ..., xk−1 ]
f [x0 , ..., xk ] =
xk − x0

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 12 / 15

 Interpolación de Newton

Interpolación de Newton

Con lo anterior2 podemos definir por recurrencia la siguiente lista de polinomios:

P0 (x) = f (x0 )

P1 (x) = f (x0 ) + f [x0 , x1 ](x − x0 ) = P0 (x) + f [x0 , x1 ](x − x0 )

P2 (x) = P1 (x) + f [x0 , x1 , x2 ](x − x0 )(x − x1 )

Pk (x) = Pk−1 + f [x0 , x1 , ..., xk ](x − x0 )(x − x1 )...(x − xk−1 )

Así siguiendo, se obtiene un polinomio de grado a lo sumo n representado por




Pn (x) = f (x0 ) + (x − x0 ) f [x0 , x1 ] + (x − x1 ) f [x0 , x1 , x2 ] + (x − x2 ) ... + f [x0 , ..., xn ](x − xn−1 )

que se conoce como el Polinomio Interpolador de Newton.

2
Consideremos que f [x0 ] = f (x0 )
Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 13 / 15
 Interpolación de Newton

Ejemplo 4

Hallar el polinomio de grado mínimo tal que P (xk ) = f (yk ):

xk -1 0 1 3
f (xk ) -9 4 13 139

Calculamos las diferencias divididas necesarias.

Primer orden

f (0) − f (−1) 4 − (−9)

f [−1, 0] = = = 13
0 − (−1) 1
f (1) − f (0) 13 − 4
f [0, 1] = = =9
1−0 1
f (3) − f (1) 139 − 13
f [1, 3] = = = 63
3−1 2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 14 / 15

 Interpolación de Newton

Ejemplo 4

Hallar el polinomio de grado mínimo tal que P (xk ) = f (yk ):

xk -1 0 1 3
f (xk ) -9 4 13 139

Calculamos las diferencias divididas necesarias.

Primer orden

f (0) − f (−1) 4 − (−9)

f [−1, 0] = = = 13
0 − (−1) 1
f (1) − f (0) 13 − 4
f [0, 1] = = =9
1−0 1
f (3) − f (1) 139 − 13
f [1, 3] = = = 63
3−1 2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 14 / 15

 Interpolación de Newton

Ejemplo 4

Hallar el polinomio de grado mínimo tal que P (xk ) = f (yk ):

xk -1 0 1 3
f (xk ) -9 4 13 139

Calculamos las diferencias divididas necesarias.

Primer orden

f (0) − f (−1) 4 − (−9)

f [−1, 0] = = = 13
0 − (−1) 1
f (1) − f (0) 13 − 4
f [0, 1] = = =9
1−0 1
f (3) − f (1) 139 − 13
f [1, 3] = = = 63
3−1 2

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 14 / 15

 Interpolación de Newton

Ejemplo 4

Segundo orden

f [0, 1] − f [−1, 0] 9 − 13
f [−1, 0, 1] = = = −2
1 − (−1) 2
f [1, 3] − f [0, 1] 63 − 9
f [0, 1, 3] = = = 18
3−0 3
Tercer orden

f [0, 1, 3] − f [−1, 0, 1] 18 − (−2)

f [−1, 0, 1, 3] = = =5
3 − (−1) 4
Aplicamos la fórmula de Newton:
P (x) = f [−1] + f [−1, 0](x + 1) + f [−1, 0, 1](x + 1)(x − 0) + f [−1, 0, 1, 3](x + 1)(x − 0)(x − 1) =

= −9 + 13(x + 1) − 2(x + 1)x + 5(x + 1)x(x − 1) = 5x3 − 2x2 + 6x + 4

que es el polinomio que cumple lo pedido.

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 15 / 15

 Interpolación de Newton

Ejemplo 4

Segundo orden

f [0, 1] − f [−1, 0] 9 − 13
f [−1, 0, 1] = = = −2
1 − (−1) 2
f [1, 3] − f [0, 1] 63 − 9
f [0, 1, 3] = = = 18
3−0 3
Tercer orden

f [0, 1, 3] − f [−1, 0, 1] 18 − (−2)

f [−1, 0, 1, 3] = = =5
3 − (−1) 4
Aplicamos la fórmula de Newton:
P (x) = f [−1] + f [−1, 0](x + 1) + f [−1, 0, 1](x + 1)(x − 0) + f [−1, 0, 1, 3](x + 1)(x − 0)(x − 1) =

= −9 + 13(x + 1) − 2(x + 1)x + 5(x + 1)x(x − 1) = 5x3 − 2x2 + 6x + 4

que es el polinomio que cumple lo pedido.

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 15 / 15

 Interpolación de Newton

Ejemplo 4

Segundo orden

f [0, 1] − f [−1, 0] 9 − 13
f [−1, 0, 1] = = = −2
1 − (−1) 2
f [1, 3] − f [0, 1] 63 − 9
f [0, 1, 3] = = = 18
3−0 3
Tercer orden

f [0, 1, 3] − f [−1, 0, 1] 18 − (−2)

f [−1, 0, 1, 3] = = =5
3 − (−1) 4
Aplicamos la fórmula de Newton:
P (x) = f [−1] + f [−1, 0](x + 1) + f [−1, 0, 1](x + 1)(x − 0) + f [−1, 0, 1, 3](x + 1)(x − 0)(x − 1) =

= −9 + 13(x + 1) − 2(x + 1)x + 5(x + 1)x(x − 1) = 5x3 − 2x2 + 6x + 4

que es el polinomio que cumple lo pedido.

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 15 / 15

 Interpolación de Newton

Ejemplo 4

Segundo orden

f [0, 1] − f [−1, 0] 9 − 13
f [−1, 0, 1] = = = −2
1 − (−1) 2
f [1, 3] − f [0, 1] 63 − 9
f [0, 1, 3] = = = 18
3−0 3
Tercer orden

f [0, 1, 3] − f [−1, 0, 1] 18 − (−2)

f [−1, 0, 1, 3] = = =5
3 − (−1) 4
Aplicamos la fórmula de Newton:
P (x) = f [−1] + f [−1, 0](x + 1) + f [−1, 0, 1](x + 1)(x − 0) + f [−1, 0, 1, 3](x + 1)(x − 0)(x − 1) =

= −9 + 13(x + 1) − 2(x + 1)x + 5(x + 1)x(x − 1) = 5x3 − 2x2 + 6x + 4

que es el polinomio que cumple lo pedido.

Carlos Alliera (calliera@[Link]) 17 de mayo de 2021 15 / 15