Interpolación polinómica de Lagrange

Interpolación polinómica de Lagrange

En análisis numérico, el polinomio de Lagrange, llamado así en honor a Joseph-Louis de Lagrange, es el polinomio que interpola un conjunto de puntos dado en la forma de Lagrange. Fue descubierto por Edward Waring en 1779 y redescubierto más tarde por Leonhard Euler en 1783.

Dado que existe un único polinomio interpolador para un determinado conjunto de puntos, resulta algo confuso llamar a este polinomio el polinomio interpolador de Lagrange. Un nombre más conciso es interpolación polinómica en la forma de Lagrange.

En esta imagen se muestran, para cuatro puntos ((−9, 5), (−4, 2), (−1, −2), (7, 9)), la interpolation polinómica (cúbica) L(x), que es la suma de la bases polinómicas escaladas y0l0(x), y1l1(x), y2l2(x) y y3l3(x). La interpolación polinómica pasa exactamente por los cuatro puntos (llamados puntos de control) y cada base polinómica escalada pasa por su respectivo punto de control y se anula cuando x corresponde a los otros puntos de control.

Contenido

Definición

Dado un conjunto de k + 1 puntos

(x_0, y_0),\ldots,(x_k, y_k)

donde todos los xj se asumen distintos, el polinomio interpolador en la forma de Lagrange es la combinación lineal

L(x) = \sum_{j=0}^{k} y_j \ell_j(x)

de bases polinómicas de Lagrange

\ell_j(x) = \prod_{i=0,\, i\neq j}^{k} \frac{x-x_i}{x_j-x_i} = \frac{x-x_0}{x_j-x_0}\cdots \frac{x-x_{j-1}}{x_j-x_{j-1}}\frac{x-x_{j+1}}{x_j-x_{j+1}}\cdots \frac{x-x_{k}}{x_j-x_{k}}

Demostración

La función que estamos buscando es una función polinómica L(x) de grado k con


El problema de interpolación puede tener tan solo una solución, pues la diferencia entre dos tales soluciones, sería otro polinomio de grado k a lo sumo, con k+1 ceros.

Por lo tanto, L(x) es el único polinomio interpolador.

Concepto

La resolución de un problema de interpolación lleva a un problema de álgebra lineal en el cual se debe resolver un sistema de ecuaciones. Usando una base monómica estándar para nuestro polinomio interpolador, llegamos a la matriz de Vandermonde. Eligiendo una base distinta, la base de Lagrange, llegamos a la forma más simple de matriz identidad = δi,j, que puede resolverse inmediatamente.

Uso

Ejemplo

La función tangente y su interpolador.

Se desea interpolar f(x) = tan(x) en los puntos

x0 = − 1.5 f(x0) = − 14.1014
x1 = − 0.75 f(x1) = − 0.931596
x2 = 0 f(x2) = 0
x3 = 0.75 f(x3) = 0.931596
x4 = 1.5 f(x4) = 14.1014

Con cinco puntos, el polinomio interpolador tendrá, como máximo, grado cuatro (es decir, la máxima potencia será cuatro), al igual que cada componente de la base polinómica.

La base polinómica es:

\ell_0(x)={x - x_1 \over x_0 - x_1}\cdot{x - x_2 \over x_0 - x_2}\cdot{x - x_3 \over x_0 - x_3}\cdot{x - x_4 \over x_0 - x_4}
             ={1\over 243} x (2x-3)(4x-3)(4x+3)
\ell_1(x)={x - x_0 \over x_1 - x_0}\cdot{x - x_2 \over x_1 - x_2}\cdot{x - x_3 \over x_1 - x_3}\cdot{x - x_4 \over x_1 - x_4}
             =-{8\over 243} x (2x-3)(2x+3)(4x-3)
\ell_2(x)={x - x_0 \over x_2 - x_0}\cdot{x - x_1 \over x_2 - x_1}\cdot{x - x_3 \over x_2 - x_3}\cdot{x - x_4 \over x_2 - x_4}
             ={1\over 243} (243-540x^2+192x^4)
\ell_3(x)={x - x_0 \over x_3 - x_0}\cdot{x - x_1 \over x_3 - x_1}\cdot{x - x_2 \over x_3 - x_2}\cdot{x - x_4 \over x_3 - x_4}
             =-{8\over 243} x (2x-3)(2x+3)(4x+3)
\ell_4(x)={x - x_0 \over x_4 - x_0}\cdot{x - x_1 \over x_4 - x_1}\cdot{x - x_2 \over x_4 - x_2}\cdot{x - x_3 \over x_4 - x_3}
             ={1\over 243} x (2x+3)(4x-3)(4x+3)

Así, el polinomio interpolador se obtiene simplemente como la combinación lineal entre los \ell_i(x) y los valores de las abscisas:

{1\over 243}\Big(f(x_0)x (2x-3)(4x-3)(4x+3)-8f(x_1)x (2x-3)(2x+3)(4x-3)
+f(x_2)(243-540x^2+192x^4)-8f(x_3)x (2x-3)(2x+3)(4x+3) \,
+f(x_4)x (2x+3)(4x-3)(4x+3)\Big)\,
=-1.47748x+4.83456x^3.\,

Desventajas de su uso

No siempre funciona correctamente con cantidades mayores de seis puntos. A medida que crece el grado del polinomio interpolador, se perciba una creciente variación entre puntos de control consecutivos, lo que produce que la aproximación entre dos puntos continuos sea muy distinta a la que se esperaría. Es complicado para cálculos manuales.

Otras aplicaciones

Aunque el polinomio interpolador de Lagrange se emplea mayormente para interpolar funciones e implementar esto fácilmente en una computadora, también tiene otras aplicaciones en el campo del álgebra exacta, lo que ha hecho más célebre a este polinomio, por ejemplo en el campo de los proyectores ortogonales:

Sea un espacio vectorial complejo de dimensión finita E en el que definimos un producto escalar (no necesariamente el usual). Sea F un operador normal, tal que gracias al teorema de la descomposición espectral es igual a \sum_{i=1}^{n} \lambda_{i}P_{i}. Donde Pi son los proyectores ortogonales y λi los autovectores de F asociados a cada proyector. Entonces:

P_i=\prod_{i\neq j}\frac{F-\lambda_jI}{\lambda_i-\lambda_j}=l_i(F)

Siendo I la matriz identidad.


Demostración:

Haciendo uso de la descomponsición espectral y aplicando las propiedades de los proyectores:

l_i(F)=l_i( \sum_{j=1}^{n} \lambda_jP_j )= \sum_{j=1}^{n} l_i( \lambda_j )P_j = \sum_{j=1}^{n} \delta_{ij}P_j = P_i


Véase también

Enlaces externos


Wikimedia foundation. 2010.

См. также в других словарях:

  • Interpolación polinómica — En análisis numérico, la interpolación polinómica es una técnica de interpolación de un conjunto de datos o de una función por un polinomio. Es decir, dado cierto número de puntos obtenidos por muestreo o a partir de un experimento se pretende… …   Wikipedia Español

  • Lagrange — Saltar a navegación, búsqueda Lagrange puede referirse a: Joseph Louis de Lagrange, matemático y físico italiano; Joseph Lagrange, militar francés de la época de las Guerras Napoleónicas; Marie Joseph Lagrange, fundador y primer director de la… …   Wikipedia Español

  • Lagrange (desambiguación) — Lagrange puede referirse a: Contenido 1 Apellido 1.1 Personajes 2 Matemática y física 3 Divisiones administrativas …   Wikipedia Español

  • Interpolación trigonométrica — En matemática, la interpolación trigonométrica es una interpolación con polinomios trigonométricos. La interpolación es el método por el cual se encuentra una función a partir de un conjunto de puntos. Para la interpolación trigonométrica, esta… …   Wikipedia Español

  • Esquema de Shamir — Saltar a navegación, búsqueda Adi Shamir, desarrollador del sistema de compartición de secretos que lleva su nombre. El sistema de compartición de secretos de Shamir es un algoritmo criptográfico. Es una forma de compartición de secretos donde un …   Wikipedia Español

  • Edward Waring — Saltar a navegación, búsqueda Edward Waring Nacimiento 1734 Old Heath, Shropshire, Inglaterra …   Wikipedia Español

  • Regla de Simpson — La función f (x) (azul) es aproximada por una función cuadrática P (x) (rojo). En análisis numérico, la regla o método de Simpson (nombrada así en honor de Thomas Simpson) y a veces llamada regla de Kepler es un método de integración numérica que …   Wikipedia Español

  • Matriz de Vandermonde — es, en álgebra lineal, una matriz que presenta una progresión geométrica en cada fila. Esta matriz recibe dicho nombre en honor al matemático francés Alexandre Théophile Vandermonde. Los índices de la matriz de tamaño n×n están descritos por para …   Wikipedia Español


Поделиться ссылкой на выделенное

Прямая ссылка:
Нажмите правой клавишей мыши и выберите «Копировать ссылку»