1

Capı́tulo 3: Interpolación y aproximación polinómica

Interpolación y polinomio de Lagrange

Considere las siguientes funciones que definen un polinomio de primer grado que
pasa por los puntos (x0 , y0 ) y (x1 , y1 )
x − x1 x − x0
L0 (x) = y L1 (x) =
x0 − x1 x1 − x0
Entonces se define un polinomio como

P (x) = L0 (x)f (x0 ) + L1 (x)f (x1 )

Como

L0 (x0 ) = 1 y L0 (x1 ) = 0
L1 (x0 ) = 0 y L1 (x1 ) = 1

Entonces se tiene que

P (x0 ) = 1· f (x0 ) + 0· f (x1 ) = f (x0 ) = y0

P (x1 ) = 0· f (x0 ) + 1· f (x1 ) = f (x1 ) = y1

Para poder generalizar este resultado, consideremos la construcción de un polino-

mio de grado máximo n que pase por los n + 1 puntos.

(x0 , f (x0 )), (x1 , f (x1 )), ..., (xn , f (xn ))

Como k = 0, 1, 2, ..., n construimos la función Lnk (x) con la propiedad que Lnk (xi ) = 0
cuando i 6= k y Lnk (xk ) = 1. Para satisfacer que Lnk (xi ) = 0 para cada i 6= k se requiere
que el numerador de Lnk (x) contenga el término

(x − x0 )(x − x1 ) · · · (x − xk−1 )(x − xk+1 ) · · · (x − xn )

 2

Para satisfacer Lnk (xk ) = 1 el denominador de Lnk (x) de coincidir con el tértino
cuando se evalúe en x = xk es decir,
(x − x0 )(x − x1 ) · · · (x − xk−1 )(x − xk+1 ) · · · (x − xn )
Lnk =
(xk − x0 )(xk − x1 ) · · · (xk − xk−1 )(xk − xk+1 ) · · · (xk − xn )

Teorema: Si x0 , x1 , ..., xn son n + 1 números distintos y si f es una función cuyos

valores están dados en sus números, entonces existen un único polinomio P (x) de
grado a lo más n, con la propiedad de que

f (xk ) = P (xk ) para cada k = 0, 1, ..., n

Este polinomio está dado por

n
X
P (x) = f (x0 )Ln,0 + · · · + f (xn )Ln,n = f (xk )Ln,k (x)
k=0

donde para cada k = 0, 1, ..., n

(x − x0 )(x − x1 ) · · · (x − xk−1 )(x − xk+1 ) · · · (x − xn )
Lnk =
(xk − x0 )(xk − x1 ) · · · (xk − xk−1 )(xk − xk+1 ) · · · (xk − xn )
n
Y (x − xi )
=
i=0
(xk − xi )
i6=k
Escribimos Lnk (x) simplemente como Lk (x) cuando no haya confusión respecto a
su grado.

Ejemplo: Si queremos utilizar los números x0 = 2, x1 = 2, 5 y x2 = 4, para obte-

1
ner el segundo polinomio interpolante para f (x) = .
x

Teorema: Supongamos que x0 , x1 , ..., xn son números distintos en el intervalo [a, b] y

que f ∈ C n+1 [a, b]. Entonces para cada x en [a, b] existe un número ξ(x) en ]a, b[ con
f n+1 (ξ(x))
f (x) = P (x) + (x − x0 )(x − x1 ) · · · (x − xn )
(n + 1)!
 3

Donde P (x) es la ecuación interpolante de Lagrange

Nótese que la cota de error para el polinomio de Lagrange esta dada por

f n+1 (ξ(x))
(x − x0 )(x − x1 ) · · · (x − xn )
(n + 1)!

Ejemplo: Considere la siguiente tabla, que muestra valores de una función en di-
versos puntos. Aproxime f (1, 5) por medio de Lagrange.
x f (x)
1,0 0,7651977
1,3 0,6200860
1,6 0,4554022
1,9 0,2818186
2,2 0,1103623

Definición: Sea f una función definida en x0 , x1 , x2 , ..., xn y supongamos que m1 , m2 , ..., mk

son k enteros distintos con 0 6 mi 6 n para cada i. El polinomio de Lagrange que
concuerda con f en los k puntos xmi , xm2 , ..., xmk se denota por Pm1 ,m2 ,...,mn (x)

Ejemplo: Si x0 = 1, x1 = 2, x2 = 3, x3 = 4, x4 = 6 y si f (x) = ex , obtenga P1,2,4 (x)

 4

Teorema: Si f está definida en x0 , x1 , ..., xk y xj y xi son dos números distintos de

este conjunto entonces
(x − xj )P0,1,...,j−1,j+1,...,k (x) − (x − xi )P0,1,...,i−1,i+1,...,k (x)
P (x) =
(xi − xj )
Describe el polinomio de grado k de Lagrange que interpola f en las k + 1 puntos
x0 , x1 , ..., xk

Como conseguencia al teorema anterior, los polinimios interpolantes pueden ge-

nerarse de manera recursiva. Los ponemos generar con la siguiente tabla:

x0 P0 = Q0,0
x1 P1 = Q1,0 P0,1 = Q1,1
x2 P2 = Q2,0 P1,2 = Q2,1 P0,1,2 = Q2,2
x3 P3 = Q3,0 P2,3 = Q3,1 P1,2,3 = Q3,2 P0,1,2,3 = Q3,3
x4 P4 = Q4,0 P3,4 = Q4,1 P2,3,4 = Q4,2 P1,2,3 = Q4,3 P0,1,2,3,4 = Q4,4

A este proceso se le conoce con el nombre de método de Neville

Ejemplo: Los valores de polinomio interpolan en x = 1, 5. Aproxime f (1, 5) usan-

do x0 = 1, 0, x1 = 1, 3, x2 = 1, 6, x3 = 1, 9, x4 = 2, 2

Ejemplo: Considere los siguientes valores de f (x) = ln(x)

i xi ln(xi )
0 2 0,6931
1 2,2 0,7885
2 2,3 0,8329

Aproximar f (2, 1)
 5

Diferencias divididas

Supongamos que Pn (x) es el n-ésimo polinimio de Lagrange que concuerda con

la función f en los números distintos x0 , x1 , ..., xn . Las diferencias divididas de f
respecto a x0 , x1 , ..., xn se usan para expresar Pn (x) en la forma

Pn (x) = a0 + a1 (x − x0 ) + a2 (x − x0 )(x − x1 ) + · · · + an (x − x0 )(x − x1 ) · · · (x − xn−1 )

Para las constantes apropiadas a0 , a1 , ..., an

Se realizará la idea esencial de como usar diferencias divididas de manera general.

Note que para determinar a0 se obtiene al sustituir Pn (x0 ) = f (x0 ) = a0 ,

Ahora de manera similar al evaluar P (x1 ) se obtiene a0 + (x1 − x0 ) = Pn (x1 ) = f (x1 )
f (x1 ) − f (x0 )
Ahora como a0 = f (x0 ) y despejando a1 se obtiene a1 =
x1 − x0

El resto de las diferencias divididas se definen de manera inductiva, es decir,

La diferencia dividida cero: f [xi ] = f [xi ]

f [xi+1 ] − f [xi ]
La primera diferencia dividida f [xi , xi+1 ] =
xi+1 − xi
f [xi+1 , xi+2 ] − f [xi , xi+1 ]
La segunda diferencia dividida f [xi , xi+1 , xi+2 ] =
xi+2 − xi
Siguiendo el proceso
La k-ésima diferencia dividida
f [xi+1 , xi+2 , ..., xi+k ] − f [xi , xi+1 , ..., xi+k−1 ]
f [xi , xi+1 , ..., xi+k−1 , xi+k ] =
xi+k − xi

Con esta notación podemos expresar a1 = f [x0 , x1 ] y ahora el polinimio interpolante

se puede expresar como

Pn (x) = f [x0 ] + f [x0 , x1 ](x − x0 ) + a2 (x − x0 )(x − x1 ) + · · · + an (x − x0 )(x − x1 ) · · · (x − xn−1 )

como cabe suponer tras evaluar a0 y a1 , las constantes requeridas son

ak = f [x0 , x1 , ..., xk ]
 6

Para cada k = 0, 1, 2..., n, por lo tanto se puede reescribir Pn (x) como

n
X
Pn (x) = f [x0 ] + f [x0 , x1 , ..., xk ](x − x0 ) · · · (x − xk−1 )
k=1

A esta ecuación se le conoce con el nombre de Fórmula de diferencias divididas

interpolantes de Newton
Nota: Para ver un esquema de las diferencias divididas ver el libro de Burden pági-
na 124

Ejemplo: Construya el polinomio de diferencias divididas según la siguiente in-

formación

i xi f [xi ]
0 1 0,7651977
1 1,3 0,6200860
2 1,6 0,4554022
3 1,9 0,2818186
4 2,2 0,1103623

Teorema: Supongamos que f ∈ C n [a, b] y x0 , x1 , ..., xn son números distintos en [a, b].
Entonces existe un número ξ en ]a, b[ con
f nξ
f [x0 , x1 , ..., xn ] =
n!

La fórmula de las diferencias divididas interpolantes de Newrton se pueden expre-

sar de forma simplificada, si introducimos h = xi+1 − xi para cada i = 0, 1, 2, ..., n − 1
y sea x = x0 + sh entonces se obtiene

Pn (x) = Pn (x0 + sh)

f [x0 ] + shf [x0 , x1 ] + s(s − 1)h2 f [x0 , x1 , x2 ] + · · · + s(s − 1)(s − n + 1)hn f [x0 , x1 , ..., xk ] =
Xn
s(s − 1)(s − k + 1)hk f [x0 , x1 , ..., xk ]
k=0
 7

Si utilizamos la notación del coeficiente binomial

 
s s(s − 1) · · · (s − k + 1)
=
k k!
Podemos expresar Pn (x) en forma más compacta como:
n  
X s
Pn (x) = Pn (x0 + sh) = f [x0 ] + k!hk f [x0 , x1 , ..., xk ]
k=1
k

A ésta ecuación se le llama fórmula de las diferecias divididas progresivas de

Newton

Definición: Dada la sucesión {pn }∞

n=0 defina la diferencia regresiva 5pn (nabla pn )

por medio de

5pn = pn − pn+1 para n > 1

Las potencias mayores se definen recursivamente por

5k pn = 5(5k−1 pn ) para k > 2

Entonces la ecuación de diferencias divididas interpolantes de Newton, se puede

expresar como

1
f [xk , xn−1 ] = 5 f (xn )
h
1
f [xn , xn−1 , xn−2 ] = 2 52 f (xn )
2h
y en general
1
f [xn , xn−1 , ..., xn−k ] = k
5k f (xk )
k!h
por lo tanto
s(s + 1) 2 s(s + 1) · · · (s + n − 1) n
Pn (x) = f [xn ] + s 5 f (xn ) + 5 f (xn ) + · · · + 5 f (xn )
2 n!
Usando la notación binomial
 
−s −s(−s − 1) · · · (−s − k + 1) s(s + 1) · · · (s + k − 1)
= = (−1)k
k k! k!
 8

Entonces la fórmula de diferencias regresiva de Newton es

n  
X
k −s
Pn = f (xn ) + (−1) 5k f (xn )
k=1
k

Ejemplo: Usando los datos de la tabla anterior aproxime f (1, 1) y f (2, 0) usando las
fórmulas progresivas y regresivas de Newton

Fórmula de diferencias centrada

Fórmula de Stirling

sh
Pn (x) = P2m+1 (x) = f [x0 ] + (f [x−1 , x0 ] + f [x0 , x1 ]) + s2 h2 f [x−1 , x0 , x1 ]
2
s(s2 − 1)h3
+ (f [x−2 , x−1 , x0 , x1 ] + f [x−1 , x0 , x1 , x2 ])
2
+ · · · + s2 (s2 − 1)(s2 − 4) · · · (s2 − (m − 1)2 )h2m f [x−m , ..., xm ]
s(s2 − 1) · · · (s2 − m2 )h2m+1
+ (f [x−m−1 , ..., xn ] + f [x−m , ..., xm+1 ])
2

Ejemplo: Con los datos anteriores aproxime f (1, 5)

 9

Interpolación de Hermite

Las polinomios oculantes representan una generalazación de los polinomios de

Taylor y de Lagrange
Definición: Sean x0 , x1 , ..., xn , n + 1 números distintos de [a, b] y mi un entero no
negativo asociado a xi para i = 0, 1, 2, ..., n. Supongase que f ∈ C m [a, b] y que m =
máx06i6n mi . El polinomio osculante que aproxima f es el polinonio P (x) de menor
grado tal que

dk (P (xi )) dk f (xi )
= para cadai = 0, 1, ..., n y k = 0, 1, ..., mi
dxk dxk
Teorema: Si f ∈ C 1 [a, b] y si x0 , ..., xn ∈ [a, b] son distintos, el polimonio único de
menor grado que concuerda con f y f 0 en x0 , ..., xn es el polinimio de Hermite de
grado a lo más 2n + 1 que esta dado por
n
X n
X
H2n+1 (x) = f (xj )Hn,j (x) + f 0 (xj )H
b n,j (x)
j=0 j=0

donde
Hn,j (x) = 1 − 2(x − xj )L0n,j (xj ) L2n,j (x)
 

y
b n,j (x) = (x − xj )L2 (x)
H n,j

Dentro de este contexto Ln,j (x) denota el j-ésimo polinomio de Lagrange de grado n.
Si f ∈ C 2n+2 [a, b] entonces para x ∈ [a, b]

(x − x0 )2 · · · (x − xn )2 (2n+2)
f (x) = H2n+1 (x) + f (ξ)
(2n + 2)!
Para alguna ξ con a < ξ < b
Ejemplo: Utilice el polinomio de Hermite, para aproximar f (1, 5)

k xk f (xk ) f 0 (xk )
0 1,3 0,6200860 -0,5220232
1 1,6 0,4554022 -0,5698959
2 1,9 0,2818186 -05811571
 10

Este es un proceso muy tedioso entonces creando una nueva sucesión z0 , z1 , ..., z2n+1
entonces
z2i = z2i+1 = xi para cada i = 0, 1, 2, ..., n

entonces
f [z2i , z2i+1 ] = f 0 (z2i ) = f 0 (xi )

entonces el polinomio de Hermite está dado por

2n+1
X
H2n+1 (x) = f [z0 ] + f [z0 , ..., zk ](x − z0 )(x − z1 ) · · · (x − zk−1 )
k=1

Ejemplo: Obtengamos H5 (x) según los valores del ejemplo anterior, para aproximar
f (1, 5)
 11

Trazador Cúbico y curvas paramétricas

Definición: Dada una función en [a, b] y un conjunto de nodos a = x0 < x1 < · · · <
xn = b un interpolante de trazador cúbico S para f es una función que cumple con
las siguientes condiciones:

1. S(x) es un polinomio cúbico, denotado Si (x) en el subintervalo [xi − xi+1 ] para

cada i = 0, 1, ..., n − 1

2. S(xi ) = f (xi ) para cada i = 0, 1, ..., n

3. Si+1 (xi+1 ) = Si (xi+1 ) para cada i = 0, 1, ..., n − 2

0
4. Si+1 (xi ) = Sj0 (xi+1 ) para cada i = 0, 1, ..., n − 2

00
5. Si+1 (xi+1 ) = Si00 (xi+1 ) para cada i = 0, 1, ..., n − 2

6. Una de las siguientes condiciones de frontera se satisface:

a) S 00 (x0 ) = S 00 (xn ) = 0(Frontera libre o natural)

b) S 0 (x0 ) = f 0 (x0 ) y S 0 (xn ) = f 0 (xn ) (frontera sujeta)

Si queremos construir el interpolante del trazador cúbico de determinada función

f , aplicamos las de la definición a los polinomios cúbicos:

Si = ai + bi (x − xi ) + ci (x − xi )2 + di (x − xi )3

Este tema se verá todo en la pizarra