Tratamento Numérico de EDOs - ROMA
Agenda do Curso
-
-
Página
-
Página
-
Arquivo
-
Tarefa
-
Tarefa
-
Tarefa
-
Tarefa
-
Tarefa
-
Arquivo
-
URL
-
Arquivo
-
Arquivo
-
Arquivo
-
Arquivo
-
Página
-
Arquivo
-
Arquivo
-
Arquivo
-
AULA #01 (2009-01-05)
- Informações gerais: objetivos da disciplina, ementa, critério de avaliação, sobre a lista de presença, bibliografia e como serão os plantões de dúvidas.
- Problema de Cauchy. Notação (e o abuso dela).
- Discretização de um problema de Cauchy.
- Alguns exemplos
- EXERCÍCIOS SUPLEMENTARES.
AULA #02 (2009-01-06)
- Exemplos de métodos numéricos: Euler, Euler Implícito e Adams-Bashforth de 4 passos. Definições de métodos numéricos de passo único (ou de 1 passo) e métodos numéricos de passo múltiplo (ou de n passos).
- Discretização da EDO y´(x)=y(x)*sin(y(x)) com yo = y(xo).
- Método de Passo Único (explícito) Genérico: yk+1 = yk + h*Ø(xk,yk,h), onde Ø(xk,yk,h) é a "função de discretização" (no Método de Euler, Ø(xk,yk,h)=f(xk,yk)).
- Primeiros conceitos fundamentais: erro local de discretização e erro global de discretização. Convergência.
- EXERCÍCIOS SUPLEMENTARES.
ERRO GLOBAL DE DISCRETIZAÇÃO: É o erro total produzido (acumulado e transformado pela Ø(xk,yk,h)) após k passos de integração.
AULA #03 (2009-01-07) -
- Exemplos de funções que não são Lipschitzianas: g(x)=x2, para x peretencente aos reais e g(x,y)=xy/(1-x2) numa vizinhança do ponto (1,0).
- Sob as hipóteses de Ø(xk,yk,h) ser Lipschitz na segunda variável e do erro local de discretização ter ordem p, mostramos que o método de passo único yk+1 = yk + h*Ø(xk,yk,h) tem seu erro global de discretização indo a zero com velocidade hp-1. Se, portanto, p>=2, então o método converge.
- Definição de consistência com a EDO e de ordem de consistência.
- Métodos da Série de Taylor. AVISO: Revejam série de Taylor em 1 e em 2 variáveis.
- EXERCÍCIOS SUPLEMENTARES.
AULA #04 (2009-01-08) - EXERCÍCIOS
- Plantões presenciais de terça-feira as 14h. Procure o monitor na sala 257A.
- Divulgação da Lista 1.
- Quem quiser saber mais sobre o Butcher envie um email para o monitor.
- Exemplo de programas em C e Scilab e uso do gnuplot.
- Exemplos: pêndulo simplificado, modelos de epidemiologia SIR e SIM.
- Transformação de EDO de ordem >=2 em sistema de EDO's de primeira ordem.
- Comentários sobre funções lipschitzianas.
- Veja um Nanotutorial de Scilab.
-
AULA #05 (2009-01-12)
Expansão assintótica do erro global de discretização em métodos de passo único (e suas aplicações na validação da implementação do método e na estimativa do erro global).
- Conceito de estabilidade absoluta e sobre sua importância na escolha do passo de integração (soluções exatas limitadas devem ser aproximadas por soluções aproximadas limitadas!).
- Método formal de obter o fator de amplificação (basta aplicar o método numérico ao problema modelo y' = lambda y, com lambda < 0.
- Como obter algo parecido com o lambda do problema modelo no caso de problemas de Cauchy genéricos, y' = f(x,y).
- Veja exemplos: Método de Euler Implícito (EXERCÍCIOS SUPLEMENTARES c.f. ex. 08 do ano passado!).
- O Método da Série de Taylor tem ordem q quando a função de discretização tem a forma Ø = f + hDf/2! + h2D2f/3! + ... + hq-1Dq-1f/q!
AULA#07 (2009-01-14)
Alguns comentários sobre estabilidade e sobre Métodos da Série de Taylor.
- Na prática, o passo de integração h>0 é escolhido analisando-se a estabilidade absoluta do método numérico considerado.
- Definição de região de estabilidade absoluta e de intervalo de estabilidade absoluta (usando-se o problema modelo y'= \lambda y).
- Importante: para o Método de Euler, a escolha do passo de integração h>0 num problema genérico y'=f(x,y) depende localmente de linearizações do problema, isto é, do problema linearizado, y'= df/dy(x,y) y (lambda, neste caso, é df/dy(x,y) suposto aproximadamente constante).
- Introdução aos métodos de Runge-Kutta (explícitos) com alguns exemplos.
AULA #08 (2009-01-15) - AULA DE EXERCÍCIOS
- Resolução de exercícios da lista
- Série de taylor de várias variáveis
- Validação de código
- Conceito de estabilidade absoluta e sobre sua importância na escolha do passo de integração (soluções exatas limitadas devem ser aproximadas por soluções aproximadas limitadas!).
-
AULA #09 (2009-01-19)
AULA #10 (2009-01-20) - AULA DE EXERCÍCIOS
- Resolução de exercícios da lista
- Sistemas de equações diferenciais ordinárias lineares
- Auto-valores e auto-vetores
- Retratos de fase
AULA #11 (2009-01-21)
AULA #12 (2009-01-22) - PROVA