MÉTODOS NUMÉRICOS PARA
INGENIEROS
Interpolación
Martha Elena Delgado Osorio, Ing. PhD.
[Link]@[Link]
INTERPOLACIÓN
Interpolar significa estimar un valor desconocido de una función en un punto, tomando una medida ponderada de
sus valores conocidos en puntos cercanos al lado.
Dados 𝑛 + 1 puntos en ℝ2 , 𝑥0 , 𝑦0 , 𝑥1 , 𝑦1 , … , 𝑥𝑛 , 𝑦𝑛 , con 𝑦𝑖 = 𝑓 𝑥𝑖 (siendo 𝑓(𝑥) no necesariamente conocida)
en los cuales 𝑥0 , 𝑥1 , … , 𝑥𝑛 son números distintos que se distribuyen en el intervalo 𝑥0 , 𝑥𝑛 , se quiere encontrar un
polinomio 𝑝𝑛 𝑥 de grado menor o igual a 𝑛 tal que:
𝑝𝑛 𝑥𝑘 = 𝑦𝑘 , 𝑘 = 0,1, … , 𝑛
Si para estimar un valor 𝑦, se emplea el polinomio 𝑝𝑛 𝑥 de grado menor o igual a 𝑛 que pasa por los puntos, la
aproximación se denomina interpolación polinomial (lineal, cuando sólo se emplean dos puntos) y a 𝑝𝑛 𝑥 se
le denomina polinomio de interpolación o polinomio interpolante.
Para un valor 𝑥 distinto de los dados:
• Si 𝑥0 < 𝑥 < 𝑥𝑛 entonces 𝑦 será un valor interpolado.
• Si 𝑥 < 𝑥0 o 𝑥𝑛 < 𝑥 entonces 𝑦 será un valor extrapolado.
INTERPOLACIÓN
INTERPOLACIÓN
INTERPOLACIÓN
INTERPOLACIÓN POLINOMIAL DE
NEWTON
INTERPOLACIÓN POLINOMIAL DE
LAGRANGE - Ejemplo
Considere la función:
𝑓 𝑥 = 2𝑥
Use el polinomio de interpolación de Lagrange que interpole 𝑓 los nodos o puntos 𝑥0 = −1, 𝑥1 = 0 𝑦 𝑥2 = 1, para
aproximar el valor de 2.
Solución:
𝑓 𝑥 = 2𝑥
𝑥−0 𝑥−1 𝑥 𝑥−1
𝑥0 = −1 𝑦0 = 1/2 𝐿0 (𝑥) = ∙ =
−1 − 0 −1 − 1 2
𝑥1 = 0 𝑦1 = 1
𝑥2 = 1 𝑦2 = 2
INTERPOLACIÓN POLINOMIAL DE
LAGRANGE - Ejemplo
Considere la función:
𝑓 𝑥 = 2𝑥
Use el polinomio de interpolación de Lagrange que interpole 𝑓 los nodos o puntos 𝑥0 = −1, 𝑥1 = 0 𝑦 𝑥2 = 1, para
aproximar el valor de 2.
Solución:
𝑓 𝑥 = 2𝑥 𝑥−0 𝑥−1 𝑥 𝑥−1
𝐿0 (𝑥) = ∙ =
−1 − 0 −1 − 1 2
𝑥0 = −1 𝑦0 = 1/2
𝑥 − −1 𝑥−1
𝑥1 = 0 𝑦1 = 1 𝐿1 (𝑥) = ∙ = 1 − 𝑥2
0 − −1 0−1
𝑥2 = 1 𝑦2 = 2
INTERPOLACIÓN POLINOMIAL DE
LAGRANGE - Ejemplo
Considere la función:
𝑓 𝑥 = 2𝑥
Use el polinomio de interpolación de Lagrange que interpole 𝑓 los nodos o puntos 𝑥0 = −1, 𝑥1 = 0 𝑦 𝑥2 = 1, para
aproximar el valor de 2.
Solución:
𝑓 𝑥 = 2𝑥 𝑥−0 𝑥−1 𝑥 𝑥−1
𝐿0 (𝑥) = ∙ =
−1 − 0 −1 − 1 2
𝑥0 = −1 𝑦0 = 1/2
𝑥 − −1 𝑥−1
𝑥1 = 0 𝑦1 = 1 𝐿1 (𝑥) = ∙ = 1 − 𝑥2
0 − −1 0−1
𝑥2 = 1 𝑦2 = 2
𝑥 − −1 𝑥−0 𝑥 𝑥+1
𝐿2 (𝑥) = ∙ =
1−0 1 − −1 2
INTERPOLACIÓN POLINOMIAL DE
LAGRANGE - Ejemplo
Considere la función:
𝑓 𝑥 = 2𝑥
Use el polinomio de interpolación de Lagrange que interpole 𝑓 los nodos o puntos 𝑥0 = −1, 𝑥1 = 0 𝑦 𝑥2 = 1, para
aproximar el valor de 2.
Solución:
𝑓 𝑥 = 2𝑥 𝑥 𝑥−1
𝐿0 𝑥 =
2
𝑥0 = −1 𝑦0 = 1/2
𝑥1 = 0 𝑦1 = 1 𝐿1 𝑥 = 1 − 𝑥 2 1 𝑥 𝑥−1 𝑥 𝑥+1
𝑃2 𝑥 = ∙ + 1 ∙ 1 − 𝑥2 + 2 ∙
𝑥2 = 1 𝑦2 = 2 2 2 2
𝑥 𝑥+1 𝒙𝟐 + 𝟑𝒙 + 𝟒
𝐿2 𝑥 = 𝑷𝟐 𝒙 =
2 𝟒
INTERPOLACIÓN POLINOMIAL DE
LAGRANGE - Ejemplo
Solución:
𝑥 2 + 3𝑥 + 4
𝑃2 𝑥 =
4
Se quiere aproximar el valor de 2 con 𝑓 𝑥 = 2𝑥 . Así:
2𝑥 = 2
1
𝑥=
2
Por lo tanto:
1 23
2 ≈ 𝑃2 = = 1.4375
2 16
INTERPOLACIÓN POLINOMIAL DE
LAGRANGE - Pseudocodigo
Función PLagrange(x0, y0, u) Paso 4: Ciclo interno para construir L_k(x)
Entradas: Para j desde 0 hasta N hacer:
x0 → Vector con las coordenadas x de los puntos conocidos Si j ≠ k entonces:
y0 → Vector con las coordenadas y = f(x) de los puntos conocidos pNum ← pNum * (x - x0[j+1])
u → Punto donde se desea evaluar el polinomio interpolante pDen ← pDen * (x0[k+1] - x0[j+1])
Salidas: Fin del ciclo interno
p → Polinomio interpolante simbólico de Lagrange
v → Valor del polinomio evaluado en u L_k(x) ← pNum / pDen // Término base de Lagrange
p(x) ← p(x) + y0[k+1] * L_k(x) // Sumar el término al polinomio total
Paso 1: Definir una variable simbólica x
Definir x como variable simbólica Paso 5: Simplificar el polinomio simbólico
Inicializar p(x) como polinomio cero Redondear resultados simbólicos a 6 cifras (opcional)
Simplificar p(x)
Paso 2: Establecer el número de puntos
N ← longitud de x0 - 1 // número de términos en la interpolación Paso 6: Evaluar el polinomio en u
v ← Evaluar p(u) y convertir a número decimal
Paso 3: Ciclo externo sobre k (cada término de Lagrange)
Para k desde 0 hasta N hacer: Paso 7: Graficar los puntos dados
Inicializar pNum ← 1 // Numerador del polinomio base L_k(x) Graficar los puntos (x0, y0) como estrellas '*'
Inicializar pDen ← 1 // Denominador del polinomio base L_k(x) (Opcional: trazar el polinomio con fplot)
Retornar p(x), v
Fin función
Referencias
Cooper, J. (2001). A MATLAB companion for multivariable calculus. Academic Press.
Chapra, S. C. (2012). Applied numerical methods. With MATLAB for Engineers and Scientists.
Chapra, S. C. (2007). Métodos Numéricos para Ingenieros.
Yang, W. Y., Cao, W., Chung, T. S., & Morris, J. (2005). Applied numerical methods using MATLAB. John Wiley &
Sons.
Nieves, A., Domínguez, F. Métodos numéricos aplicados a la ingeniería
MathWorks Documentation (2017) Primer MATLAB 2017a [Link]
[Link]
13