|
发表于 2017-1-15 09:06:28
|
显示全部楼层
抱歉,鄙人回覆怠慢了。
個人在天文計算也經歷幾個階段。
1. NOAA 的 excel 模板,在某些計算仍然是近似值,例如: EOT 的計算即是。可惜,NOAA 的計算主要在太陽的黃赤、地平、仰角。不能解決中國七政四餘推步、天星計算及回溯古代的天象驗證。
2. Jean Meeus 先生的天文計算書籍 Astronomical Algorithms 是必讀的經典,該書從1982年初版、1985、1988、1990, 已經到第4版。但鄙人發現其中有錯誤,直到 1991 出版,才更正書中的一些錯誤。若您個人想看此書,請選購 1991 年以後版本。鄙人是購買其英文原版書及電子書。國內專家許劍偉先生有把該書做重點翻譯《天文算法》。
* Jean Meeus 先生的天文計算( Astronomical Algorithms ) 書籍提供行星模型及計算公式,是非常好的進階書籍。鄙人乃逐一按其書中算法,撰寫數萬行的計算機程式,將其整本書的算法全部用軟件完成。這是一個苦功夫,但也非常值得。試想,在中國術數界,有人會花時間去閱讀Jean Meeus 先生英文原版的《天文算法》? 除了許劍偉先生外,在中國術數界是否有人將Jean Meeus 先生英文原版的《天文算法》用軟件實現?
* Jean Meeus 先生的天文計算( Astronomical Algorithms ) 書中的行星模型及公式,是採用 VSOP 行星理論。VSOP是"Variations Seculaires des Orbites Planetaires"的縮寫,由巴黎的P.Bretagnon 於1982年,發表了他的行星理論VSOP82。VSOP82由大星行(水星到海王星)的長長的週期項序列組成。給定一個行星及一個時間,對它的序列取和計算,即可獲得密切軌道參數。
Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中採用 1987 年的模型, 即 VSOP87。
VSOP82 行星理論是採用複雜的微積分及上千個多項式週期項,當不需要完全精度時,可以把計算公式在適當處截斷。幸運的是,1987年 Bretagnon 和 Francou 創建了VSOP87,它提供了直接計算行星日心座標的週期序列項。也就是說可以直接計算出:
L, 日心黃經
B, 日心黃緯
R, 行星到太陽的距離
Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中附錄中有列出七政四餘的 VSOP 87 的簡化週期項數值。
但是如果您按 Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中附錄中列出七政四餘的 VSOP 87 的簡化週期項數值, 去計算七政四餘三王的天文位置,其精度仍然是不夠的。因為原本 VSOP87 的行星模式,是把微積分公式,用數千個週期項展開,而Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中附錄中列出七政四餘的 VSOP 87 的簡化週期項數量, 最多只有 64-100 之間。
以下列出个人所使用 VSOP87 行星模型中【地球】(也就是計算太陽位置) 的完整周期項的函式及公式:
(而這僅是計算太陽行星位置的最基礎的數據。)
|
|