Polinômio de Newton

Definição

O polinômio de Newton é definido como:

P(x)=y0+xx0x1x0(y1y0)+(xx0)(xx1)(x2x0)(x1x0)(y2y0x2x0x1x0(y1y0))+...

onde xi são os pontos de interpolação e yi são as respectivas funções.

Cálculo Do Polinômio

Para calcular o polinômio de Newton, é necessário seguir os seguintes passos:

  1. Defina os pontos de interpolação: Escolha os pontos x0,x1,,xn onde deseja realizar a interpolação.
  2. Calcule as diferenças divididas: Calcule as diferenças divididas xixjxkxj para cada par de pontos (xi,xj) e (xk,xj).
  3. Calcule os coeficientes do polinômio: Utilize as diferenças divididas para calcular os coeficientes do polinômio P(x).

Exemplos

Interpolação de Funções

Suponha que deseja interpolar a função f(x)=x2+1 nos pontos x0=1,x1=0,x2=1. O polinômio de Newton pode ser calculado como:

P(x)=f(1)+x+11(1)(f(0)f(1))+(x+1)(x1)1(1)2(f(1)f(1)1(1)1(1)(f(0)f(1)))

Interpolação de Dados

Suponha que deseja interpolar os dados (x0,y0)=(2,3),(x1,y1)=(4,5),(x2,y2)=(6,7). O polinômio de Newton pode ser calculado como:

P(x)=y0+xx0x1x0(y1y0)+(xx0)(xx1)(x2x0)(x1x0)(y2y0x2x0x1x0(y1y0))

Propriedades Do Polinômio de Newton

Aplicações Do Polinômio de Newton

Limitações Do Polinômio de Newton

Erro na Interpolação Polinomial de Newton

Dado um conjunto de pontos x0,x1,,xn distintos, e uma função f interpolada por um polinômio de Newton Pn(x), o erro da interpolação em um ponto x (não necessariamente um dos xi) é dado por:

f(x)Pn(x)=f(n+1)(ξ)(n+1)!(xx0)(xx1)(xxn)

com:

Derivação Do Erro em Newton

A fórmula de Newton baseia-se em diferenças divididas:

Pn(x)=f[x0]+f[x0,x1](xx0)++f[x0,,xn](xx0)(xxn1)

Mas existe um resultado teórico importante:

A diferença dividida de ordem n+1, f[x0,x1,,xn,x], está relacionada com a derivada de ordem (n+1) da função:

f[x0,x1,,xn,x]=f(n+1)(ξ)(n+1)!

Assim, adicionando esse termo ao polinômio de grau n, temos a forma completa com erro:

f(x)=Pn(x)+f(n+1)(ξ)(n+1)!(xx0)(xx1)(xxn)

Estimativa Do Erro (limitante superior)

Se |f(n+1)(t)|M em todo t[a,b], então o erro máximo possível em x é:

|f(x)Pn(x)|M(n+1)!|(xx0)(xx1)(xxn)|

Exemplo

Interpole f(x)=cos(x) em [0,π/2] com os pontos x0=0, x1=π/4, x2=π/2.

Derivadas de f(x):

Logo, |f(3)(x)|1 para todo x[0,π/2], então M=1.

Erro em x=π/3:

ω(x)=(π30)(π3π4)(π3π2)=π3π12(π6)

Então:

|f(x)P2(x)|13!|π33126|=16|π3216|=π31296