模拟轨道

时间:2017-07-07 19:48:57

标签: python-3.x vpython orbital-mechanics

所以我试图模拟绕太阳运行的地球,地球的速度由它与原点和水平面的角度决定。我通过创建一个使用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

非常感谢一些帮助, 干杯

1 个答案:

答案 0 :(得分:0)

我不懂python,但我知道物理。

在每一步中,您沿着切线将地球移动到轨道,而不是沿着轨道本身移动。这会给你一个向外的螺旋(当你外出时实际上会变得不那么严重)。

尝试使时间增量更小(例如,将位置调整除以100),螺旋效果会变得更小。

如果你想做得更好,你需要一个不同的公式。你可以施加一个圆形轨道,或者根据守恒量做一些事情(这需要对基本物理学有相当的了解)。