02 Notas - Interpolacao PDF
02 Notas - Interpolacao PDF
02 Notas - Interpolacao PDF
Cálculo Numérico
Notas de aulas
Interpolação Polinomial
Ouro Preto
2013
(Última revisão em novembro de 2013)
Depto de Computação – Instituto de Ciências Exatas e Biológicas – Universidade Federal de Ouro Preto
Sumário
1 - Introdução......................................................................................................................... 3
2 - Existência e unicidade do polinômio interpolador ........................................................... 4
3 - Erro de truncamento ......................................................................................................... 5
4 - Métodos de obtenção do polinômio interpolador............................................................. 7
4.1 – Método de Lagrange ................................................................................................. 7
4.2 – Método das diferenças divididas ............................................................................ 11
4.2.1 – O operador diferença dividida ......................................................................... 11
4.2.2 – O polinômio interpolador com diferenças divididas ....................................... 12
4.3 – Método das diferenças finitas ascendentes ............................................................. 16
4.3.1 – O Operador Diferença Finita Ascendente ....................................................... 16
4.3.2 – O polinômio interpolador com diferenças finitas ascendentes ........................ 18
5 – Complexidade dos métodos de interpolação ................................................................. 23
6 – Considerações finais ...................................................................................................... 23
Anexos ................................................................................................................................. 25
a) Teorema do Valor Médio ............................................................................................ 25
b) Operador linear ............................................................................................................ 26
Interpolação polinomial
1 - Introdução
Em geral, dispõe-se de dados que são fornecidos em um conjunto discreto de valores, den-
tro de um contínuo de possibilidades. Entretanto, pode ser necessário fazer estimativas em
pontos que estão entre os valores discretos, ou seja, não constam do conjunto. Ocorre, tam-
bém, a situação na qual se faz necessária uma versão simplificada de uma função compli-
cada. Ambas as aplicações são conhecidas como ajuste de curvas. Há duas abordagens
gerais para o ajuste de curvas, as quais se distinguem com base na quantidade de erro asso-
ciada com os dados.
Primeiro, quando os dados exibirem um grau significativo de erro, a estratégia será deter-
minar uma única curva que represente a tendência geral dos dados. Como cada ponto indi-
vidual poderá estar incorreto, não será feito qualquer esforço para passar a curva por todos
os pontos. Em vez disto, a curva é escolhida para seguir o padrão dos pontos considerados
como um grupo. Uma abordagem desta natureza é chamada de regressão por mínimos
quadrados.
Segundo, quando se souber que os dados são muito precisos, a abordagem básica é ajustar
uma curva ou uma série de curvas que passam diretamente por cada um dos pontos. Este
tipo de abordagem, que é o objeto deste texto, é chamada de interpolação.
Interpolar uma função, y = f(x), em um intervalo finito (a, b), consiste em substituí-la, ou
aproximá-la, por outra função, y = g(x). A necessidade de se utilizar este procedimento
ocorre, basicamente, quando a função:
a) não é conhecida na sua forma analítica, mas, apenas por meio de um conjunto de pontos
(xi, yi), i = 0, 1, ..., n; esta situação ocorre com muita freqüência, na prática, quando se
trabalha com dados obtidos de forma experimental;
b) é conhecida analiticamente, mas operações como a diferenciação e a integração são difí-
ceis (ou mesmo impossíveis) de realizar, ou seja, a função é de difícil tratamento.
Teoricamente, a função y = g(x) pode ser qualquer, mas o caso mais comumente conside-
rado é aquele em que pertence à classe das funções polinomiais.
A aproximação de funções por polinômios é uma das idéias mais antigas da análise numé-
rica, e ainda das mais utilizadas. É fácil entender a razão. Os polinômios são facilmente
computáveis, suas derivadas e integrais são, novamente, polinômios, seus zeros podem ser
determinados com facilidade, etc. O uso de polinômios interpoladores é importante, por
exemplo, para a obtenção de valores intermediários em tabelas, na integração numérica, no
cálculo de raízes de equações e na resolução de equações diferenciais ordinárias.
Prof. José Álvaro Tadeu Ferreira - Notas de aulas de Cálculo Numérico 3
Depto de Computação – Instituto de Ciências Exatas e Biológicas – Universidade Federal de Ouro Preto
As funções interpolantes polinomiais são as mais populares não só por suas propriedades
algébricas, mas, sobretudo, pela justificativa fornecida pelo teorema de aproximação de
Weierstrass que, de fato, garante a existência de um polinômio capaz de aproximar uma
função f tão bem quanto se queira.
Teorema (Weierstrass)
Se f é uma função contínua em um intervalo fechado [a, b], então, dado > 0, existe algu-
ma função polinomial, p, de ordem n = n(), tal que
|f(x) – p(x)| < , para x [a, b]
Objetivo
Sendo (xi, yi), i = 0, 1, ..., n; pontos, com abscissas distintas, de uma função y = f(x), obter
o polinômio, y = p(x) tal que:
p( x 0 ) a n x 0n a n - 1x 0n 1 ... a1x 0 a 0 y 0
Que é um sistema de (n + 1) equações lineares com (n + 1) incógnitas a0, a1, a2, ..., an. A
sua matriz dos coeficientes é dada por:
x n x n 1 ... x 1
0 0 0
x n x n 1 ... x 1
X 1 1 1
............................
x nn x nn 1 ... x n 1 (2.1)
Trata-se de uma ma matriz de Vandermonde. O seu determinante é calculado da seguinte
maneira
det(X) = (x0 – x1) (x0 – x2) ... (x0 – xn) (x1 – x2) (x1 – x3) ... (x1 – xn) ... (xn - 1 – xn)
Como, por condição, x0, x1, ..., xn são valores distintos, então tem-se que o determinante de
X é não nulo e o sistema linear admite solução única. Portanto, existe um único polinômio,
y = p(x), tal que p(xi) = f(xi) = yi, i = 0, 1, ..., n. Conclui-se, ainda, que o polinômio tem
grau máximo n, uma vez que os coeficientes, ai, i = 0, 1, ..., n; podem assumir qualquer
valor real, zero inclusive.
3 - Erro de truncamento
Teorema 3.1
Sejam:
(i) (xi, yi), i = 0,1, ..., n pontos com abscissas distintas de uma função y = f(x);
(ii) y = f(x) uma função com (n + 1) derivadas contínuas no intervalo [x0, xn].
Então, para cada x [x0, xn], existe um número ξ (x0, xn), que depende de x, tal que
f n 1 (( x ))
f ( x ) - p(x) E t ( x ) (x - x 0 ).(x - x1 ) ... (x - x n ).
(n 1)! (3.1)
A importância do Teorema 3.1 é mais teórica do que prática, uma vez que não é possível
determinar o ponto ξ. Na prática, para estimar o erro cometido, quando a função é dada
analiticamente, é utilizado o corolário a seguir.
Corolário 3.1
Se f(x) e suas derivadas até a ordem (n + 1) são contínuas no intervalo [x0, xn], então:
M
E t ( x ) (x - x 0 ).(x - x1 ) ... (x - x n ) .
(n 1)! (3.2)
Exemplo 3.1
Sabendo-se que os pontos a seguir são da função f(x) = x.e3.x, calcular um limitante superi-
or para o erro de truncamento quando se avalia y para x = 0,25.
i 0 1 2
xi 0,2 0,3 0,4
f(xi) 1,8221 2,4596 3,3201
Solução
De (3.2) tem-se que
M
E t ( x ) (x - x 0 ).(x - x1 ) ... (x - x n ) .
(n 1)!
Onde M = max f ' ' ' ( x ) no intervalo [0,2; 0,4]. Como f(x) = x.e3.x, segue que:
No intervalo [0,2; 0,4], f ’’’(x) é máxima para x = 0,4. Logo M = f ’’’(0,4) = 125,4998. Sen-
do assim:
125,4998
E t (0,25) (0,25 - 0,2).(0,25 - 0,3).(0,25 - 0,4) .
3!
E t (0,25) 0,0078
Note-se que y = p(x) não necessariamente converge para y = f(x) em [a, b] à medida que se
aumenta o número de pontos de interpolação. Polinômios interpoladores de grau elevado
podem produzir grandes oscilações nos extremos do intervalo, é o Fenômeno de Runge.
Este fenômeno demonstra que polinômios de grau elevado são normalmente pouco reco-
mendáveis para a interpolação porque aumentam o erro em valores próximos aos extremos
do intervalo de interpolação e melhoram a aproximação em valores próximos ao centro.
O problema pode ser evitado usando interpolação polinomial por partes com polinômios de
grau moderado. Desta forma, pode-se tentar diminuir o erro de interpolação aumentando o
número de peças de polinômios usadas, em vez de aumentar o grau do polinômio.
Exemplos típicos: interpolação linear por partes (uma reta para cada par de pontos) e inter-
polação quadrática por partes (uma parábola para cada três pontos), curvas spline.
Para determinar cada Li(x), i = 0,1, ..., n; basta considerar que todo xj, j = 0, 1, ..., n; é um
zero de Li(x) quando i j.
Seja a determinação de L0(x). Tem-se, por condição, que:
L0(x0) = 1
L0(xj) = 0; j = 1, 2, ..., n
1
c0
( x 0 x1 )(x 0 x 2 ) ( x 0 x n )
( x x1 )(x x 2 ) ( x x n )
L 0 (x)
( x 0 x1 )(x 0 x 2 ) ( x 0 x n ) (4.2)
L1(x1) = 1
L1(xj) = 0; j = 0, 2, ..., n
De modo análogo ao que foi feito anteriormente, para determinar o coeficiente c1 basta
considerar o valor numérico de L1(x) em x = x1 que, por condição, é igual a 1, obtendo-se
então
1
c1
( x1 x 0 )( x1 x 2 ) ( x1 x n )
( x x 0 )(x x 2 ) ( x x n )
L1 ( x )
( x1 x 0 )(x1 x 2 ) ( x1 x n ) (4.3)
( x x 0 )(x x1 ) ( x x i 1 )(x x i 1 ) ( x x n )
Li (x) , i = 0, 1, ..., n (4.4)
( x i x 0 )(x i x1 ) ( x i x i 1 )(x i x i 1 ) ( x i x n )
Exemplo 4.1
Seja y = f(x) uma função dada nos pontos a seguir. Utilizando interpolação polinomial,
método de Lagrange, determinar o polinômio que a interpola.
i 0 1 2 3
xi 0 1 2 4
yi 4 11 20 44
Solução
O polinômio interpolador é:
L(x) = y0.L0(x) + y1.L1(x) + y2.L2(x) + y3.L3(x)
Seja, então, a obtenção de Li(x), i = 0, 1, 2, 3
L(x) = x2 + 6.x + 4
Exemplo 4.2
Sendo y = f(x) uma função conhecida nos pontos:
i 0 1 2
xi 0,9 1 1,1
yi 0,6216 0,5403 0,4536
Pede-se:
(i) Utilizando interpolação polinomial, método de Lagrange, estimar o valor de y para
x = 1,07.
Solução
O polinômio interpolador é:
L(x) = y0.L0(x) + y1.L1(x) + y2.L2(x)
Neste item, pede-se para calcular L(1,07) que é dado por:
(x - x 0 )(x - x 1 ) (x - 0,9).(x - 1)
L 2 (x) L 2 (1,07) 0,5950
(x 2 - x 0 )(x 2 - x 1 ) (1,1 - 0,9).(1,1 - 1)
Portanto
L(1,07) = (0,6216).(-0,1050) + (0,5403).(0,5100) + (0,4536).(0,5950) ⇒ L(1,07) = 0,4802
(ii) Sabendo-se que os pontos dados são relativos à função y = cos(x), estimar o erro de
truncamento máximo cometido no item (i).
Solução
Sabe-se que o erro de truncamento máximo cometido é dado por:
M
E t ( x ) (x - x 0 ).(x - x1 ) ... (x - x n ) .
(n 1)!
onde M = máx|f n + 1(x)| no intervalo [x0, xn].
Tem-se que f ´´´(x) = sen(x), cujo módulo é máximo no intervalo [0,9; 1,1] para x = 1,1 e
f ´´´(1,1) = 0,8912 = M. Sendo assim,
0,8912
E t (1,07) | (1,07 - 0,9).(1,07 - 1).(1,07 - 1,1) | . E t (1,07) 5,3 x 10 - 5 0,0001
3!
Prof. José Álvaro Tadeu Ferreira - Notas de aulas de Cálculo Numérico 10
Depto de Computação – Instituto de Ciências Exatas e Biológicas – Universidade Federal de Ouro Preto
f ( x h) - f(x)
f ' (x) lim (4.5)
h 0 h
f ( x i h) - f(x i )
f ' (x i ) lim
h 0 h
Seja
xi + h = xi + 1 h = xi + 1 - xi
Sendo assim
f ( x i 1 ) - f(x i )
f ' (x i ) lim (4.6)
x i x i1 xi 1 - xi
Definição 4.2
Sendo (xi, yi), i = 0, 1, ..., n; um conjunto de pontos, com abscissas distintas, de uma fun-
ção y = f(x), define-se o operador diferença dividida de primeira ordem como:
f ( x i 1 ) - f(x i ) y i 1 - y i
Dyi , i = 0, 1, ..., n – 1 (4.7)
xi 1 - xi xi 1 - xi
Observe-se que este operador nada mais é do que uma aproximação do valor numérico da
primeira derivada de uma função em um ponto.
Pode ser demonstrado que as diferenças divididas de ordem superior são aproximações
para as derivadas de ordem superior.
A diferença dividida de segunda ordem é definida como:
Dy i 1 - Dy i
D 2 yi , i = 0, 1, ..., n – 2 (4.8)
xi 2 - xi
D 2 yi 1 - D 2 yi
D3 yi , i = 0, 1, ..., n – 3 (4.9)
xi 3 - xi
Prof. José Álvaro Tadeu Ferreira - Notas de aulas de Cálculo Numérico 11
Depto de Computação – Instituto de Ciências Exatas e Biológicas – Universidade Federal de Ouro Preto
Considerando as definições (4.7), (4.8) e (4.9), tem-se que a diferença dividida de ordem k,
é definida como:
D k - 1y i 1 - D k - 1y i k 1, 2, ..., n
D k yi , (4.10)
xi k - xi i 0, 1, ..., n - k
p(x) = a0 + a1.(x – x0) + a2.(x – x0)(x – x1) + ... + an.(x – x0)(x – x1) ... (x – xn - 1) (4.12)
Tendo em vista que y = p(x) deve ser tal que p(xi) = f(xi) = yi, i = 0, 1, ..., n
Então
Sabe-se que:
y1 - y 0
D Y0 Y1-Y0 = Dy0(x1-x0)
x1 - x 0
y -y
D Y1 2 1 Y2-Y1 = Dy1(x2-x1)
x 2 - x1
1
Contribuição do Professor José Américo Trivellato Messias
Prof. José Álvaro Tadeu Ferreira - Notas de aulas de Cálculo Numérico 12
Depto de Computação – Instituto de Ciências Exatas e Biológicas – Universidade Federal de Ouro Preto
Y2 - Y0 - Dy0 ( x 2 x 0 )
a2
( x 2 x 0 )( x 2 x1)
Dy 0 .x 1 Dy 0 .x 0 Dy1 ( x 2 x 1 ) - Dy 0 x 2 Dy 0 .x 0 Dy1 ( x 2 x 1 ) - Dy 0 x 2 Dy 0 .x 1
a2
( x 2 x 0 )(x 2 x1) ( x 2 x 0 )(x 2 x1)
Dy1 ( x 2 x 1 ) - Dy 0 ( x 2 x 1 )
a2
( x 2 x 0 )(x 2 x1)
Portanto
Dy1 - Dy 0
a2 (4.18)
x2 - x0
Com base na definição 4.8, conclui-se que 4.18 é a diferença dividida de segunda ordem.
Sendo assim
a2 = D2y0 (4.19)
p(x) = y0 + (x – x0) .Dy0 + (x – x0)(x – x1) .D2y0 + ... + (x – x0)(x – x1) ... (x – xn - 1).Dny0
(4.20)
f n ( )
Dn y0
n! (4.21)
Demonstração
Seja
e(x) = f(x) – p(x)
Onde p(x) é o polinômio que interpola f(x) nos pontos dados. Assim sendo, a função e(x)
tem n + 1 zeros distintos, o que implica, pelo Teorema de Rolle Generalizado, que e`(x)
tem n zeros em [x0, xn] e, assim, sucessivamente. Assim, conclui-se que existe um ξ ∈ [a,b]
tal que en(ξ) = 0. Ou seja
f n (x)
D n f (x)
n! (4.22)
Corolário 4.2
Se y = f(x) e suas derivadas até a ordem (n + 1) são contínuas no intervalo [x0, xn], então:
ET(x) ≤ |(x – x0). (x – x1). ... . (x – xn)|.máx|Dn + 1f(x)|
Tendo em vista o teorema 4.1 e o corolário 4.2, na ausência de informação sobre f n + 1(x),
uma estimativa para o erro de truncamento máximo pode ser obtida utilizando-se uma dife-
rença dividida de ordem (n + 1), caso estas não variem muito.
Exemplo 4.3
A tabela a seguir apresenta valores da voltagem, V, em função da corrente elétrica, I. Utili-
zando interpolação polinomial, método das diferenças divididas, estimar o valor de V
quando I = 3A.
i 0 1 2 3
I = xi 1 2 4 8
V = yi 120 94 75 62
Solução
Inicialmente, são determinados os valores das diferenças divididas.
Tem-se, então:
p(x) = y0 + (x - x0).Dy0 + (x - x0).(x – x1).D2y0 + (x - x0).(x – x1).(x – x2).D3y0
p(3) = 120 + (3 - 1).(- 26) + (3 - 1).(3 – 2).(5,5) + (3 - 1).(3 – 2).(3 – 4).(- 0,64)
p(3) = 80,28V
Exemplo 4.4
Uma barra de metal está presa em duas paredes separadas pela distância de 12m. A 5m da
parede A, um corpo apoiado sobre a barra faz com que esta toque no solo. Os pontos de
engate nas duas paredes estão a 8m (parede A) e 3m (parede B) do solo, conforme mostra a
figura a seguir. Usando interpolação polinomial, Método das Diferenças Divididas, pede-se
estimar:
3m
SOLO
d=12m
Solução
a) Os pontos a considerar são os da tabela a seguir.
i xi V = yi Dyi D2yi
0 0 8 - 1,6 0,169
1 5 0 0,429
2 12 3
Prof. José Álvaro Tadeu Ferreira - Notas de aulas de Cálculo Numérico 15
Depto de Computação – Instituto de Ciências Exatas e Biológicas – Universidade Federal de Ouro Preto
i xi yi Dyi D2yi
0 0 8
1 2 y
2 5 0
Para que este trecho seja representado por um polinômio de grau um, é necessário que a
diferença dividida de segunda ordem seja nula. Então, fazendo:
y = 4,8m
Da definição (4.24), verifica-se que o operador ∆(.) é linear (ver anexo), sendo assim, as
diferenças finitas ascendentes de ordem superior são definidas, por recorrência, da seguinte
maneira.
Segunda ordem.
Terceira ordem.
Prof. José Álvaro Tadeu Ferreira - Notas de aulas de Cálculo Numérico 16
Depto de Computação – Instituto de Ciências Exatas e Biológicas – Universidade Federal de Ouro Preto
k 1, 2, ..., n
∆kyi = ∆k - 1yi + 1 – ∆k - 1yi (4.28)
i 0, 1, ..., n - k
Teorema 4.3
Sendo y = f(x) uma função com derivadas contínuas até a ordem k, tem-se que:
Hipótese de indução
Admita-se que a relação vale para k – 1.
∆k – 1f(x) = hk – 1.fk – 1(ξk – 1), ξk – 1 ∈ (x, x + (k – 1).h)
Passagem de indução
Provar que a relação é válida para k.
∆k[f(x)] = ∆k - 1[∆[f(x)]] = ∆k - 1[f(x + h) − f(x)]
= ∆k - 1[f(x + h)] − ∆k - 1[f(x)]
Corolário 4.3
[∆kf(x) / hk].é uma aproximação para f (k)(x) e o erro cometido tende a zero quando h tende
a zero.
k y i
D k yi , k = 0, 1, 2, ..., n; i = 0, 1, 2, ..., n – k (4.31)
h k .k!
Demonstração:
A demonstração é feita por meio de indução finita em k.
Hipótese de indução
Admita-se que o argumento é válido para a ordem k – 1.
k 1 y i
D k -1 y i , i = 0,1, 2, … , n – k + 1
h k 1.(k 1)!
Passagem de indução
Provar que é válido para ordem k. Por definição
D k 1y i 1 D k 1y i
D k yi , i = 0,1, 2, … , n – k
xi k xi
k 1 y i 1 k 1 y i
h k 1.(k 1)! h k 1.(k 1)!
D k yi , i = 0,1, 2, … , n – k
k.h
k 1y i 1 - k 1 y i
D k yi , i = 0,1, 2, … , n - k
h.h k 1.k.(k 1)!
Portanto
k k y i
D yi , k = 0, 1, 2, ..., n; i = 0, 1, 2, ..., n – k
h k .k!
c.q.d.
Seja a variável
x - x0
z (4.32)
h
De onde vem que
x = x0 + h.z
x - x0 = h.z
x – x1 = x – (x0 + h) = x – x0 – h = h.z – h = h.(z – 1)
x – x2 = x – (x0 + 2.h) = x – x0 – 2.h = h.z – 2.h = h.(z – 2)
.
.
.
x – xn - 1 = h.[z – (n - 1)]
Exemplo 4.5
Os pontos a seguir relacionam a solubilidade, S, da água no óleo mineral, em partes por
milhão, com a temperatura, t, em graus centígrados. Utilizando interpolação polinomial,
método das diferenças finitas ascendentes, estime o valor de t quando S = 200ppm.
S t
i xi yi Δyi Δ2yi Δ3yi
0 30 15 35 -19 13
1 130 50 16 -6
2 230 66 10
3 330 76
Sabe-se que
x - x0 200 - 30
z z 1,7
h 100
Logo
o
Tem-se, então, que p(200) = 62,4 C
Exemplo 4.6
Uma hidroelétrica tem capacidade máxima de 60MW, que é determinada por três gerado-
res de 30MW, 15MW e 15MW, respectivamente. A demanda de energia varia num ciclo
de 24h, sendo que a demanda mínima ocorre entre 2h e 5h e a máxima entre 14h e 17h.
Utilizando interpolação polinomial, método das diferenças finitas ascendentes, estime a
demanda mínima e a máxima e o horário em que cada uma ocorre, considerando os dados
a seguir.
i 0 1 2 3
Hora (xi) 2 3 4 5
Demanda (yi) 16,4 15,2 14,9 16,0
i 0 1 2 3
Hora (xi) 14 15 16 17
Demanda (yi) 36,5 43,0 34,0 31,2
Solução
Demanda mínima
Inicialmente, são calculados os valores das diferenças finitas ascendentes.
i xi yi yi 2yi 3yi Sendo
0 2 16,4 - 1,2 0,9 0,5 x - x0
z
1 3 15,2 - 0,3 1,4 h
2 4 14,9 1,1
3 5 16,0 então z = x – 2 e x = 2 + z
Assim,
z(z 1) z(z 1)(z 2)
p(2 z) 16,4 z.(1,2) (0,9) (0,5)
2! 3!
Para estimar a demanda mínima basta obter a primeira derivada de p(2 + z) e determinar os
seus zeros. Tem-se, então:
Trata-se de uma função do segundo grau. Seus zeros são – 3,46, que não tem sentido para
este problema, e 1,79. A questão, agora, é verificar que z = 1,79 é abscissa de ponto de
mínimo. Para isto toma-se a segunda derivada de p(2 + z) e verifica-se, facilmente, que:
Demanda máxima
Cálculo das diferenças finitas ascendentes.
i xi yi yi 2yi 3yi Sendo
0 14 36,5 6,5 - 15,5 21,7 x - x0
z
1 15 43,0 -9 6,2 h
2 16 34,0 - 2,8
3 17 31,2 então z = x – 14 e x = 14 + z
Assim,
z(z 1) z(z 1)(z 2)
p(14 z) 36,5 z.(6,5) (15,5) (21,7)
2! 3!
Cujos zeros são z = 0,74 e z = 2,69. Basta, agora, calcular o valor numérico da segunda
derivada de p(14 + z) em cada um destes pontos para verificar qual deles é abscissa de
ponto de máximo.
Sendo
p’’(14 + z) = 21,72.z - 37,2
Para z = 0,74, tem-se que p’’(14,74) = - 21,13 e, para z = 2,69, p’’(16,69) = 21,23. Portan-
to, z = 0,74 é abscissa de ponto de máximo e, calculando o valor numérico do polinômio
interpolador neste ponto, tem-se a estimativa para a demanda máxima que é
p(14,74) = 43,7MW
e verifica-se que ela ocorre às 14h44min, aproximadamente.
Diferenças finitas n2 n 3 2 5
.n .n 1
n2 + n + 1 n
ascendentes 2 2 2 2
Tabela 5.1: Complexidade dos métodos de interpolação (n é o grau do polinômio)
Tomando como exemplo um polinômio interpolador de grau dez verifica-se que o número
total de operações efetuadas pelo Método de Lagrange é igual a 450, pelo Método das Di-
ferenças Divididas 285 e, pelo Método das Diferenças Finitas Ascendentes, 176. O que
leva a verificar que o Método das Diferenças Finitas Ascendentes apresenta maior eficiên-
cia quando comparado com os outros dois métodos estudados.
6 – Considerações finais
(a) Os métodos que utilizam diferenças (divididas ou finitas ascendentes) são eficientes
quando se deseja aumentar (ou diminuir) o grau do polinômio obtido, pois basta, sim-
plesmente, acrescentar (ou retirar) termos. Logo, para cálculos exploratórios, estes mé-
todos, em geral, são preferíveis.
(b) No método de Lagrange a alteração do grau do polinômio exige que os cálculos sejam,
todos, refeitos.
(c) O método de Lagrange ocupa menos memória, uma vez que não é necessário o cálculo
e o armazenamento de uma tabela de diferenças divididas ou finitas ascendentes.
Anexos
a) Teorema do Valor Médio
Se y = f(x) é uma função que satisfaz as condições:
(i) é contínua no intervalo fechado [ a, b ]
(ii) é derivável no intervalo aberto ( a, b )
Então, existe pelo menos um número c em (a, b), tal que:
f(b) - f(a)
f ' ( c)
b-a
Geometricamente, o teorema do valor médio diz que se f é uma função "suave" que liga os
pontos A= ( a ,f( a )) e B = ( b ,f( b )) existe pelo menos um ponto c, entre a e b , tal que a
reta tangente ao gráfico de f em c é paralela a reta secante que passa por A e B. A figura a
seguir ilustra o teorema.
O teorema do valor médio é a tradução matemática para um fato que aparece de forma cor-
riqueira em muitas situações cotidianas. Por exemplo, se a média de velocidade, em uma
viagem de carro é de 80 km/h, então, em algum momento da viagem, o velocímetro do
carro deve ter marcado 80km/h.
Para traduzir a afirmação em termos matemáticos, considere-se que s(t) é a posição do car-
ro em um instante t. Se a viagem começa em t = a (horas) e termina em t = b (horas), a ve-
locidade média é dada por:
s(b) - s(a)
vm
b-a
A afirmação de que em algum momento da viagem a velocidade instantânea deve ser igual
à velocidade média, significa que em algum tempo c tem-se:
s(b) - s(a)
vm v(c) s ' (c)
b-a
O Teorema do Valor Médio estabelece as condições mínimas que uma função s deve satis-
fazer para que a igualdade acima seja verdadeira.
b) Operador linear
Um operador é linear se, e somente se:
(i) .(u ± w) = .u ± .w
(ii) .(k.w) =k.( .w), onde k é uma constante real