Exercícios Resolvidos de Método de Gauss-Jacobi

Enunciado

Considere o sistema:

x - 3 y + z = 5     2 y + 5 z = 9 2.1 x + y + z = 3

Calcule, pelo método de Gauss-Jacobi uma primeira aproximação X 1 , partindo de X 0 = 0,0 , 0 t e usando o sistema de forma tal que a convergência do método de Gauss-Jacobi esteja garantida.

Passo 1

Opaaa!! O objetivo do método de Gauss-Jacobi é encontrar uma aproximação desejada para a solução de sistemas lineares. Porém, para poder aplicar esse método, precisamos primeiro falar sobre convergência, pois só poderemos aplicar esse método se o sistema linear convergir.

Show de bola?!

Tá... e como funciona esse critério de convergência?

Uma condição suficiente para a convergência do método de Gauss-Jacobi é o critério das linhas. Nesse critério precisamos calcular um α para cada linha da matriz. Se o maior desses α n for menor do que 1, então o nosso critério foi satisfeito e o método irá convergir!!

Vamos fazer as continhas pra ficar mais claro!

Aqui como o nosso sistema tem 3 equações, calcularemos α 1 ,   α 2 e α 3 :

α 1 = a 12 + a 13 a 11

α 2 = a 21 + a 23 a 22

e

α 3 = a 31 + a 32 a 33

Beleza?! MAS PERA! Quem são esses a 12 , a 23 ...?

Podemos reescrever o nosso sistema na forma de matriz, olha só:

1 - 3 1 0 2 5 2.1 1 1 x y z = 5 9 3

Os elementos da matriz dos coeficientes serão os elementos

a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33

Que usaremos para calcular cada α

Show de bola?!

Passo 2

Já que para que a convergência do método esteja garantida precisamos atender ao critério das linhas, podemos trocar as linhas de forma que os maiores elementos de cada linha fiquem na diagonal principal. Assim, eles estarão no denominador de cada α e a chance deles serem menores do que 1 aumenta!!

Assim, podemos fazer as trocas L 1 ↔ L 2 :

    2 y + 5 z = 9 x - 3 y + z = 5 2.1 x + y + z = 3

E em seguida a troca L 1 ↔ L 3 :

2.1 x + y + z = 3 x - 3 y + z = 5     2 y + 5 z = 9

Fizemos isso para que

Dessa forma temos:

α 1 = 2 2.1 < 1

α 2 = 2 3 < 1

α 3 = 2 5 < 1

Então o critério das linhas é atendido e a convergência é garantida. Show de bola!? Agora vamos achar as funções de iteração e fazer as continhas!!

Passo 3

Para encontrar as funções de iteração vamos isolar o x na primeira equação, o y na segunda equação e o z na terceira equação. Assim, as funções de iteração são:

x k + 1 = 3 - y k - z k 2.1

y k + 1 = - 5 + x k + z k 3

z k + 1 = 9 - 2 y k 5

E temos do enunciado que X 0 = x 0 , y 0 , z 0 = 0,0 , 0

Passo 4

Iteração 1 : Vamos substituir os valores de X 0 = x 0 , y 0 , z 0 = 0,0 , 0 nas funções de iteração para achar X 1 :

x 1 = 3 - y 0 - z 0 2.1

x 1 = 3 2.1 = 1.429

y 1 = - 5 + x 0 + z 0 3

y 1 = - 5 3 = - 1.667

z 1 = 9 - 2 y 0 5

z 1 = 9 5 = 1.8

Pronto, obtemos nossa aproximação X 1 = 1.429 , - 1.667,1.8 .

Resposta

X 1 = 1.429 , - 1.667,1.8 .