Erro na Interpolação

Teoria Completa

Introdução Erro Exato

Fala meus amigos, tudo certinho? Veja bem, quando nós estudamos anteriormente os polinômios interpoladores nós vimos que são aproximações de uma função real, e essas aproximações podem ser tãaaaaaaao próximas quanto nós quisermos. O lance é, como toda aproximação, nossa função terá um erro (tava demorando né?).

Então, vamos dar uma estudada nesse erro da interpolação! Vamos aprender um negócio que basicamente é apresentar um limitante pro erro envolvido na aproximação de uma função feita por um polinômio interpolador.

Imagina a seguinte situação:

Seja f x = l n   x   + 3 x 2 . Usando a tabela abaixo encontre um polinômio interpolador de grau 2 e calcule:

(a) uma aproximação para f 2.5

(b) o erro desta aproximação

(c) um limitante para o erro

(d) uma aproximação para o erro assumindo que f x não é conhecida.

Usaremos esse exemplo para explicar a teoria, vem comigo!

Primeiro o item (a).

O exercício está pedindo uma aproximação para f 2.5 . Já aprendemos vários métodos para calcular uma aproximação e podemos usar qualquer um deles porém, como em seguida o exercício pede o erro, faça pelo método de Newton (a menos que o enunciado mande você fazer por outro método, é claro). Isso porque para calcular a aproximação do erro precisamos da tabela de diferenças divididas! Então já que vamos ter que montá-la de qualquer jeito, aproveitamos e já fazemos a aproximação pelo método de Newton.

Uma dica: se você quer estimar o erro em um polinômio de grau n , calcule os operadores até grau n + 1 já que vamos precisamos do maior operador diferenças divididas de ordem n + 1 ! (não se preocupe logo fará sentido!). Nesse caso, o polinômio que o exercício quer é de grau 2 então vamos fazer a tabela até d 3 ,

Passo 1. Vamos montar nossa tabela! Ela fica assim:

Repare que com 5 nós poderíamos ter calculado até a Ordem 4 , mas como não vamos precisar dela melhor evitarmos trabalho haha.

Passo 2. Como queremos uma aproximação para f 2.5 , vamos usar x 0 = 1 , x 1 = 2 e x 2 = 3 . Da tabela, temos:

d 0 = 3

d 1 = 9.6931

d 2 = 2.8562

Nosso polinômio de grau 2 é:

p 2 x = d 0 + d 1 x - x 0 + d 2 x - x 0 x - x 1

p 2 x = 3 + 9.6931 x - 1 + 2.8562 x - 1 x - 2

p 2 x = 2.8562 x 2 + 1.1245 x - 0.9807

Passo 3. Nossa aproximação será:

p 2 2.5 = 19.6818

Erro Exato

Para resolver o item (b), precisamos saber o que é o erro exato.

Como já falamos anteriormente, sempre que realizamos uma aproximação não obtemos o valor verdadeiro da função certo? (justamente por se tratar de uma aproximação). Existe um erro entre o polinômio e a função, e a diferença entre os dois nós chamamos de Erro Exato.

Tendo uma função f x e sua função interpoladora p n x , calculamos o erro como:

E x = f x - p n x

Pronto! já sabemos responder o item (b) do nosso exemplo:

O erro exato é:

E x = f x - p n x

E 2.5 = p 2 2.5 - f 2.5 = 19.6818 - 19.6663 = 0.0155

Limitante para o erro

O item (c) pede uma limitante para o erro. Ou seja, ele quer saber o erro máximo que seu polinômio interpolador vai ter naquele intervalo que você está aproximando f x . E como fazemos isso?

Seja f x uma função definida e (n+1) vezes diferenciável num intervalo [a,b]. Nós temos os pontos x 0 , x 1 , ..., x n , (n+1), que são pontos igualmente espaçados.

Ou seja, se x 1 - x 0 = x 2 - x 1 = … = x n - x n - 1 = h deste intervalo. Se P x interpola f x nestes pontos, então o limitante superior para o erro é dado:

E x ≤ h n + 1 M n + 1 4 n + 1

E o M n + 1 o valor máximo do módulo das derivadas de ordem n + 1 de f x no intervalo x 0 , x n , ou seja, o maior valor em módulo das derivadas da função aplicadas no intervalo x 0 , x n .

M n + 1 = f n + 1 x  

Mas, o que são pontos não serem equidistantes? Então:

E x ≤ x - x 0 x - x 1 … x - x n M n + 1 n + 1 !

Vamos tentar resolver o item (c) agora?

Para encontrarmos um limitante para o erro precisamos da derivada terceira de f x para calcular o valor Máximo M n + 1 :

f 3 x = 2 x 3

Como essa função é contínua, decrescente e positiva no intervalo 1,3 , ela terá seu maior valor em módulo em x = 1 , onde:

M n + 1 = f n + 1 x   = f 3 1 = 2 1 3 = 2

Como os valores de x são igualmente espaçados de h = 1 :

E x ≤ h n + 1 M n + 1 4 n + 1 = 1 3 ∙ 2 4 3 = 0.1667

Repare que usando essa forma, o resultado independe do valor x , ou seja, obtemos o valor máximo do erro no intervalo 1,3 .

Estimativa para o erro

No último item do nosso exemplo pede uma aproximação para o erro assumindo que f x não é conhecida. Veja bem, o problema vem agora. Geralmente não conhecemos nossa f x (daí a necessidade de encontrarmos uma aproximação), e o nosso ξ x nunca é conhecido, na maioria dos casos não podemos calcular o erro e nem um limitante para ele! Além disso, nem sempre a aproximação baseada em todos os pontos é a que se aproxima mais do valor verdadeiro.

Lembra que lá no começo, na hora de calcular a aproximação pelo método de Newton, eu falei para calcularmos d n + 1 ? Uma ordem acima da necessária para o polinômio? Então agora que ele entra em cena.

Bom, por sorte algum maluco percebeu que:

M n + 1 n + 1 ! ≈ d n + 1  

Onde d n + 1 são os operadores diferenças divididas de ordem n + 1 ! Nesse caso, dizemos que:

E x ≈ x - x 0 x - x 1 … x - x n d n + 1  

Voltando agora para o nosso exemplo:

Imaginamos que não conhecemos a f x e vamos calcular uma aproximação para o erro. O maior operador diferenças divididas de ordem 3 é:

d 3 = 0,0283

Aplicando agora a fórmula:

E x ≈ x - x 0 x - x 1 … x - x n d n + 1  

E 2.5 ≈ 2.5 - 1 2.5 - 2 2.5 - 3 0.0283 = 0.0106

Que é uma aproximação razoável para 0.0155 .

Pronto galera! Concluímos o exemplo. Partiu botar em prática agora!

Expandir

Exercício Resolvido #1

Márcia A. Gomes Ruggiero e Vera Lúcia da Rocha Lopes, Cálculo Numérico – Aspectos Teóricos e Computacionais, 2ª ed. São Paulo: MAKRON Books, 1996, pag 256, exercício 2

Verifique que na interpolação linear

E x ≤ h 2 M 2 8

onde h = x 1 - x 0 .

MOSTRAR SOLUÇÃO COMPLETA

Passo 1

-

Como na interpolação linear nós só temos dois pontos, sempre podemos dizer que os pontos são igualmente espaçados haha. Logo, sempre é válida a fórmula:

E x ≤ h n + 1 M n + 1 4 n + 1

Substituindo n = 1 já que estamos usando interpolação linear,

E x ≤ h 2 M 2 8

Resposta

Exercício Resolvido #2

Márcia A. Gomes Ruggiero e Vera Lúcia da Rocha Lopes, Cálculo Numérico – Aspectos Teóricos e Computacionais, 2ª ed. São Paulo: MAKRON Books, 1996, pag 257, exercício 5

Queremos construir uma tabela que contenha valores de cos ⁡ x para pontos igualmente espaçados no intervalo I = 1,2 .

Qual deve ser o menor número de pontos desta tabela para se obter, a partir dela, o cos ⁡ x , usando interpolação linear com erro menor que 10 - 6 para qualquer x no intervalo 1,2 .

MOSTRAR SOLUÇÃO COMPLETA

Passo 1

Nesse caso usaremos o número de pontos como uma forma de “dimensionar o h ”. Em uma interpolação linear temos apenas dois pontos separados por uma distância h . Logo, sempre podemos dizer que:

E x ≤ h 2 M 2 8

Vamos começar encontrando M 2 . Precisamos da derivada segunda de cos ⁡ x :

f x = cos ⁡ x

f ' x = - sin ⁡ x

f ' ' x = - cos ⁡ x

Como o erro deve ser menor que 10 - 6 no intervalo 1,2 , usamos o valor máximo de - cos ⁡ x nesse intervalo como M 2 . E esse valor é em 1 , onde - cos ⁡ 1 ≅ 0.54 .

Passo 2

Agora precisamos escrever h em função do número de pontos. Então vamos montar uma tabelinha:

Com essa relação podemos substituir tudo na fórmula:

E x ≤ 1 n - 1 2 0.54 8 = 0.54 8 n 2 - 2 n + 1 = 10 - 6

n 2 - 2 n + 1 = 67.5 × 10 3

Resolvendo a equação de segundo grau, temos:

n = - 258.808   o u   260.808

Como não faz sentido termos um número de pontos negativo e nem quebrado, n = 261 pontos.

Resposta

n = 261

Exercício Resolvido #3

Lista 4 – Professor Mauro – UFRJ – Questão 5

Deseja-se construir uma tabela da função f x = e x no intervalo 0,1 com pontos x i igualmente espaçados. Seja h o passo. Qual o valor máximo de h para que o erro da interpolação linear em qualquer ponto do intervalo seja menor ou igual a ϵ = 10 - 2 .

MOSTRAR SOLUÇÃO COMPLETA

Passo 1

Em uma interpolação linear temos apenas dois pontos separados por uma distância h . Logo, sempre podemos dizer que:

E x ≤ h 2 M 2 8

Vamos começar encontrando M 2 . A derivada segunda de