11.08.2013 Views

Apunte

Apunte

Apunte

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

94 5. INTERPOLACI ÓN<br />

Si se eligen 10 puntos equiespaciados se obtiene un polinomio como muestra la Figura 5.2.<br />

Si consideramos el error numérico, que es el que se obtiene como diferencia máxima entre f<br />

y el polinomio evaluados en una malla suficientemente fina (puntos equiespaciados con paso<br />

h = 0,01) se tiene un error de 0.4303. . . Tan sólo al considerar 25 puntos equiespaciados (tomados<br />

a intervalos de longitud 0.25) se obtiene un error numérico menor que 10 −6 . En este caso, en<br />

una figura como la anterior, los gráficos del polinomio y la función se confunden.<br />

4<br />

3<br />

2<br />

1<br />

0<br />

−1<br />

−2<br />

−3<br />

−4<br />

−5<br />

−6<br />

−1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1<br />

Figura 5.2. Interpolación de f(x) = cos(x) 3 en [−3, 3], (10 puntos equiespaciados)<br />

3. Forma de Newton<br />

La forma de Newton es conveniente para calcular el polinomio, en Pn, que interpola a una<br />

función f en x0, . . . , xn−1, xn una vez conocido el polinomio interpolador de f en x0, . . . , xn−1.<br />

La forma de Newton del polinomio interpolador puede verse como una generalización del polinomio<br />

de Taylor asociado a una función. En esta construcción aparecen las diferencias divididas<br />

que presentamos a continuación.<br />

Primera diferencia dividida<br />

Segunda diferencia dividida<br />

f[x0, x1] = f(x1) − f(x0)<br />

.<br />

x1 − x0<br />

f[x0, x1, x2] = f[x1, x2] − f[x0, x1]<br />

.<br />

x2 − x0<br />

Así sucesivamente se define la diferencia de orden k asociada a los puntos x0, . . . , xk,

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

Saved successfully!

Ooh no, something went wrong!