Métodos Numéricos para Derivadas: Aumentando a Acurácia

Calcular a derivada de uma função numericamente é uma tarefa comum em muitas áreas da ciência e engenharia. Uma das primeiras e mais intuitivas abordagens é usar a definição da derivada com um pequeno passo $h$. Podemos, por exemplo, usar a fórmula

$$ f'(x) \approx \frac{f(x+h) - f(x-h)}{2h} $$

que é um método eficiente e amplamente utilizado. Vamos detalhar essa e outras técnicas para obter maior acurácia.

A Diferença Central de Segunda Ordem

A expressão $f'(x) \approx \frac{f(x+h)-f(x-h)}{2h}$ é conhecida como a **diferença central de segunda ordem**. Ela é chamada assim porque utiliza pontos simetricamente posicionados ($x+h$ e $x-h$) em torno do ponto $x$.

$$ f'(x) \approx \frac{f(x+h) - f(x-h)}{2h} $$

Esta fórmula tem a vantagem de ser exata para polinômios de até grau 2. O principal termo de erro associado a esta aproximação é proporcional a $h^2$ (o erro é $O(h^2)$). Isso significa que, se você reduzir o passo $h$ pela metade, o erro de truncamento será dividido por 4.

Outras Expressões para Maior Acurácia

Existem várias maneiras de melhorar a precisão do cálculo numérico da derivada, geralmente utilizando mais pontos de avaliação ou combinando resultados de diferentes passos.

1. Diferença Central de Quarta Ordem

Para obter uma acurácia ainda maior, podemos usar uma fórmula que emprega mais pontos. A diferença central de quarta ordem é exata para polinômios de até grau 4. O termo de erro principal aqui é proporcional a $h^4$.

$$ f'(x) \approx \frac{1}{12h} \left[ -f(x+2h) + 8f(x+h) - 8f(x-h) + f(x-2h) \right] $$

Esta fórmula utiliza 5 pontos: $f(x-2h), f(x-h), f(x), f(x+h), f(x+2h)$. A redução do erro é significativa quando comparada com a fórmula de segunda ordem, especialmente para funções suaves.

2. Diferença Central de Sexta Ordem

A Diferença Central de Sexta Ordem é um método numérico avançado para aproximar a primeira derivada de uma função $f(x)$. Ela se destaca pela sua alta precisão, sendo exata para polinômios de até grau 5 e com um erro de truncamento da ordem de $O(h^6)$.

Para aplicar esta fórmula, são necessários 6 pontos da função, distribuídos simetricamente em torno do ponto $x$.

A expressão para a primeira derivada $f'(x)$ utilizando a Diferença Central de Sexta Ordem é a seguinte:

$$ f'(x) \approx \frac{1}{60h} \left[ -f(x+3h) + 9f(x+2h) - 45f(x+h) + 45f(x-h) - 9f(x-2h) + f(x-3h) \right] $$

Alternativamente, a fórmula pode ser escrita de forma mais simétrica, agrupando os termos que envolvem pontos em direções opostas:

$$ f'(x) \approx \frac{1}{60h} \left[ (f(x+3h) - f(x-3h)) - 9(f(x+2h) - f(x-2h)) + 45(f(x+h) - f(x-h)) \right] $$

Componentes da Fórmula

• $f(x)$: A função cujos pontos são avaliados.
• $x$: O ponto onde a derivada está sendo calculada.
• $h$: Um pequeno passo, positivo e real. O tamanho de $h$ afeta tanto o erro de truncamento quanto o erro de arredondamento.
• $f(x+kh)$ e $f(x-kh)$: Avaliações da função em pontos adiante e atrás de $x$, com múltiplos $k$ do passo $h$.
• Os coeficientes (1, -9, 45 no termo agrupado simetricamente) são determinados para maximizar a ordem de acurácia.

Esta fórmula é particularmente útil quando se busca alta precisão em cálculos de derivadas, desde que a função seja suave o suficiente e os pontos necessários estejam disponíveis para avaliação.

3. Fórmulas de Maior Ordem (Diferenças Diretas/Retrocesso)

Em situações onde não é possível usar pontos de ambos os lados de $x$ (por exemplo, nos limites de um intervalo), fórmulas diretas (forward differences) ou de retrocesso (backward differences) de ordem superior podem ser empregadas. Embora geralmente menos precisas que as diferenças centrais para o mesmo número de pontos, elas são úteis em contextos específicos.

Diferença Direta de Terceira Ordem (usando 4 pontos):

$$ f'(x) \approx \frac{1}{6h} \left[ -11f(x) + 18f(x+h) - 9f(x+2h) + 2f(x+3h) \right] $$ O erro é $O(h^3)$.

Diferença de Retrocesso de Terceira Ordem (usando 4 pontos):

$$ f'(x) \approx \frac{1}{6h} \left[ -2f(x-3h) + 9f(x-2h) - 18f(x-h) + 11f(x) \right] $$ O erro é $O(h^3)$.

4. Extrapolação de Richardson

A Extrapolação de Richardson é uma técnica poderosa para melhorar a acurácia de métodos de aproximação numérica. A ideia é calcular a aproximação da derivada com um passo $h$ e, em seguida, com um passo menor (por exemplo, $h/2$), e combinar esses dois resultados para "eliminar" os termos de erro de ordem inferior.

Para a diferença central de segunda ordem, que tem erro $O(h^2)$, podemos obter uma aproximação de quarta ordem. Seja $D(h) = \frac{f(x+h)-f(x-h)}{2h}$.

Uma aproximação melhorada é dada por:

$$ f'(x) \approx \frac{4 \cdot D(h/2) - D(h)}{3} $$

Esta fórmula, obtida pela combinação de duas aproximações com passos diferentes, resulta em uma precisão $O(h^4)$, similar à da diferença central de quarta ordem, mas com a conveniência de usar apenas a fórmula de segunda ordem como base. Podemos obter um valor ainda melhor para a derivada $f'(x)$ usando a aproximação de quarta ordem ($D(h) \approx \frac{1}{12h}(-f(x+2h) + 8f(x+h) - 8f(x-h) + f(x-2h)) $) e a Extrapolação de Richardson. A expressão se torna:

$$ f'(x) \approx \frac{16 \cdot D(h/2) - D(h)}{15} $$

Qual Método Escolher?

A escolha do método depende dos requisitos de acurácia e das restrições computacionais:

• Para alta acurácia em funções suaves, a **diferença central de quarta ordem** ou a **extrapolação de Richardson** (partindo da diferença central de segunda ordem) são geralmente as melhores opções.
• Se você precisa calcular a derivada perto das bordas de um domínio onde pontos simétricos não estão disponíveis, fórmulas diretas ou de retrocesso de maior ordem podem ser necessárias.

É importante notar que, independentemente do método, a escolha do passo $h$ é crucial. Um $h$ muito grande resulta em um erro de truncamento significativo, enquanto um $h$ muito pequeno pode levar a erros de arredondamento devido à subtração de números muito próximos. Geralmente, um $h$ entre $10^{-3}$ e $10^{-6}$ é um bom ponto de partida para muitas aplicações.