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 74<br />

Ordem de Convergência<br />

Daremos aqui a ordem de convergência do método das secantes.<br />

Teorema 3.7 - A ordem de convergência do método das secantes é p = (1 + √ 5)/2 1.618.<br />

Prova: A prova deste teorema pode ser encontrada em [Ostrowski,1966].<br />

Observe que apesar da ordem de convergência do método das secantes ser inferior à do método de<br />

Newton, ele fornece uma alternativa viável, desde que requer somente um cálculo da função f por passo,<br />

enquanto que dois cálculos (f(xk) e f ′ (xk)) são necessários para o método de Newton.<br />

Exercícios<br />

3.12 - Considere a equação dada no exemplo 3.10. Obtenha a raiz positiva com quatro casas decimais<br />

corretas. Usando (3.5) confirme que a ordem de convergência do método das secantes é p 1.618.<br />

3.13 - Determinar, pelo método das secantes, uma raiz de cada uma das equações:<br />

a) x = 2.7 ln x,<br />

b) log x − cos x = 0,<br />

c) e −x − log x = 0.<br />

3.14 - A equação x = tgx tem uma raiz entre π 2 e 3π 2<br />

. Determiná-la pelo método das secantes com<br />

erro inferior a 10−3 .<br />

3.5 Método Regula Falsi<br />

O método Regula Falsi, é uma variação do método das secantes. Ele consiste em tomar duas<br />

aproximações iniciais x0 e x1 tais que f(x0) e f(x1) tenham sinais opostos, isto é:<br />

f (x0) × f (x1) < 0.<br />

Uma nova aproximação é determinada usando o método das secantes, ou seja:<br />

Se <br />

x2 = x0 f (x1) − x1 f (x0)<br />

.<br />

f (x1) − f (x0)<br />

x2 − x0<br />

x2<br />

<br />

<br />

<br />

< ɛ ou <br />

<br />

x2 − x1<br />

x2<br />

<br />

<br />

<br />

< ɛ ,<br />

para um ɛ pré-fixado, então x2 é a raiz procurada. Caso contrário, calculamos f(x2) e escolhemos entre x0<br />

e x1 aquele cuja f tenha sinal oposto ao da f(x2). Com x2 e esse ponto calculamos x3 usando a fórmula<br />

das secantes, isto é, usando (3.9) e assim sucessivamente. O processo iterativo deve ser continuado até<br />

que se obtenha a raiz com a precisão pré-fixada.<br />

Uma interpretação geométrica do método regula falsi é dada na Figura 3.13. Observe que, na Figura<br />

3.13, xk+1 é o ponto de interseção da corda que une os pontos (xk−1, f(xk−1)) e (xk, f(xk)) com o eixo<br />

dos x. Nesse caso o novo intervalo contendo a raiz será (xk−1, xk+1). A aproximação xk+2 será o ponto<br />

de interseção da corda que une os pontos (xk−1, f(xk−1)) e (xk+1, f(xk+1)) com o eixo dos x. Observe<br />

ainda que a aplicação do método regula falsi sempre mantém a raiz procurada entre as aproximações<br />

mais recentes.

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

Saved successfully!

Ooh no, something went wrong!