显示奇异行为的二维轨道物理模型

Question

我试图创建一个二维环境，它模拟基本的轨道物理使用Pygame。然而，当我运行它时，它会表现出一些奇怪的行为（编辑：以及一些更奇怪的行为），比如卫星从x轴和y轴上弹出，以及卫星沿着太阳轴以波浪状的模式移动。下面是计算由中央太阳引起的物体加速度的代码。BM为体质量，BP为体位[x，y]，SM为太阳质量，SP为太阳位[x，y]。

def CalcGrav(BM,BP,SM,SP):  
    Dist = SubList(BP,SP)  
    Mass = BM*SM/100  
    for i in range(0,2):  
        if(Dist[i] == 0):  
            Dist[i] += 0.01  
    AV = [-(Mass/Dist[0]),-(Mass/Dist[1])]  
return AV  

然后在以下上下文中使用该函数：

TraceAcc = CalcGrav(BallMass,TracePosition,SunMass,SunPosition)
TraceVector = SumList(TraceVector, TraceAcc)
TracePosition = SumList(TracePosition,TraceVector)

SumList（A,B）只是将一个列表的项相加在一起：（1,3)+(2,6)=(3,9)
当前速度向量和加速度向量存储为列表:[x,y]

Answer 1

问题是你没有构建一个现实的轨道运动物理模型。

你需要做的是建立一个简化的开普勒轨道模型来进行一个比较真实的模拟。您还需要在极坐标（半径，角度）中工作，并转换为x，y坐标。

你需要做的是用代码构造椭圆轨道的开普勒轨道方程。