Aula 08

Interpolação polinomial
Método de Lagrange
Método de Newton
ERRO de Interpolação
GEX 114 – Cálculo Numérico
Profª Evelise
 Motivação:
• Considere a tabela abaixo sobre o calor específico da
água relacionado a algumas temperaturas:

Às vezes é necessário encontrar o calor específico em

uma temperatura que não foi medida.
Como proceder???

2
 Definição de interpolação
Interpolar uma função f(x) consiste em
aproximar essa função por uma outra função
g(x), escolhida entre uma classe de funções pré-
definidas, e que satisfaça a algumas
propriedades.

3
 Conceito de interpolação
Considere como:
• Nós de interpolação: (n+1) pontos distintos.
• Os valores de f(x) nesses pontos
 p( x0 )  f ( x0 )
 p( x )  f ( x1 )
 1
 p( x2 )  f ( x2 )


 p( xn )  f ( xn )
4
 Graficamente, para n=5

Note que para n=5 teremos 6 pontos de interpolação.

5
 Interpolação Polinomial
Dados os pontos (x0,f(x0)), (x1,f(x1)), ..., (xn,f(xn)), queremos
encontrar um polinômio da forma:
f ( xk )  pn ( xk ) k=0,1,2,...,n
A partir dessa condição, levando em conta que
pn ( x)  a0  a1 x  a2 x  ...  an x
2 n

Montamos o seguinte sistema linear:

6
 Matriz das incógnitas:

• Note que V é uma matriz de Vandermonde, e logo, Det(A)≠0,

e a solução do sistema é única. 7
 Formas de se obter os polinômios

1. Solução de sistema linear

2. Polinômio de Lagrange
3. Polinômio de Newton

8
 Exemplo 01- Resolução do sistema
linear
• Encontre o polinômio de grau ≤ 2, que
interpola os pontos da tabela

x -1 0 2
F(x) 4 1 -1

9
 Forma de Lagrange

• Ao invés de utilizarmos um polinômio

interpolador da forma:

pn ( x)  a0  a1 x  a2 x  ...  an x
2 n

• Lagrange optou por efetuar a interpolação

com um polinômio da forma:

10
 Como obter Lk(x) ?

• Os termos Lk(x) são chamados de fatores de Lagrange

• Fatores de Lagrange são obtidos através da

formulação:

Lk(xi) =

11
 Forma final do polinômio
interpolador
• Polinômio:
n
pn ( x)   yk Lk
k 0

• Fatores de Lagrange:
n

 (x  x )
j 0
j

j k
Lk ( x)  n

 (x
j 0
k  xj )
j k
12
 Exemplo 02

• Encontre o polinômio de Lagrange, que

interpola os pontos da tabela

x -1 0 2
F(x) 4 1 -1

13
 Exemplo 03 – Tarefa

• Dada a função f(x)=5x+ln(x+1), encontre o

polinômio de Lagrange que passa pelos
pontos x=1 e x=2. Calcule também P(1.1)

14
 Exemplo 04
• Os dados da tabela a seguir representam a
velocidade de um móvel em relação ao tempo.
Use-os para calcular aproximadamente a
velocidade do móvel em t=15.
t (Tempo 0 10 20 30 40
em s)
v 0 30 75 90 100
(velocidad
e em m/s)

Obs: Escolha o grau do polinômio

15
 Exemplo 05: Interpolação Linear - Tarefa

Faça a interpolação genérica entre dois pontos

distintos (x0,f(x0)), (x1,f(x1)). Como temos dois
pontos, n=1, e esta interpolação é chamada
de interpolação linear.

*lousa

16
 Exemplo 05 - Tarefa

• Utilize o método de Lagrange para encontrar

um polinômio capaz de estimar a taxa de
metabolismo de um ser vivo em função de sua
massa, baseada nos seguintes dados:

Animal Massa, kg Metabolismo,

watts
Vaca 400 270
Homem 70 82
Ovelha 45 50

17
 Forma de Newton
• Ao invés de utilizarmos um polinômio
interpolador da forma:

pn ( x)  a0  a1 x  a2 x  ...  an x
2 n

• A idéia de Newton foi interpolar dados pré

definidos aplicando um polinômio da forma:

Operador das diferenças divididas 18

 Como obter os operadores de diferencias
divididas?
Para obter os operadores, basta fazer:

19
Forma genérica dos operadores:

20
Outra forma de visualizar:

21
Exemplo 06 – Treino de montagem

Monte uma tabela de diferenças divididas a

partir dos seguintes dados:

x -1 0 1 2 3
F(x) 1 1 0 -1 -2

22
Exemplo 07 – Execução do método

Utilize a forma de Newton para encontrar o

polinômio P2(x) que interpola os pontos dados
na tabela:
x -1 0 2
F(x) 4 1 -1

23
 Exemplo 08 - Aplicação

• Utilize a forma de Newton para encontrar um

polinômio capaz de estimar a taxa de
metabolismo de um ser vivo em função de sua
massa, baseada nos seguintes dados:

Animal Massa, kg Metabolismo,

watts
Vaca 400 270
Homem 70 82
Ovelha 45 50

24
 Formulação Matemática do Erro

• Se o erro é a diferença entre a função original e o

polinômio, podemos escrever:

En ( x)  f ( x)  pn ( x)

Para todo x no intervalo [x0 , xn]

Este erro é o erro verdadeiro
e para calculá-lo precisamos
da especificação de f(x).

25
 TEOREMA
Sejam x0 < x1 < x2 <...< xn com n+1 pontos.
Seja f(x) com derivadas até ordem (n+1) para todo x
pertencente ao intervalo [x0 , xn] e pn(x) o polinômio
interpolador de f(x) nos nós.
Então, em qualquer ponto x do intervalo, o erro é dado
por:

26
 Problemas com o teorema:
• Uso muito limitado na prática, pois:
1. A derivada não é conhecida na maior parte dos
casos.
2. O ponto ξ (arbitrário para x) nunca é conhecido

Formulação serve somente para um estudo teórico de

embasamento para o cálculo do limitante do erro.
SOLUÇÃO: Encontrar um limitante para o erro

27
 Corolário
Sob as hipóteses do Teorema, se as derivadas de f forem
contínuas, teremos a seguinte relação:

M n 1
En ( x)  f ( x)  pn ( x)  ( x  x0 )( x  x1 )...( x  xn )
(n  1)!
( n 1)
M n 1  max f ( x)
xI Derivada da função

OBS IMPORTANTE: Corolário se aplica somente para o caso

em que conhecemos a função f(x), pois precisamos
calcular a derivada.
28
 Exemplo 09
Seja a função f(x)= ex + x -1 tabelada abaixo.
x 0 0.5 1 1.5 2
f(x) 0.0 1.1487 2.7183 4.9811 8.3890

Obtenha f(0.7) aplicando método de Newton de

primeira ordem e faça uma análise do erro
cometido da seguinte forma:
a) Calcule o erro verdadeiro cometido, com quatro
casas decimais (É possível, pois temos a função
f).
b) Calcule um limitante para o erro.
29
 Estimativa para o erro
Aplica-se quando:
• A função f(x) não é conhecida e os valores estão
especificados somente numa tabela.
• Não dá pra calcular o valor de Mn+1
SOLUÇÃO: Aproximar M pelo maior valor em módulo das
diferenças divididas até ordem n+1 (ordem das derivadas)

• A estimativa do erro será dada por:

En ( x)  ( x  x0 )( x  x1 )...( x  xn ) max(diferenças | divididas )

30
 Exemplo 10

Considere o conjunto de dados especificados na

tabela:
x 0.2 0.34 0.4 0.52 0.6 0.72

f(x) 0.16 0.22 0.27 0.29 0.32 0.37

a) Obtenha f(0.47) usando a Forma de Newton um

polinômio de grau 2 e quatro casas decimais.
b) Calcule uma estimativa para o erro.

31
 Exemplo 11
Considere o conjunto de dados especificados na
tabela:
x 0 0.2618 0.5234 0.7854 1.0472

f(x) 0 1,0353 2 2.8284 3.4641

a) Obtenha f(0.6) usando a forma de Newton de

grau 2, grau 3 e quatro casas decimais.
b) Calcule uma estimativa para o erro para cada
uma das formas encontradas na parte a. Qual é
mais vantajosa?
32