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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

CAPÍTULO 3. EQUAÇÕES NÃO LINEARES 68<br />

3.3 Método de Newton<br />

O método de Newton é uma das técnicas mais populares para se determinar raízes de equações não<br />

lineares. Existem várias maneiras de deduzir o método de Newton, a que apresentaremos aqui é baseada<br />

no método de iteração linear. Assim, para descrever tal método, consideremos a equação (3.2), isto é:<br />

ψ(x) = x + A(x)f(x) , com f ′ (x) = 0 , (3.6)<br />

onde a função A(x) deve ser escolhida de tal forma que A(¯x) = 0.<br />

Vimos pelo Teorema 3.4 que temos garantia de convergência se max|ψ ′ (x)| < 1 para x ∈ I. Assim se<br />

escolhermos A(x) tal que ψ ′ (¯x) = 0, teremos que para x ∈ I, (I suficientemente pequeno), max |ψ ′ (x)| < 1,<br />

garantindo então a convergência do método.<br />

Derivando (3.6) em relação a x, obtemos:<br />

Fazendo x = ¯x, segue que:<br />

e colocando:<br />

ψ ′ (x) = 1 + A ′ (x)f(x) + A(x)f ′ (x).<br />

ψ ′ (¯x) = 1 + A(¯x)f ′ (¯x) , pois f(¯x) = 0 ,<br />

ψ ′ (¯x) = 0, teremos A(¯x) = − 1<br />

f ′ (¯x) = 0 desde que f ′ (¯x) = 0 .<br />

Tomando então: A(x) = − 1<br />

f ′ ,<br />

(x)<br />

obtemos ψ(x) = x −<br />

f(x)<br />

f ′ definido por:<br />

(x)<br />

xk+1 = xk −<br />

f (xk)<br />

f ′ (xk)<br />

e o processo iterativo então<br />

é chamado Método de Newton, que converge sempre que |x0 − ¯x| for suficientemente pequeno.<br />

Uma interpretação geométrica do método de Newton é dada na Figura 3.9.<br />

f(x) ✻<br />

¯x xk+1<br />

Figura 3.9<br />

Dado xk, o valor xk+1 pode ser obtido graficamente traçando-se pelo ponto (xk, f(xk)) a tangente à curva<br />

y = f(x). O ponto de interseção da tangente com o eixo dos x determina xk+1.<br />

α<br />

xk<br />

f(xk)<br />

✲<br />

x<br />

(3.7)

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

Saved successfully!

Ooh no, something went wrong!