INTERPOLACIÓN DE LAGRANGE
El método de Interpolación de Lagrange tiene el mismo fundamento que el método de Newton.
Sólo que ahora replantea las fórmulas de la siguiente manera:
𝑛𝑛
𝑓𝑓𝑛𝑛 (𝑥𝑥) = � 𝐿𝐿𝑖𝑖 (𝑥𝑥) ∙ 𝑓𝑓(𝑥𝑥𝑖𝑖 )
𝑖𝑖=0
Donde los 𝐿𝐿𝑖𝑖 (𝑥𝑥) conocidos como los langrangianos de la ecuación, se calculan mediante:
𝑛𝑛
𝑥𝑥 − 𝑥𝑥𝑗𝑗
𝐿𝐿𝑖𝑖 (𝑥𝑥) = �
𝑥𝑥𝑖𝑖 − 𝑥𝑥𝑗𝑗
𝑗𝑗=0
𝑗𝑗≠𝑖𝑖
Así, para la interpolación lineal tenemos:
𝑓𝑓1 (𝑥𝑥) = 𝐿𝐿0 (𝑥𝑥)𝑓𝑓(𝑥𝑥0 ) + 𝐿𝐿1 (𝑥𝑥)𝑓𝑓(𝑥𝑥1 )
Donde:
𝑥𝑥−𝑥𝑥1
𝐿𝐿0 (𝑥𝑥) = Lagrangiano 0 (la excepción es cuando j=0)
𝑥𝑥0 −𝑥𝑥1
𝑥𝑥−𝑥𝑥0
𝐿𝐿1 (𝑥𝑥) = Lagrangiano 1 (la excepción es cuando j=1)
𝑥𝑥1 −𝑥𝑥0
Como se puede observar, el símbolo 𝜋𝜋 de la fórmula de los lagrangianos, significa una sucesión
de productos desde j=0 hasta al grado de interpolación n, excepto cuando j=i. Ya que eso
significaría una indeterminación del término, debido al cero que se produciría en el
denominador.
Por lo tanto, el polinomio de interpolación lineal de Lagrange es:
𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥0
𝑓𝑓1 (𝑥𝑥) = 𝑓𝑓(𝑥𝑥0 ) + 𝑓𝑓(𝑥𝑥1 )
𝑥𝑥0 − 𝑥𝑥1 𝑥𝑥1 − 𝑥𝑥0
Ahora, observe como resulta el polinomio de interpolación cuadrada o grado 2 de Lagrange:
𝑓𝑓1 (𝑥𝑥) = 𝐿𝐿0 (𝑥𝑥)𝑓𝑓(𝑥𝑥0 ) + 𝐿𝐿1 (𝑥𝑥)𝑓𝑓(𝑥𝑥1 ) + 𝐿𝐿2 (𝑥𝑥)𝑓𝑓(𝑥𝑥2 )
Beatriz Adriana Jaime Fonseca
Donde:
𝑥𝑥−𝑥𝑥1 𝑥𝑥−𝑥𝑥2
𝐿𝐿0 (𝑥𝑥) = ∙ Lagrangiano 0 (la excepción es cuando j=0)
𝑥𝑥0 −𝑥𝑥1 𝑥𝑥0 −𝑥𝑥2
𝑥𝑥−𝑥𝑥0 𝑥𝑥−𝑥𝑥2
𝐿𝐿1 (𝑥𝑥) = ∙ Lagrangiano 1 (la excepción es cuando j=1)
𝑥𝑥1 −𝑥𝑥0 𝑥𝑥1 −𝑥𝑥2
𝑥𝑥−𝑥𝑥0 𝑥𝑥−𝑥𝑥1
𝐿𝐿2 (𝑥𝑥) = ∙ Lagrangiano 2 (la excepción es cuando j=2)
𝑥𝑥2 −𝑥𝑥0 𝑥𝑥2 −𝑥𝑥1
Por lo tanto, el polinomio de interpolación cuadrática o grado 2 de Lagrange es:
𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥1
𝑓𝑓2 (𝑥𝑥) = ∙ 𝑓𝑓(𝑥𝑥0 ) + ∙ 𝑓𝑓(𝑥𝑥1 ) + ∙ 𝑓𝑓(𝑥𝑥2 )
𝑥𝑥0 − 𝑥𝑥1 𝑥𝑥0 − 𝑥𝑥2 𝑥𝑥1 − 𝑥𝑥0 𝑥𝑥1 − 𝑥𝑥2 𝑥𝑥2 − 𝑥𝑥0 𝑥𝑥2 − 𝑥𝑥1
Observe que los lagrangianos de primer orden o de interpolación lineal, tienen un solo término,
mientras que los de interpolación cuadrática, tienen dos. Por lo que es esperado que los
lagrangianos de la interpolación cúbica tengan tres términos y así sucesivamente. Otra cosa que
podemos observar es que los términos que se restan en el numerador y denominador, son
iguales; es decir:
En el numerador siempre será 𝑥𝑥 al que se le va a restar y en el denominador se tendrá la
𝑥𝑥𝑠𝑠𝑠𝑠𝑠𝑠í𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛 𝑑𝑑𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 menos algún término. Esto es:
Ya sólo deben de cuidar las excepciones para que no se indeterminen los términos de los
lagrangianos.
Beatriz Adriana Jaime Fonseca
Conforme a lo anterior, los lagrangianos del polinomio de interpolación cúbica resultan:
𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥3
𝐿𝐿0 = ∙ ∙
𝑥𝑥0 − 𝑥𝑥1 𝑥𝑥0 − 𝑥𝑥2 𝑥𝑥0 − 𝑥𝑥3
𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥3
𝐿𝐿1 = ∙ ∙
𝑥𝑥1 − 𝑥𝑥0 𝑥𝑥1 − 𝑥𝑥2 𝑥𝑥1 − 𝑥𝑥3
𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥3
𝐿𝐿2 = ∙ ∙
𝑥𝑥2 − 𝑥𝑥0 𝑥𝑥2 − 𝑥𝑥1 𝑥𝑥2 − 𝑥𝑥3
𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥2
𝐿𝐿3 = ∙ ∙
𝑥𝑥3 − 𝑥𝑥0 𝑥𝑥3 − 𝑥𝑥1 𝑥𝑥3 − 𝑥𝑥2
Y el polinomio se calculará con:
𝑓𝑓3 (𝑥𝑥) = 𝐿𝐿0 𝑓𝑓(𝑥𝑥0 ) + 𝐿𝐿1 𝑓𝑓(𝑥𝑥1 ) + 𝐿𝐿2 𝑓𝑓(𝑥𝑥2 ) + 𝐿𝐿3 𝑓𝑓(𝑥𝑥3 )
Ejemplo:
a) Estime el ln 2 empleando la interpolación lineal de Lagrange, considerando los siguientes datos:
𝑥𝑥0 = 1 𝑓𝑓(𝑥𝑥0 ) = 0
𝑥𝑥1 = 6 𝑓𝑓(𝑥𝑥1 ) = 1.791759
Aplicando la siguiente fórmula con 𝑥𝑥 = 2:
𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥0 2−6 2−1
𝑓𝑓1 (𝑥𝑥) = 𝑓𝑓(𝑥𝑥0 ) + 𝑓𝑓(𝑥𝑥1 ) = (0) + (1.791759) = 0.358351
𝑥𝑥0 − 𝑥𝑥1 𝑥𝑥1 − 𝑥𝑥0 1−6 6−1
Que es el mismo resultado que nos arrojó la interpolación lineal de Newton.
Sabiendo que el valor real es ln 2 = 0.693147 calculamos el error verdadero:
0.693147 − 0.358351
𝜀𝜀𝑡𝑡 = � � ∗ 100% = 48.3%
0.693147
b) Estime el ln 2 empleando la interpolación cuadrática de Lagrange, considerando los siguientes datos:
𝑥𝑥0 = 1 𝑓𝑓(𝑥𝑥0 ) = 0
𝑥𝑥1 = 4 𝑓𝑓(𝑥𝑥1 ) = 1.386294
𝑥𝑥2 = 6 𝑓𝑓(𝑥𝑥2 ) = 1.791759
𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥1
𝑓𝑓2 (𝑥𝑥) = ∙ 𝑓𝑓(𝑥𝑥0 ) + ∙ 𝑓𝑓(𝑥𝑥1 ) + ∙ 𝑓𝑓(𝑥𝑥2 )
𝑥𝑥0 − 𝑥𝑥1 𝑥𝑥0 − 𝑥𝑥2 𝑥𝑥1 − 𝑥𝑥0 𝑥𝑥1 − 𝑥𝑥2 𝑥𝑥2 − 𝑥𝑥0 𝑥𝑥2 − 𝑥𝑥1
2−4 2−6 2−1 2−6 2−1 2−4
= ∙ (0) + ∙ (1.386294) + ∙ (1.791759) = 0.565844
1−4 1−6 4−1 4−6 6−1 6−4
Beatriz Adriana Jaime Fonseca
Lo que representa un error real del:
0.693147 − 0.565844
𝜀𝜀𝑡𝑡 = � � ∗ 100% = 18.36%
0.693147
c) Estime el ln 2 empleando la interpolación cúbica de Lagrange, considerando los siguientes datos:
𝑥𝑥0 = 1 𝑓𝑓(𝑥𝑥0 ) = 0
𝑥𝑥1 = 4 𝑓𝑓(𝑥𝑥1 ) = 1.386294
𝑥𝑥2 = 6 𝑓𝑓(𝑥𝑥2 ) = 1.791759
𝑥𝑥3 = 5 𝑓𝑓(𝑥𝑥3 ) = 1.609437
Calculamos primero los lagrangianos:
𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥3 2−4 2−6 2−5
𝐿𝐿0 (2) = ∙ ∙ = ∙ ∙ = 0.4
𝑥𝑥0 − 𝑥𝑥1 𝑥𝑥0 − 𝑥𝑥2 𝑥𝑥0 − 𝑥𝑥3 1 − 4 1 − 6 1 − 5
𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥2 𝑥𝑥 − 𝑥𝑥3 2−1 2−6 2−5
𝐿𝐿1 (2) = ∙ ∙ = ∙ ∙ =2
𝑥𝑥1 − 𝑥𝑥0 𝑥𝑥1 − 𝑥𝑥2 𝑥𝑥1 − 𝑥𝑥3 4 − 1 4 − 6 4 − 5
𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥3 2−1 2−4 2−5
𝐿𝐿2 (2) = ∙ ∙ = ∙ ∙ = 0.6
𝑥𝑥2 − 𝑥𝑥0 𝑥𝑥2 − 𝑥𝑥1 𝑥𝑥2 − 𝑥𝑥3 6 − 1 6 − 4 6 − 5
𝑥𝑥 − 𝑥𝑥0 𝑥𝑥 − 𝑥𝑥1 𝑥𝑥 − 𝑥𝑥2 2−1 2−4 2−6
𝐿𝐿3 (2) = ∙ ∙ = ∙ ∙ = −2
𝑥𝑥3 − 𝑥𝑥0 𝑥𝑥3 − 𝑥𝑥1 𝑥𝑥3 − 𝑥𝑥2 5 − 1 5 − 4 5 − 6
Sustituimos en la ecuación:
𝑓𝑓3 (𝑥𝑥) = 𝐿𝐿0 𝑓𝑓(𝑥𝑥0 ) + 𝐿𝐿1 𝑓𝑓(𝑥𝑥1 ) + 𝐿𝐿2 𝑓𝑓(𝑥𝑥2 ) + 𝐿𝐿3 𝑓𝑓(𝑥𝑥3 )
𝑓𝑓3 (2) = (0.4)(0) + (2)(1.386294) + (0.6)(1.791759) + (−2)(1.609437) = 0.628769
Lo que representa un error real del:
0.693147 − 0.628769
𝜀𝜀𝑡𝑡 = � � ∗ 100% = 9.28%
0.693147
Que como se puede observar, coincide con los resultados del Método Interpolación de Newton.
Beatriz Adriana Jaime Fonseca
