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

Enunciado

Resolva o sistema linear:

5 x + y + z =   5 3 x + 4 y + z =       6   3 x + 3 y + 6 z =     0

Pelo método de Gauss-Seidel com x 0 y 0 z 0 = 0 0 0 e ε = 5 × 10 - 2 .

Obs.: O critério de Sassenfeld é satisfeito!!

Passo 1

Fala Aí!!! O objetivo do método de Gauss-Seidel é encontrar uma aproximação desejada para a solução de sistemas lineares. A primeira coisa que precisaríamos fazer aqui é verificar se o critério de Sassenfeld é satisfeito, mas como o enunciado me garante isso, vamos pular essa parte, beleza?!

Então o próximo passo é encontrar as funções de iterações pra fazermos as continhas!

Passo 2

Podemos reescrever o caso do enunciado na forma de matriz para nos ajudar a encontrar as funções de iteração. Ficaria assim ó,

5 x + y + z =   5 3 x + 4 y + z =       6   3 x + 3 y + 6 z =     0

Assim, para encontrar as funções de iteração, vamos isolar o x 1 na primeira equação, o x 2 na segunda equação e o x 3 na terceira equação! (Aqui é como se o x 1 = x , x 2 = y e z = x 3 , beleza?)

Portanto, as funções de iteração são:

x k + 1 = 1 - 0.2 y k - 0.2 z k

y k + 1 = 1.5 - 0.75 x k + 1 - 0.25 z k

z k + 1 = 0 - 0.5 x k + 1 - 0.5 y k + 1

Como o critério de Sassenfeld é satisfeito podemos usar qualquer chute inicial que quisermos que o método vai convergir! Quanto melhor o chute, mais rápido vamos convergir, mas podemos escolher qualquer vetor como aproximação inicial. Aqui vamos usar como aproximação inicial x 0 = 0 ,   0 ,   0 .

Show de bola?!

Passo 3

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

x 1 = 1 - 0.2 * 0 - 0.2 * 0 = 1

y 1 = 1.5 - 0.75 * 1 - 0.25 * 0 = 0.75

z 1 = 0 - 0.5 * 1 - 0.5 * 0.75 = - 0.875

x 1 y 1 z 1 = 1 0.75 - 0.875

Será que já podemos parar?

Critério de Parada:

| x 1 - x 0 | = 1 - 0 = 1 > ε

Falhou logo no primeiro. Então precisamos de mais iterações.

Passo 4

Iteração 2:

x 2 = 1 - 0.2 * 0.75 - 0.2 * - 0.875 = 1.025

y 2 = 1.5 - 0.75 * 1.025 - 0.25 * - 0.875 = 0.95

z 2 = 0 - 0.5 * 1.025 - 0.5 * 0.95 = - 0.9875

x 2 y 2 z 2 = 1.025 0.95 - 0.9875

Será que já podemos parar?

Critério de Parada:

| x 2 - x 1 | = | 1.025 - 1 | = 0.025 < ε

| y 2 - y 1 | = | 0.95 - 0.75 | = 0.2 > ε

Bora fazer mais iterações então!

Passo 5

Iteração 3:

x 3 = 1 - 0.2 * 0.95 - 0.2 * - 0.9875 = 1.0075

y 3 = 1.5 - 0.75 * 1.0075 - 0.25 * - 0.9875 = 0.9913

z 3 = 0 - 0.5 * 1.0075 - 0.5 * 0.9913 = - 0.9994

x 3 y 3 z 3 = 1.0075 0.9913 - 0.9994

Será que já podemos parar?

Critério de Parada:

| x 3 - x 2 | = | 1.0075 - 1.025 | = 0.0175 < ε

| y 3 - y 2 | = | 0.9913 - 0.95 | = 0.0413 < ε

| z 3 - z 2 | = | - 0.9994 + 0.9875 | = 0.0119 < ε

Alô, Alô Graças a DEUS!!! Já pode parar!!!

Resposta

x = 1.0075 ,   0.9913 , - 0.9994