Aproximação da Linha Tangente
A Ideia Geométrica
Uma equação diferencial ordinária dy/dx = f(x,y) atribui uma inclinação a cada ponto no plano (x,y) — um campo de direções. A solução verdadeira y(x) é uma curva que segue sempre aquelas inclinações atribuídas.
O método de Euler converte o campo de direção contínuo em uma caminhada discreta:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
De um ponto (xₙ, yₙ), mover uma distância h ao longo da linha tangente. Chegar em um próximo ponto aproximado. Repetir.
Erro geométrico: a linha tangente em (xₙ, yₙ) tem inclinação f(xₙ, yₙ), mas a curva verdadeira tem uma inclinação diferente em todos os pontos ao longo do intervalo [xₙ, xₙ + h]. O passo de Euler usa a inclinação no ponto extremo esquerdo em todo o tempo - 'a inclinação que foi'. O erro por passo cresce como h².
Erro Acumulado
Ao longo de N passos para atingir um ponto fixo x = a, com h = a/N:
- Erro local de truncamento por passo: O(h²)
- Número de passos: N = a/h
- Erro global: O(h²) × (a/h) = O(h) - precisão de primeira ordem
O método de Euler é de primeira ordem: dividir h por dois reduz o erro global por dois.
Executando o Método de Euler
Considere dy/dx = y, com condição inicial y(0) = 1. Solução verdadeira: y(x) = eˣ, então y(1) = e ≈ 2.71828.
Aplicar o método de Euler com h = 0.5, de x = 0 para x = 1 (2 passos):
Passo 1: y₁ = y₀ + h·f(x₀, y₀) = 1 + 0.5·(1) = 1.5. Novo ponto: (0.5, 1.5).
Passo 2: y₂ = y₁ + h·f(x₁, y₁) = 1.5 + 0.5·(1.5) = 2.25. Novo ponto: (1.0, 2.25).
Euler dá 2.25 contra o valor verdadeiro 2.71828. Erro: 0.468. Taxa de erro relativa: ~17%.
Derivando a Região de Estabilidade de Euler
Para a equação de teste dy/dx = λy (onde λ é um número complexo), o método de Euler fornece:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
O fator de amplificação por passo: z = 1 + hλ.
Condição de estabilidade: a solução computada permanece limitada se e somente se |z| ≤ 1, ou seja, |1 + hλ| ≤ 1.
Isso é uma condição geométrica no plano complexo hλ: o ponto hλ deve estar dentro do círculo de raio 1 centrado em (-1, 0).
Região de estabilidade de Euler: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
Para valores reais e negativos de λ (uma ODE decaindo como dy/dx = -2y): hλ deve estar no intervalo (-2, 0] no eixo real. Com λ = -2 e h = 0.5: hλ = -1. Isso está exatamente na fronteira de estabilidade - o método é marginalmente estável, o que explica o fracasso qualitativo no exemplo anterior.
Com h = 1 e λ = -2: hλ = -2, colocando-nos fora da região de estabilidade. A solução oscila com amplitude crescente.
Encontrando a Fronteira de Estabilidade
Runge-Kutta 4 (RK4) tem uma região de estabilidade maior do que o método de Euler, o que é uma das razões pelas quais é preferido para a maioria dos problemas.
Para valores negativos reais de λ, o RK4 permite hλ até aproximadamente -2.785 no eixo real (em comparação com o limite de Euler de -2).
Para equações rígidas com autovalores λ com magnitudes muito diferentes - digamos λ₁ = -1 e λ₂ = -1000 - a estabilidade exige que hλ₂ fique dentro da região. Para o RK4 no eixo real: h·(-1000) ≥ -2.785, então h ≤ 0.002785.
Este pequeno passo, ditado pelo autovalor rígido λ₂, torna a simulação cara, mesmo que o componente lento λ₁ pudesse usar h = 2.
Pontos Fixos & Bacias de Atração
O método de Euler aplicado a dy/dx = f(y) define uma mapa discreto: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
Um ponto fixo deste mapa: y tal que g(y) = y. Para o Euler em dy/dx = f(y), os pontos fixos satisfazem f(y) = 0 - equilíbrios da ODE.
Estabilidade de um ponto fixo: se |g'(y)| < 1, as iterações próximas convergem para y. Se |g'(y*)| > 1, elas divergem.
g'(y) = 1 + h·f'(y). Em um ponto fixo y: |1 + h·f'(y)| < 1 para estabilidade.
Isso é exatamente a condição de estabilidade de Euler com λ = f'(y*) - a linearização da ODE no equilíbrio.
Bacia de atração: o conjunto de condições iniciais que convergem para y* sob o mapa de Euler. Para sistemas não lineares, a fronteira da bacia define onde a simulação seguirá confiavelmente o equilíbrio da ODE ou divergirá para outro atrator.
O loop de simulação é um sistema dinâmico discreto. Sua comportamento qualitativo - convergência, oscilação, divergência - depende do tamanho do passo h em relação à geometria do campo de direções da ODE.
Conectando a Geometria ao Design de Simulação
A geometria da simulação numérica reduz-se a três questões:
1. Onde está a região de estabilidade? Para Euler: o disco |1 + hλ| ≤ 1. Maiores para RK4. Ilimitada (todo o semiplano esquerdo) para métodos implícitos.
2. Onde estão os autovalores das EDOs? Os autovalores λ da Jacobiano de f em cada ponto determinam qual a região de estabilidade que deve conter hλ.
3. Qual o h que mantém todos os hλ dentro da região? O h permitido máximo = (raio da região de estabilidade) / max|λ|.
Para sistemas rígidos: max|λ| é enorme, forçando h minúsculo para métodos explícitos. Métodos implícitos são caros por passo, mas permitem grandes h — eles trocam custo por passo por estabilidade.
A intuição de Hamming se traduz: a escolha do método numérico encoda uma aposta sobre a geometria do espectro de autovalores das EDOs. Tornar essa aposta explícita é a primeira decisão de design em qualquer simulação.