PRÁCTICA 3: INTERPOLACIÓN
En esta sesión interpolaremos usando polinomios y también segmentaria-
mente (con splines).
1) Polinomios en Octave/MatLab: Recuerda que para representar, por
ejemplo, la función y = x2 , se usa un array x=linspace(-1,1,100) para las
abcisas y otro array y=x.^2 para las ordenadas, y luego plot(x,y). Los
polinomios se representan en MatLab/Octave como vectores de coeficientes.
Se multiplican con la convolución conv, y se evaluan con polyval (y las raı́ces
con roots).
Objetivo 1: Construye el polinomio f (x) = (x − 2)9 a partir de x − 2 usando
conv. Haz la gráfica de f (x) en un entorno de 2. ¿Qué se aprecia si el entorno
es pequeño?
Objetivo 2: Recuerda que el polinomio de Taylor de la función f (x) en el
punto a hasta el grado k es:
0 f 00 (a) 2 f (k) (a)
f (x) = f (a) + f (a)(x − a) + (x − a) + · · · + (x − a)k
2! k!
Representa la función sen x junto con sus polinomios de Taylor en 0 hasta
grado 7. Lo mismo con cos x. Representa la función ex junto con sus polino-
mios de Taylor en 1 hasta grado 5.
2) Polinomios de interpolación
Objetivo 1: Desarrolla una función que calcule el polinomio de interpola-
ción, ya sea por Lagrange o por Newton. Testarlo con la función polyfit.
Objetivo 2: Comparar las instrucciones anteriores usando la función f (x) =
x sen(x2 ) en el intervalo [−π, π] tomando 8 puntos equiespaciados.
1
2 PRÁCTICA INTERPOLACIÓN
Objetivo 3: Calcular el polinomio de interpolación de la función f (x) =
−x2
e utilizando un soporte de 9 puntos equiespaciados en el intervalo [−2, 2].
Representar la función y el polinomio de interpolación.
3) Interpolación segmentaria: Para evitar los fenómenos que surgen
al aumentar el grado del polinomio, se usa la interpolación por segmentos
spline.
Objetivo 1: Representar los splines de interpolación de las funciones ante-
riores.
4) Interpolación de curvas: En el tema de interpolación polinómica se
ha visto como interpolar los valores de una tabla:
xi x0 x1 · · · xn
yi y0 y1 · · · yn
produciendo una función polinómica y = p(x) que pasa por esos puntos.
Pero este esquema solo vale para funciones, es decir, si los valores xi son
distintos entre sı́. No se puede hallar la función que interpola los puntos
A = (0, 5), B = (5, 5), C = (5, 0), D = (0, 0), o sea la tabla
PRÁCTICA INTERPOLACIÓN 3
x 0 5 5 0
y 5 5 0 0
porque no hay una función que en el punto x = 5 valga a la vez y = 0 (por
C) y y = 5 (por B).
Pero sı́ se puede definir una curva parametrizada, de parámetro t, cuyas
funciones x = x(t); y = y(t) sean cada una un polinomio de interpolación, de
t 1 2 3 4 t 1 2 3 4
la tablas y
x 0 5 5 0 y 5 5 0 0
Objetivo 1: Hallar los polinomios de interpolación x(t), y(t) de modo que la curva
resultante pase por los puntos A, B, C y D.
Objetivo 2: Representar esa curva con Octave/Matlab.
Objetivo 3: ¿Qué sucede si los cuatro puntos se escogen en otro orden? Por ejemplo
en zig-zag.
[ Pedro J. Pascual Broncano, Dept. de Ingenierı́a Informática, U. Autóno-
ma de Madrid
[Link]
es/~pedro/graficos/teoria/CurvasIntro/[Link] ]
4 PRÁCTICA INTERPOLACIÓN
Referencias
1. Curso de Cálculo Numérico, Miguel Pasadas Fernández, Departamento
de Matemática Aplicada, Universidad de Granada, 2008-09.
2. Prácticas de Métodos Matemáticos, Damián Ginestar Peiró, Departa-
mento de Matemática Aplicada, Universitat Politècnica de València,
2018-2019.
