M Etodo de Interpolaci On Mediante Diferencias Finitas

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

Método de Interpolación mediante Diferencias Finitas

Edgar Ramos Flores

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:

Definición 1 Se tiene las siguinetes definiciones:

1. Las primeras diferencias finitas están dadas por:

∆1 yi = yi+1 − yi i = 0, 1, 2, . . .

2. Las segundas diferencias finitas están dadas por:

∆2 yi = ∆1 yi+1 − ∆1 yi i = 0, 1, 2, . . .

3. Las terceras diferencias finitas están dadas por:

∆3 yi = ∆2 yi+1 − ∆2 yi i = 0, 1, 2, . . .

4. Las k−ésima diferencias finitas están dadas por:

∆k yi = ∆k−1 yi+1 − ∆k−1 yi i = 0, 1, 2, . . . , k = 1, 2, 3, . . .

1
Observación 1 Las diferencias finitas se pueden expresar con los puntos dados:

∆1 y 0 = y1 − y0

∆2 y 0 = ∆1 y1 − ∆1 y0 = (y2 − y1 ) − (y1 − y0 ) = y2 − 2y1 + 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

(1.0, 5), (1.5, 7), (2.0, 10), (2.5, 8)

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)

sustituyendo los puntos y despejando a0 y a1



 y = a + a (x − x )
0 0 1 0 0
 y1 = a0 + a1 (x1 − x0 )

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 )

Sustituyendo los tres puntos y despejando a0 , a1 y a2 se obtienen, incluyendo la notacion usual:

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

Definición 2 Dado el conjunto de puntos (xi , yi ), i = 0, 1, . . . , n. El polinomio de diferencias finitas o


polinomio de Newton esta dado por la siguinete formula

∆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

Polinomio de tercer grado de diferencias finitas:


∆1 y 0 ∆2 y0 ∆3 y 0
P3 (x) = y0 + (x − x0 ) + (x − x 0 )(x − x 1 ) + (x − x0 )(x − x1 )(x − x2 )
h 2!h2 3!h3
1 −4 6
= 5 + (x − 2) + (x − 2)(x − 3) + (x − 2)(x − 3)(x − 4)
1 2!(1)2 3!(1)3
= x3 − 11x2 + 37x − 33

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

4. Construir la tabla de diferencias finitas para el conjunto de puntos siguientes

x 0.125 0.250 0.375 0.500 0.625 0.750


y 0.79168 0.77334 0.74371 0.70413 0.65632 0.60228

Usar la formula de interpolacion por diferencias finitas para estimar f (0.158) y f (0.636)

También podría gustarte