Entradas

Mostrando entradas de junio, 2021

Clase del día 29 de Junio de 2021 parte 2

Imagen
import  scipy.integrate  as  integrate import  matplotlib.pyplot  as  plt import  numpy  as  np Las condiciones iniciales que se puede presentar en el sistema serian,  t = 0 , x = x 0 ,  y  x = 0 . Por otra parte la otra combinacion seria  t = 0 , x = 0  y  x ˙ = v 0 Aplicamos la primera combinacion de condiciones iniciales en la ecuacion de  x  teniendo: x 0 = A ∗ ( 0 ) + B ∗ ( 1 ) B = x 0 Ahora sustituimos en  x ˙ 0 = A ∗ w n ∗ ( 1 ) A = 0 Definimos entonces que: x 1 = x 0 c o s ( w n t ) Continuamos sustituyendo las condiciones iniciales restantes 0 = A ∗ ( 0 ) + B ∗ ( 1 ) B = 0 Sustituyendo en  x ˙ v 0 = A ∗ w n ∗ ( 0 ) A = v 0 w n De ahi tenemos que x 2 = v 0 w n s i n ( w n t ) Entonces x = x 1 + x 2 Desarrollando x = v 0 w n s i n ( w n t ) + x 0 c o s ( w n t ) pi=np.pi sqrt=np.sqrt cos=np.cos sin=np.sin def   deriv_x ( d2x ,  phi ):  x, xdot=d2x   return [xdot...