Descargue como PDF, TXT o lea en línea desde Scribd
Descargar como pdf o txt
Está en la página 1/ 3
INTERPOLACIÓN POLINÓMICA.
- FORMA DE LAGRANGE
Vamos a recordar, en primer lugar, cómo se define el problema de interpolación polinómica.
Supongamos que se conocen, a priori, los valores que toma una función f (que probablemente no sepamos cómo se define explı́citamente) en n + 1 puntos diferentes, x0 , x1 , . . . , xn , con n ≥ 0 (denominados nodos). Es decir, se conocen los valores f (x0 ), f (x1 ), ..., f (xn ). El problema de interpolación polinómica consiste en determinar un polinomio pn , de grado menor o igual que n, tal que
De esta forma, tendrı́amos una aproximación a la desconocida función f dada a través del polinomio pn . Esto nos permite, a continuación, dar una estimación del valor de f en un punto x̄ diferente de los nodos, es decir, dar una aproximación del valor f (x̄), mediante pn (x̄). Determinar esta aproximación a través de un polinomio es una ventaja, pues los polinomios son fáciles de evaluar. Al polinomio pn anterior se le denomina polinomio de interpolación o polinomio inter- polador. Si consideramos que la expresión general del mismo viene dada por
pn (x) = a0 + a1 x + a2 x2 + . . . + an xn , (1)
con a0 , a1 , . . . , an números reales, entonces se satisface que
a0 + a1 x0 + a2 x20 + . . . + an xn0 = f (x0 ) a0 + a1 x1 + a2 x2 + . . . + an xn = f (x1 ) 1 1 .. . . = .. a + a x + a x2 + . . . + a xn = f (x ) 0 1 n 2 n n n n
Como ya hemos estudiado en el texto base, el sistema anterior tiene solución única aunque la resolución del mismo puede ser una tarea complicada. La forma de Newton mediante diferencias divididas facilita enormemente la obtención del polinomio interpolador. Este método se estudia en el texto base. A continuación, se considera otro método, la forma de Lagrange, que permite también determinar el polinomio pn de una forma menos laboriosa y, por tanto, con un menor coste computacional. En primer lugar, hay que introducir los polinomios li (x), i = 0, 1, . . . , n, definidos del siguiente modo: (x − x0 )(x − x1 ) · · · (x − xi−1 )(x − xi+1 ) · · · (x − xn ) li (x) = . (xi − x0 )(xi − x1 ) · · · (xi − xi−1 )(xi − xi+1 ) · · · (xi − xn )
Observamos que en la definición de li aparecen todos los factores (x − xj ) menos el corres-
pondiente al nodo xi . Por tanto, en particular, el denominador no se anula nunca (recordamos que los nodos además son todos distintos). Hay que notar también que el denominador es un número, y que li es un polinomio de grado n. Por otro lado, es fácil comprobar que li (xi ) = 1, para todo i = 0, 1, . . . , n, y li evaluado en cualquier nodo que no sea xi nos da el valor cero. En la forma de Lagrange, el polinomio interpolador de grado menor o igual que n, pn , se expresa del siguiente modo: pn (x) = b0 l0 (x) + b1 l1 (x) + · · · + bn ln (x),
1 con b0 , b1 , . . . , bn números reales a determinar. Por tanto, en este método, pn se expresa como una combinación lineal de los polinomios li (x), i = 0, 1, . . . , n, en vez de expresarse en potencias de x, como en (1). Ahora bien, con esta expresión, observamos que
se deduce que b0 = f (x0 ), b1 = f (x1 ), . . . , bn = f (xn ), por lo que pn (x) = f (x0 )l0 (x) + f (x1 )l1 (x) + . . . + f (xn )ln (x).
Vamos a ver un ejemplo para comprender mejor las notaciones y la aplicación del método. ————————————————————————————————————————— Ejemplo. Consideremos los datos de los Ejemplos 3.17 y 3.18 del texto base. En estos ejem- plos, se conoce el valor de la temperatura de una habitación en cuatro instantes de tiempo representados mediante los nodos x0 = 0, x1 = 1, x2 = 2 y x3 = 4. La temperatura en esos instantes viene dada, respectivamente, por los valores f (x0 ) = 15, f (x1 ) = 18, f (x2 ) = 20 y f (x3 ) = 23. Vamos a determinar el polinomio interpolador p3 para los datos anteriores mediante la forma de Lagrange. Para ello, construimos, en primer lugar, los polinomios l0 (x), l1 (x), l2 (x) y l3 (x). Según la definición, se tiene que
2 Utilizando, por ejemplo, el comando expand de Maxima, vemos que la expresión anterior coin- cide con la del polinomio obtenido en el Ejemplo 3.18 del texto base.
