A Aproximação da Reta 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 em todos os lugares segue essas inclinações atribuídas.
O método de Euler converte o campo de direções contínuo em um passeio discreto:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
Do ponto (xₙ, yₙ), mova uma distância h ao longo da reta tangente. Chegue a um próximo ponto aproximado. Repita.
Erro geométrico: a reta tangente em (xₙ, yₙ) tem inclinação f(xₙ, yₙ), mas a curva verdadeira tem uma inclinação diferente em cada ponto ao longo do intervalo [xₙ, xₙ + h]. O passo de Euler usa a inclinação no ponto final esquerdo em todo o intervalo — 'a inclinação que era.' O erro por passo cresce como h².
Erro Acumulado
Ao longo de N passos para alcançar um ponto final fixo x = a, com h = a/N:
- Erro de truncamento local 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: reduzir h pela metade reduz o erro global pela metade.
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.
Aplique o método de Euler com h = 0.5, de x = 0 a 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 vs valor verdadeiro 2.71828. Erro: 0.468. Erro relativo: ~17%.
Derivando a Região de Estabilidade de Euler
Para a equação teste dy/dx = λy (onde λ é um número complexo), o método de Euler dá:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
O fator de amplificação por passo: z = 1 + hλ.
Condição de estabilidade: a solução calculada permanece limitada se & somente se |z| ≤ 1, ou seja, |1 + hλ| ≤ 1.
Esta é 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 λ real & negativo (uma equação diferencial decadente como dy/dx = -2y): hλ deve estar no intervalo (-2, 0] no eixo real. Com λ = -2 & h = 0.5: hλ = -1. Isso está exatamente na fronteira de estabilidade — o método é marginalmente estável, o que explica a falha qualitativa no exemplo anterior.
Com h = 1 & λ = -2: hλ = -2, nos colocando 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 que Euler, o que é uma razão pela qual é preferida para a maioria dos problemas.
Para λ real negativo, RK4 permite hλ até aproximadamente -2.785 no eixo real (vs limite de -2 para Euler).
Para equações stiff com autovalores λ em magnitudes muito diferentes — digamos λ₁ = -1 & λ₂ = -1000 — a estabilidade requer que hλ₂ permaneça dentro da região. Para RK4 no eixo real: h·(-1000) ≥ -2.785, então h ≤ 0.002785.
Este tamanho de passo minúsculo, ditado pelo autovalor stiff λ₂, 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 um mapa discreto: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
Um ponto fixo deste mapa: y tal que g(y) = y. Para Euler em dy/dx = f(y), pontos fixos satisfazem f(y) = 0 — equilíbrios da equação diferencial.
Estabilidade de um ponto fixo: se |g'(y)| < 1, 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.
Esta é exatamente a condição de estabilidade de Euler com λ = f'(y*) — a linearização da equação diferencial 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 rastreará com confiabilidade o equilíbrio da equação diferencial vs divergir para outro atrator.
O loop de simulação é um sistema dinâmico discreto. Seu comportamento qualitativo — convergência, oscilação, divergência — depende do tamanho do passo h relativo à geometria do campo de direções da equação diferencial.
Conectando Geometria ao Design de Simulação
A geometria da simulação numérica se reduz a três questões:
1. Onde está a região de estabilidade? Para Euler: o disco |1 + hλ| ≤ 1. Maior para RK4. Ilimitado (meio-plano esquerdo inteiro) para métodos implícitos.
2. Onde estão os autovalores da equação diferencial? Os autovalores λ do Jacobiano de f em cada ponto determinam qual região de estabilidade deve conter hλ.
3. Qual h mantém todo hλ dentro da região? O h máximo permitido = (raio da região de estabilidade) / max|λ|.
Para sistemas stiff: max|λ| é enorme, forçando h pequeno para métodos explícitos. Métodos implícitos são caros por passo, mas permitem h grande — eles trocam custo por passo por estabilidade.
A percepção de Hamming se traduz: a escolha do método numérico codifica uma aposta sobre a geometria do espectro de autovalores da equação diferencial. Tornar essa aposta explícita é a primeira decisão de design em qualquer simulação.