0

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

非常感谢一些帮助,干杯

4

1 回答 1

0

我不懂python,但我懂物理。

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

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

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

于 2017-07-07T20:01:05.733 回答