Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra...
Transcript of Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra...
![Page 1: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/1.jpg)
Met Comp
Metodos Computacionais em Fısica
Erica PolycarpoSandra Amato
Instituto de FısicaUniversidade Federal do Rio de Janeiro
Segundo Semestre de 2008
![Page 2: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/2.jpg)
Met Comp
Metodos Computacionais em Fısica
1 Determinacao de RaızesMetodo da BissecaoMetodo de Newton-Raphson
![Page 3: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/3.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Determinacao de Raızes de Funcoes
Em muitos problemas e necessario encontrar a solucao daequacao f (x) = 0 e sao muito frequentes os casos em quea solucao so pode ser obtida numericamente.
Exemplo: Encontre as raızes da seguinte equacaotranscendental:
f (x) = e−x − sin(πx2
)= 0.
![Page 4: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/4.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Determinacao de Raızes de Funcoes
Em muitos problemas e necessario encontrar a solucao daequacao f (x) = 0 e sao muito frequentes os casos em quea solucao so pode ser obtida numericamente.Exemplo: Encontre as raızes da seguinte equacaotranscendental:
f (x) = e−x − sin(πx2
)= 0.
![Page 5: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/5.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Determinacao de Raızes de Funcoes
Em muitos problemas e necessario encontrar a solucao daequacao f (x) = 0 e sao muito frequentes os casos em quea solucao so pode ser obtida numericamente.Exemplo: Encontre as raızes da seguinte equacaotranscendental:
f (x) = e−x − sin(πx2
)= 0.
![Page 6: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/6.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao
Ü Sabemos visualmente onde estao as raızes pelamudanca de sinal da funcao.
Ü Escolhemos um intervalo onde ha uma unica raizÜ A raiz esta entre entre x = a e x = b tais quef (a)f (b) < 0.Ü O ponto medio e uma primeira aproximacao para aposicao da raiz. Iterando ...
![Page 7: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/7.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao
Ü Sabemos visualmente onde estao as raızes pelamudanca de sinal da funcao.Ü Escolhemos um intervalo onde ha uma unica raiz
Ü A raiz esta entre entre x = a e x = b tais quef (a)f (b) < 0.Ü O ponto medio e uma primeira aproximacao para aposicao da raiz. Iterando ...
![Page 8: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/8.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao
Ü Sabemos visualmente onde estao as raızes pelamudanca de sinal da funcao.Ü Escolhemos um intervalo onde ha uma unica raizÜ A raiz esta entre entre x = a e x = b tais quef (a)f (b) < 0.
Ü O ponto medio e uma primeira aproximacao para aposicao da raiz. Iterando ...
![Page 9: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/9.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao
Ü Sabemos visualmente onde estao as raızes pelamudanca de sinal da funcao.Ü Escolhemos um intervalo onde ha uma unica raizÜ A raiz esta entre entre x = a e x = b tais quef (a)f (b) < 0.Ü O ponto medio e uma primeira aproximacao para aposicao da raiz.
Iterando ...
![Page 10: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/10.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao
Ü Sabemos visualmente onde estao as raızes pelamudanca de sinal da funcao.Ü Escolhemos um intervalo onde ha uma unica raizÜ A raiz esta entre entre x = a e x = b tais quef (a)f (b) < 0.Ü O ponto medio e uma primeira aproximacao para aposicao da raiz. Iterando ...
![Page 11: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/11.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Encerrando um Algoritmo
Ü Se a solucao do problema e xraiz e a precisao desejadae ε, estamos procurando um ponto xN tal que
|xN − xraiz | < ε
Ü Mas nao conhecemos xraiz !Ü Se a sequencia {xi} converge para xraiz devemos ter:
|xN+1 − xN | < ε
Ü Esse e um criterio de convergencia muito usado
![Page 12: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/12.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Encerrando um Algoritmo
Ü Se a solucao do problema e xraiz e a precisao desejadae ε, estamos procurando um ponto xN tal que
|xN − xraiz | < ε
Ü Mas nao conhecemos xraiz !
Ü Se a sequencia {xi} converge para xraiz devemos ter:
|xN+1 − xN | < ε
Ü Esse e um criterio de convergencia muito usado
![Page 13: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/13.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Encerrando um Algoritmo
Ü Se a solucao do problema e xraiz e a precisao desejadae ε, estamos procurando um ponto xN tal que
|xN − xraiz | < ε
Ü Mas nao conhecemos xraiz !Ü Se a sequencia {xi} converge para xraiz devemos ter:
|xN+1 − xN | < ε
Ü Esse e um criterio de convergencia muito usado
![Page 14: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/14.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Encerrando um Algoritmo
Ü Se a solucao do problema e xraiz e a precisao desejadae ε, estamos procurando um ponto xN tal que
|xN − xraiz | < ε
Ü Mas nao conhecemos xraiz !Ü Se a sequencia {xi} converge para xraiz devemos ter:
|xN+1 − xN | < ε
Ü Esse e um criterio de convergencia muito usado
![Page 15: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/15.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Encerrando um Algoritmo
Ü Se a solucao do problema e xraiz e a precisao desejadae ε, estamos procurando um ponto xN tal que
|xN − xraiz | < ε
Ü Mas nao conhecemos xraiz !Ü Se a sequencia {xi} converge para xraiz devemos ter:
|xN+1 − xN | < ε
Ü Esse e um criterio de convergencia muito usado
![Page 16: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/16.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao - Observacoes
O metodo sempre converge para uma raiz caso ela exista eesteja no intervalo inicial, mas nao funciona em algumassituacoes:
8 varias raızes8 singularidades8 pontos extremos8 pode ser muito
lento
![Page 17: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/17.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao - Observacoes
O metodo sempre converge para uma raiz caso ela exista eesteja no intervalo inicial, mas nao funciona em algumassituacoes:
8 varias raızes
8 singularidades8 pontos extremos8 pode ser muito
lento
![Page 18: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/18.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao - Observacoes
O metodo sempre converge para uma raiz caso ela exista eesteja no intervalo inicial, mas nao funciona em algumassituacoes:
8 varias raızes8 singularidades
8 pontos extremos8 pode ser muito
lento
![Page 19: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/19.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao - Observacoes
O metodo sempre converge para uma raiz caso ela exista eesteja no intervalo inicial, mas nao funciona em algumassituacoes:
8 varias raızes8 singularidades8 pontos extremos
8 pode ser muitolento
![Page 20: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/20.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo da Bissecao - Observacoes
O metodo sempre converge para uma raiz caso ela exista eesteja no intervalo inicial, mas nao funciona em algumassituacoes:
8 varias raızes8 singularidades8 pontos extremos8 pode ser muito
lento
![Page 21: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/21.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Algoritmo
m Escolher visualmente xmin e xmax
m Solicitar xmin e xmax e verificar se esses limites delimitamuma raiz (f (xmin)f (xmax) < 0)
m se delimitam, repetir enquanto xmax − xmin > precisao e onumero de iteracoes < limite :
Ü contar o numero de iteracoesÜ calcular o valor medio (xmedio)Ü calcular a funcao no ponto medio: f (xmedio)Ü escolher o intervalo entre xmin e (xmedio) ou entre (xmedio) e xmax :
Se f (xmin)f (xmedio) < 0 (intervalo da esquerda): xmax = xmedio
Senao (intervalo da direita): xmin = xmedio
m escrever numero de iteracoes e raiz=xmedio
![Page 22: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/22.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Exercıcio
8 Escreva um programa que implemente o algoritmo dometodo da bissecao. Aplique-o para encontrar as raızes daequacao:
f (x) = e−x − sin(πx
2
)= 0.
8 O programa deve imprimir, alem da raiz, o numero deiteracoes utilizadas8 Usando o gnuplot faca um grafico da funcao indicando ointervalo inicial.
![Page 23: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/23.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Um metodo do seculo 17
8 O metodo foi publicado pela primeira vez em 1690 porJoseph Raphson (1648-1715). Entre suas obras ha umatraducao de uma obra de Newton (Arithmetica Universalis)para o ingles.8 Isaac Newton (1643-1727) desenvolveu o mesmometodo, de forma independente, em 1671 mas nao opublicou em vida.8 E o metodo mais usado para se encontrar as raızes deuma equacao.
![Page 24: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/24.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson
Ü Em vez de escolher um intervalo
⇒ use a reta tangente.Ü O ponto em que a reta corta o eixo e o candidato a raizÜ Iterando . . .
![Page 25: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/25.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson
Ü Em vez de escolher um intervalo ⇒ use a reta tangente.
Ü O ponto em que a reta corta o eixo e o candidato a raizÜ Iterando . . .
![Page 26: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/26.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson
Ü Em vez de escolher um intervalo ⇒ use a reta tangente.Ü O ponto em que a reta corta o eixo e o candidato a raiz
Ü Iterando . . .
![Page 27: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/27.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson
Ü Em vez de escolher um intervalo ⇒ use a reta tangente.Ü O ponto em que a reta corta o eixo e o candidato a raizÜ Iterando . . .
![Page 28: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/28.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Ilustracao
![Page 29: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/29.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Ilustracao
![Page 30: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/30.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Ilustracao
![Page 31: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/31.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Ilustracao
![Page 32: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/32.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Ilustracao
![Page 33: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/33.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Cuidado: O ponto inicial nao pode ser qualquer um.
![Page 34: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/34.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Cuidado: O ponto inicial nao pode ser qualquer um.8 Equacao da reta tangente a curva f (x) no ponto x0:
r(x) = f (x0) + (x − x0) f ′(x0)
![Page 35: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/35.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Equacao da reta tangente a curva f (x) no ponto x0:
r(x) = f (x0) + (x − x0) f ′(x0)
8 Calculo do ponto xi+1 (mudando ligeiramente anotacao:
0 = f (xi) + (xi+1 − xi) f ′(xi)
−f (xi) = (xi+1 − xi) f ′(xi)
− f (xi)
f ′(xi)= (xi+1 − xi)
xi+1 = xi −f (xi)
f ′(xi)
![Page 36: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/36.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Equacao da reta tangente a curva f (x) no ponto x0:
r(x) = f (x0) + (x − x0) f ′(x0)
8 Calculo do ponto xi+1 (mudando ligeiramente anotacao:
0 = f (xi) + (xi+1 − xi) f ′(xi)
−f (xi) = (xi+1 − xi) f ′(xi)
− f (xi)
f ′(xi)= (xi+1 − xi)
xi+1 = xi −f (xi)
f ′(xi)
![Page 37: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/37.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Equacao da reta tangente a curva f (x) no ponto x0:
r(x) = f (x0) + (x − x0) f ′(x0)
8 Calculo do ponto xi+1 (mudando ligeiramente anotacao:
0 = f (xi) + (xi+1 − xi) f ′(xi)
−f (xi) = (xi+1 − xi) f ′(xi)
− f (xi)
f ′(xi)= (xi+1 − xi)
xi+1 = xi −f (xi)
f ′(xi)
![Page 38: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/38.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Equacao da reta tangente a curva f (x) no ponto x0:
r(x) = f (x0) + (x − x0) f ′(x0)
8 Calculo do ponto xi+1 (mudando ligeiramente anotacao:
0 = f (xi) + (xi+1 − xi) f ′(xi)
−f (xi) = (xi+1 − xi) f ′(xi)
− f (xi)
f ′(xi)= (xi+1 − xi)
xi+1 = xi −f (xi)
f ′(xi)
![Page 39: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/39.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Equacao da reta tangente a curva f (x) no ponto x0:
r(x) = f (x0) + (x − x0) f ′(x0)
8 Calculo do ponto xi+1 (mudando ligeiramente anotacao:
0 = f (xi) + (xi+1 − xi) f ′(xi)
−f (xi) = (xi+1 − xi) f ′(xi)
− f (xi)
f ′(xi)= (xi+1 − xi)
xi+1 = xi −f (xi)
f ′(xi)
![Page 40: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/40.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Metodo de Newton-Raphson - Detalhamento
8 Equacao da reta tangente a curva f (x) no ponto x0:
r(x) = f (x0) + (x − x0) f ′(x0)
8 Calculo do ponto xi+1 (mudando ligeiramente anotacao:
0 = f (xi) + (xi+1 − xi) f ′(xi)
−f (xi) = (xi+1 − xi) f ′(xi)
− f (xi)
f ′(xi)= (xi+1 − xi)
xi+1 = xi −f (xi)
f ′(xi)
![Page 41: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/41.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Exercıcio
8 Escreva um algoritmo que procure a raiz de uma equacaopelo metodo de Newton Raphson.
8 Escreva um programa que implemente o algoritmo dometodo de Newton Raphson. Aplique-o para encontrar asraızes da equacao:
f (x) = e−x − sin(πx
2
)= 0.
8 O programa deve imprimir, alem da raiz, o numero deiteracoes utilizadas8 Compare o resultado com o obtido pelo metodo dabissecao.
![Page 42: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/42.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Exercıcio
8 Escreva um algoritmo que procure a raiz de uma equacaopelo metodo de Newton Raphson.8 Escreva um programa que implemente o algoritmo dometodo de Newton Raphson. Aplique-o para encontrar asraızes da equacao:
f (x) = e−x − sin(πx
2
)= 0.
8 O programa deve imprimir, alem da raiz, o numero deiteracoes utilizadas8 Compare o resultado com o obtido pelo metodo dabissecao.
![Page 43: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/43.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Observacoes
Ü O metodo de Newton Raphson pode resolver casos emque o metodo da bissecao falha.
![Page 44: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/44.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Observacoes
Ü O metodo de Newton Raphson pode resolver casos emque o metodo da bissecao falha.
Ü Nem sempre converge!
![Page 45: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/45.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Observacoes
Ü O metodo de Newton Raphson pode resolver casos emque o metodo da bissecao falha.
Ü Nem sempre converge!Ü O calculo da derivada pode ser um problema!
![Page 46: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/46.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
Observacoes
Ü O metodo de Newton Raphson pode resolver casos emque o metodo da bissecao falha.
Ü Nem sempre converge!Ü O calculo da derivada pode ser um problema!Ü Em geral o metodo de Newton-Raphson e mais
eficiente mas ha casos em que e melhor usar outrosmetodos.
![Page 47: Metodos Computacionais em F´ısica´ Erica Polycarpo Sandra ...sandra/MetComp/2008-2/doc/Aula6.pdf · Met Comp´ Metodos Computacionais em F´ısica´ 1 Determinac¸ao de Ra´ızes˜](https://reader031.fdocument.pub/reader031/viewer/2022021808/5bea0d3209d3f260758b8486/html5/thumbnails/47.jpg)
Met Comp
RaızesBissecao
Newton-Raphson
BOA PROVA