轻帆模拟 光帆物理学 在此模拟中,我们假设最简单的航行方式是:可以控制其姿态的单块平板。 我们将$ \ vec {r} $定义为位置,将$ \ vec {v} $定义为速度,将$ \ hat {p} $定义为帆的指向。 让我们还将在$ r = 1AU $处面向太阳$(\ theta = 0)$的帆的加速度定义为$ a_E $。 我们可以看到帆的加速度为:$$ \ frac {d ^ 2 \ vec {r}} {dt ^ 2} = \ vec {a} = \ left(\ frac {a_Er_E ^ 2} {r ^ 2 } \ cos {\ theta} \ right)\ hat {p} $$ 在哪里 $$ \ cos {\ theta} = \ frac {\ vec {r} \ cdot \ hat {p}} {r} $$ 然后,我们可以使用scipy.integrate.