Homework_2018_1_04_1032478036.ipynb

Kernel: Python 3

TAREA: Para el polinomio de grado 2, $$\begin{equation} P (t) = -4.9t^2 + 14.7t \end{equation}$$

que interpola la caída libre de 3 segundos de vuelo de lanzamiento vertical hacia arriba, grafique:

a) P(t)

b) P'(t)

c) P''(t)

Interprete fisicamente los resultados

import numpy as np
from scipy.misc import derivative
import matplotlib.pyplot as plt
%pylab inline

Populating the interactive namespace from numpy and matplotlib

def derivate(func,x0,dx=1,n=1,args=(),order=3):
    try:
        nn=np.array(x0).shape[0]
    except IndexError:
        nn=-1
    if nn>-1:
        y=[]
        for xx in x0:
             y.append(derivative(func,xx,dx=dx,n=n, order=order))
                
        y=np.array(y)   
    else:
        y=derivative(func,xx, dx=dx, n=n, order=order)
        
    return y

t=np.linspace(0,3,50)
pol=lambda t: -4.9*t**2+14.7*t

plt.figure( figsize = (16,4) )

plt.subplot(1,3,1)
plt.plot(t,pol(t), color="blue", linewidth=3, label="P(t)")
plt.ylabel("Posición(m)")
plt.xlabel("Tiempo(s)")
plt.legend( loc="upper right", fontsize=8 )
plt.grid()

plt.subplot(1,3,2)
plt.plot(t, derivate(pol,t,dx=1E-5,n=1,order=3),'g:', linewidth=3, label="1ra derivada")
plt.plot(t, derivate(pol,t,dx=1E-5,n=2,order=5),'r:', linewidth=3, label="2da derivada")
plt.xlabel("Tiempo(s)")
plt.legend( loc="upper right", fontsize=8 )
plt.grid()

La función P(t) se obtuvo como producto de la interpolación de puntos que describían el movimiento de caída libre. De esta manera, la gráfica que se muestra en azul es dicha representación e indica la posición de la partícula desde que arranca el movimiento, alcanza su altura máxima en la mitad del tiempo (1.5s) y regresa al punto de partida a los 3s.

Por otro lado la primera derivada de P(t), representada por la línea verde, corresponde a la velocidad y su comportamiento muestra que al ser lanzamiento vertical hacia arriba, su valor inicial es distinto de cero y en este caso corresponde a $14.7 \frac{m}{s}$; luego cuando la altura es máxima su valor es cero y después de esto los valores se vuelven negativos, indicando un cambio en la dirección del movimiento (regresa cayendo a la altura inicial), hasta volver a alcanzar los $14.7 \frac{m}{s}$ pero en dicho sentido.

Por último, la segunda derivada corresponde a la velocidad y permite verificar que en este tipo de movimiento la aceleración se mantiene constante pues solo actúa la gravedad y por ende su representación gráfica es una recta de pendiente cero en un valor constante igual a $g=9.8 \frac{m}{s²}$ y con signo de acuerdo al sistema de coordenadas.