DEPARTAMENTO DE CIENCIAS EXACTAS

ANÁLISIS NUMÉRICO

GRUPO N°03

Interpolación por Newton

Introducción

El método de interpolación de Newton es un poco más complicado que el de LaGrange, pero es más
preciso. Hay ocasiones en las que resulta útil construir varios polinomios aproximados, P1(x), P2(x),…,
PN(x) y, después, elegir el más adecuado a las necesidades de interpolación. Si se utiliza los polinomios
de interpolación de LaGrange, uno de los inconvenientes es que no se pueden utilizar los cálculos
realizados en la construcción de PN-1(x) para la de PN(x); cada polinomio debe construirse
individualmente y para calcular polinomios degrado elevado es necesario hacer muchas operaciones.
Por lo que el método de Newton es una mejor alternativa a utilizar en el proceso de interpolación.

Fundamento teórico

CARACTERISTICAS DE METODO

 Menor cantidad de cálculos necesaria para la interpolación

 Evaluación más sencilla del error
 Pueden utilizarse los resultados de los cálculos previos cuando el número de datos tiene que
incrementarse o decrementarse
 La interpolación de Newton es más adecuada que la interpolación de Lagrange para obtener otros
modelos numéricos (por ejemplo las aproximaciones de derivada por diferencias) o para
desarrollar una interpolación por medio de series de potencias
 La interpolación de Newton puede realizarse hacia adelante y hacia atrás, ambas son
matemáticamente equivalentes pero representan expresiones distintas.
 Dependiendo de cómo se aplique la formula una de las interpolaciones puede ser más conveniente
que la otra. Para ajustar los datos en una tabla, es preferible la interpolación hacia adelante
DEDUCCION DE LA FORMULA
Los polinomios de interpolación se calculan mediante un esquema recursivo:

Afirmación Justificación
𝑃1 (𝑥 ) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 )
𝑃2 (𝑥 ) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 ) Polinomio de
𝑃3 (𝑥 ) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) + 𝑎3 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 ) Newton, añadido
… termino a termino
𝑃𝑁 (𝑥 ) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 ) + ⋯ + 𝑎𝑁 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) … (𝑥 − 𝑥𝑁 )

:
𝑃𝑁 = 𝑃𝑁−1 (𝑥 ) + 𝑎𝑁 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 ) … (𝑥 − 𝑥𝑁−1 ) El polinomio PN(x) se
obtiene a partir de PN-1(x)
usando la recurrencia

El polinomio PN(x) calculado así es el polinomio de interpolación de Newton.

A. Multiplicación encajada
Para evaluar el polinomio PN(x), lo más eficiente (menos operaciones) es usar el esquema de
multiplicaciones encajadas. Para P3(x), por ejemplo:
𝑆3 = 𝑎3
𝑆2 = 𝑆3 (𝑥 − 𝑥2 ) + 𝑎2
𝑆1 = 𝑆2 (𝑥 − 𝑥1 ) + 𝑎1
𝑆0 = 𝑆1 (𝑥 − 𝑥0 ) + 𝑎0

Esta última cantidad S0 es P3(x).

B. Cálculo del polinomio de interpolación de Newton

Para interpolar una función dada f(x) se debe encontrar los coeficientes ak de todos los polinomios
P1(x), P2(x),…, PN(x). Entonces cada Pk(x) es el polinomio de Newton que tiene como nodos x 0, x1,…, xk.
Para el polinomio P1(x), los coeficientes a0 y a1 tienen un significado familiar:

Afirmación Justificación

𝑓 (𝑥0 ) = 𝑃1 (𝑥0 ) = 𝑎0 + 𝑎1 (𝑥0 − 𝑥0 ) = 𝑎0 Primer polinomio de Newton

evaluado en el punto 𝑥0
𝑎0 = 𝑓(𝑥0 )
𝑓 (𝑥1 ) = 𝑃1 (𝑥1 ) = 𝑎0 + 𝑎1 (𝑥1 − 𝑥0 ) = 𝑓 (𝑥0 ) + 𝑎1 (𝑥1 − 𝑥0 ) Primer polinomio de Newton
evaluado en el punto 𝑥1

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

Es decir, a1 es la pendiente de la recta que pasa por los puntos (𝑥0 , 𝑓 (𝑥0 )) y (𝑥1 , 𝑓(𝑥1 )). Los coeficientes
a0 y a1 son los mismos para P1(x) y P2(x).
Afirmación Justificación

𝑓 (𝑥2 ) = 𝑃2 (𝑥2 ) = 𝑎0 + 𝑎1 (𝑥2 − 𝑥0 ) + 𝑎2 (𝑥2 − 𝑥0 )(𝑥2 − 𝑥1 ) Ahora evaluamos la expresión en el

nodo x2

𝑓(𝑥2 )−𝑎0−𝑎1(𝑥2 −𝑥0)

𝑎2 = Coeficiente 𝑎2 despejado
(𝑥2 −𝑥0 )(𝑥2 −𝑥0)

𝑓(𝑥2 )−𝑓(𝑥1 ) 𝑓(𝑥1 )−𝑓(𝑥0 )

( − )
𝑥2 −𝑥1 𝑥1 −𝑥0
𝑎2 = Otra manera de visualizar
(𝑥2 −𝑥0 )
𝑎2

El cálculo de los coeficientes se puede realizar de forma más rápida y sencilla utilizando la notación de
las diferencias divididas, estas se definen como:

Afirmación Justificación

𝑓 [𝑥𝑘 ] = 𝑓(𝑥𝑘 ) La diferencia dividida de orden cero

𝑓[𝑥𝑘 ]−𝑓[𝑥𝑘−1 ]
𝑓 [𝑥𝑘−1 , 𝑥𝑘 ] = La diferencia dividida de primer orden
𝑥𝑘 −𝑥𝑘−1
𝑓[𝑥𝑘−1 ,𝑥𝑘 ]−𝑓[𝑥𝑘−2 ,𝑥𝑘−1 ]
𝑓 [𝑥𝑘−2 , 𝑥𝑘−1 , 𝑥𝑘 ] = 𝑥𝑘 −𝑥𝑘−2
La diferencia dividida de segundo orden

𝑓[𝑥𝑘−2 ,𝑥𝑘−1 ,𝑥𝑘 ]−𝑓[𝑥𝑘−3 ,𝑥𝑘−2 ,𝑥𝑘−1 ]

𝑓 [𝑥𝑘−3 , 𝑥𝑘−2 , 𝑥𝑘−1 , 𝑥𝑘 ] = La diferencia dividida de tercer orden
𝑥𝑘 −𝑥𝑘−3
𝑓[𝑥𝑘−𝑗+1 ,…,𝑥𝑘 ]−𝑓[𝑥𝑘−𝑗 ,…,𝑥𝑘−1 ]
𝑓[𝑥𝑘−𝑗 , 𝑥𝑘−𝑗+1 , … , 𝑥𝑘 ] = Aplicando el método de recursión
𝑥𝑘 −𝑥𝑘−𝑗

Fig 1. Polinomio interpolador de Newton

Esta fórmula se utiliza para calcular la Tabla de Diferencias Divididas:

𝑥 𝑓(𝑥) Primera diferencia Segunda diferencia Tercera diferencia

𝑓 [𝑥0 , 𝑥1 ] = 𝑎1 𝑓 [𝑥0 , 𝑥1 , 𝑥2 ] = 𝑎2 𝑓[𝑥3 , 𝑥2 , 𝑥1 , 𝑥0 ] = 𝑎3
𝑥0 𝑓 (𝑥0 ) 𝑓 (𝑥1 ) − 𝑓(𝑥0 ) 𝑓[𝑥1 𝑥2 ] − 𝑓[𝑥0 𝑥1 ] 𝑓 [𝑥3 𝑥2 𝑥1 ] − 𝑓[𝑥2 𝑥1 𝑥0 ]
= 𝑎0 = = =
𝑥1 − 𝑥0 𝑥2 − 𝑥0 𝑥3 − 𝑥0
𝑓 [𝑥1 , 𝑥2 ] 𝑓[𝑥1 , 𝑥2 , 𝑥3 ]
𝑥1 𝑓(𝑥1 ) 𝑓 (𝑥2 ) − 𝑓(𝑥1 ) 𝑓 [𝑥2 𝑥3 ] − 𝑓[𝑥1 𝑥2 ]
= =
𝑥2 − 𝑥1 𝑥3 − 𝑥1
[
𝑓 𝑥2 , 𝑥3 ]
𝑥2 𝑓(𝑥2 ) 𝑓 (𝑥3 ) − 𝑓(𝑥2 )
=
𝑥3 − 𝑥2
𝑥3 𝑓(𝑥3 )

C. Propiedades de las diferencias divididas

 La diferencia dividida de orden K es el coeficiente de 𝑥 𝑘 en 𝑃𝑘 (𝑥 )

 La diferencia dividida de cualquier orden es independiente del orden en el que se tomen los
nodos
 La diferencia dividida de orden K se calcula recursivamente a partir de dos diferencias divididas
de orden K-1

D. Elección de los puntos de interpolación

El polinomio interpolador puede presentar comportamientos patológicos en el sentido que cuando se

aumenta el número de puntos de interpolación en un intervalo, la desviación máxima entre el polinomio
y la función en dicho intervalo aumente con el grado del polinomio
1
Esto sucede frecuentemente con las funciones 𝐵(𝑥 ) = |𝑥 | y 𝑅(𝑥 ) = 1+25𝑥2 en el intervalo [-1; 1] que
corresponden a las funciones de Bernstein y Runge. Cuando se aumenta el grado del polinomio
interpolados de estas funciones con puntos igualmente espaciados, sucede que:
lim ( max |𝑓 (𝑥 ) − 𝑃𝑛 (𝑥)|) = ∞
𝑛→∞ −1≤𝑥≤1
Con un polinomio 𝑃𝑛 (𝑥) de grado ≤ 𝑛 , la interpolación resultante oscila hacia los extremos del
intervalo, es decir, cerca de -1 y 1.
Por lo tanto se deben escoger los mejores puntos para la interpolación en caso de que se tenga libertad
para elegirlos. Para que el error se minimice se eligen los puntos de interpolación como ceros de los
polinomios de Chebychev, definidos por:
𝑇𝑛 (𝑥 ) = cos(𝑛arccos(𝑥))
De su definición los polinomios 𝑇𝑛 (𝑥 ) solo toman valores en el intervalo [-1; 1], pues fuera de este
intervalo la función arccos(𝑥) no está definida. Como los ceros sus están en el intervalo [-1; 1], para una
función f(x) definida en el intervalo |a, b| se tiene que realizar un cambio de variables a la variable t
definida por:
𝑏−𝑎 𝑏+𝑎
𝑥= 𝑡+
2 2
De forma que f(x) tome valores en [-1; 1]. Los valores de los ceros de los polinomios de Chebychev
vienen dados por los valores de x que satisfacen cos(𝑛𝜃 ) = 0, donde 𝜃 = arccos 𝑥. Por lo tanto se tiene

(2𝑘 + 1)𝜋
𝑛𝜃 =
2
(2𝑘 + 1)𝜋
𝜃=
2𝑛

De esta manera los ceros de 𝑇𝑛 (𝑥 ) vienen dados por la expresión

(2𝑘 + 1)𝜋
𝑥𝑛 = cos ( )
2𝑛

Fig 2. Elección de los puntos de interpolación

ERROR DEL METODO

Sean 𝑓 ∈ 𝐶 𝑛 [𝑎, 𝑏] y 𝑥0 , 𝑥1 , … , 𝑥𝑛 n+1 puntos distintos en [a,b]. Entonces ∃𝜉 ∈ (𝑎, 𝑏) tal que:

Afirmación Justificación
𝑛) (
𝑓 𝜉)
𝑓 [𝑥0 , 𝑥1 , … , 𝑥𝑛 ] =
𝑛!

𝑔(𝑥 ) = 𝑓 (𝑥 ) − 𝑃𝑛 (𝑥 ) = 𝑓 (𝑥 ) − (𝑓 (𝑥0 ) + 𝑓[𝑥0 𝑥1 ](𝑥 − 𝑥0 ) + ⋯ Sea 𝑔(𝑥 ) = 𝑓 (𝑥 ) − 𝑃𝑛 (𝑥), siendo

𝑃𝑛 (𝑥) el polinomio de
interpolación de Newton de f(x),
en los n+1 puntos
+𝑓 [𝑥0 𝑥1 … 𝑥𝑛 ](𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) … (𝑥 − 𝑥𝑛−1 )

𝑛) 𝑓𝑛) (𝜉)
𝑓 𝑛 (𝜉 ) − 𝑃𝑁 (𝜉 ) = 0 → 𝑓 𝑛) (𝜉 ) = 𝑛! 𝑓 [𝑥0 𝑥1 … 𝑥𝑛 ] → 𝑓 [𝑥0 𝑥1 … 𝑥𝑛 ] = 𝑔 (𝑥 ) Se anula en
𝑛!
𝑥0 , 𝑥1 , … , 𝑥𝑛 , y aplicando el teorema de Rolle generalizado existirá un 𝜉 ∈ (𝑎, 𝑏) tal que 𝑔𝑛) (𝜉 ) = 0
 𝑓𝑛+1) (𝜉) 𝑛
𝐸𝑁 (𝑥 ) = (𝑛+1)!
∏𝑖=0(𝑥 − 𝑥𝑖 ) == 𝑓[𝑥0 , 𝑥1 , … , 𝑥𝑛 , 𝑧] ∏𝑛𝑖=0(𝑥 − 𝑥𝑖 ) Utilizando esta relación,
obtenemos el término de error

Siendo 𝑧 ∈ (𝑎, 𝑏).

Esta expresión permite estimar el error del polinomio de interpolación de grado N cuando no se conoce
la función, siempre que sea posible añadir un nodo más y obtener la diferencia dividida de orden N+1.
Aplicación en Electrónica
Calcule el polinomio de interpolación para la función f(x)=|x| en los nodos xi= -4. -1, 2, 5 y 7. Grafique la
función y la aproximación de la interpolación, además estime el valor de y para x=4.5.

Grafica

Polinomio de interpolación:
𝑥4 43𝑥 3 149𝑥 2 59𝑥 875
− − + + +
7128 3564 792 324 891

Valor interpolado = 4.45

Valor evaluado en la función = 4.50
Error generado al evaluar el valor ingresado comparado entre la función y el polinomio
generado es: 0.04581404
  Opción 1: Ingresar puntos de manera directa con el ratón.
Interpolación de
 Opción 2: Ingresar la función.
Newton
 Opción 3: Ingresar (n+1) puntos y el grado que se desea el
INICIO
Polinomio de Newton

Seleccione una opción:

Desea evaluar
Se despliega el plano sí 1 el polinomio
interpolado ?

no
sí
Seleccionar los un Ingresar
mínimo de 2 puntos con 2 sí la Ingresar el sí
el ratón función valor de X a
evaluar
no

3
x=Téminos a
interpolar?

Ingresar el grado de
Polinomio de no
no
Newton deseado (n)
Se calcula el Error sí

Ingresar (n+1)
Términos (x.y) Se imprime el Error

Se imprime el polinomio
Pn(x) evaluado
Desea ver el proceso
para hallar la matriz de
diferencias divididas
Se imprime el tiempo de
convergencia del método
implementado
sí
Empezar a tomar el tiempo del programa
implementado

Desea repetir
Se imprime el proceso para hallar la matriz de el proceso?
diferencias divididas

no
Se imprime la matriz de diferencias divididas
FIN

Se imprime el polinomio interpolador de

Newton
 · Opción 1: Ingresar puntos de manera directa con el ratón.
Interpolación de
· Opción 2: Ingresar la función.
Newton
· Opción 3: Ingresar (n+1) puntos y el grado que se desea el
INICIO
Polinomio de Newton

Seleccione una opción:

Desea evaluar
Se despliega el plano sí 1 el polinomio
interpolado ?

no
sí
Seleccionar los un Ingresar
mínimo de 2 puntos con 2 sí la Ingresar el sí
el ratón función valor de X a
evaluar
no

3
x=Téminos a
interpolar?

Ingresar el grado de
Polinomio de no no
Newton deseado (n)
Se calcula el Error sí
Ingresar (n+1)
Términos (x,y) Se imprime el Error

Se imprime el polinomio
Pn(x) evaluado
Desea ver el proceso
para hallar la matriz de
diferencias divididas
Se imprime el tiempo de
convergencia del método
implementado
sí
Empezar a tomar el tiempo del programa
implementado

Desea repetir
Se imprime el proceso para hallar la matriz de el proceso?
diferencias divididas

no
Se imprime la matriz de diferencias divididas

FIN

Se imprime el polinomio interpolador de

Newton

MATLAB 2015
02-01-17 Cerda, Landázuri
02-01-17 Vélez
03-01-17 Bernis DIAGRAMA DE FLUJO NA
INTERPOLACIÓN DE NEWTON

Norma:"ANSI"
DS1-001

Core i7 x 64