Interpol Ac I On

Interpolación
3 / 43
Interpolación
Supongamos que tenemos un conjunto de

n + 1 datos que relacionan a dos variables
xi , y i . x y
x0 y0
En principio no existe ninguna restricción x1 y1
sobre las diferencias Δx = xi+1 − xi . x2 y2
.. ..
. .
xi yi
.. ..
. .
xn yn
4 / 43
Interpolación lineal
Aproxima la función entre dos puntos por una línea recta:
y = mx + b
Sean dos puntos (x0 , y0 ) y (x1 , y1 ). Para encontrar m y b se imponen las

condiciones:
y0 = mx0 + b
y1 = mx1 + b
Se resuelve el sistema de ecuaciones de 2 × 2, obteniendo:

y1 − y 0
m =
x1 − x 0
y1 − y 0
b = y0 − x0
x1 − x 0
5 / 43
Interpolación lineal
La ecuación de la recta es:
y1 − y 0 y1 − y 0
y= x + y0 − x0
x1 − x 0 x1 − x 0
Y la re-escribimos separando lo que depende de y0 y lo que depende de y1 :
x1 − x x − x0
y = y0 + y1 (1)
x0 − x 1 x1 − x 0
Ejemplo: interpLineal_0.m
6 / 43
Interpolación cuadrática
Aproxima la función entre tres puntos por una parábola. Se propone:
f (x ) = Ax 2 + Bx + C
y se resuelve el sistema de 3 × 3:
f (x0 ) = Ax02 + Bx0 + C

f (x1 ) = Ax12 + Bx1 + C
f (x2 ) = Ax22 + Bx2 + C
Se deben calcular A, B y C para cada conjunto de 3 puntos a lo largo del

segmento de función estudiado.
7 / 43
Interpolación cúbica
Aproxima la función entre cuatro puntos por una cúbica. Se propone:
f (x ) = Ax 3 + Bx 2 + Cx + D
y se resuelve el sistema de 3 × 3:
f (x0 ) = Ax03 + Bx02 + Cx0 + D

f (x1 ) = Ax13 + Bx12 + Cx1 + D
f (x2 ) = Ax23 + Bx22 + Cx2 + D
f (x3 ) = Ax33 + Bx32 + Cx3 + D
Se deben calcular A, B, C y D para cada conjunto de 4 puntos a lo largo

del segmento de función estudiado.
8 / 43
Interpolación de Lagrange
Construcción del polinomio interpolador (o interpolante).

Para un polinomio de grado 1, se definen las siguientes funciones:
x1 − x x − x0
L0 (x ) = ; L1 (x ) =
x1 − x 0 x1 − x 0
Dicho polinomio es:
P(x ) = L0 (x )f (x0 ) + L1 (x )f (x1 )
Con esta definición se recupera la ecuación (1)

x − x1 x − x0
P(x ) = f (x0 ) + f (x1 )
x0 − x 1 x1 − x 0
9 / 43
Generalizando esta idea, se considera la construcción de un polinomio de

grado n que pase por los n + 1 puntos conocidos:
Para cada k = 0, 1, 2, ..., n se construye la función Ln,k (x ) con la
propiedad:

 0 si
 i �= k
Ln,k (xi ) =

 1 si i =k
Para satisfacer Ln,k (xk ) = 1, el numerador debe ser el mismo que el

denominador:
(x − x0 )(x − x1 ) · · · (x − xk−1 )(x − xk+1 ) · · · (x − xn )

Ln,k (x ) =
(xk − x0 )(xk − x1 ) · · · (xk − xk−1 )(xk − xk+1 ) · · · (xk − xn )
10 / 43
Teorema: Si x0 , x1 , ..., xn son n + 1 números distintos y si f es una
función con valores bien definidos en esos puntos, entonces ∃! polinomio
p(x ) de grado a lo más n con la propiedad p(xk ) = f (xk ) para cada
k = 0, 1, 2, ..., n.
El polinomio p(x ) esta dado por:
p(x ) = f (x0 )Ln,0 (x ) + · · · + f (xn )Ln,n (x )
es decir:
n
�
p(x ) = f (xk )Ln,k (x )
k=0
con
n
� x − xi
Ln,k (x ) =
i=0, i�=k
xk − x i
Ejemplo: interplagr_0.m
11 / 43
Polinomio interpolante
La forma del polinomio interpolante es la siguiente
Pn (x ) = a0 + a1 (x − x0 ) + a2 (x − x0 )(x − x1 ) +
a3 (x − x0 )(x − x1 )(x − x2 ) + · · · +
an (x − x0 )(x − x1 )(x − x1 ) · · · (x − xn−1 )
¿Cómo se calculan los coeficientes a0 , a1 , a2 , ..., an ?
12 / 43
Sistema de ecuaciones a resolver
y0 = a 0
y1 = a0 + a1 (x1 − x0 )
y2 = a0 + a1 (x2 − x0 ) + a2 (x2 − x0 )(x2 − x1 )
..
.
yn = a0 + aa (xn − x0 ) + a2 (xn − x0 )(x − x1 ) +
a3 (xn − x0 )(xn − x1 )(xn − x2 ) + · · · +
an (xn − x0 )(xn − x1 )(xn − x1 ) · · · (xn − xn−1 )
13 / 43
Solución
Teniendo el valor del coeficiente a0 = y0 , se calcula a1 despejándola de la

siguiente ecuación:
y1 − y 0
a1 = .
x1 − x 0
Después, se sustituyen a0 y a1 en la siguiente ecuación y se despeja a2 .
Sumando un cero de la forma 0 = y1 − y1 y después de algunas
manipulaciones algebraicas, se tiene
y2 −y1 y1 −y0
x2 −x1 − x1 −x0
a2 =
x2 − x 0
Y así sucesivamente hasta calcular an .
14 / 43
Diferencias divididas
Definición
Diferencias divididas de orden 0
dd(xi , 0) = y0

dd(xi+1 , 0) − dd(xi , 0)
dd(xi , 1) =
xi+1 − xi
dd(xi+1 , 1) − dd(xi , 1)
dd(xi , 2) =
xi+2 − xi
dd(xi+1 , 2) − dd(xi , 2)
dd(xi , 3) =
xi+3 − xi
etc.
15 / 43
Coeficientes
Los coeficientes del polinomio interpolador de Newton son:
a0 = dd(x0 , 0)
a1 = dd(x0 , 1)
a2 = dd(x0 , 2)
..
.
an = dd(x0 , n)
dd(xi , 0) dd(xi , 1) dd(xi , 2) dd(xi , 3) dd(xi , 4)

y1 = a 0 a1 a2 a3 a4
y2 ... ... ... 0
y3 ... ... 0 0
y4 ... 0 0 0
y5 0 0 0 0
16 / 43
Ejemplo
Se tienen los siguientes datos sobre 6 propiedades del agua a diferentes
temperaturas:
Figura: Tomada de http://materias.fi.uba.ar/6731/Tablas/Tabla9.pdf
ρ: densidad,
Cp : calor específico,
k: conductividad térmica,
κ: coeficiente de difusión térmica,
µ: viscosidad dinámica,
ν: viscosidad cinemática.
Ejemplo: interpNewton.m
17 / 43
Incrementos constantes
Si xi+1 − xi = Δx =cte., el polinomio interpolante de Newton se puede

expresar como: � �
n
� s
pn (x ) = f [x0 ] + Δk f (x0 )
k=1
k
usando diferencias progresivas, y
n
� �
� −s
k
pn (x ) = f [xn ] + (−1) ∇k f (xn )
k=1
k
donde � �
−s s!
= (−1)k
k k!(s − k)!
18 / 43
Diagrama de rombos
19 / 43
Interpolación bilineal
Si
a = f (x0 , y0 ),
b = f (x1 , y0 ),
c = f (x0 , y1 ),
d = f (x1 , y1 ),
entonces
(x1 − x )(y1 − y ) (x − x0 )(y1 − y )

f (x , y ) = a+ b
(x1 − x0 )(y1 − y0 ) (x1 − x0 )(y1 − y0 )
(x1 − x )(y − y0 ) (x − x0 )(y − y0 )
+ c+ d
(x1 − x0 )(y1 − y0 ) (x1 − x0 )(y1 − y0 )
20 / 43
Ejemplo
� � � �
2 x 2 y
f (x , y ) = sen + cos
7 4
21 / 43

Interpol Ac I On

Cargado por

Copyright:

Formatos disponibles

Interpol Ac I On

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Interpol Ac I On

Cargado por

Copyright:

Formatos disponibles

Interpolación

Supongamos que tenemos un conjunto de

Aproxima la función entre dos puntos por una línea recta:

Sean dos puntos (x0 , y0 ) y (x1 , y1 ). Para encontrar m y b se imponen las

Se resuelve el sistema de ecuaciones de 2 × 2, obteniendo:

Aproxima la función entre tres puntos por una parábola. Se propone:

f (x0 ) = Ax02 + Bx0 + C

Se deben calcular A, B y C para cada conjunto de 3 puntos a lo largo del

Aproxima la función entre cuatro puntos por una cúbica. Se propone:

f (x0 ) = Ax03 + Bx02 + Cx0 + D

Se deben calcular A, B, C y D para cada conjunto de 4 puntos a lo largo

Construcción del polinomio interpolador (o interpolante).

P(x ) = L0 (x )f (x0 ) + L1 (x )f (x1 )

Con esta deﬁnición se recupera la ecuación (1)

Generalizando esta idea, se considera la construcción de un polinomio de

Para satisfacer Ln,k (xk ) = 1, el numerador debe ser el mismo que el

(x − x0 )(x − x1 ) · · · (x − xk−1 )(x − xk+1 ) · · · (x − xn )

La forma del polinomio interpolante es la siguiente

¿Cómo se calculan los coeﬁcientes a0 , a1 , a2 , ..., an ?

Teniendo el valor del coeﬁciente a0 = y0 , se calcula a1 despejándola de la

Diferencias divididas de orden 1

dd(xi , 0) dd(xi , 1) dd(xi , 2) dd(xi , 3) dd(xi , 4)

Figura: Tomada de http://materias.fi.uba.ar/6731/Tablas/Tabla9.pdf

Si xi+1 − xi = Δx =cte., el polinomio interpolante de Newton se puede

(x1 − x )(y1 − y ) (x − x0 )(y1 − y )

También podría gustarte