You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2. ERROR DE INTERPOLACIÓN 93<br />
reemplazo tenga alguna validez numérica es importante conocer una estimación del error que<br />
se comete. Para esto será necesario suponer que la función f verifica algunas condiciones de<br />
suavidad. Llamemos a este error:<br />
En(x) = f(x) − pn(x), x ∈ [a, b].<br />
Con el siguiente teorema, damos el primer paso para poder estimar el error cometido; es decir,<br />
damos una expresión para En(x).<br />
Dados los puntos x0, . . . , xn, utilizaremos la notación Wn+1 para designar al polinomio mónico<br />
de grado n + 1 que se anula en esos puntos. Es decir,<br />
Wn+1(x) = (x − x0) · · · (x − xn)<br />
Teorema 5.4. Sean f ∈ C n+1 [a, b] y pn ∈ Pn el polinomio interpolador de f en x0, . . . , xn<br />
puntos del intervalo [a, b]. Para cada x ∈ [a, b], existe ξ ∈ [a, b], ξ = ξ(x), tal que<br />
En(x) = f(x) − pn(x) = f (n+1) (ξ)<br />
(n + 1)! Wn+1(x).<br />
Demostración. Notar que En(xj) = 0 y Wn+1(xj) = 0 para todo j. Por lo tanto, podemos<br />
suponer x = xj. Fijado x definimos la siguiente función de t,<br />
F (t) = f(t) − pn(t) − αWn+1(t)<br />
donde α se elige de modo que F (x) = 0. O sea, α = f(x)−pn(x)<br />
Wn+1(x) , que está bien definida pues<br />
Wn+1(x) = 0. Observemos que para todo j,<br />
F (xj) = f(xj) − pn(xj) − αWn+1(xj) = 0.<br />
Entonces F se anula en los n + 2 puntos x0, . . . , xn, x. En consecuencia, por el teorema de Rolle,<br />
F ′ tiene al menos n + 1 ceros, F ′′ al menos n ceros y así siguiendo se tiene que existe un punto<br />
ξ ∈ (a, b) tal que F (n+1) (ξ) = 0. Como<br />
Se obtiene,<br />
lo que concluye la demostración.<br />
F (n+1) (t) = f (n+1) (t) − (n + 1)!α<br />
f (n+1) (ξ)<br />
(n + 1)!<br />
= f(x) − pn(x)<br />
Wn+1(x)<br />
Ejemplo 5.5. Analicemos qué sucede si se quiere interpolar la función f(x) = cos(x) 3 en el<br />
intervalo [−3, 3] por un polinomio.