Este documento apresenta 4 exercícios que envolvem a resolução numérica de equações diferenciais ordinárias (EDOs) usando os métodos de Euler e Picard. O primeiro exercício compara as soluções do oscilador harmônico usando esses dois métodos. O segundo resolva a equação de Bessel numericamente. O terceiro estuda o caos na equação de van der Pol. E o quarto resolve diversas EDOs de valor inicial.
Este documento apresenta 4 exercícios que envolvem a resolução numérica de equações diferenciais ordinárias (EDOs) usando os métodos de Euler e Picard. O primeiro exercício compara as soluções do oscilador harmônico usando esses dois métodos. O segundo resolva a equação de Bessel numericamente. O terceiro estuda o caos na equação de van der Pol. E o quarto resolve diversas EDOs de valor inicial.
Este documento apresenta 4 exercícios que envolvem a resolução numérica de equações diferenciais ordinárias (EDOs) usando os métodos de Euler e Picard. O primeiro exercício compara as soluções do oscilador harmônico usando esses dois métodos. O segundo resolva a equação de Bessel numericamente. O terceiro estuda o caos na equação de van der Pol. E o quarto resolve diversas EDOs de valor inicial.
Este documento apresenta 4 exercícios que envolvem a resolução numérica de equações diferenciais ordinárias (EDOs) usando os métodos de Euler e Picard. O primeiro exercício compara as soluções do oscilador harmônico usando esses dois métodos. O segundo resolva a equação de Bessel numericamente. O terceiro estuda o caos na equação de van der Pol. E o quarto resolve diversas EDOs de valor inicial.
Baixe no formato PDF, TXT ou leia online no Scribd
Fazer download em pdf ou txt
Você está na página 1de 2
Lista 3 - Métodos de Euler e Picard
Fı́sica Computacional II (092444)
Prof. Matheus Paes Lima
Exercı́cio 1 - Oscilador Harmônico
Utilizando o método de Euler, no qual sem perda de generalidade uma EDO pode ser escrita na forma: d~y = ~g(y, t) (1) dt (a) Escreva um programa para o oscilador harmônico (my 00 + ky = 0). utilize m = k = 1, T max = 20., itmax = 200 (número de intervalos de tempo). Considere como condições iniciais a partı́cula parada em y = 1. Apresente o gráfico da posição em função do tempo (b) Faça a correção Preditora-corretora baseada no método de Picard nas mesmas condições do item (a), e compare os resultados. Apresente o gráfico, e comente o que você observou. (c) Aumente o itmax para 2000 e faça um gráfico comparando as soluções obtidas com o método de Euler, e preditora-corretora-Picard. Compare as quatro curvas obtidas, e comente o que foi observado. (d) Estime o valor mı́nimo de itmax para que o método preditor-corretor- Picard forneça resultados razoáveis no intervalo adotado. Justifique os critérios utilizados.
Exercı́cio 2 - Solução numérica da equação de Bessel
As funções de Bessel são amplamente utilizadas na fı́sica, tendo aplicações em Ondas eletromagnéticas, Condução de calor, Vibração, Difusão, Proces- samento de sinais (filtro Bessel), Soluções para a equação de Schrödinger radial (em coordenadas esféricas e cilı́ndricas) de uma partı́cula livre, etc. A EDO que define uma função de bessel é: x2 y 00 + xy 0 + (x2 − p2 )y = 0 (2) Na qual y = y(x) é uma função da variável independente x, e p é uma constante. (a) Coloque a EDO de Bessel no formato d~y = ~g(y, t), (3) dt 1 obtendo uma função ~g(y, t) válida. Apresente a função obtida. (b) Obtenha a solução da equação acima para p = 0 com o método preditor- corretor de Picard, tendo como condições iniciais x0 = 0.0001, y(x0 ) = 1, y 0 (x0 ) = −0.0001, no intervalo x0 < x < 20.0, considerando um grid com 200 subdivisões neste intervalo. Aumente o grid para 2000 subdivisões. Apre- sente um gráfico comparativo e comente seus resultados. (c) Obtenha ao menos uma solução ya (x) para p=1. Observe que as condições de contorno podem variar.
Exercı́cio 3 - Caos: A equação de van der Pol
As coordenadas x e y para um certo sistema hipotético são determinadas pelas seguintes EDOs acopladas: dx = y − fµ (x) (4) dt dy = −x (5) dt onde,
fµ (x) = x3 − µx. (6)
µ ∈ [−1, +1]. Dependendo do valor de µ, o sistema pode ser atraı́do para um
polo. Determine a dinâmica do sistema (gráfico de x .vs. y) para diversos valores da constante real µ, e responda para que valores de µ a solução é atraı́da para um polo. (dica: utilize o método de Picard, tendo como condições iniciais (x0 , y0 )=(0, 1), (1, 0), (0, −1), e (−1, 0)). Neste exercı́cio, os parâmetros de simulação devem ser determinados por você. Exercı́cio 4 - Soluções de EDOs diversas. Para cada um dos problemas de valor inicial abaixo coloque as EDOs no formato d~ y dt = ~g(y, t), e faça um gráfico com a solução. Para cada item, envie o gráfico e o programa utilizado. (Observe que há mais de uma condição inicial) (a) 2y 0 (t) + ty(t) = 2, y(t = 0) = 1 e y(t = 0) = −1, e 0 ≤ t ≤ 6 (b) y 0 + (1/2)y = 10 − 5 sin (2t), y(t = 0) = 0, e y(t = 0) = 10, 0 ≤ t ≤ 20 (c) y 00 + 5y 0 + 6y = 0, y(0) = 2, y 0 (0) = 3, 0 ≤ t ≤ 2.5 (d) y 0000 − y = 0, y(0) = 7/2, y 0 (0) = −4, y 00 (0) = 5/2, y 000 (0) = −2, 0 ≤ t ≤ 15