Descargue como PDF, TXT o lea en línea desde Scribd
Descargar como pdf o txt
Está en la página 1/ 10
Para resolver este tipo de problemas deberemos tener en cuenta la construc-
ción de un polinomio llamado de Lagrange, de grado n, nulo en todos sus
puntos xk salvo en uno, que será el xj , cuyo valor será 1. Éste será el polinomio general, aunque cada problema particular tendrá un grado k, de tal forma que existirán k + 1 puntos de interpolación. El intervalo necesariamente debe estar acotado de tal forma que verifique, siendo los puntos de interpolación, todos distintos:
{x0 , x1 , ..., xk } / xi = xj
∀ζ ∈ [a, b] , (a = x0 , b = xk ) ∃ P (x) ≃ f(x) / |E(x)| |Emáx |
dado que fuera de ese intervalo no existirá interpolación y el error absoluto
E(x) para esos puntos puede incrementarse por encima de su error máximo Emáx . La interpolación de Lagrange significa matemáticamente aproximar una fun- ción f(x) en un intervalo definido por los puntos donde se va a interpolar, como ya se ha visto, de tal forma que P (x) sea aproximadamente igual a la función en dicho intervalo.
1 He creado un programa informático denominado Lagrange para las cal- culadoras TI 92 Plus / Voyage 200 de Texas Instruments tras estudiar la asignatura AF.3 del curso modular MEF. Voy a realizar el cálculo paso a paso con las pantallas de este programa que es didáctico y presenta el cálculo se- cuencialmente. Vale para calcular cualquier interpolación de Lagrange para aproximar funciones en un intervalo de puntos. Realmente el intervalo en sí lo deciden los puntos, pero aún así se ha puesto de esta forma. El no de puntos puede ser de 2 a 6.
En el intervalo de validez, podemos en lugar de establecer el intervalo en-
tre corchetes, cualquier caracter/es alfanumérico/s sin los corchetes. De esta forma, el programa asumirá que el intervalo estará comprendido en el cerrado correspondiente al punto primero más bajo y al punto último más alto, tal y como está en la 1a de las pantallas de presentación de arriba. Podemos elegir entre resolver el problema paso a paso o ir directo al resultado del polinomio aproximador de Lagrange. A continuación se muestra un cálculo secuencial, paso a paso.
2 Primeramente definimos la función f(x) introducida en los datos y obte- nemos los polinomios Li (x), siendo i = 0 a 2.
Una vez hallados los polinomios de Lagrange, obtenemos el valor de la fun-
ción en cada uno de los 3 puntos dados.
3 Como se sabe el polinomio de Lagrange aproximador tiene la forma siguiente (para 3 puntos).
Su cálculo es el siguiente, observando que la solución del polinomio la da en
varios formatos, primero factorizado, segundo expandido en fracciones simples, y por último en decimales.
En el texto Cálculo Numérico AF.3 del Curso Modular "Teoría y apli-
cación práctica del método de los elementos finitos y simulación", en su página III.15, viene solo el polinomio de interpolación omitiéndose su cálculo, de la forma:
4 √ √ 2/2 1− 2 π P2 (x) = π/4 x + π2 /8 x x − 4
Haciendo operaciones vemos como la expresión anterior es la misma que la
mostrada en la calculadora en su formato expandido:
√ √ 2 2 8(1− 2) 2 π 8 2 √ √ P2 (x) = π x + π2 x − 4x = π2 x − π2 x − 8 π22 x2 + 4 2 π x
Ahora el programa Lagrange calculará el error máximo Emáx . Si quere-
mos hallar una cota máxima del error de interpolación E(x) en cada punto, tendremos:
n
f (n+1) (ξ) E(x) = (n+1)! (x − xi ) , ξ ∈ [a, b] i=0
De tal forma que hemos logrado un mayorando para E(x), con lo que lo- gramos cuantificar el valor máximo del error al aproximar la función mediante la interpolación de Lagrange:
E(x) ≤ Emáx , ∀ξ ∈ [a, b]
5 Calculemos el error Emáx , con el programa. Nuevamente lo calculamos paso a paso, aunque se podría ir directo al resultado del error.
6 El programa a continuación muestra un menú para introducir un punto de aproximación que debe estar en el intervalo de puntos, pues como se sabe, el polinomio solo interpola en dicho intervalo. Elegimos π6 que como se ve, cumple:
π
6 ∈ 0, π2
A continuación el programa muestra el error absoluto y el error relativo en
% del punto seleccionado anteriormente (en nuestro caso π6 ). Ambos errores no se dan en valores absolutos en un principio. Es por ello, que se han estimado tal y como aparecen en la siguiente pantalla, ofreciéndose primero como un valor algebraico y después con sus cifras decimales aproximadas:
7 Podríamos haber elegido un punto que cayera fuera del intervalo, por ejem- plo, π. El programa detecta que se está extrapolando a los valores del intervalo, por lo que aparece un mensaje de advertencia y no deja calcular los errores, dado que el polinomio no aproxima la función en valores fuera del intervalo, como es sabido.
Una vez obtenido los datos de errores del punto se muestra una pantalla de gráficos en 2D, con eje Y de ordenadas y eje X de abcisas. Ahora debemos establecer un dominio visual para apreciar la gráfica. Para ello, tenemos que definir las esquinas del rectángulo sobre el cual aparecerán las gráficas f(x) y P (x). Una vez establecido el rectángulo, se pulsa nuevamente ENTER en la calculadora. Vemos como podemos hacer el rectángulo como queramos, pero siempre pensando en cómo es la función f(x) para ello.
8 Entonces, aparecen las dos gráficas correspondientes, donde se aprecia fácil- mente que el polinomio efectivamente se aproxima a la función en los valores del intervalo dado.
Una vez obtenidas las gráficas podemos acceder al menú de la pantalla si- guiente, donde podemos elegir varios tipos de zoom: normal, acercado y alejado. Además si seleccionamos la opción 2:Otro punto de aproximación, podremos volver a seleccionar un rectángulo a nuestra elección.
A continuación se observa un pantallazo de las gráficas con el zoom normal.
9 Otros pantallazos con el zoom acercado.
Pongamos como importante aclaración que si introducimos un intervalo de
puntos en el primer menú que no se corresponde con los puntos introducidos, el programa calculará los errores absolutos y relativos muy elevados sobre esos puntos. Esto habrá que tenerlo en cuenta y no es un fallo del programa.
En definitiva, el programa Lagrange es una solución académica portable
en una calculadora para estudiar polinomios aproximadores de Lagrange para cualquier tipo de función de una variable desde los 2 a los 6 puntos, es decir para interpolaciones de P(x) que van de 1 a 5.
