Universidade Federal do Rio Grande do Sul
Instituto de Física - Departamento de Física
FIS01206 - UNIDADE I - Lista de Problemas


1.No problema do decaimento radioativo,

$\displaystyle \frac{dx}{dt} = -\,\alpha\, x \; ; \; x(0)\,=\,x_0\,\,,
$

determine o limite superior do passo de tempo $ \Delta t$ para o qual o método de Euler converge para uma solução assintótica nula.


2.No problema do decaimento radioativo,

$\displaystyle \frac{dx}{dt} = -\,\alpha\, x \; ; \; x(0)\,=\,x_0\,\,,
$

determine o limite superior do passo de tempo $ \Delta t$ para o qual o método do Ponto Médio (Runge-Kutta em segunda ordem) converge para uma solução assintótica nula.


3.No problema do decaimento radioativo,

$\displaystyle \frac{dx}{dt} = -\,\alpha\, x \; ; \; x(0)\,=\,x_0\,\,,
$

determine o limite superior do passo de tempo $ \Delta t$ para o qual o método de Heun converge para uma solução assintótica nula.

O método de Heun usa a média entre as derivadas no início e final do intervalo $ \Delta t$:

    $\displaystyle x(t+\Delta t)=x(t)+\Delta t\times$ (1)
    $\displaystyle \{f(t,x(t))+f[t+\Delta t, x(t)+\Delta t f(t,x(t))]\}/2.$ (2)


4.No problema do decaimento radioativo,

$\displaystyle \frac{dx}{dt} = -\,\alpha\, x \; ; \; x(0)\,=\,x_0\,\,,
$

determine o limite superior do passo de tempo $ \Delta t$ para o qual o método de Ralston converge para uma solução assintótica nula.

O método de Ralston usa uma média ponderada entre as derivadas no início e final do intervalo $ \Delta t$:

    $\displaystyle x(t+\Delta t)=x(t)+\Delta t \{f(t,x(t))+$ (3)
    $\displaystyle 2f[t+\frac{3}{4}\Delta t, x(t)+\frac{3}{4}\Delta t f(t,x(t))]\}/3.$ (4)


5.Usando o método de Adams-Bashforth

$\displaystyle x_{n+1} = x_{n} + \frac{3}{2}\, \Delta t \, f(x_{n}) \,-\, \frac{1}{2}\, \Delta t \, f(x_{n-1})
$

e Adams-Moulton

$\displaystyle x_{n+1} = x_{n} + \frac{1}{2}\, \Delta t \, f(x_{n+1}) \,+\, \frac{1}{2}\, \Delta t \, f(x_{n})
$

a)
Escreva um algoritmo do tipo predição-correção para o problema do sistema abaixo

$\displaystyle \frac{d^2 x}{d t^2} = - k\,x .
$

b)
Na linguagem que lhe for familiar, escreva o programa correspondente.

c)
Com base na solução analítica, faça uma estimativa para um valor inicial de $ \Delta t$. Justifique a sua resposta.

6.Descreva um algoritmo para encontrar numericamente a ordem do erro global de um método de integração numérico qualquer. Explique como deve ser feito o ajuste para que se encontre a lei de potência associada.

7.Na aproximação de ângulos pequenos, a energia total de um pêndulo simples é

$\displaystyle E = \frac{1}{2} m L^2 \omega^2 + \frac{1}{2} m g L \theta^2 - m\,g\,L \,.
$

Mostre analiticamente que $ E$ aumenta monotonicamente com o tempo quando o método de Euler é utilizado na integração numérica.

O que acontece quando utilizamos o método de Euler-Cromer?

8.Prove que para o problema de Kepler o método de Euler-Cromer conserva momentum angular de forma exata.

O que ocorre quando se utiliza o método de Euler?




Leonardo Gregory Brunnet 2017-05-01