Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Capítulo 5<br />
Interpolación<br />
El objetivo de este capítulo es estudiar cómo puede aproximarse una función por polinomios.<br />
Una forma de hacer esto es construir los polinomios de manera que coincidan con la función<br />
dada en algunos puntos predeterminados, lo que recibe el nombre de interpolación polinomial.<br />
Analizaremos distintos métodos para resolver este problema y estudiaremos el error que se<br />
comete al reemplazar una función por un polinomio interpolante.<br />
Hay diversos motivos para estudiar este problema. Por un lado, el polinomio interpolante puede<br />
utilizarse para reconstruir una función f a partir de una tabla de valores. Por otra parte, es<br />
una herramienta fundamental para integración y diferenciación numérica, como veremos más<br />
adelante.<br />
1. Interpolación de Lagrange<br />
En lo que sigue, si n ∈ IN0, llamaremos Pn al conjunto de polinomios de grado menor o igual<br />
que n, incluyendo el polinomio nulo.<br />
Supongamos que se sabe que la tabla de valores<br />
(xj): x0 x1 x2 . . . xn<br />
(yj): y0 y1 y2 . . . yn<br />
corresponde con datos de una función continua que se desconoce. Queremos poder modelizar<br />
dicha función por medio de un polinomio. Es decir, queremos encontrar un polinomio tal que<br />
p(xj) = yj, ∀j = 0, 1, . . . , n. (5.1)<br />
Nuestro primer paso será dar un resultado básico que establece que esto es posible. Mostraremos<br />
una forma concreta de hallar un polinomio p que verifique (5.1) y además veremos que si el<br />
polinomio es de grado menor o igual que n, éste es único. Vamos a basar nuestra demostración<br />
en la Base de Lagrange que es una base de polinomios que construimos a continuación.<br />
Base de Lagrange: Para cada punto xj, j = 0, . . . , n, buscamos un polinomio de grado n que<br />
se anule en todos los xi salvo xj donde queremos que valga 1. Por ejemplo, ℓ0 será un polinomio<br />
en Pn tal que se anula en x1, . . . , xn y ℓ0(x0) = 1.