M Etodo de Interpolaci On Mediante Diferencias Finitas
M Etodo de Interpolaci On Mediante Diferencias Finitas
M Etodo de Interpolaci On Mediante Diferencias Finitas
12 de mayo de 2022
1. Diferencias finitas
Supongamos que se tiene los puntos (xi , yi ), i = 0, 1, 2, . . . , n, tales que las abscisas estan espacidas
regularmente en una distancia h:
xi+1 − xi = h, i = 0, 1, . . . , n − 1
Figura 1:
∆1 yi = yi+1 − yi i = 0, 1, 2, . . .
∆2 yi = ∆1 yi+1 − ∆1 yi i = 0, 1, 2, . . .
∆3 yi = ∆2 yi+1 − ∆2 yi i = 0, 1, 2, . . .
1
Observación 1 Las diferencias finitas se pueden expresar con los puntos dados:
∆1 y 0 = y1 − y0
∆3 y 0 = ∆2 y1 − ∆2 y0
= (∆1 y2 − ∆1 y1 ) − (∆1 y1 − ∆1 y0 )
= ∆1 y2 − 2∆1 y1 + ∆1 y0
= y3 − y2 − 2(y2 − y1 ) + y1 − y0
= y3 − 3y2 + 3y1 − y0
Nota 1 Es util tabular las diferencias finitas en un cuadro como se muestra a continuación:
i xi yi ∆ 1 yi ∆2 y i ∆3 y i
0 x0 y0 ∆1 y0 ∆2 y0 ∆ 3 y0
1 x1 y1 ∆1 y1 ∆2 y1 ∆ 3 y1
2 x2 y2 ∆1 y2 ∆2 y2 ...
1
3 x3 y3 ∆ y3 ... ...
... ... ... ... ... ...
Cada diferencia finita se obtiene restando los dos valores consecutivos de la columna anterior.
Ejemplo 1 Hallar la tabla de las diferencias finitas correspondientes a los siguientes datos
Solución
Tenemos que
1
∆ y0 = y1 − y0 = 7 − 5 = 2
∆1 y1 = y2 − y1 = 10 − 7 = 3
∆1 y2 = f3 − y2 = 8 − 10 = −2
∆2 y = y − 2y + y = 1
0 2 1 0
∆2 y1 = y3 − 2y2 + y1 = −5
∆3 y0 = y3 − 3y2 + 3y1 − y0 = −6
i xi yi ∆ 1 yi ∆2 yi ∆3 y i
0 1 5 2 1 −6
1 1.5 7 3 −5
2 2.0 10 −2
3 2.5 8
2
2. Polinomio de interpolacion de diferencias finitas
Dado un conjunto de puntos (xi , yi ), i = 0, 1, 2, . . . , n espaciados en forma regular en una distancia
h y que provienen de una funcion desconocida f (x), pero supuestamente difernciable, se desea obtener
el polinomio de interpolación. Si se tienen tabuladas las diferencias finitas se puede obtener el polinomio
de interpolacion mediante un procedimiento de recurrencia y generalizacion.
Suponer que se tiene dos puntos (x0 , y0 ) y (x1 , y1 ) con los cuales se debe obtener el polinomio de
primer grado, es decir, la ecuación de la recta
y = P1 (x) = a0 + a1 (x − x0 ) (1)
a0 = y0
y1 − y0 ∆1 y 0
a1 = =
x1 − x0 h
Ahora, reemplazando en (1) y agregando la notacion usual
∆1 y0
P1 (x) = y0 + (x − x0 )
h
En el caso de tener tres puntos (x0 , y0 ), (x1 , y1 ), (x2 , y2 ) se debe obtener el polinomio de segundo
grado. Se propone la siguiente forma para este polinomio:
P2 (x) = a0 + a1 (x − x0 ) + a2 (x − x0 )(x − x1 )
a0 = y0
y1 − y0 ∆1 y 0
a1 = =
x1 − x0 h
y2 − y1 y1 − y0 ∆1 y1 ∆1 y0
− − 1 1 2
x − x1
a2 = 2
x1 − x0
= h h = ∆ y1 − ∆ y0 = ∆ y0
x2 − x0 2h 2h2 2h2
Por lo tanto,
∆1 y0 ∆2 y 0
P2 (x) = y0 + (x − x0 ) + (x − x0 )(x − x1 )
h 2h2
∆1 y0 ∆2 y 0 ∆3 y 0
Pn (x) = y0 + (x − x0 ) + (x − x 0 )(x − x 1 ) + (x − x0 )(x − x1 )(x − x2 )
h n 2!h2 3!h3
∆ y0
+··· + (x − x0 )(x − x1 ) . . . (x − xn−1 )
n!hn
Ejemplo 2 Dados los siguientes puntos (2, 5), (3, 6), (4, 3), (5, 2), encuentre el polinomio de interpolación
que incluye a los cuatro datos usando el método de diferencias finitas.
3
Solución
El metodo de diferencias finitas es aplicable pues h es constante e igual 1. Tabla de diferencias finitas:
i xi yi ∆1 y i ∆2 yi ∆3 y i
0 2 5 1 −4 6
1 3 6 −3 2
2 4 3 −1
3 5 2
Figura 2:
3. Ejercicios propuestos
1. Con la formulacion de diferencias finitas construya un polinomio interpolador que pasa por los
puntos dados en la siguiente tabla:
x 8 16 24
y 5 6 4
2. Con la formulacion de diferencias finitas construya un polinomio interpolador que pasa por los
puntos dados en la siguiente tabla:
x 5 6 7 8
y 1 3 14 15
4
3. Con la formulacion de diferencias finitas construya un polinomio interpolador que pasa por los
puntos dados en la siguiente tabla:
x 1 2 3 4 5
y 9 5 7 13 26
Usar la formula de interpolacion por diferencias finitas para estimar f (0.158) y f (0.636)