Metodos Jacobi e Gauss
Calculo Numérico
Prof. Fumachi
Demerval dos Santos Izidoro RA. 6659369829
Jacarei - 2013
Diferenças entre os metodos de Jacobi e Gauss-Seidel
A forma mais simples de se determinar a matriz, a partir do sistema Ax = b
Seja A a matriz do sistema, da forma
A =
Vamos supor que A foi reordenada de modo que todos os seus elementos da diagonal sejam não-nulos:
.
Vamos então tirar o valor de cada xi na i-ésima equação (i = 1, 2, . . ., n). Como assumimos que aii é não nulo, podemos escrever:
Se considerarmos o lado esquerdo do sistema como os elementos de um novo passo de iteração (k+1) e os elementos do lado direito como elementos do passo anterior (k), teremos:
J =
Vamos resolver o sistema :
2.x1 + x2 = 5 x1+ 2.x2 = 4
Tiramos inicialmente o valor de x1 na primeira equação e de x2 na segunda equação: x1 = (5/2) - (1/2) x2 x1= 0.x1 - (1/2).x2 + 5
{ ou { (3.6) x2 = 2 - (1/2) x1 x2 = - (1/2).x1 + 0.x2 + 2
X = , J = , C =
Agora façamos o seguinte:
1. Chamamos de e as aproximações iniciais (arbitrárias, como vamos ver posteriormente) das componentes de X, ou seja, definimos um vetor :
2. Aplicamos do lado direito do sistema obtendo um novo valor para x1 e x2. Digamos que escolhemos = = 0; assim obtemos os valores:
3. Usamos estes valores novamente no sistema (3.6) obtendo os valores:
4. O próximo passo será:
5. Para os demais: etc... Como vemos, o valor das componentes de X(i) vão se aproximando da solução exata, x1 = 2 e x2 = 1, na medida em que vamos calculando novas iterações. Como já dissemos anteriormente, esse método é chamado Método Iterativo de Jacobi e a matriz J é a sua matriz de iteração.
Podemos, entretanto, introduzir uma variação na escolha dos índices (k) e (k+1), caracterizando um novo processo iterativo. Com o intuito de aproveitar os valores já encontrados em