Hello 天文

自動翻譯,請以簡體中文原文校核。

日曆上換一天,計算機怎麼知道行星繞到哪裡了?

把一連串有節奏的小變化加起來,得到該時刻真實計算的位置。

行星與太陽系小天體

行星不只沿著一條固定橢圓運動,其他天體的擾動也會影響位置。原書用一組週期項,把這些運動規律整理成可計算的日心經度、緯度和距離。

開啟互動課程 →

先猜一下

日心位置的原點在誰那裡?

查看答案與解釋

太陽 — 日心座標以太陽為原點,還不是望遠鏡要指向的地平座標。

再把體驗和科學連結起來

先明確原點在太陽

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 章

行星與太陽系小天體 · 全部互動課程