INTERPOLACIÓN POLINÓMICA

DE NEWTON EN DIFERENCIAS DIVIDIDAS

 Interpolación Polinómica de Newton
en Diferencias Divididas
Con frecuencia se encontrará que se tiene que estimar
valores intermedios entre datos definidos por puntos.
El método más común es la Interpolación Polinomial
de Newton en Diferencia Divididas. Para un
polinomio de n-enésimo grado.

(ec. 1)
 Interpolación Polinómica de Newton
en Diferencias Divididas
Dados n+1 puntos , existe uno y solo un polinomio
de grado n que pasa por todos los puntos.

Ejemplos de Interpolación Polinomial:

(a) De primer grado (lineal) que une dos puntos.
(b) de segundo grado (cuadrática o parabólica) que une tres puntos.
(c) De tercer grado (cúbica) que une cuatro puntos
 Interpolación Polinómica de Newton
en Diferencias Divididas
 Interpolación Lineal: forma más simple de interpolación que consiste
en unir dos puntos es mediante una línea recta.

𝑓 1 ( 𝑥 ) − 𝑓 ( 𝑥 0 ) 𝑓 ( 𝑥 1) − 𝑓 ( 𝑥 0 )
=
𝑥 − 𝑥0 𝑥1 − 𝑥0

𝑓 ( 𝑥 1) − 𝑓 ( 𝑥 0 )
𝑓 1 ( 𝑥 )= 𝑓 ( 𝑥 0 ) + ( )
𝑥 − 𝑥 0ec.2
𝑥1− 𝑥0

La notación designa que este es un polinomio de interpolación de

primer grado.
Representa la pendiente de la línea que une los puntos. Es una
aproximación en diferencia dividida finita a la primera derivada
 Interpolación Polinómica de Newton
en Diferencias Divididas
 Interpolación Cuadrática: Si se tienen tres puntos como datos, estos pueden
ajustarse en un polinomio de segundo grado. La forma general queda
expresada:
 ec.3
 Esta ecuación es similar a la primera

 Agrupando términos:

 Donde
 Interpolación Polinómica de Newton
en Diferencias Divididas
 Interpolación Cuadrática:
Para encontrar los valores de los coeficientes se evalúa en la
ecuación (3) en su orden , , para obtener:

ec.4

ec.5

ec.6
 Interpolación Polinómica de Newton
en Diferencias Divididas
 Forma general de los polinomios de interpolación de Newton:
Se puede generalizar para ajustar a un polinomio de n-enésimo grado a n+1
datos:
ec.7
Para un polinomio de n-enésimo grado se requieren n + 1 puntos.
Los coeficientes resultantes serán:

ec.8
 .
ec.9
ec.10
 .
ec.11
 Interpolación Polinómica de Newton
en Diferencias Divididas
 Forma general de los polinomios de interpolación de Newton:
Donde la evaluaciones de la función colocadas entre paréntesis son
diferencias divididas finitas. La primera diferencia dividida finita se
representa por:

ec.12
La segunda diferencia dividida finita , que representa la diferencia de las dos
primeras diferencias divididas, se expresa

ec.13
 Interpolación Polinómica de Newton
en Diferencias Divididas
 Forma general de los polinomios de interpolación de Newton:
En forma similar, la enésima diferencia dividida es
ec.14
Estas diferencias sirven para evaluar los coeficientes en las ecuaciones (ec.8) a
(ec.11), los cuales se sustituirán en la ecuación (ec.7) para obtener el polinomio
de interpolación de Newton en diferencias divididas

ec.15

Representación gráfica de la naturaleza recursiva de las diferencias divididas

finitas.

Problema
1_ Dado los datos:
x 1 2 3 3.5 5 6
f(x) 1 5 7.5 8 5 1 𝑏 4= 𝑓 ( 𝑥 5 , 𝑥 4 , 𝑥 3 , 𝑥 2 , 𝑥1 , 𝑥 0 )
Calcule f(3.4) mediante polinomio de interpolación de Newton de orden 2 a 4. Escoja la
secuencia de puntos para su estimación con el fin de obtener la mejor exactitud posible.
Construya en papel milimetrado la gráfica correspondiente a los datos y la curva de cada
polinomio obtenido. Emplee colores diferentes para cada una.
i x f(x) Primero Segundo tercero cuarto quinto
0 1 1 𝑓 ( 𝑥1 ) − 𝑓 ( 𝑥 0 ) 𝑓 ( 𝑥 2 ) − 𝑓 ( 𝑥1 ) − 𝑓 ( 𝑥 1 ) − 𝑓 ( 𝑥 0 )𝑏 = 𝑓 𝑥 , 𝑥 , 𝑥 ,𝑏𝑥5 = , 𝑥0 )( 5 4 3 2 1
𝑓 𝑥 , 𝑥 ,𝑥 ,𝑥 ,𝑥 , 𝑥
1 2 5 𝟒 −𝟎. 𝟕𝟓 −𝟎 . 𝟏 −𝟎.𝟎𝟏𝟔𝟔𝟔𝟔𝟔
𝑏1 = 𝑓 ( 𝑥 1 , 𝑥 0 )=
𝑏 2= 𝑓 ( 𝑥 𝑥 2 , 1𝑥 −1 ,𝑥𝑥00 ) = 𝑓 ( 𝑥 ) − 𝑓 1( 𝑥3 ) (𝑓 (3𝑥 1) −
𝑥 2 − 𝑥𝑏 = 𝑓 𝑥 ,
𝑥 𝑥 −
𝑥2 2 − 𝑥0 2 𝑓 ( 𝑥 1 )
,
2 01𝑥
𝑥 , 𝑥 0)
4 ( 4 3 2 1

−𝟏 −𝟎.𝟏𝟔𝟔𝟔𝟔𝟕𝟎 . 𝟏 𝑏 4= 𝑓 ( 𝑥 5 , 𝑥 4 , 𝑥 3 , 𝑥 2 , 𝑥1 )
3
𝑓 ( 𝑥 2 ) − 𝑓 ( 𝑥1 ) −
𝑏1 = 𝑓 ( 𝑥 2 , 2
𝑥𝑏1 ) =
= 𝑓 ( 𝑥 , 𝑥 , 𝑥 )= 𝑥 3 − 𝑥𝑏 2 =𝑓 𝑥 𝑥
( , 2 −,𝑥
𝑥 𝑥 1 ,𝑥
1)
2 3 7.5 2 𝑥3 2 −2𝑥1 1 3
𝑥 − 𝑥
4 3 2
𝑓 ( 𝑥3 ) − 𝑓 ( 𝑥 2 ) 𝑓 ( 𝑥 ) − 𝑓 ( 𝑥 ) 𝑓 ( 𝑥 ) − 𝑓 ( 𝑥 2)

𝟏 . 𝟎 −𝟏 . 𝟓
4 33 3 3
𝑏1 = 𝑓 ( 𝑥 3 , 𝑥 2 )= −
𝑥 − 𝑥𝑏 = 𝑓 ( 𝑥 5 ,𝑥𝑥3 4−, 𝑥 𝑥23 , 𝑥 2 )
𝑏 2= 𝑓 ( 𝑥𝑥43,− 𝑥 3𝑥,2𝑥 2 ) = 𝑓 𝑥 4 − 𝑓 3 3𝑥
3 3.5 8 ( 5) ( 𝑥4 4) − 𝑥𝑓 2( 𝑥 4 ) − 𝑓 ( 𝑥 3 )
−𝟐 −𝟎 . 𝟖
𝑏1 = 𝑓 ( 𝑥 4 , 𝑥 𝑏
3 )= = 𝑓
𝑓 ( 𝑥4) − 𝑓 ( 𝑥3)
( 𝑥𝑥 , 𝑥
5 −4𝑥 3 , 𝑥 ) =
𝑥5− 𝑥4
−
𝑥4 − 𝑥3
4 5 5 2
4 3
𝑓 ( 𝑥5 ) − 𝑓 ( 𝑥 4 ) 𝑥5− 𝑥3
−𝟒
𝑏1 = 𝑓 ( 𝑥 5 , 𝑥 4 )=
𝑥5− 𝑥4
5 6 1
 Continuación del problema

Polinomios generados de :

Segundo orden:

f2(x) = f(x0) +f[x1x0](x-x0) + f[x2x1x0](x-x0)(x-x1)

f2(x) = 1 + 4(x - 1) - 0.75(x- 1)(x - 2)

f2(3.4) = 1 + 4(3.4 - 1) - 0.75(3.4 - 1)(3.4 - 2) = 8.08

Tercer orden:

f3(x) = f(x0) +f[x1x0](x-x0) + f[x2x1x0](x-x0)(x-x1) + f[x3x2x1x0](x-x0)(x-x1)(x-x2)

f3(x) = 1 + 4(x - 1) - 0.75(x - 1)(x - 2) - 0.1(x - 1)(x – 2)(x - 3)

f3(3.4) = 1 + 4(3.4 - 1) - 0.75(3.4 - 1)(3.4 - 2) - 0.1(3.4 - 1)(3.4 – 2)(3.4 - 3) = 7.9456

Cuarto orden:

f4(x) = f(x0) +f[x1x0](x-x0) + f[x2x1x0](x-x0)(x-x1) + f[x3x2x1x0](x-x0)(x-x1)(x-x2) + f[x4x3x2x1x0](x-x0)(x-x1)(x-x2)(x-x3)

f4(x) = 1 + 4(x - 1) - 0.75(x- 1)(x - 2) - 0.1(x - 1)(x – 2)(x - 3) - 0.01666666(x - 1)(x – 2)(x - 3)(x - 3.5)

f4(3.4) = 1 + 4(3.4 - 1) - 0.75(3.4- 1)(3.4 - 2) - 0.1(3.4 - 1)(3.4 – 2)(3.4 - 3)

- 0.01666666(3.4 - 1)(3.4 – 2)(3.4 - 3)(3.4 - 3.5) = 7.94784

Quinto orden:
f5(x) = f(x0) +f[x1x0](x-x0) + f[x2x1x0](x-x0)(x-x1) + f[x3x2x1x0](x-x0)(x-x1)(x-x2) + f[x4x3x2x1x0](x-x0)(x-x1)(x-x2)(x-x3) +
f[x5x4x3x2x1x0](x-x0)(x-x1)(x-x2)(x-x3)(x-x4)

f5(x) = 1 + 4(x- 1) - 0.75(x - 1)(x - 2) - 0.1(x - 1)(x – 2)(x - 3) - 0.0166666(3.4 - 1)(3.4 – 2)(3.4 - 3)(3.4 - 3.5) +
0.02333333(3.4 - 1)(3.4 – 2)(3.4 - 3)(3.4 - 3.5)(3.4 - 5)

f5(3.4) = 1 + 4(3.4 - 1) - 0.75(3.4 - 1)(3.4 - 2) - 0.1(3.4 - 1)(3.4 – 2)(3.4 - 3)- 0.0166666(3.4 - 1)(3.4 – 2)(3.4 - 3)(3.4 - 3.5) +
0.02333333(3.4 - 1)(3.4 – 2)(3.4 - 3)(3.4 - 3.5)(3.4 - 5)
= 7.9528576
 Interpolación Polinómica de Newton
en Diferencias Divididas
Datos ajustados a un polinomio de Segundo orden ( 3 primeros puntos)
f2(x) = 1 +4(x-1) -0.75(x-1)(x-2)
10
Datos
9 Polinomio segundo orden
Estimación 2do orden f(3.4) = 8.08
8

3
f(x)

-1

-2

-3

-4

-5
0 0.5 1 1.5 2 2.5 3 3.4
3.5 4 4.5 5 5.5 6 6.5 7
x
 Interpolación Polinómica de Newton
en Diferencias Divididas
Datos ajustados a un polinomio de Tercer orden ( 4 primeros puntos)
f3(x) = 1 + 4(x-1) - 0.75(x-1)(x-2) - 0.1(x-1)(x-2)(x-3)
10

9
Datos
8 Polinomio tercer orden
Estimación 3er orden f(3.4) = 7.9456
7

3
f(x)

-1

-2

-3

-4

-5 3.4
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7
x
 Interpolación Polinómica de Newton
en Diferencias Divididas
Datos ajustados a un polinomio de cuarto orden ( 5 primeros puntos)
f4(x) = 1 + 4(x-1) - 0.75(x-1)(x-2) - 0.1(x-1)(x-2)(x-3) - 0.01666666(x-1)(x-2)(x-3)(x-3.5)
10
9
Datos
8 Polinomio cuarto orden
7 Estimación 4to orden f(3.4) = 7.94784
6
5
4
3
2
1
f(x)

0
-1
-2
-3
-4
-5
-6
-7
-8
-9
-10 3.4
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7
x
 Interpolación Polinómica de Newton
en Diferencias Divididas
Datos ajustados con polinomios de segundo, tercero, cuarto y quinto
orden
10

2
f(x)

-2 Datos
Polinomio segundo orden
polinomio tercer orden
-4 polinomio cuarto orden
polinomio quinto orden
Estimación 2do f(3.4) = 8.08
-6 Estimacion 3er f(3.4) = 7.9456
Estimación 4to f(3.4) = 7.94784
Estimación 5to f(3.4) = 7.9528576
-8

-10
3.4
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7
x