文件名称:N-body-Gravity-Simulator:用于模拟天体轨道的Python脚本
文件大小:138KB
文件格式:ZIP
更新时间:2024-06-07 02:17:34
Python
N体重力模拟器 这是计算物理研究生院的一个数值方法项目。 我的python脚本根据引力对天体的轨道进行建模。 行星的初始条件来自NASA。 方法 为了计算行星的轨迹,我使用了Leapfrog积分方法。 Leapfrog是一种二阶方法,因为它具有时间可逆性,因此对机械师很有用。 如果您将n步进行积分,然后再进行n步向后积分,那么您将到达初始位置。 蛙跳也很辛苦,因为它可以节省能量。 Runge-Kutta是一种四阶方法,也很流行,但是它不节省能量,并且系统会随着时间的流逝而漂移。 我探查了太阳系,看它是否装满了。 在模拟中,您可以在火星和木星之间添加一个地球大小的星球。 考虑到小行星带的质量,这是有道理的。 结果 我们期望总能量和总角动量得以守恒。 能量和动量变化的舍入误差会随着时间的流逝而抵消并很好地振荡。 在过去的40年中,请注意能量和动量变化的幅度很小。 我将冥王星包含在模拟中,
【文件预览】:
N-body-Gravity-Simulator-master
----img()
--------NBodyOrbit40Energy.png(42KB)
--------NBodyOrbit40.png(49KB)
--------NBodyOrbit40Momentum.png(47KB)
--------readme(1B)
----starData.py(7KB)
----NBody.py(15KB)
----README.md(2KB)