# Problemas (Física)

## Problema (Movimento de partícula)

A posição de uma partícula à medida que se move ao longo do eixo $y$ é dada por 
$y(t) = 2\sin\left(\frac{\pi t}{4}\right),$
com $t$ em segundos e $y$ em centímetros. Calcule: 

a. A velocidade média da partícula entre $t = 0$ e $t=2.0 \, s$.

b. A velocidade instantânea da partícula em $t = 0, \, 1.0$ e $2.0 \, s$.

c. A aceleração média da partícula entre $t = 0$ e $t=2.0 \, s$.

d. A aceleração instantânea da partícula em $t = 0, \, 1.0$ e $2.0 \, s$.

(Adaptado de Halliday & Resnick)

### Resolução


```python
import numpy as np
import sympy as sym
import matplotlib.pyplot as plt

def y(t):
    return 2*np.sin( (np.pi*t)/4 )

# calculo de média
def m(tf,t0,f):
    return (f(tf) - f(t0))/(tf - t0)
```


```python
# a.
print(f'Velocidade média = {np.round(m(0,2,y),2)} cm/s') 
```

    Velocidade média = 1.0 cm/s



```python
# b. c. e d.

# Cria expressão simbólica
tsym = sym.Symbol('t')
def ysym(tsym):
    return 2*sym.sin( (sym.pi*tsym)/4 )
print(f'Velocidade média = {np.round(m(0,2,y),2)} cm/s') 

# Calcula velocidade instantânea
for t0 in [0.,1.0,2.0]:    
    dydt = sym.diff(ysym(tsym),tsym) # dy/dt  
    vinst = dydt.subs({'t':t0})    
    print(f'Velocidade instantânea em t = {t0} é {vinst.evalf(3)} cm/s') 

for t0 in [0.,1.0,2.0]: 
    if t0 == 2.0:
        amed = (dydt.subs({'t':t0}) - dydt.subs({'t':0}))/2.0
        print(f'Aceleração média = {amed.evalf(2)} cm/s^2')     
                 
    d2ydt2 = sym.diff(ysym(tsym),tsym,2) # d2y/dt2
    ainst = d2ydt2.subs({'t':t0})    
    print(f'Aceleração instantânea em t = {t0} é {ainst.evalf(3)} cm/s^2') 
```

    Velocidade média = 1.0 cm/s
    Velocidade instantânea em t = 0.0 é 1.57 cm/s
    Velocidade instantânea em t = 1.0 é 1.11 cm/s
    Velocidade instantânea em t = 2.0 é 0 cm/s
    Aceleração instantânea em t = 0.0 é 0 cm/s^2
    Aceleração instantânea em t = 1.0 é -0.872 cm/s^2
    Aceleração média = -0.79 cm/s^2
    Aceleração instantânea em t = 2.0 é -1.23 cm/s^2



```python
tnum = np.linspace(0,5)
plt.plot(tnum,y(tnum));
```

583e734c95b85e884f71bd1b6d20a6713af4f8ec
ipynb
Jupyter Notebook
_build/html/_sources/ipynb/problems/fisica.ipynb
gcpeixoto/FMECD
9bca72574c6630d1594396fffef31cfb8d58dec2
CC0-1.0
_build/html/_sources/ipynb/problems/fisica.ipynb
gcpeixoto/FMECD
9bca72574c6630d1594396fffef31cfb8d58dec2
CC0-1.0
_build/html/_sources/ipynb/problems/fisica.ipynb
gcpeixoto/FMECD
9bca72574c6630d1594396fffef31cfb8d58dec2
CC0-1.0
Qwen/Qwen-72B
1. YES
2. YES