Manual - MN - 2015 - Cap V PDF

5.
RESOLUO NUMRICA DE SISTEMAS DE EQUAES
Estima-se que em 75% dos problemas cientficos, a soluo de um sistema linear de

equaes aparece em algum estgio da soluo. Podemos, entre outros, citar os
seguintes problemas:
Interpolao, ajuste de curvas, soluo de sistemas de equaes no lineares, soluo
de equaes diferenciais usando diferenas finitas e clculo de autovalores e
autovetores.
Um sistema de equaes lineares, s vezes simplesmente dito sistema linear, uma
equao do tipo Ax = b . Usando os valores de A , x e b temos:
Um sistema linear pode tambm ser representado na seguinte forma:
_ _ _ _
Resolver o sistema de equaes, significa determinar um vetor x = ( x1 ; x 2 ;....; x n ) que
satisfaa todas equaes do sistema.

Existem duas classes de mtodos para soluo de sistemas lineares. Os denominados
mtodos directos e os denominados mtodos iterativos. A escolha de uma ou outra
classe para a soluo do sistema Ax = b vai depender de caractersticas e
propriedades da matriz A .
5.1 MTODOS DIRECTOS
Os mtodos diretos so aqueles em que aps um nmero finito de etapas e no se

considerando os erros de arredondamento encontramos a soluo do sistema.
Um sistema linear da forma
onde a ii 0 para i = 1, 2, . . . , n dito um sistema triangular. Sua soluo dada por
retro substituio atravs da frmula de recorrncia.
bn
Onde x n = . Este algoritmo pode ser incorporado ao mtodo de eliminao de
a nn
Gauss para encontrar a soluo de um sistema.
5.1.1 MTODO DE ELIMINAO DE GAUSS

Este um esquema eficiente para resolver um sistema Ax = b de n equaes a n
incgnitas de pequeno tamanho (ordem 10). O passo crucial construir um sistema
triangular superior A' x = b' equivalente que pode ser resolvido por retro substituio.
Dois sistemas lineares de dimenso n x n so equivalentes desde que os seus conjuntos
de solues sejam os mesmos. Teoremas de lgebra linear mostram que quando certas
operaes so aplicadas a um dado sistema, os conjuntos solues no muda.
As seguintes operaes, quando aplicadas a um sistema linear, resultam num sistema
equivalente:
Mudana de ordem de duas equaes.
Multiplicao de uma equao por uma constante no nula.
Adio de um mltiplo de uma equao a uma outra equao.
Descreve-se, a seguir, o mtodo de eliminao de Gauss para um sistema de ordem 3,

sendo que o mesmo processo pode ser aplicado a sistemas de qualquer ordem.
Assim, considere um sistema do tipo
que pode ser representado pela matriz aumentada
O objectivo obter um sistema triangular da forma

Para melhor perceber o mtodo, a seguir apresenta-se um exemplo:
Considere o sistema dado por
Passo 1
Podemos eliminar a incgnita x1 nas linha l2 e na linha l3 fazendo repectivamente as

seguintes operaes:
teremos ento o seguinte sistema equivalente
Passo 2
Vamos agora eliminar a incgnita x2 na linha l3 fazendo a seguinte operao:
Teremos ento o seguinte sistema equivalente
Temos agora um sistema triangular que pode ser resolvido com o algoritmo da retro
substituio.
5.2 MTODOS ITERATIVOS
Nesta seco descreve-se os mtodos iterativos de Jacobi e Gauss-Seidel. Tais

mtodos so eficazes para sistemas de grande porte, economizando tempo
computacional, principalmente os esparsos, na qual os mtodos directos so
ineficientes. A tcnica iterativa para resolver um sistema linear Ax = b de nn
elementos, ser proposto da seguinte forma: parte-se de uma aproximao inicial X0 e
ento constri-se consecutivamente os vetores x1, x2, ..., at que a condio de
convergncia seja satisfeita.
5.2.1 MTODO DE JACOBI
Considere o sistema:
E exigindo que .
Passo 1: Isola-se o vector X mediante a separao do elemento diagonal, conforme:

Passo 2: Extrai-se a matriz G
Passo 3: Verifica-se a convergncia da transformao
G <1
Passo 4: Determina-se o vector inicial:
X 0 = ( x1( 0 ) ; x 2( 0) ;...; x n( 0) )
Caso o problema no indique de forma explcita, nesta disciplina iremos adoptar a

seguinte soluo inicial: X 0 = (0;0;...;0)
Passo 5: Preenche-se a tabela recorrendo ao mtodo de Jacobi que corresponde a

resolver
(k ) 1 ( k 1) ( k 1)
x1 = a (b1 a12 x 2 ... a1n x n )
11
(k ) 1
x2 = (b2 a 21 x1( k 1) .... a 2 n x n( k 1) )
a 22
....
(k ) 1 ( k 1) ( k 1)
x n = a (bn a n1 x1 ... a nn 1 x n 1 )
nn
At que o critrio de convergncia seja satisfeito.
k x1( k ) x 2( k ) . x n(k ) X ( k ) X ( k 1)
5.2.2 MTODO DE GAUSS-SEIDEL
Seja dado o seguinte sistema:
E exigindo que . O Mtodo de Gauss Seidel para

determinao da soluo aproximada do sistema dado, compreende os seguintes
passos:
Passo 1: Isola-se o vector X mediante a separao do elemento diagonal, conforme:

Passo 2: Extrai-se a matriz G
Passo 3: Verifica-se a convergncia da transformao
G <1
Passo 4: Determina-se o vector inicial:
X 0 = ( x1( 0 ) ; x 2( 0) ;...; x n( 0) )
Caso o problema no indique de forma explcita, nesta disciplina iremos adoptar a

seguinte soluo inicial: X 0 = (0;0;...;0)
Passo 5: Preenche-se a tabela recorrendo ao mtodo de Gauss-Seidel que

corresponde a resolver
(k ) 1 ( k 1) ( k 1)
x1 = a (b1 a12 x 2 ... a1n x n )
11
(k ) 1
x2 = (b2 a 21 x1( k ) .... a 2 n x n( k 1) )
a 22
....

(k ) 1
x n = a (bn a n1 x1 ... a nn 1 x n 1 )
(k ) (k )
nn
At que o critrio de convergncia seja satisfeito.
k x1( k ) x 2( k ) . x n(k ) X ( k ) X ( k 1)
Passo 6: R% X = ( x1 x1 ; x 2 x 2 ;...; x n x n )
Exemplo:
Pelo mtodo de Gauss-Seidel, determine a soluo do seguinte sistema de equaes:
EXERCCIOS
1. Resolva o sistema linear abaixo utilizando o mtodo de eliminao de Gauss

2. Calcule a norma de A para as seguintes matrizes:
3. Determine a inversa da matriz A
4. Determine a soluo do sistema Ax = b, sabendo que:
5. Dado o sistema Ax = b, mostre que o mtodo iterativo de Jacobi converge e

determine a sua soluo aproximada:
6. Resolver o sistema linear abaixo, pelo mtodo de Jacobi com chute inicial
x ( 0) = {1,1,1,1} , tolerncia = 10 3 :
7. Resolver o sistema linear abaixo, pelo mtodo de Gauss Seidel com chute inicial
x ( 0) = {0,0,0, 0}, tolerncia = 10 3 :
8. Considere o sistema linear Ax = b, onde:
3
Pelo mtodo de Gauss Seidel, determine a soluo do sistema com erro = 10

Manual - MN - 2015 - Cap V PDF

Enviado por

Direitos autorais:

Formatos disponíveis

Manual - MN - 2015 - Cap V PDF

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Manual - MN - 2015 - Cap V PDF

Enviado por

Direitos autorais:

Formatos disponíveis

5.

RESOLUO NUMRICA DE SISTEMAS DE EQUAES

Estima-se que em 75% dos problemas cientficos, a soluo de um sistema linear de

Um sistema linear pode tambm ser representado na seguinte forma:

satisfaa todas equaes do sistema.

5.1 MTODOS DIRECTOS

Os mtodos diretos so aqueles em que aps um nmero finito de etapas e no se

Um sistema linear da forma

onde a ii 0 para i = 1, 2, . . . , n dito um sistema triangular. Sua soluo dada por

retro substituio atravs da frmula de recorrncia.

5.1.1 MTODO DE ELIMINAO DE GAUSS

Descreve-se, a seguir, o mtodo de eliminao de Gauss para um sistema de ordem 3,

que pode ser representado pela matriz aumentada

O objectivo obter um sistema triangular da forma

Considere o sistema dado por

Podemos eliminar a incgnita x1 nas linha l2 e na linha l3 fazendo repectivamente as

teremos ento o seguinte sistema equivalente

Vamos agora eliminar a incgnita x2 na linha l3 fazendo a seguinte operao:

Teremos ento o seguinte sistema equivalente

Nesta seco descreve-se os mtodos iterativos de Jacobi e Gauss-Seidel. Tais

5.2.1 MTODO DE JACOBI

Passo 1: Isola-se o vector X mediante a separao do elemento diagonal, conforme:

Passo 3: Verifica-se a convergncia da transformao

Passo 4: Determina-se o vector inicial:

Caso o problema no indique de forma explcita, nesta disciplina iremos adoptar a

Passo 5: Preenche-se a tabela recorrendo ao mtodo de Jacobi que corresponde a

At que o critrio de convergncia seja satisfeito.

5.2.2 MTODO DE GAUSS-SEIDEL

Seja dado o seguinte sistema:

E exigindo que . O Mtodo de Gauss Seidel para

Passo 1: Isola-se o vector X mediante a separao do elemento diagonal, conforme:

Passo 3: Verifica-se a convergncia da transformao

Passo 4: Determina-se o vector inicial:

Caso o problema no indique de forma explcita, nesta disciplina iremos adoptar a

Passo 5: Preenche-se a tabela recorrendo ao mtodo de Gauss-Seidel que

At que o critrio de convergncia seja satisfeito.

Pelo mtodo de Gauss-Seidel, determine a soluo do seguinte sistema de equaes:

1. Resolva o sistema linear abaixo utilizando o mtodo de eliminao de Gauss

3. Determine a inversa da matriz A

4. Determine a soluo do sistema Ax = b, sabendo que:

5. Dado o sistema Ax = b, mostre que o mtodo iterativo de Jacobi converge e

8. Considere o sistema linear Ax = b, onde:

Você também pode gostar