Método do Ponto Fixo para Sistemas de 3 Equações e 3 Variáveis
O método do ponto fixo é uma técnica iterativa utilizada para resolver sistemas de equações não lineares. O método pode ser eficaz para encontrar soluções aproximadas.
Formulação do Problema
Considere o sistema de equações:
f₁(x, y, z) = 0 f₂(x, y, z) = 0 f₃(x, y, z) = 0
Para aplicar o método do ponto fixo, precisamos reescrever o sistema na forma:
x = g₁(x, y, z) y = g₂(x, y, z) z = g₃(x, y, z)
Algoritmo do Método do Ponto Fixo
- Escolha uma aproximação inicial (x₀, y₀, z₀)
- Para k = 0, 1, 2, …, até convergir, calcule:
- xk+1 = g₁(xk, yk, zk)
- yk+1 = g₂(xk, yk, zk)
- zk+1 = g₃(xk, yk, zk)
- Pare quando |xk+1 – xk| < ε, |yk+1 – yk| < ε e |zk+1 – zk| < ε
Condições de Convergência
O método do ponto fixo converge se as funções g₁, g₂ e g₃ forem contrativas em uma vizinhança da solução. Uma condição suficiente é que a norma da matriz Jacobiana das funções g seja menor que 1 nessa vizinhança.
Exemplo Numérico
Considere o sistema:
x = 0.5 cos(y) + 0.5 sin(z) y = 0.5 sin(x) + 0.5 cos(z) z = 0.5 cos(x) + 0.5 sin(y)
Com aproximação inicial (0.5, 0.5, 0.5) e ε = 0.0001, o método produzirá uma sequência convergente para a solução.
Implementação em Python
def fixed_point_3d(g1, g2, g3, x0, y0, z0, tol=1e-6, max_iter=100):
x, y, z = x0, y0, z0
for i in range(max_iter):
x_new = g1(x, y, z)
y_new = g2(x, y, z)
z_new = g3(x, y, z)
if abs(x_new - x) < tol and abs(y_new - y) < tol and abs(z_new - z) < tol:
return x_new, y_new, z_new
x, y, z = x_new, y_new, z_new
return x, y, zVantagens e Limitações
- Vantagens: Simplicidade conceitual, fácil implementação
- Limitações: Convergência nem sempre garantida, pode ser lento
O método do ponto fixo é útil quando as funções g podem ser expressas de forma que a convergência seja assegurada, sendo uma alternativa aos métodos de Newton para sistemas não lineares.
Veja e teste o algoritmo em JSou leia mais sobre este assunto clicando aqui.