Tema III Interpolación, Derivación e Integración Numéricas

4. Interpolación, Derivación e Integración Numéricas

Objetivo: El alumno analizará y comparará algunos de los métodos numéricos para interpolar,
derivar e integrar funciones.

INTRODUCCIÓN

El análisis numérico proporciona diversos métodos para interpolar, derivar e integrar

funciones definidas en forma tabular, las cuales pueden ser resultado de experimentos o la
simple tabulación de alguna función que se desea integrar o derivar.

DIFERENCIAS FINITAS HACIA ADELANTE O DE AVANCE

Sea una función y = f ( x) , continua y diferenciable en el intervalo cerrado x0,xn, la cual es

expresada en forma tabular como sigue:

x y = f ( x)
x0 y0
x1 y1
. .
. .
xn yn

Se definen como primeras diferencias hacia adelante o de avance de y = f ( x) a:

y 0 = y1 − y 0
y1 = y 2 − y1

yi = yi +1 − yi i = 0,1,2,3, , n − 1
donde  es el operador diferencia.

Las segundas diferencias hacia adelante son las diferencias de las primeras diferencias, las
terceras diferencias hacia adelante son las diferencias de las segundas diferencias; de manera
general las k-ésimas diferencias hacia adelante de una función se pueden obtener como:

i = 0,1,2, , n − k
k yi = k −1 yi +1 − k −1 yi
k = 1,2,3, , n

-1-
 Tema III Interpolación, Derivación e Integración Numéricas
3.2 TABLA DE DIFERENCIAS

La función expresada en forma tabular y sus diferencias se pueden agrupar en una tabla de
diferencias como sigue:

TABLA DE DIFERENCIAS HACIA ADELANTE DE y=f(x)

x y=f(x) yi   yi   yi   yi yi
x0 y0
y0
x1 y1  y0
y1  y0
x2 y2  y1
  y0
y2  y1

x3 y3  y2  y1  y0
y3  y2

x4 y4  y3  y2  y1
y4  y3

. .  y4  y3  y2
.  y4
  .
. . .  yn-4
.  yn-5
.  yn-3
. .  yn-2


yn-1
xn yn

El proceso de obtención de diferencias es finito y deja de aplicarse cuando una de estas

diferencias se vuelve constante o aproximadamente constante, sin anularse; por ejemplo en un
polinomio de grado enésimo, puede comprobarse que es posible obtener hasta las n-ésimas
diferencias, si el incremento en la variable independiente es constante; por ejemplo: Sean la s
funciones, y = x , y = x 2 , y = x 3 ; tabulando algunos pares de puntos y obteniendo sus
diferencias hacia adelante:

x y=x y x y = x2 y 2y x y = x3 y 2y y

0 0 0 0 0 0
1 1 1
1 1 1 1 2 1 1 6
1 3 7 6
2 2 2 4 2 2 8 12
1 5 19 6
3 3 3 9 2 3 27 18
1 7 37
4 4 4 16 4 64

-2-
 Tema III Interpolación, Derivación e Integración Numéricas
INTERPOLACIÓN CON INCREMENTO CONSTANTE EN LA VARIABLE INDEPENDIENTE

El problema de interpolación consiste en determinar el valor de una función f(x) para un valor de
la variable independiente x que se encuentren entre dos valores consecutivos conocidos es
decir x0: < xk <x1.

Se considera que la función se aproxima a un polinomio de grado n que coincide con todos los
puntos de la función ( expresada en forma tabular ) y que el incremento de la variable
independiente es constante e igual a h:

x y = f ( x)
x0 y0
x1=x0+h y1
x2= x0+2h y2
x3= x0+3h y3
. .
. .
xn=x0+nh yn

h = x1 − x 0 = x 2 − x1 = = x n − x n −1

3.3.1 POLINOMIO DE INTERPOLACIÓN DE AVANCE DE NEWTON

Si se construye una tabla de diferencias para la función antes presentada, será similar a la
mostrada en el apartado 3.2 y de acuerdo con dicha tabla de diferencias y la definición de
diferencias se puede expresar cualquier valor yi de la función en términos de y0 y de las
diferencias hacia adelante de y0 de la siguiente manera:
y1 = y 0 + y0 (1)
y 2 = y1 + y1 = y 0 + y 0 + y1
y si :

y1 = y 0 + 2 y 0
2 y 0 = y1 − y 0
sustituyendo en y2 :

y 2 = y 0 + 2 y 0 + 2 y 0 (2)

para y3 se sigue un procedimiento similar obteniéndose:

y 3 = y 0 + 3y 0 + 32 y 0 + 3 y 0 (3)

Las ecuaciones (1), (2) y (3) permiten ver que los coeficientes que multiplican a las diferencias
hacia adelante de y0 corresponden a los del binomio de Newton, por lo que se puede escribir
para un valor cualquiera yk de la función:

-3-
 Tema III Interpolación, Derivación e Integración Numéricas

 k  k  k  k
y k = y0 +   y0 +   2 y0 +   3 y0 ++  n y0 (4)
 1  2  3  n
desarrollando los números combinatorios:

k ( k − 1) 2 k ( k − 1)( k − 2) 3 k ( k − 1)( k − 2)( k − 3) 4

y k = y 0 + ky 0 +  y0 +  y0 +  y 0 + (5)
2! 3! 4!

La ecuación (5) es el Polinomio de Interpolación de avance de Newton (algunos autores lo

llaman de Newton-Gregory); donde:

yk valor aproximado de la función que desea obtenerse para un valor x= xk

y0 valor de la función, incluido en la tabulación, inmediato anterior al valor que se va a calcular,

correspondiente a x0

y0, y0, y0,..., diferencias hacia adelante de y0

k fracción del incremento constante h que debe agregarse a y0 para obtener yk puede calcularse
de considerar que x k = x 0 + kh ; por lo tanto:
x − x0
k= k (6)
h

El polinomio de avance de Newton además de aplicarse a interpolación, permite también

realizar extrapolaciones, es decir, determinar valores de la función que no se encuentre dentro
del intervalo de la tabulación.

Cuando las diferencias hacia adelante, de orden n, de la función se hacen constantes e iguales
la función se puede aproximar a un polinomio de avance de Newton de grado n y obtener
prácticamente sin error el resultado de la interpolación; si las n-ésimas diferencias se hacen
aproximadamente constantes o no se hacen constantes, la función no se comportará
exactamente como un polinomio de grado n y entonces se presentan errores de discretización,
redondeo y truncamiento.

POLINOMIOS INTERPOLANTES Y DIAGRAMA DE ROMBOS PARA LA INTERPOLACIÓN

Si se toma en cuenta que los valores de x0 y y0 (inmediato anterior al que se desea interpolar)
no necesariamente se encuentran al inicio del intervalo de valores de la función expresada en
forma tabular, sino que pueden ubicarse en cualquier posición; es posible construir una tabla de
diferencias de manera que los valores de yo y sus diferencias hacia adelante pueden ubicarse
en una posición intermedia; además de estas diferencias es factible incluir números
combinatorios por arriba y por debajo de las diferencias y unir las columnas de la función y sus
diferencias por medio de diagonales; a dicho arreglo se le conoce como diagrama de rombos
(por la forma tan particular que toma) y dicho diagrama permite, al seguir determinada
trayectoria, la construcción de nuevos polinomios que permiten la interpolación cuando existe
incremento constante en la variable independiente; el diagrama de rombos es simplemente una
tabla de diferencias ampliada, que incluye los coeficientes que afectan a las diferencias, tiene la
forma siguiente:
-4-
 Tema III Interpolación, Derivación e Integración Numéricas

x-4 y-4 y-5 y-6

 k + 4  k + 5
   
 2   4 
1 y-4 y-5 y-6
 k + 3  k + 4  k + 4
     
x-3 y-3  1  y-4  3  y-5  5 

1 y-3  k + 3
  y-4  k + 4
  y-5
 2   4 

 k + 2  k + 3  k + 4
x-2 y-2  
y-3    
 1   3  y-4  5 
 k + 3
 
1 y-2  4 
 k + 2
  y-3 y-4
 2 
 k + 1  k + 3
 
x-1 y-1  1  y-2  k + 2 y-3 
 5 

 
 3 

1 y-1  k + 1
  y-2  k + 2
  y-3
 2   4 
 k
 
 1
x0 y0  k + 1
 
 k + 2
 
y-1  3  y-2  5 

1 y0
 k
  y-1  k + 1
  y-2
 k − 1  2  4 
x1 y1  
 1  y0  k  k + 1
   
 3 y-1  5 
1 y1  k − 1
  y0  k
 
 k − 2  2   4
x2 y2

 1 
 y-1
 k − 1  k
 
y1  3 
y0
 
 5
1 y2
 k − 2
  y1  k − 1
  y0
 k − 3  2   4 
x3 y3  
 1  y2  k − 2 y1  k − 1
   
 3   5 
1
y3  k − 3
  y2  k − 2
  y1
 k − 4  2   4 
x4 y4  
 1 
 k − 3  k − 2
 
y3  3 
y2

 5 

1 y4
 k − 4
  y3  k − 3
  y2
x5 y5  2 
 4 
y4 y3

para poder generar polinomios de interpolación similares al de avance de Newton porque su

primer término es y0 y los sumandos restantes están formados por un coeficiente expresado
como número combinatorio multiplicando a una diferencia en particular, a partir de un diagrama
de rombos se procede de la manera siguiente:

a) El primer término se localiza en la columna de los valores de yi, llamándose y0 (valor

inmediato anterior al que se va a interpolar.

-5-