所以我试图模拟绕太阳运行的地球,地球的速度由它与原点和水平面的角度决定。我通过创建一个使用tanh(相反/相邻)规则的三角形O_correction(x,y)
来完成此操作。问题是,它不是圆形轨道而是螺旋式旋转,我不知道为什么。
scene = canvas()
scene.background = color.white
O = 0
ball = sphere(pos=vector(10,0,0), radius=0.1, color=color.blue)
x = ball.pos.x
y = ball.pos.y
def O_correction(x,y):
O = math.atan((((y)**2)**0.5)/(((x)**2)**0.5))
answer = O
if x >= 0 and y >= 0:
answer = O
if x < 0 and y >= 0:
answer = pi - O
if x <= 0 and y < 0:
answer = O + pi
if x > 0 and y < 0:
answer =pi*2 - O
return answer
t =0
while t < 100:
x = ball.pos.x
y = ball.pos.y
print = (float(O_correction(x,y))
print = ((x**2) + (y**2))**0.5)
ball.pos.x -= sin(O_correction(x,y))
ball.pos.y += cos(O_correction(x,y))
print(" ")
t += 1
非常感谢一些帮助, 干杯
答案 0 :(得分:0)
我不懂python,但我知道物理。
在每一步中,您沿着切线将地球移动到轨道,而不是沿着轨道本身移动。这会给你一个向外的螺旋(当你外出时实际上会变得不那么严重)。
尝试使时间增量更小(例如,将位置调整除以100),螺旋效果会变得更小。
如果你想做得更好,你需要一个不同的公式。你可以施加一个圆形轨道,或者根据守恒量做一些事情(这需要对基本物理学有相当的了解)。