11.08.2013 Views

Apunte

Apunte

Apunte

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!