UNIVERSIDAD DE CONCEPCIÓN
FACULTAD DE CIENCIAS FÍSICAS Y MATEMÁTICAS
DEPARTAMENTO DE INGENIERÍA MATEMÁTICA
Guía N◦ 4: Interpolación - Parte I
Cálculo Numérico 521230, 2020-2
Nota: El comando polyfit de Octave permite encontrar los coeficientes, en la base canónica del espacio correspondiente,
del polinomio que interpola un conjunto de datos dado.
Por otro lado, el comando polyval es útil para evaluar dicho polinomio.
Para mayor información sobre ambos, escriba help polyfit o help polyval en el terminal de Octave.
1. Determine, si es posible, el polinomio que interpola a los siguientes pares ordenados. Represéntelo como combinación
lineal de los polinomios en la base canónica del espacio correspondiente.
a) (0, 1), (2, 3), (3, 0), c) (−1, 0), (2, 1), (3, 1), (5, 2),
b) (−1, 1), (0, 0), (1, 1), d ) (0, 1), (1, 2), (1, −1).
Compruebe los resultados obtenidos utilizando el comando polyfit de Octave. Grafique los puntos y el polinomio
obtenido.
2. Considere los mismos puntos del ejercicio anterior. Escriba, en cada caso en que fue posible encontrar el polinomio
de interpolación, escriba al polinomio de interpolación utilizando los polinomios de Lagrange.
3. El Cuadro 1 muestra datos de temperatura de una sala a partir de las 6:00 hrs y cada 20 minutos.
Minutos después de 6:00 hrs temperatura (0 C).
0 10
20 20
40 30
Cuadro 1. Datos para ejercicio 3
a) Encuentre el polinomio que interpola a los datos de la tabla.
b) Deduzca la temperatura de la sala a las 6:05 y 6:35 hrs.
4. Determine el polinomio que interpola a la función sin |[0,π] en los siguientes puntos:
a) x0 = 0, x1 = π/2 y x2 = π.
b) x0 = 0, x1 = π/4, x2 = π/2 y x3 = π.
c) x0 = 0, x1 = π/4, x2 = π/2 , x3 = 3π/4 y x4 = π.
Compruebe los resultados obtenidos en cada caso utilizando el comando polyfit de Octave, grafique la función,
el polinomio obtenido y los puntos de interpolación.
5. Obtenga una cota para el error de interpolación en cada uno de los casos anteriores.
6. Considere la función f = ln(x), x ∈ [1, 3].
a) Escriba, como combinación lineal de polinomios en una base de Lagrange, el polinomio p que interpola a ln en
los puntos (1, ln 1), (2, ln 2) y (3, ln 3).
b) Grafique f , p y los puntos de interpolación.
c) Utilice el polinomio calculado para aproximar los valores de ln(1,5) y ln(2,4).
d ) Determine una cota para el error de interpolación cometido en cada caso.
7. Considere los siguientes pares ordenados (−1, 5), (0, 1), (1, 1), (2, 11).
a) Muestre que los polinomios p(x) = x3 + 2x2 − 3x + 1 y q(x) = 18 x4 + 43 x3 + 15 2 11
8 x − 4 x + 1 los interpolan.
b) Explique por qué esto no contradice el teorema visto en clases sobre unicidad del polinomio de interpolación.
1
2
8. Suponga se quiere determinar un polinomio de grado menor o igual que 1 que interpole a la función f : [−1, 1] −→ R.
La manera en que se escojan los puntos x0 , x1 ∈ [−1, 1] en los que el polinomio interpola a f influyen en el error de
interpolación
(x − x0 )(x − x1 ) 00
f (x) − p(x) = f (ξx ).
2
Suponga que máxx∈[−1,1] |f 00 (x)| = 1.
a) Determine una cota para |f (x) − p(x)| en el caso en que x0 = −1√y x1 = 1. √
b) Determine una cota para |f (x) − p(x)| en el caso en que x0 = − 2/2 y x1 = 2/2.
c) Determine una cota para |f (x) − p(x)| en el caso en que x0 = − 12 y x1 = 21 .
¿En cuál de los casos obtiene la menor cota para el error de interpolación?
Observación: Dado que en este ejemplo máxx∈[−1,1] |f 00 (x)| = 1
1 1
|f (x) − p(x)| = |(x − x0 )(x − x1 )||f 00 (ξx )| ⇒ |f (x) − p(x)| ≤ máx |(x − x0 )(x − x1 )|.
2 2 x∈[−1,1]
1
Una cota para |f (x) − p(x)| con x ∈ [−1, 1] es 2 máxx∈[−1,1] |(x − x0 )(x − x1 )|.
Concepción, 3 de abril de 2021