UNIDAD 1.

INTERPOLACIN

UNIDAD 1

INTERPOLACIN

1. Introduccin.
La idea bsica de la interpolacin es hallar un polinomio o funcin que cumpla con pasar
por todos los puntos de datos (1 , 1 ), (2 , 2 ), , ( , ), y poder estimar los valores
entre ellos por medio de un polinomio.
2. Interpolacin de Lagrange.
Para un conjunto de puntos, el polinomio de Lagrange es:

( )
= () =
= =
( )

donde
( , ) = Coordenadas del punto interpolado.
= Nmero de puntos empleados en el proceso. El grado del polinomio es 1 y lo ms
usual es emplear entre tres y cinco puntos.
, = Sirven para enumerar los puntos conocidos y los trminos en la sumatoria y en el
producto.
= = Sumatoria que representa la suma de los trminos colocados a su derecha.

= = Indica multiplicacin.

La expresin anterior se puede reescribir como

() = ()
=0

donde


() =

=0

es llamado polinomios fundamentales de Lagrange.

M. GUERRERO RODGUEZ 1
 UNIDAD 1. INTERPOLACIN

Actividad 1. Escriba los polinomios de Lagrange para dos puntos, tres puntos y cuatro
puntos.
Ejemplo 1. Dados los puntos (9.0, 2.1972), (9.5, 2.2513), (11.0, 2.3979), obtener un
polinomio cuadrtico, 2 , al introducir el vector u de los valores de , el vector v de los
valores de . Calcular el valor de 2 , cuando = 9.2 .
Solucin.
La respuesta se presenta en Excel.

x=u y =v
9 2.1972
9.5 2.2513
11 2.3979

La grfica de los valores anteriores, se muestra a continuacin, ntese que esta grfica no
es totalmente una lnea recta.

2.45

2.4

2.35

2.3
y

2.25

2.2

2.15
8.5 9 9.5 10 10.5 11 11.5
x

Un polinomio de segundo orden tiene la forma:

= () = 0 + 1 + 2 2
Tomando como referencia los puntos proporcionados, se pueden obtener las ecuaciones
siguientes:

M. GUERRERO RODGUEZ 2
 UNIDAD 1. INTERPOLACIN

(9,2.1972) 2.1972 = 0 + 1 (9) + 2 (9)2 = 0 + 91 + 812

(9.5,2.2513) 2.2513 = 0 + 1 (9.5) + 2 (9.5)2 = 0 + 9.51 + 90.252
(11,2.3979) 2.3979 = 0 + 1 (11) + 2 (11)2 = 0 + 111 + 1212

La solucin del sistema de ecuaciones anterior se presenta enseguida, el cual es resuelto en

Excel.

1 9 81 1 9 81
= 1 9.5 90.25 1 9.5 90.25
1 11 121 1 11 121

104.5 -132 28.5

1 = -20.5 26.6666667 -6.16666667
1 -1.33333333 0.33333333

2.1972 2.1972
= 2.2513 2.2513
2.3979 2.3979

0 0.77595
= 1 0.20501667
2 -0.00523333

El polinomio de segundo orden queda de la forma siguiente:

() = . + . .

La grfica del polinomio, en Geogebra, se muestra enseguida:

M. GUERRERO RODGUEZ 3
 UNIDAD 1. INTERPOLACIN

Evaluando el polinomio para = 9.2:

(9.2) = 0.77595 + 0.20501667(9.2) 0.00523333(9.2)2

p(9.2) = 2.219154

Ejemplo 2. Dados los valores de la siguiente tabla

1 2 3 5
1.54 1.5 1.42 0.66

Encontrar el polinomio de interpolacin de Lagrange.

Solucin.
El diagrama de dispersin de estos puntos se presenta a continuacin en Geogebra

M. GUERRERO RODGUEZ 4
 UNIDAD 1. INTERPOLACIN

Polinomios fundamentales de Lagrange.

( 1 )( 2 )( 3 ) ( 2)( 3)( 5)
0 () = =
(0 1 )(0 2 )(0 3 ) (1 2)(1 3)(1 5)
1
= ( 3 10 2 + 31 30)
8
( 0 )( 3 )( 4 ) ( 1)( 3)( 5)
1 () = =
(1 0 )(1 3 )(1 4 ) (2 1)(2 3)(2 5)
1
= ( 3 9 2 + 23 15)
3
( 0 )( 1 )( 4 ) ( 1)( 2)( 5)
2 () = =
(3 0 )(3 1 )(3 4 ) (3 1)(3 2)(3 5)
1
= ( 3 8 2 + 17 10)
4
( 1 )( 2 )( 3 ) ( 1)( 2)( 3)
3 () = =
(4 1 )(4 2 )(4 3 ) (5 1)(5 2)(5 3)
1 3
= ( 6 2 + 11 6)
24

El polinomio de interpolacin de Lagrange es

() = 0 0 () + 1 1 () + 2 2 () + 3 3 ()

Sustituyendo en la ecuacin anterior cada uno de los polinomios fundamentales de

Lagrange y los valores de (proporcionados en la tabla), resulta que el polinomio tiene la
forma siguiente:

() = . + . . + .

La grfica del polinomio anterior, en Geogebra, es:

M. GUERRERO RODGUEZ 5
 UNIDAD 1. INTERPOLACIN

A continuacin se presenta la solucin en wxMaxima, se obtiene el diagrama de dispersin,

el polinomio de Lagrange y su representacin grfica.

M. GUERRERO RODGUEZ 6
 UNIDAD 1. INTERPOLACIN

Diagrama de dispersin.

Grfica del polinomio.

M. GUERRERO RODGUEZ 7
 UNIDAD 1. INTERPOLACIN

Actividad 2. Resolver el ejemplo 1 haciendo uso de una interpolacin de Lagrange.

Actividad 3. Un ingeniero debe resolver un problema de fsica, para lo cual requiere conocer
la densidad del agua a 43.70 C con la mayor precisin posible y recurre a una tabla de datos
de un texto:

t, 0C Densidad, g/cm3 Volumen de 1 g/cm3

0 0.99987 1.00013
2 0.99997 1.00003
4 1.00000 1.00000
6 0.99997 1.00003
10 0.99973 1.00027
20 0.99823 1.00177
50 0.98807 1.01207
75 0.97489 1.02576
100 0.95838 1.04343

Emple:
a) Interpolacin grfica.
b) Un polinomio de segundo grado: = 0 + 1 + 2 2
c) Un polinomio de Lagrange de tres puntos.

3. Mtodo de diferencias divididas (Polinomio de Interpolacin de Newton).

Definicin.

Dado 0, se tiene [ ] = , = 0,1, , , al que se denomina diferencia dividida de orden

de en el punto , y se representa por [ , +1 , , + 1 , + ], al cociente:

[ , +1 , , + 1] [+1 , , + 1 , + ]
[ , +1 , , + 1 , + ] =
+

A partir de las diferencias divididas, el polinomio de interpolacin se determina como

() = [0 ] + ([0 , 1 , , ] ( ))
=1 =0

M. GUERRERO RODGUEZ 8
 UNIDAD 1. INTERPOLACIN

el cual se conoce como polinomio de interpolacin de Newton.

Al agregar un nuevo par de interpolacin (+1 , +1 ), el nuevo polinomio es

+1 () = () + [0 , 1 , , +1 ] ( )
=0

Ejemplo 3. Dada la siguiente tabla, obtener el polinomio de interpolacin de Newton.

2 4 6 8
4 8 14 16

Solucin.

El ejemplo es resuelto en Excel.

x f
2 4
2,4 = 2
4 8 2,4,6 = 0.25
4,6 = 3 2,4,6,8 = -0.125
6 14 4,6,8 = -0.5
6,8 = 1
8 16

Polinomio de interpolacin de Newton.

3 = 4 + 2 2 + 0.25 2 4 0.125 2 4 6

= . + . +

x f p3(x)
2 4 4
4 8 8
6 14 14
8 16 16
-10 358

M. GUERRERO RODGUEZ 9
 UNIDAD 1. INTERPOLACIN

-5 92.375
0 8
10 8
11 -1.625
20 -392

A continuacin se presenta un programa en wxMaxima para resolver este ejemplo.

M. GUERRERO RODGUEZ 10
 UNIDAD 1. INTERPOLACIN

4. Interpolacin con splines cbicos.

Una funcin spline est formada por varios polinomios, cada uno definido en un intervalo y que se
unen bajo ciertas condiciones de continuidad. Se ha observado que las splines cbicas son las ms
adecuadas.

Dada una particin

= { = 0 < 1 < < = }

se interpola con funciones de clase dos, restringidas a los intervalos [ , +1 ] con polinomios de
grado tres. Si () es la restriccin del spline () en el intervalo [ , +1 ], = 0,1, , 1,
entonces

() = 3 + 2 + + = 0,1, , 1

El spline cumple con las condiciones siguientes:

a) Condiciones de interpolacin.
( ) = ( ), = 0,1, , 1, 1 ( ) = ( )

M. GUERRERO RODGUEZ 11
 UNIDAD 1. INTERPOLACIN

b) Condiciones de continuidad (en nodos interiores).

(+1 ) = +1 (+1 ), = 0,1, , 2
c) Condiciones de suavidad (en nodos interiores).
(+1 ) = +1 (+1 ), = 0,1, , 2
(+1 ) = +1 (+1 ), = 0,1,2, , 2

Ntese que el nmero de ecuaciones es 4 2, lo cual nos indica que para determinar el spline
() de forma nica es necesario tener dos condiciones adicionales. Las condiciones ms habituales,
impuestas sobre los extremos del intervalo, son:

0 () = 0, 1 () = 0, (spline cbico natural),

0 () = (), 1 () = (), (spline cbico sujeto)

Ejemplo 4. Construir un spline cbico sujeto que se ajuste a los datos de la siguiente tabla

0 1 2 3
0 1/2 2 3/2

y que compruebe las condiciones (0) = 1/5 y (3) = 1.

Para ello indique como se realicen los clculos manualmente, escriba un programa en wxMaxima y
elabore la grfica respectiva.

Solucin.

Usando Excel se obtiene el diagrama de dispersin de los valores de la tabla anterior.

x f
0 0
1 0.5
2 2
3 1.5

M. GUERRERO RODGUEZ 12
 UNIDAD 1. INTERPOLACIN

Modelos matemticos del spline cbico para la tabla de datos proporcionada.

1 3 + 1 2 + 1 + 1 si [0,1]
() = {2 3 + 2 2 + 2 + 2 si [1,2]
3 3 + 3 2 + 3 + 3 si [2,3]

Las ecuaciones que se obtienen al sustituir los valores de cada intervalo, en base al modelo anterior,
son:

1 = 0

1 + 1 + 1 + 1 = 0.5

2 + 2 + 2 + 2 = 0.5

82 + 42 + 22 + 2 = 2

83 + 43 + 23 + 3 = 2

273 + 93 + 33 + 3 = 1.5

Ecuaciones de la primer derivada de ().

31 2 + 21 + 1 si [0,1]
() = {32 2 + 22 + 2 si [1,2]
33 2 + 23 + 3 si [2,3]

M. GUERRERO RODGUEZ 13
 UNIDAD 1. INTERPOLACIN

De los intervalos se observa que las posibles discontinuidades de la derivada estn en = 1 y en

= 2.

Las siguientes ecuaciones permiten hacer que la primera derivada sea continua:

31 + 21 + 1 = 32 + 22 + 2

122 + 42 + 2 = 123 + 43 + 3

Ecuaciones de la segunda derivada de ().

61 + 21 si [0,1]
() = {62 + 22 si [1,2]
63 + 23 si [2,3]

Tambin, se observa que las posibles discontinuidades de la segunda derivada estn en = 1 y en

= 2.

Para que la segunda derivada sea continua se obtienen las ecuaciones siguientes:

61 + 21 = 62 + 22 31 + 1 = 32 + 2

122 + 22 = 123 + 23 62 + 2 = 63 + 3

Condiciones sobre la primera derivada:

(0) = 0.2 1 = 0.2

(3) = 1 273 + 63 + 3 = 1

Resumen de ecuaciones:

1 = 0

1 + 1 + 1 + 1 = 0.5

2 + 2 + 2 + 2 = 0.5

82 + 42 + 22 + 2 = 2

83 + 43 + 23 + 3 = 2

273 + 93 + 33 + 3 = 1.5

M. GUERRERO RODGUEZ 14
 UNIDAD 1. INTERPOLACIN

31 + 21 + 1 = 32 + 22 + 2

122 + 42 + 2 = 123 + 43 + 3

31 + 1 = 32 + 2

62 + 2 = 63 + 3

1 = 0.2

273 + 63 + 3 = 1

Enseguida se tiene el sistema de ecuaciones en forma matricial.

Matriz de coeficientes.

0 0 0 1 0 0 0 0 0 0 0 0
1 1 1 1 0 0 0 0 0 0 0 0
0 0 0 0 1 1 1 1 0 0 0 0
0 0 0 0 8 4 2 1 0 0 0 0
0 0 0 0 0 0 0 0 8 4 2 1
= 0 0 0 0 0 0 0 0 27 9 3 1
3 2 1 0 -3 -2 -1 0 0 0 0 0
0 0 0 0 12 4 1 0 -12 -4 -1 0
3 1 0 0 -3 -1 0 0 0 0 0 0
0 0 0 0 6 1 0 0 -6 -1 0 0
0 0 1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 27 6 1 0

Matriz inversa.

1.2 -1.2 -0.6 0.6 0.2 -0.2 0.46666667 -0.13333333 0.26666667 -0.06666667 0.73333333 0.06666667
-2.2 2.2 0.6 -0.6 -0.2 0.2 -0.46666667 0.13333333 -0.26666667 0.06666667 -1.73333333 -0.06666667
0 0 0 0 0 0 0 0 0 0 1 0
1 0 0 0 0 0 0 0 0 0 0 0
-0.6 0.6 0.8 -0.8 -0.6 0.6 -0.4 0.4 0.2 0.2 -0.2 -0.2
1=
3.2 -3.2 -3.6 3.6 2.2 -2.2 2.13333333 -1.46666667 -1.06666667 -0.73333333 1.06666667 0.73333333
-5.4 5.4 4.2 -4.2 -2.4 2.4 -3.6 1.6 1.8 0.8 -1.8 -0.8
2.8 -2.8 -0.4 1.4 0.8 -0.8 1.86666667 -0.53333333 -0.93333333 -0.26666667 0.93333333 0.26666667
0.2 -0.2 -0.6 0.6 1.2 -1.2 0.13333333 -0.46666667 -0.06666667 0.26666667 0.06666667 0.73333333
-1.6 1.6 4.8 -4.8 -8.6 8.6 -1.06666667 3.73333333 0.53333333 -2.13333333 -0.53333333 -4.86666667
4.2 -4.2 -12.6 12.6 19.2 -19.2 2.8 -9.8 -1.4 5.6 1.4 10.4
-3.6 3.6 10.8 -10.8 -12.6 13.6 -2.4 8.4 1.2 -4.8 -1.2 -7.2

M. GUERRERO RODGUEZ 15
 UNIDAD 1. INTERPOLACIN

Vector de trminos independientes.

0
0.5
0.5
2
2
= 1.5
0
0
0
0
0.2
-1

Coeficientes del polinomio.

a1 = 0.48
b1 = -0.18
c1 = 0.2
d1 = 0
a2 = -1.04
b2= 4.38
c2 = -4.36
d2 = 1.52
a3 = 0.68
b3 = -5.94
c3 = 16.28
d3 = -12.24

Spline cbicos.

Grfica de los spline cbicos.

M. GUERRERO RODGUEZ 16
 UNIDAD 1. INTERPOLACIN

x s0(x) s1(x) s2(x)

0 0
0.2 0.03664
0.4 0.08192
0.6 0.15888
0.8 0.29056
1 0.5 0.5
1.2 0.79808
1.4 1.14704
1.6 1.49696
1.8 1.79792
2 2 2
2.2 2.06704
2.4 2.01792
2.6 1.88528
2.8 1.70176
3 1.5

Nota. Escribir el programa respectivo para el problema anterior en wxMaxima.

4. Interpolacin lineal.

Consiste en trazar una recta entre cada par de puntos de datos, razn por la cual tambin es llamada
interpolacin por trazadores lineales o splines de primer orden.

M. GUERRERO RODGUEZ 17
 UNIDAD 1. INTERPOLACIN

Para el siguiente conjunto de datos (1 , 1 ), (2 , 2 ), , ( , ), entre dos puntos consecutivos del

conjunto de datos se tiene una recta cuya pendiente es

+1
=
+1

Y que pasa por el conjunto inicial ( , ), entonces la ecuacin de la recta que interpola entre ese
par de puntos es

+1
= ( ) +
+1

Ntese que la interpolacin lineal se hace por pedazos y no entrega un solo polinomio para todo el
conjunto de datos como en el caso de los polinomios de Lagrange.

Ejemplo 5. Dados los puntos (9.0, 2.1972), (9.5, 2.2513), (11.0, 2.3979), tomados del ejemplo 1.
Calcular el valor (mediante una interpolacin lineal) cuando = 9.2 .

Solucin.

El valor que se desea interpolar se encuentra entre los dos primeros puntos por lo que los valores
de los parmetros son:

= 9

= 2.1972

+1 = 9.5

+1 = 2.2513

= 9.2

Sustituyendo en la ecuacin anterior resulta:

2.2513 2.1972
= (9.2 9) + 2.1972
9.5 9

= 2.21884

M. GUERRERO RODGUEZ 18
 UNIDAD 1. INTERPOLACIN

Actividad 4. Empleando interpolacin lineal calcule el valor del logaritmo de 4.7 . Use los datos
siguientes:

Nm. Logaritmo10
3 0.477121
4.5 0.653212
5 0.698970
a) Interpolar entre 3 y 5.
b) Interpolar entre 4.5 y 5.
Calcular el error relativo porcentual para ambos incisos, basado en el valor verdadero,
log10(4.7) = 0.672097857936

Actividad 5. La siguiente tabla muestra los valores de entalpa de vapor sobrecalentado a diferentes
temperaturas y una presin de 60 bar. Utilice polinomios de Lagrange de tercero y cuarto orden
para determinar la entalpa a 110 0C, realiza una interpolacin lineal sobre la zona ms prxima y
compara resultados.

Temperatura,
50 75 100 150 200 250
0
C
Entalpa,
214.4 318.7 423.5 635.6 854.2 1085.8
Kj / kg

M. GUERRERO RODGUEZ 19