INTERPOLACIÓN POLINOMIAL.
30/10/2023
—
Análisis Numérico
—
Álvaro Espinoza Pérez
1. Dada la función halla su polinomio de Taylor de orden en torno a
y aproxima . Encuentre un límite superior para el error
utilizando la formula de error y comparela con el error real. Grafique y
SOLUCIÓN
Debemos hallar las derivadas de la función
Se evaluan en
Se aplica el teorema:
INTERPOLACIÓN POLINOMIAL PÁGINA 2
Acotamos el error
Comparación entre el polinomio interpolado y la curva original
Evaluamos en 0.5
El error es similar, es una buena aproximación
INTERPOLACIÓN POLINOMIAL PÁGINA 3
2. Compara tu respuesta con el método anterior y obtén alguna Considere la función de Bessel
Tenemos la siguiente información
a. Obtener la forma de Lagrange del polinomio interpolante
b. Interpolar J0(0.25). Calcule el error.
Se calcula el polinomio de lagrange como se muestra a continuación
Interpolamos el punto si gráficamos a ya tenemos
INTERPOLACIÓN POLINOMIAL PÁGINA 4
El error para esta aproximación viene dado por
Las derivadas de la función son:
El valor para ya que en este punto la tercera derivada en máxima
INTERPOLACIÓN POLINOMIAL PÁGINA 5
3. Verifique directamente que si 𝑥1, 𝑥2, 𝑥3 son puntos distintos, entonces:
𝑓 [𝑥1, 𝑥2, 𝑥3] = 𝑓 [𝑥3, 𝑥2, 𝑥1] = 𝑓 [𝑥2, 𝑥3, 𝑥1]
Calculamos
Si calculamos esta será:
Probamos que:
Las diferencias y son iguales al estar los puntos equiespaciados
Se ha llegado a una tautología, por tanto si se cumple que: , de
manera que ahora probamos que se cumple así:
Esta igualdad no puede reducirse ya que las diferencias y no son iguales en
magnitud, por tanto no es cierto que
INTERPOLACIÓN POLINOMIAL PÁGINA 6
4. En la siguiente tabla:
donde 𝑖 es la corriente y 𝑣 el voltaje consumido por un arco magnético. Aproxime el valor de 𝑣
𝑖 = 3.5 por un polinomio de Newton en diferencias divididas y calcule el error. Compare con el valor
para
dado por la fórmula empírica:
SOLUCIÓN
Se calculan las diferencias divididas con la ayuda de OCTAVE
Por inducción se puede establecer que, en general, para un polinomio de grado 𝑛, el polinomio
El método de Newton de diferencias divididas es otra forma de obtener el polinomio interpolador.
interpolador de Newton se escribe de la forma:
Los coeficientes son:
Calculamos en
INTERPOLACIÓN POLINOMIAL PÁGINA 7
Calculamos el error
Se gráfica el polinomio
INTERPOLACIÓN POLINOMIAL PÁGINA 8
5. De la siguiente tabla de valores de 𝑙𝑛 (𝑥), interpole para obtener 𝑙𝑛 (2.352) y 𝑙𝑛 (2.387) usando
la forma de diferencia hacia adelante de Newton del polinomio de interpolación:
Usando los valores correctamente redondeados 𝑙𝑛 2.352 ≈ 0.85527 y 𝑙𝑛 2.387 ≈ 0.87004, muestre
que la fórmula de diferencia hacia adelante es más exacta cerca de la parte superior de la tabla que
lo que está cerca de la parte inferior.
SOLUCIÓN
Se calcula el polinomio de newton con la diagonal superior
Se calcula el polinomio de newton con la diagonal inferior
INTERPOLACIÓN POLINOMIAL PÁGINA 9
Al observar el error, los errores mas bajos se cometieron con la diagonal superior, por lo tanto el
polinomio fue mas exacto
INTERPOLACIÓN POLINOMIAL PÁGINA 10
6. INTERPOLACIÓN DE HERMITE:
Usando el polinomio de Hermite que interpola 𝑓(𝑥)
en los puntos 𝑥0 = −2, 𝑥1 = −1, 𝑥2 = 0 𝑦 𝑥3 = 1. Halle el valor de 𝑓(0,5). Calcule el error.
a. Consideremos la función
La derivada de la función es
Por tanto se puede construir la siguiente tabla de datos
Se tiene la igualdad y se procede a calcular las diferencias siguiendo el mismo
procedimiento del punto anterior:
Construimos el polinomio de Hermite con los valores de la fila principal
INTERPOLACIÓN POLINOMIAL PÁGINA 11
Interpolemos el valor de la función en
Como se dispone de la función original
El erros es:
INTERPOLACIÓN POLINOMIAL PÁGINA 12
b. Use el polinomio de Hermite que concuerda con los datos listados en la tabla para encontrar una
aproximación de f(1.5).
Se construye la tabla de fierencias divididas como se muestra a continuación
Usando los valores de la primera fila empezando en la columna tres se construye el polinomio de
Hermite
Al simplificar los productos
Interpolamos el valor que la función toma en
INTERPOLACIÓN POLINOMIAL PÁGINA 13
7. INTERPOLACIÓN POR SPLINES CÚBICOS:
a. Interpolar los siguientes datos de la tabla utilizando splines cúbicos
Cada subintervalo será interpolado por una función cubica de la forma:
para un
Es decir se tiene 3 sub intervalos
Para en intervalo [-1, 1]
Para en intervalo [1, 2]
Para en intervalo [2, 4]
Los valores de las constantes coinciden con , ademas se calculan los valores
Calculamos ahora las diferencias divividas
INTERPOLACIÓN POLINOMIAL PÁGINA 14
Para determinar las constantes se usa la matriz:
Los resultados son
Se calculan las constantes mediante la relación:
Esto arroja que
INTERPOLACIÓN POLINOMIAL PÁGINA 15
Los valores ara las constantes se obtienen con la formula:
Los resultados son
De manera que ya se disponen de todas las constantes, ahora se construyen los polinomios
INTERPOLACIÓN POLINOMIAL PÁGINA 16
Grafica
INTERPOLACIÓN POLINOMIAL PÁGINA 17
b. Interpolar por splines cúbicos la función en el intervalo 0 ≤ 𝑥 ≤ 1 tomando los seis
puntos de abscisas , 𝑘 = 0, 1, 2, 3, 4, 5.
Se van a interpolar los valores mostrados en la siguiente tabla
Se tienen 5 sub intervalos, por lo tanto se requieren de 5 funciones cúbicas para interpolar la curva
Para en intervalo [0, 1/5]
Para en intervalo [1/5, 2/5]
Para en intervalo [2/5, 3/5]
Para en intervalo [3/5, 4/5]
Para en intervalo [4/5, 1]
Las constantes se hallan mediante por tanto:
INTERPOLACIÓN POLINOMIAL PÁGINA 18
Los espaciamientos son iguales en cada intervalo, es decir
Calculamos ahora las diferencias divididas requeridas en la matriz
La matriz queda así:
INTERPOLACIÓN POLINOMIAL PÁGINA 19
Si despejamos mediante las constantes c tenemos el resultado:
Se tiene finalmente
Se calculan las constantes mediante la relación:
Esto arroja que
INTERPOLACIÓN POLINOMIAL PÁGINA 20
Los valores aran las constantes se obtienen con la formula:
Los resultados son
Ya que se tienen las constantes se construyen las funciones cúbicas de cada subintervalo
INTERPOLACIÓN POLINOMIAL PÁGINA 21
La representación gráfica muestra el resultado de la interpolación:
INTERPOLACIÓN POLINOMIAL PÁGINA 22
interpolador en dichos nodos y un valor aproximado de 𝑓(0.5) mediante el polinomio de Lagrange.
8. Dada la función , determina los cuatro nodos de Chebyshev en [1,5], el polinomio
Realiza los cálculos con 4 decimales
Calculamos el polinomio de Chebyshev de grado 4
Calculamos sus raices mediante:
con y
De manera que las raices son:
Se calculan los nodos en el intervalo [1, 5] con:
con y
Entonces se tiene los nodos
Determinamos el polinomio de Lagrange en los puntos , , y
-
+
-
Tras simplificar se tiene
De igual forma, usando los nodos de Chebyshev en el Polinomio de Lagrange
Tras simplificar se tiene
INTERPOLACIÓN POLINOMIAL PÁGINA 23
Si calculamos la gráfica de la función y los dos polinomios interpolantes tendremos:
Sí aproximamos el valor de mediante el polinomio de lagrange tenemos
Mediante el polinomio de Chebyshev
En comparación con el valor real de
INTERPOLACIÓN POLINOMIAL PÁGINA 24
9. APLICACIONES
a. Las torres de soporte de una línea de alta tensión tienen un vano de 400 m. Se ha medido la
altura a que se encuentra el cable en distintas partes del recorrido obteniéndose los siguientes
Obtener la forma de Newton del polinomio de interpolación y utilizarla para obtener la flecha
máxima (punto medio de la catenaria).
Se calculan las diferencias divididas con la siguiente tabla, usando OCTAVE:
Se introducen los siguientes datos:
INTERPOLACIÓN POLINOMIAL PÁGINA 25
El polinomio de Newton tiene la forma:
Los valores para los coeficientes es:
Se calcula la altura en el centro de la catenaria:
Calculamos el vaor del parametro f
INTERPOLACIÓN POLINOMIAL PÁGINA 26
b. En la tabla siguiente se presentan los alargamientos de un resorte correspondientes a fuerzas
de
diferente magnitud que lo deforman
Encuentre un polinomio interpolante que
permita calcular la longitud del resorte en
4.5 𝑘𝑔𝑓.
SOLUCIÓN
metros, si la fuerza que se le aplica es de
Calculamos las diferencias divididas para el polinomio de Newton, mediante el siguiente scrip:
El polinomio de Newton tiene la forma:
Los valores para los coeficientes es:
INTERPOLACIÓN POLINOMIAL PÁGINA 27
Rectificamos en Octave
INTERPOLACIÓN POLINOMIAL PÁGINA 28
a. Ahora representamos la función de 𝑅𝑢𝑛𝑔𝑒 𝑓(𝑥) y el polinomio interpolador 𝑃(𝑥) para
10. OCTAVE:
comparar
ambas gráficas en OCTAVE:
% 𝐹𝑢𝑛𝑐𝑖ó𝑛 𝑑𝑒 𝑅𝑢𝑛𝑔𝑒
>> 𝑓 = @(𝑥) 1./(1 + 25 ∗ 𝑥. ^2); % 𝐹𝑢𝑛𝑐𝑖ó𝑛 𝑞𝑢𝑒 𝑒𝑣𝑎𝑙ú𝑎 𝑒𝑙 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑜
𝑖𝑛𝑡𝑒𝑟𝑝𝑜𝑙𝑎𝑑𝑜𝑟 𝑑𝑒 𝑔𝑟𝑎𝑑𝑜 𝑛.
>> 𝑓𝑢𝑛𝑐𝑡𝑖𝑜𝑛 𝑣 = 𝑝𝑧(𝑧, 𝑛)
>> 𝑥 = 𝑙𝑖𝑛𝑠𝑝𝑎𝑐𝑒(−1,1, 𝑛 + 1); % 𝑛 + 1 𝑛𝑜𝑑𝑜𝑠 𝑒𝑞𝑢𝑖𝑒𝑠𝑝𝑎𝑐𝑖𝑎𝑑𝑜𝑠
>> 𝑦 = 1./(1 + 25 ∗ 𝑥. ^2);
>> 𝑝 = 𝑝𝑜𝑙𝑦𝑓𝑖𝑡(𝑥, 𝑦, 𝑛);
>> 𝑣 = 𝑝𝑜𝑙𝑦𝑣𝑎𝑙(𝑝, 𝑧);
>> 𝑒𝑛𝑑𝑓𝑢𝑛𝑐𝑡𝑖𝑜𝑛
SOLUCIÓN
La gráfica generada por el scrip, es la representación de la función y dos
polinomios interpolante de 8 y 14 nodos igualmente espaciados para evitar el fenomeno de RUNGE
función 𝑓.
En general, un polinomio interpolante se podría ver afectado por el conjunto (𝑥1,...,𝑥𝑛) y por la
INTERPOLACIÓN POLINOMIAL PÁGINA 29
de vapor sobrecalentado a una presión de 0.1 𝑀𝑃a
b. En ingeniería se trabaja frecuentemente con datos tabulados. Se va considerar la siguiente tabla
Aplicar, con OCTAVE, la función de interpolación lineal para determinar la energía interna a 225 [C].
De igual modo, determinar la temperatura si la energía interna es 2735 [kJ/kg] usando el siguiente
>> 𝑇 = [100, 160, 200, 240, 300, 400, 500];
script:
>> 𝑢 = [2506.7, 2597.8, 2658.1, 2718.5, 2810.4, 2967.9, 3131.6];
>> 𝑇_𝑖𝑛𝑡𝑒𝑟 = 225;
>> 𝑢_𝑖𝑛𝑡𝑒𝑟 = 2735;
>> 𝐸𝑛𝑒𝑟_𝑖𝑛𝑡𝑒𝑟 = 𝑖𝑛𝑡𝑒𝑟𝑝1(𝑇, 𝑢, 𝑇_𝑖𝑛𝑡𝑒𝑟)
>> 𝑇𝑒𝑚𝑝_𝑖𝑛𝑡𝑒𝑟 = 𝑖𝑛𝑡𝑒𝑟𝑝1(𝑢, 𝑇, 𝑢_𝑖𝑛𝑡𝑒𝑟)
¿Cuál es la solución?
La solución es:
La energia interna a 225°C es. y la temperatura interna con una energia de
es 250.77°C
INTERPOLACIÓN POLINOMIAL PÁGINA 30
>> 𝑥 = [1 3 5 7 9 11];
c. Use esta rutina en la consola de OCTAVE y comente los resultados:
>> 𝑦 = [1 2 1.5, 1 3 5];
>> 𝑥𝑝 = 1: 0.25: 11;
>> 𝑖𝑛𝑡 = 𝑖𝑛𝑡𝑒𝑟𝑝1(𝑥, 𝑦, 𝑥𝑝, ′𝑠𝑝𝑙𝑖𝑛𝑒′);
>> 𝑝𝑙𝑜𝑡(𝑥, 𝑦, ′𝑜′, 𝑥𝑝,𝑖𝑛𝑡, ′: . ′);
>> 𝑥𝑙𝑖𝑚([0 11]);
>> 𝑡𝑖𝑡𝑙𝑒(′𝐼𝑛𝑡𝑒𝑟𝑝𝑜𝑙𝑎𝑐𝑖ó𝑛 𝑚𝑒𝑑𝑖𝑎𝑛𝑡𝑒 𝑠𝑝𝑙𝑖𝑛𝑒𝑠 𝑐ú𝑏𝑖𝑐𝑜𝑠 ′);
SOLUCIÓN
La gráfica muestra una interpolación mediante splines, de los puntos x=[1 3 5 7 9 11], y=[1 2 1.5 1 3
5]; con resultados correctos, ya que la curva resultante intersecta cada nodo.
INTERPOLACIÓN POLINOMIAL PÁGINA 31
función para calcular el polinomio interpolante de 𝑓(𝑥) = , en los puntos 𝑥𝑖 = 𝑖, con 𝑖 = 0, … ,6.
d. Programe en OCTAVE una función que implemente la interpolación de Lagrange. Utilice la
Grafique la función 𝑓(𝑥) y el polinomio interpolante 𝑃𝑛(𝑥) y aproxime cuando 𝑥 = 3.85 .
SOLUCIÓN
Se muestra la función creada
En la ventana de comandos se llama a la función mediante:
La gráfica generada es:
Se aproxima el valor x=3.85 con
INTERPOLACIÓN POLINOMIAL PÁGINA 32
e. Considere la siguiente tabla de datos:
de Newton. Pruebe el programa que desarrolló para calcular 𝑓(4). Grafique el interpolador juntos a
Desarrolle, depure y pruebe un programa en OCTAVE, para implantar la interpolación de polinomios
los puntos dados. .
SOLUCIÓN
Se muestra el scrip
Se muestra la ventana de comandos:
INTERPOLACIÓN POLINOMIAL PÁGINA 33
La gráfica generada es:
Al usar la función, se interpola f(4)
INTERPOLACIÓN POLINOMIAL PÁGINA 34
INTERPOLACIÓN POLINOMIAL PÁGINA 35
