Polinomio de interpolación de Lagrange:
Concepto general:
La interpolación de Lagrange es un método utilizado para aproximar una función
desconocida a partir de un conjunto de puntos conocidos. Es un método de interpolación
polinómica donde busca hallar un polinomio que pase justamente por los puntos dados. El
método de interpolación de Lagrange utiliza un polinomio llamado polinomio de Lagrange,
que se construye de la siguiente manera:
donde:
Primer Grado
Se busca un polinomio que sigua por cada uno de los puntos establecidos en la tabla. Si hay
n puntos el polinomio será de grado n-1 o menor.
n=1 Interpolación lineal Polinomio que interpola los datos (x0, y0), (x1, y1)
Formamos los polinomios de Lagrange
n=2 Interpolación cuadrática Polinomio que interpola los datos (x0, y0), (x1, y1),(x2, y2)
Formamos los polinomios de Lagrange:
y el polinomio p (x) =
n=3 Interpolación cubica Polinomio que interpola los datos (x0, y0), (x1, y1),(x2, y2)
y (x3, y3)
Formamos los polinomios de Lagrange
y el polinomio p (x) =
Ejemplo:
Segundo Grado
Polinomio de grado menor o igual que 2 que interpola a la función f (x) = ex en
los puntos x0 = −1, x1 = 0, x2 = 1
El polinomio es
Primer Grado:
El polinomio de Lagrange de grado 1 P(x)P(x) es:
Sustituyendo los valores:
Simplificando:
Tercer Grado:
Para la función tabular. Obtener el valor de y para x = 2.
La ecuación (1) expresada en forma extendida para un total de n = 4 es la siguiente:
Sustituyendo valores se obtiene:
Para este mismo ejemplo, si se desea obtener el polinomio que pasa por todos los puntos
que
conforman la función tabular, debe repetirse el cálculo anterior dejando al valor x como
incógnita:
Método de Euler (Ecuaciones Diferenciales):
El método más simple para resolver ecuaciones de primer orden, el Método de Euler. Si bien
no es el método más eficiente, sí nos proporciona una imagen de cómo se procede y se puede
mejorar introduciendo mejores técnicas, que normalmente se cubren en un texto de análisis
numérico.
Consideremos la clase de problemas de valor inicial de primer orden de la forma
Nos interesa encontrar la solución y(x) de esta ecuación que pasa por el punto inicial(x0,y0)
en el xy -plano para valores de x en el intervalo [a,b], donde a=x0. Buscaremos
aproximaciones de la solución en N puntos, etiquetadas xn paran=1,…,N. Para puntos
igualmente espaciados tenemos Δx=x1−x0=x2−x1, etc. Podemos escribir estos como
La idea es simple. Aproximamos la derivada en la ecuación diferencial por su cociente de
diferencia:
Dado que la pendiente de la tangente a la curva en(x0,y0) es y′(x0)=f(x0,y0), podemos
escribir:
Resolviendo esta ecuación para y1, obtenemos:
Esto da y1, En cuanto a cantidades que conocemos.
Ahora procedemos a aproximar y(x2). La línea tangente correspondiente se muestra pasando
aunque (x1,y1) y luego podemos obtener el valor de y2 desde la intersección de la línea
tangente con una línea vertical, x=x2. Siguiendo los argumentos anteriores, encontramos que:
Continuando con este procedimiento para todos xn, n=1,…N, llegamos al siguiente esquema
para determinar una solución numérica al problema del valor inicial:
Ejemplo:
Iniciamos con la condición inicial: x0=0 y y0=1.
Usamos la ecuación diferencial para encontrar la pendiente en (x0,y0):f(x0,y0)=x0+y0.
Luego, avanzamos un pequeño paso h desde x0 para encontrar y1: y1=y0+h⋅f(x0,y0).
Repetimos este proceso hasta llegar al final del intervalo.
Método de Euler Mejorado
Se basa en aproximar la curva integral de la Ecuación a(xi,y(xi)) por la línea a
través(xi,y(xi)) con pendiente.
Es decir, mi es el promedio de las pendientes de las tangentes a la curva integral en los
puntos finales de[xi,xi+1]. Por lo tanto, la ecuación de la línea de aproximación es:
Ajuste x=xi+1=xi+h
como una aproximación a y(xi+1). Al igual que en nuestra derivación del método de Euler,
reemplazamos y(xi) por su valor aproximado yi; luego la Ecuación se convierte en
No obstante, esto todavía no va a funcionar, porque desconocemosy(xi+1), que aparece a la
derecha. Superamos esto reemplazando y(xi+1) por yi+hf (xi,yi), el valor al que asignaría el
método de Euleryi+1. Así, el método mejorado de Euler inicia con el valor conocido y(x0)=y0
y calcula y1,y2,,...,yn sucesivamente con la fórmula:
El cálculo indicado aquí se puede organizar convenientemente de la siguiente manera: dado
yi, computar
Ejemplo:
Retomando el mismo ejemplo de la Ecuación Diferencial y misma condición inicial:
Para aplicar el método, necesitamos un paso h. Asumamos que h=0.1 para obtener una tabla
de iteraciones:
Donde f(x,y)=x+y.
Ahora, vamos a calcular iterativamente estos valores:
x0=0, y0=1, h=0.1
Método de Runge – Kutta
Primer orden
Aunque se desconoce la posición de la curva y = y(x), se sabe que en todo punto de esta la
pendiente es igual a f(x, y): esto es simplemente la interpretación geométrica de la ecuación
diferencial. De esta forma, la pendiente de la tangente en P es y0' = f(x0, y0); esto puede
calcularse ya que tanto x0 como y0 se conocen. Si h es razonablemente pequeña, la tangente
PT no debe desviarse demasiado de la curva PQ; de esta forma la altura MT, la cual por
geometría es igual a hy0', debe ser una aproximación a la altura buscada MQ. En otras
palabras, una primera aproximación a k está dada por:
Se hace evidente que se podría obtener una mejor aproximación si se utilizara el valor y' en
algún punto intermedio de la curva PQ en vez de el valor en un extremo, pero ya que se no
se conoce ninguno de estos puntos intermedios sobre la curva, no es posible hacer esto. Sin
embargo, podría tomarse:
en donde α y β son fracciones adecuadas como segunda aproximación a k.
Segundo Orden:
Uno de los tipos más simples del método de Runge-Kutta se obtiene al definir k1 y k2 y
aproximar a k por medio de:
El problema consiste en obtener valores adecuados para α, β, W1 y W2 de manera que la
aproximación a k sea buena en términos generales.
Y utilizando la serie de Taylor con dos variables, es posible escribir:
Por otro lado, dado que: y1 - y0 = k y a partir de la serie de Taylor, podemos escribir:
por lo que: Ahora bien:
por lo que: además:
De manera que: y entonces se tiene:
A partir de lo anterior, puede observarse que la aproximación, equivale a:
Mientras que el valor correcto de k está dado por la ecuación (6.10). Se desea un
aproximación válida para cualquier ecuación diferencial, esto es para cualquier f(x, y). Con
esto nos queda:
Se tienen cuatro parámetros desconocidos y sólo tres ecuaciones que satisfacer, de manera
que existe cierta liberta de elección. Puede tomarse a como cualquier número que se desee
(excepto cero) y entonces se tiene:
Por ejemplo, puede tomarse
a = b = 1, y se obtiene W1 = W2 = ½.
Obteniéndose lo siguiente:
Tercer Orden:
El método de Runge-Kutta de tercer orden esta dado por:
Cuarto Orden:
Ejemplos:
Primer Orden:
Ecuación:
Ejemplo:
Calcular f(t0,y0):
Actualizar y:
Segundo Orden:
Ecuaciones:
Calcular k1:
Calcular k2:
Actualizar y:
Tercer Orden:
Ecuaciones:
Calcular k1:
Calcular k2:
Calcular k3:
Actualizar y:
Cuarto Orden:
Calcular k1:
Calcular k2:
Calcular k3:
Calcular k4:
Actualizar y:
Fuentes
CHAPRA, S. C. (05 de 07 de 2023). rpubs. Obtenido de
[Link]
Cortés Rosas Jesús Javier, G. C. (2019). UNAM. Obtenido de
[Link]
Granada, U. d. (s.f.). UGR. Obtenido de [Link]
Herman, R. (30 de 10 de 2022). LibreTexts. Obtenido de University of North Carolina
Wilmington:
[Link]
so_en_Ecuaciones_Diferenciales_para_Cient%C3%ADficos_e_Ingenieros_(Herman)/
03%3A_Soluciones_num%C3%A9ricas/3.01%3A_M%C3%A9todo_de_Euler
Trench, W. F. (30 de 10 de 2022). LibreTexts. Obtenido de Trinity University:
[Link]
aciones_diferenciales_elementales_con_problemas_de_valor_l%C3%ADmite_(trinch
era)/03%3A_M%C3%A9todos_num%C3%A9ricos/3.02%3A_El_M%C3%A9todo_Euler
_Mejorado_y_M%C3%A9todos_Relacionados
Veracruzana, U. (s.f.). Universidad Veracruzana. Obtenido de
[Link]
Ward Cheney, D. K. (s.f.). Métodos Numéricos y Computación. Obtenido de Cengage
Learning: [Link]
