PDF de programación - Fórmulas recursivas de Neville para el polinomio interpolante

Fórmulas recursivas de Neville
para el polinomio interpolante

Objetivos. Demostrar las fórmulas recursivas de Neville para el polinomio interpolante
y practicarla con un algunos ejemplos.

Requisitos. Concepto de interpolación polinomial, existencia y unicidad del polinomio
interpolante.

1. Notación para el método de Neville. Sean x1, . . . , xn ∈ R algunos números diferen-
tes a pares, y1, . . . , yn ∈ R algunos números, y sean i, j ∈ {1, . . . , n}, i ≤ j. Denotemos por
Pi,j al polinomio interpolante asociado a las abscisas xi, . . . , xj y las ordenadas yi, . . . , yj.
En otras palabras, Pi,j es un polinomio de grado ≤ j − i que tiene valores yi, . . . , yj en los
puntos xi, . . . , xj:

deg(Pi,j) ≤ j − i,

∀k ∈ {i, . . . , j} Pi,j(xk) = yk.

2. Ejemplo. Sean x1 = −1, x2 = 2, x3 = 3, x4 = 5, f(x) = ex. Construir P1,2.

3. Teorema (fórmula recurrente de Neville para los polinomios interpolantes).
Sean x1, . . . , xn ∈ R algunos números diferentes a pares, y sean y1, . . . , yn ∈ R algunos
números. Entonces para cualesquiera i, j ∈ {1, . . . , n} con i ≤ j,

Pi,j(x) =

(x − xi)Pi+1,j(x) − (x − xj)Pi,j−1(x)

xj − xi

.

Idea de demostración. Sea Q(x) el polinomio en el segundo miembro de la fórmula. En-
tonces es fácil ver que deg(Q) ≤ j − i y Q(xk) = yk para todo k ∈ {i, . . . , j}. Notamos que
los polinomios Pi,j y Q son de grados ≤ j − i y toman los mismos valores en j − i + 1 puntos
xi, . . . , xj. Por el teorema principal sobre la interpolación polinomial, los polinomios Q y
Pi,j coinciden.

4. El método de Neville para n = 2.

P1,1(x) = y1,

P2,2(x) = y2,

P1,2(x) = (x−x1)P2,2(x)−(x−x2)P1,1(x)

x2−x1

P3,3(x) = y3,

P2,3(x) = (x−x2)P3,3(x)−(x−x3)P2,2(x)

x3−x2

,

,

P1,3(x) = (x−x1)P2,3(x)−(x−x3)P1,2(x)

x3−x1

.

Fórmulas de Neville para el polinomio interpolante, página 1 de 2

5. Ejemplo. Aplique el método de Neville para construir el polinomio interpolante P tal
que P(−1) = 15, P(4) = 5, P(5) = 9.

x1 = −1,

P1,1(x) = 15,

x2 = 4,

x3 = 5

P2,2(x) = 5,

P1,2(x) = −2x + 13,

P3,3(x) = 9,

P2,3(x) = 4x − 11,

P1,3(x) = x2 − 5x + 9.

Respuesta: P(x) = x2 − 5x + 9.

6. Ejemplo. Aplique el método de Neville para aproximar
y los puntos x1 = 1, x2 = 4, x3 = 9.

√

2 con la función f(x) =

√

x

7. Ejercicio. Aplique el método de Neville para construir el polinomio interpolante P tal
que P(−2) = −7, P(−1) = −6, P(3) = 18. Respuesta: P(x) = x2 + 4x − 3.

8. Ejemplo. Aplique el método de Neville para aproximar
y los puntos x1 = 0, x2 = 1, x3 = 2.

√
2 con la función f(x) = 2x

9. Problema de programación. En un lenguaje de programación escribir una función
que calcule los coeficientes del polinomio interpolante usando la fórmula recursiva de
Neville. Se recomienda guardar los polinomios auxiliares en una tabla, de la misma manera
como se guardan las diferencias divididas en algorirmo de Newton. Calcular el número de
operaciones de multiplicación y división que se ejecutan en este algoritmo.

Fórmulas de Neville para el polinomio interpolante, página 2 de 2