本文整理匯總了Python中planeta.Planeta.verlet_2方法的典型用法代碼示例。如果您正苦於以下問題:Python Planeta.verlet_2方法的具體用法?Python Planeta.verlet_2怎麽用?Python Planeta.verlet_2使用的例子?那麽, 這裏精選的方法代碼示例或許可以為您提供幫助。您也可以進一步了解該方法所在類planeta.Planeta
的用法示例。
在下文中一共展示了Planeta.verlet_2方法的1個代碼示例,這些例子默認根據受歡迎程度排序。您可以為喜歡或者感覺有用的代碼點讚,您的評價將有助於係統推薦出更棒的Python代碼示例。
示例1: rk4
# 需要導入模塊: from planeta import Planeta [as 別名]
# 或者: from planeta.Planeta import verlet_2 [as 別名]
vy = sp.zeros(pasos)
energia=sp.zeros(pasos)
[x[0],y[0],vx[0],vy[0]] = condicion_inicial
energia[0]=p.energia_total()
''' Se define el segundo valor del vector yactual mediante rk4 (en este caso) puesto que verlet2 lo necesita para la iteracion'''
p.avanza_rk4(dt)
resultados = p.y_actual
x[1] = resultados[0]
y[1] = resultados[1]
vx[1] = resultados[2]
vy[1] = resultados[3]
''' Se obtienen los arreglos con todos sus valores mediante el metodo de verlet 2'''
for i in range (2,pasos):
p.verlet_2(dt,x[i-2],y[i-2])
resultados = p.y_actual
x[i] = resultados[0]
y[i] = resultados[1]
vx[i] = resultados[2]
vy[i] = resultados[3]
energia[i] = p.energia_total()
''' Se grafican los resultados, en este caso energia en funcion del tiempo y la trayectoria para el metodo de verlet 2 con alpha=0'''
plt.figure(2)
plt.subplot(2, 1, 1)
plt.subplots_adjust(hspace=.5)
plt.scatter(x,y,6,c=u'darkred',label = "trayectoria")
plt.title(" Solucion verlet 2 alfa = 0 ")
plt.xlabel("x [m]")
plt.ylabel("y [m]")