UNIVERSIDAD NACIONAL AUTONOMA DE MEXICO
FACULTAD DE ESTUDIOS SUPERIORES CUAUTITLAN
METODOS NUMERICOS
PROFESORA: MARTHA LILIA URRUTIA VARGAS
INTERPOLACION
El problema de interpolación consisten en encontrar el valor de la función F(x), de la cual sólo se
conocen algunos puntos, para un valor de x que se encuentre entre dos valores consecutivos
conocidos.
Interpolación con espacios iguales
Considere el primer intervalo de la tabla:
x y=F(x)
x0 y0
x 1=x 0+ h y1
Figura #1
Se trata de evaluar la función F(x) en un punto intermedio xk
x y=F(x)
x0 y0
x k =x 0 +kh y k =?
x 1=x 0+ h y1
Figura #2
Sin embargo, debido a que la función F(x) se conoce únicamente en dos puntos, estos se pueden
unir mediante una recta, como se muestra en la figura siguiente
Figura #3
Si se evalúa la recta en el punto x k se tendrá una aproximación de F(x). la ecuación de la recta se
puede obtener mediante el polinomio de Newton, haciendo j=1. Se obtiene:
y k = y 0+ k y 0 + k ( 0)+ k (0)+…+ k ( 0)
1 () () ()
2 3 k ()
Al simplificar:
y k = y 0+ k ∆ y 0
x k −x 0
Al sustituir k por
h
x k −x 0
y k = y 0+ ∆ y0
h
El anterior es un polinomio de primer grado que se aproxima a F(x). Debe tomarse en cuenta que
el polinomio de Newton fue desarrollado considerando a k como un numero entero positivo, y en la
interpolación k es normalmente un valor fraccionario. Esto, sin embargo, no limita la aplicación del
polinomio.
y k = y 0+ k ∆ y 0 + k ∆2 y 0+ .. . k ∆ j y 0
1 () 1 () ()
k
Debido a que si se define la combinación ( ki) como
( k−0 )( k −1 ) ( k−2 ) . .. ( k −i+ 1 )
( ki)= i!
Donde i= 1, 2, ….j
No hay inconveniente en que k sea un numero fraccional. En general, si se cuenta con una mayor
cantidad de puntos, se puede utilizar el polinomio de Newton que más se aproxime a la función,
considerando que el grado del polinomio puede ser hasta de n-1, donde n es el número de datos.
Cuando se desea obtener la interpolación entre dos puntos de F(x), no es indispensable determinar
el polinomio de Newton en función de x k. Basta con calcular el valor de k en la ecuación y
sustituirla:
Ejemplo
Dada la tabla encontrar el valor de F(x) para x= 1.5
Solución
Al obtener las diferencias hacia adelante se forma la siguiente tabla:
Como las terceras diferencias son constantes, los datos corresponden a un polinomio de tercer
grado, éste se obtiene con la ecuación:
y k = y 0+ k ∆ y 0 + k ∆ 2 y 0 + k ∆ 3 y 0
1() 2 () 3 ()
Para x0 =1 y y0 = 1 se tiene que
∆ y 0 =8 ∆ 2 y 0 =8 ∆ 3 y 0 =6
Se sustituye
y k = y 0+ k 8+ k 8+ k 6
() () ()
1 2 3
Donde el valor de k está dado por la ecuación:
x −x0 1.5−1
k= = =0.5
h 1
y 0.5 =1+ 0.5 8+ 0.5 8+ 0.5 6
( ) ( ) ( )
1 2 3
0.5(0.5−1) 0.5(0.5−1)(0.5−2)
y 0.5 =1+ 0.5 ( 8 ) + ( 8 )+ (6)
2 6
y 0.5 =1+ 4−1+0.375=4.375
Por lo tanto
F(1.5) = 4.375
Si se desea obtener la expresión de la función F(x) definida en forma tabular, se puede interpolar
para un valor x, al sustituir éste en la ecuación y se obtiene:
x −1
k= =x−1
1
A su vez, se sustituye k en la ecuación:
y=F ( x ) =1+ x−1 8+ x −1 8+ x−1 6
1 ( ) ( ) ( )
2 3
Se efectúan operaciones
( x−1 ) ( x−2 ) ( x−1 )( x−2 )( x−3 )
F ( x )=1+8 ( x−1 ) + 8+ 6
2 6
F ( x )=1+8 x +4 x2−8 x −4 x +8 x 3−6 x 2+ 11 x−6
Y se llega a que
F ( x )=8 x 3−2 x 2+7 x −5
Ejemplo
De la siguiente función tabulada
x y=F(x)
-5 0
-2 15
1 18
4 15
7 12
1 15
0
Obtener:
a) El valor de y para x= 7.1
b) El valor de y para x= 10.4
Solución
a) El valor de y para x= 7.1
La tabla de diferencias de la función es:
En este caso no es conveniente considerar x0 = -5 y y0 =0, debido a que si se utiliza un polinomio
de tercer grado, los puntos que se tomarían en cuenta para ajustar el polinomio de Newton serían
los primeros cuatro; es decir, el ajuste cubriría el intervalo -5 a 4 y no el valor que se desea
interpolar. Por esta razón es más conveniente considerar x0=1 y y0=18; con lo cual:
F ( 7.1 )= y k =18+ k (−3)+ k 0+ k (6)
()
1 2 () ()
3
x k −x 0 7.1−1
k= = =2.033
h 3
Se sustituye
y 2.033=18+ 2.033 (−3)+ 2.033 0+ 2.033 (6)
1 ( ) 2 ( 3) ( )
Se efectúan operaciones y se obtiene:
y 2.033=15.967
b) El valor de y para x= 10.4
Se ha calculado el valor de y para una x dada entre dos valores de x k de la tabla que define a la
función. Ahora se trata de calcular el valor de y para una x fuera del rango de los valores de x k de la
tabla
A este problema se le conoce como extrapolación. Dado que para x=13 no se tienen definidas las
diferencias de la función, habrá que determinarlas aceptando que las terceras diferencias son
constantes e iguales a 6. A partir de la tercera diferencia se puede determinar la segunda; con la
segunda la primera, y con esta ultima el valor de la función para x= 13. En la tabla anterior se han
marcado con un asterisco los valores encontrados con este procedimiento.
Al interpolar
10.4−13
K= =0.867
−3
y 0.867=30+ 0.867 (−15)+ 0.867 12+ 0.867 (−6)
1 ( ) 2 ( ) (
3 )
Al efectuar operaciones:
Y= 30 -13.005 - 0.692 – 0.131 = 16.172
Por lo tanto:
F(10.4) = 16.172
