Cálculo Numérico
Resolução Numérica de Equações
Prof.: Sandrerley Ramos Pires
Autores: Bruno Correia da Nóbrega Queiroz
José Eustáquio Rangel de Queiroz
Marcelo Alves de Barros
Cálculo Numérico – Métodos
Métodos Iterativos para a Obtenção de
Zeros Reais de Funções
Bissecção
Falsa Posição
Falsa Posição Modificado
Ponto Fixo
Newton-Raphson
Secante
2
Cálculo Numérico – Newton-Raphson
Método de Newton-Raphson
Dada uma função f(x) contínua no intervalo [a,b]
onde existe uma raiz única, é possível determinar
uma aproximação de tal raiz a partir da interseção
da tangente à curva em um ponto x0 com o eixo
das abscissas.
x0 - atribuído em função da geometria do método e
do comportamento da curva da equação nas
proximidades da raiz.
3
Cálculo Numérico – Newton-Raphson
f(x)
Análise Gráfica
1a iteração
2a iteração
3a iteração
4a iteração
x2
x0 x3 x1
x
Repete-se o processo até que o
valor de x atenda às condições
de parada.
4
Cálculo Numérico – Newton-Raphson
Motivação Geométrica:
Dado o ponto (xk , f(xk))
Traça-se a reta Lk(x) tangente à curva
neste ponto:
Lk(x) = f(xk) + f’(xk)(x-xk)
Determina-se o zero de Lk(x), um modelo
linear que aproxima f(x) em uma
vizinhança xk
Lk(x) = 0 x = xk - f(xk)/f’(xk)
Faz-se xk +1 = x
5
Cálculo Numérico – Newton-Raphson
Considerações Iniciais
Deste modo, escolhido x0 , a seqüência
{xk} será determinada por
f (x k )
x k 1 x k ,
f ( x k )
onde k = 0, 1, 2, ...
6
Cálculo Numérico – Newton-Raphson
Estudo da Convergência
TEOREMA 3:
Sendo f(x), f’(x) e f”(x) contínuas em um
intervalo I que contém uma raiz x = de f(x) e
supondo f’() 0, existirá um intervalo Ī I
contendo a raiz , tal que se x0 I e x0 x1 Ī,
então a seqüência {xk} gerada pela fórmula
recursiva
f (x k )
x k 1 x k
f ( x k )
convergirá para a raiz. 7
Cálculo Numérico – Newton-Raphson
Testes de Parada
A cada iteração, testa-se se a
aproximação encontrada poderá ser
considerada como a solução do problema.
|f(xk)| tolerância
|((xk+1 – xk)/xk+1 )| tolerância
8
Cálculo Numérico – Newton-Raphson
Algoritmo
k := 0; x0 := x;
while critério de interrupção não satisfeito and k L
k := k +1;
xk+1 := xk – f(xk)/f’(xk)
endwhile
9
Cálculo Numérico – Newton-Raphson
Exemplo 1: Dado a Equação
x2 + x – 6 = 0 :
Seja a raiz 2 = 2 e x0 = 1,5
Assim:
g(x) = x - f(x)/f’(x) = x – (x 2 + x – 6)/(2x + 1)
x1 = g(x0) = 1,5 – (1,52 + 1,5 – 6)/(2.1,5 + 1)
x1 = 2,062500000
x2 = g(x1) = 2,000762195
x3 = g(x2) = 2,000000116
10
Cálculo Numérico – Newton-Raphson
Comentários:
A parada poderá ocorrer na 3a iteração
(x=2,000000116), caso a precisão do cálculo
com 6 casas decimais for satisfatória para o
contexto do trabalho
Utilizando, por exemplo, o Método do Ponto
Fixo com g(x)=√6 - x só veio a produzir
x = 2,000476818 na 5a iteração
11
Cálculo Numérico – Newton-Raphson
Exemplo 2 - Considere-se a função:
f(x)=x3 - x - 1 , e tol = 0,002
cujos zeros encontram-se nos intervalos:
1 I1 = (-1, 0), 2 I2 = (1, 2)
Seja x0 = 1
xk+1 = xk - f(xk)/f’(xk)
e g(x) = x – (x3 - x - 1)/(3x2 – 1)
12
Cálculo Numérico – Newton-Raphson
Exemplo 2:
Cálculo da 1ª aproximação
g(x0) = 1 – [ (1)³ – 1 – 1 ] = 1,5
[ 3*(1)² – 1 ]
Teste de Parada
|f(x0)| =| 0,875 | = 0,875 >
13
Cálculo Numérico – Newton-Raphson
Exemplo 2:
Cálculo da 2ª aproximação
g(x1) = 1.5 – [ (1.5)³ – 1.5 – 1 ] = 1,3478261
[ 3*(1.5)² – 1 ]
Teste de Parada
|f(x1)| =| 0,100682 | = 0,100682 >
14
Cálculo Numérico – Newton-Raphson
Exemplo 2:
Cálculo da 3ª aproximação
g(x2) = 1,3478261 - [ (1,3478261)³ - 1,3478261 - 1 ]
[ 3*(1,3478261)² - 1 ]
g(x2) = 1,3252004
Teste de Parada
|f(x2)| =| 0,0020584 | = 0,0020584 >
15
Cálculo Numérico – Newton-Raphson
Exemplo 2:
A seqüência {xk} gerada pelo método de
Newton será:
Iteração x F(x)
1 1,5 0,875
2 1,3478261 0,1006822
3 1,3252004 0,0020584
4 1,3247182 9,24378.10 -7
5 1,3247178 1,86517.10-13
= 0,002
16
Cálculo Numérico – Newton-Raphson
Vantagens:
Rapidez processo de convergência;
Desempenho elevado.
17
Cálculo Numérico – Newton-Raphson
Desvantagens:
Necessidade da obtenção de f’(x) , o que
pode ser impossível em determinados casos;
O cálculo do valor numérico de f’(x) a cada
iteração;
Difícil implementação.
18
