15.05.2013 Views

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

CAPÍTULO 5. SOLUÇÃO DE SISTEMAS LINEARES: MÉTODOS ITERATIVOS 181<br />

Agora, desde que:<br />

Ap (3) =<br />

⎛<br />

⎝ −0.3532<br />

0.3685<br />

0.0398<br />

⎞<br />

⎠ ,<br />

(Ap (3) , p (3) ) = 0.0139 + 0.0150 − 0.0000 = 0.0289 ,<br />

q3 = (r(2) , r (2)<br />

(Ap (3) , p (3) )<br />

v (k+1) − v (k) ∞<br />

v (k+1) ∞<br />

= 0.0033<br />

0.0289<br />

v (3) = v (2) + q3p (3)<br />

=<br />

⎛<br />

⎝ 0.9999<br />

1.0001<br />

0.0008<br />

temos que v3 é solução do sistema dado, com ɛ < 10 −2 .<br />

Exercícios<br />

⎞<br />

⎠ .<br />

= 0.1142 ,<br />

= 0.0047<br />

1.0001 0.0047 < 10−2 ,<br />

5.7 - Mostre que no método dos gradientes conjugados o resíduo em cada passo é ortogonal ao resíduo<br />

anterior, à direção de relaxação do passo e à direção de relaxação do passo anterior.<br />

5.8 - Usando o método dos gradientes conjugados resolver os sistemas dados nos exercícios 5.5 e 5.6.<br />

5.4 Exercícios Complementares<br />

5.9 - Supomos que o sistema:<br />

⎧<br />

⎨<br />

seja resolvido iterativamente pelas fórmulas:<br />

⎩<br />

x1 − αx2 = c1<br />

−αx1 + x2 − αx3 = c2<br />

− αx2 + x3 = c3<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

x (k+1)<br />

1<br />

x (k+1)<br />

2<br />

x (k+1)<br />

3<br />

= αx (k)<br />

2<br />

= α(x (k)<br />

1<br />

= αx (k)<br />

2<br />

+ c1<br />

+ x(k) 3 ) + c2<br />

Para que valores de α a convergência do método definido acima é garantida? Justifique.<br />

+ c3

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!