Interpolación Hermite Python (Jupyter)
Publicado por Iván (5 intervenciones) el 05/10/2018 18:52:54
Buenas tardes,
En un laboratorio de clase nos han propuesto modificar la función de interpolación de Lagrange, para que calcule el polinomio interpolador de Hermite, llevo un rato intentándolo a partir de mi función de Lagrange, pero no consigo sacarlo. Dejo el enunciado del ejercicio así como la función de diferencias divididas de Lagrange. Muchas gracias!
"Modificar la función dd_interp para que calcule los coeficientes del polinomio interpolador de Hermite (en la forma de Newton).
Utilizar dicha función para conseguir el polinomio interpolador de Hermite de grado 7 para la función cos(x)
con nodos igualmente espaciados en el intervalo [0,3]. Evalúa este polinomio en los puntos 0.2i
con i=0,…,15
. Dibuja la gráfica del polinomio interpolador obtenido y la de la función cos(x)
¿Qué se observa?"
En un laboratorio de clase nos han propuesto modificar la función de interpolación de Lagrange, para que calcule el polinomio interpolador de Hermite, llevo un rato intentándolo a partir de mi función de Lagrange, pero no consigo sacarlo. Dejo el enunciado del ejercicio así como la función de diferencias divididas de Lagrange. Muchas gracias!
"Modificar la función dd_interp para que calcule los coeficientes del polinomio interpolador de Hermite (en la forma de Newton).
Utilizar dicha función para conseguir el polinomio interpolador de Hermite de grado 7 para la función cos(x)
con nodos igualmente espaciados en el intervalo [0,3]. Evalúa este polinomio en los puntos 0.2i
con i=0,…,15
. Dibuja la gráfica del polinomio interpolador obtenido y la de la función cos(x)
¿Qué se observa?"
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
def dd_interp(x,y):
import numpy as np
m=x.size
t=np.zeros((m,m))
t[:,0]=y
for i in range (1,m):
for j in range(1,i+1):
t[i,j]=(t[i,j-1]-t[i-1,j-1])/(x[i]-x[i-j])
return t.diagonal()
Valora esta pregunta
0