UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico...

22
Histórico Equações de Lorenz Conclusão UM SISTEMA CAÓTICO Cristina Teruko Ota Junho, 2013 Cristina Teruko Ota UM SISTEMA CAÓTICO

Transcript of UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico...

Page 1: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

UM SISTEMA CAÓTICO

Cristina Teruko Ota

Junho, 2013

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 2: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Conteúdo

1 Histórico

2 Equações de LorenzTolerância de erroPonto inicialParâmetro ρ

3 Conclusão

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 3: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

1 Edward Lorenz (1963)2 Teoria do caos

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 4: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

y’ = f (y) =

σ(y2 − y1)ρy1 − y2 − y1y3

y1y2 − by3

em que σ é o número de Prandtl e ρ é o número de Rayleigh,com σ, ρ,b > 0.

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 5: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 6: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Função ode23

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Neste caso: 3871 passos, 12562 avaliações da função e1.957266 s.

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 7: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Utilizando a função ode45

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Neste caso: 1425 passos, 10069 avaliações de função e1.836942 s.

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 8: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 9: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Utilizando tolerância de erro relativo de 10−6

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 10: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Utilizando tolerância de erro relativo de 10−7

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 11: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Utilizando tolerância de erro absoluto de 10−7

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 12: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Utilizando tolerância de erro relativo de 10−6

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 13: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Utilizando tolerância de erro relativo de 10−7

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 14: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Utilizando tolerância de erro absoluto de 10−7

com ρ = 28, σ = 10 e b = 83 e ponto inicial y(0) = (0,1,0)T .

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 15: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Pontos iniciais distintos

com ρ = 28, σ = 10 e b = 83 .

Pontos iniciais: y01 = (0,1,0)T , y02 = (1e − 5,1,0)T

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 16: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Pontos iniciais distintos

com ρ = 28, σ = 10 e b = 83 .

Pontos iniciais: y01 = (0,1,0)T , y02 = (1e − 5,1,0)T

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 17: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Valores de ρ distintos

com ρ = 9, σ = 10 e b = 83 .

Pontos iniciais: y01 = (0,1,0)T , y02 = (1e − 5,1,0)T

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 18: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Valores de ρ distintos

com ρ = 14, σ = 10 e b = 83 .

Pontos iniciais: y01 = (0,1,0)T , y02 = (1e − 5,1,0)T

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 19: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Valores de ρ distintos

com ρ = 14, σ = 10 e b = 83 .

Pontos iniciais: y01 = (0,1,0)T , y02 = (1e − 5,1,0)T

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 20: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Tolerância de erroPonto inicialParâmetro ρ

Valores de ρ distintos

com ρ = 23, σ = 10 e b = 83 .

Pontos iniciais: y01 = (0,1,0)T , y02 = (1e − 5,1,0)T

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 21: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

1 Teoria do caos2 Instabilidade dos resultados3 Imprevisibilidade a longo prazo de sistemas

determinísticos e dinâmicos não-lineares

Cristina Teruko Ota UM SISTEMA CAÓTICO

Page 22: UM SISTEMA CAÓTICO - ime.unicamp.brsandra/MS612/handouts/Projeto2Problema4.pdf · Histórico Equações de Lorenz Conclusão BibliografiaI E. Ott. Chaos in Dynamical Systems. 2

HistóricoEquações de Lorenz

Conclusão

Bibliografia I

E. Ott.Chaos in Dynamical Systems. 2 Edition.Cambridge University Press, Cambridge, 2002.

U. Ascher & L. Petzold.Computer Methods for Ordinary Differential Equations andDifferential-Algebraic Equations.Siam, Philadelphia, 1999.

S. H. Strogatz.Nonlinear Dynamics and Chaos.Addison-Wesley Publishing Company, Boston, 1994.

Cristina Teruko Ota UM SISTEMA CAÓTICO