日历上换一天，计算机怎样知道行星绕到哪儿了？
把一系列有节奏的小变化加起来，得到该时刻真实计算的位置。
行星不只沿一条固定椭圆运动，其他天体的扰动也会影响位置。原书用一组周期项，把这些运动规律整理成可计算的日心经度、纬度和距离。
打开互动课程 →
先猜一下
日心位置的原点在谁那里？
查看答案与解释
太阳 — 日心坐标以太阳为原点，还不是望远镜要指向的地平坐标。
再把体验和科学连起来
先明确原点在太阳
L、B、R分别是日心黄经、黄纬与日心距离。与地球上看到的方位角不是一回事。
太阳中心为原点，参考日期平均黄道
每个周期项都随时间变化
把原书附录III保留的各项相加。这里八大行星都使用书中截断表，不是完整VSOP87数千项。
每项 A cos(B+Cτ)，角参数为弧度
不同阶次乘不同的时间幂
先求L0、L1等各阶和，再乘τ的幂，最后按10⁸缩放。这里τ用儒略千年，不能误用儒略世纪。
τ=(JDE−2451545)/365250
需要精度时说明参考系统
可切换式32.3的FK5小修正。地心方向、光行时和观测条件属于下一步，不是日心L、B、R本身。
动态黄道 → FK5：加入小角修正
把计算链展开
τ = (JDE − 2451545)/365250 Sⱼ(τ) = Σᵢ Aⱼᵢ cos(Bⱼᵢ + Cⱼᵢτ) L弧度 = (S₀ + S₁τ + S₂τ² + S₃τ³ + S₄τ⁴ + S₅τ⁵)/10⁸ L度 = L弧度 × 180/π
B与R按各自保留的阶数作同样求和；B也是弧度，R为AU。A、B、C逐项取本课附录III表，内部相角用弧度。这里展开求和结构，没有逐行列出2430个系数项，也不是完整VSOP87。
试一试
先使用金星1992年书例，再换到其他行星；画面上的点与下面数值来自相同的周期项求和。
动手看一遍 →
符号不用背
- L / B
- 日心黄经与黄纬，度
- R
- 日心距离，AU
- τ
- J2000起儒略千年
算法、精度与适用范围
完整求和原书附录III中八大行星所有保留L/B/R周期项，合计2430项。L/B内部弧度，输出度；R为AU，幅系数除10⁸。默认参考日期平均动力学黄道与分点；FK5模式使用式32.3修正。只在1900–2100的教学日期范围使用。没有把书中截断表说成完整VSOP87，没有自动调用首页JPL模型；本课不含地心换算、光行时、光行差、章动、大气折射或实际观测排程。
我的理解检查
这里的τ单位是？
儒略千年 — 每365250天为一个τ单位。
书中附录III是完整VSOP87理论吗？
不是，是截断项 — 本课完整使用书中保留项，但并非原始完整VSOP87系列。
日心黄经可以直接当本地罗盘方位角吗？
不可以 — 原点和参考平面不同，需要额外坐标转换。
日期输入要正确，级数单位要正确，坐标原点也要正确；三个条件一起，位置数字才有意义。
Jean Meeus · Astronomical Algorithms, 2nd ed. · 第 32 章