老師請收一下email，另外格里曆與儒略曆問題

馮略

老師請收一下email，有占星問題想問，另外格里曆與儒略曆問題

在寫程式時，只要日期1582年後的日期，計算出的儒略日就沒問題

只要在1582年前就會有問題....這是因為公曆被教皇調整

可是網路上好像是1582年10月5號後去除了10天...

但我看到的程式有人用1582年10月25號12點來判斷耶

為何是如此啊？

謝謝

很想上老師的歷法課....不過要看運氣了....看有沒有放假

星宿老師

1582年10月5號後去除了10天是沒錯的
5日，6日，7日，8日，9日，10日，11日，12日，13日，14日
以上這十天被跳過
10月4日後接10月15日

>程式有人用1582年10月25號12點來判斷
不知意指？
如果是接4日，應該是弄錯了
（哪個程式呢？）


另外，提醒，還有公元零年和元年、-1年，與儒略日的替換的問題要注意
這個爭議比較大
因為數學上是有公元0年的
但是歷史習慣用法卻不太用公元0年（公元一年之前一年，就是公元前一年）

星宿老師

回應來信：

四、修正月球的攝動

首先計算：（參考軌道要件）
    Ms                 太陽之意義均近點角
    Mm                 月亮之意義均近點角
    Nm                 月球升交點之黃道經度
    ws                 太陽之近日點差額
    wm                 月球之近地點差額

並沒有相關資料，所以算不下去，請老師幫忙

可能是我沒有表明清楚吧
Ms就是太陽sun的M（平均近點角）
Mm就是月亮moon的M（平均近點角）

這些從前面的公式就可算出
↓
太陽之軌道要件：
    N = 0.0
    i = 0.0
    w = 282.9404 + 4.70935E-5 * d
    a = 1.000000  (AU)
    e = 0.016709 - 1.151E-9 * d
    M = 356.0470 + 0.9856002585 * d

月亮之軌道要件：
    N = 125.1228 - 0.0529538083 * d
    i = 5.1454
    w = 318.0634 + 0.1643573223 * d
    a = 60.2666  (地球半徑倍數)
    e = 0.054900
    M = 115.3654 + 13.0649929509 * d

只要算出d就會得出數值
就能帶入
    Ls = Ms + ws （+ Ns）       太陽平均經度  (太陽升交點永等於0，Ns=0)
    Lm = Mm + wm + Nm           月球平均經度
    D = Lm - Ls                 月相平均伸縮值
    F = Lm - Nm                 月球緯度差額


加上這些項目數值到月亮經度上（以degrees為單位）
    -1.274 * sin(Mm - 2*D)          (月球出差the Evection)
    +0.658 * sin(2*D)               (運行變差the Variation)
    -0.186 * sin(Ms)                (年均差the Yearly Equation)
    -0.059 * sin(2*Mm - 2*D)
    -0.057 * sin(Mm - 2*D + Ms)
    +0.053 * sin(Mm + 2*D)
    +0.046 * sin(2*D - Ms)
    +0.041 * sin(Mm - Ms)
    -0.035 * sin(D)                 (視差the Parallactic Equation)
    -0.031 * sin(Mm + Ms)
    -0.015 * sin(2*F - 2*D)
    +0.011 * sin(Mm - 4*D)

加上這些項目數值到月亮緯度上（以degrees為單位）
    -0.173 * sin(F - 2*D)
    -0.055 * sin(Mm - F - 2*D)
    -0.046 * sin(Mm + F - 2*D)
    +0.033 * sin(F + 2*D)
    +0.017 * sin(2*Mm + F)


加上這些項目數值到月亮距離上（以地球半徑為單位Earth radii）
    -0.58 * cos(Mm - 2*D)
    -0.46 * cos(2*D)


所有的攝動項目，經緯度小於0.01 degrees，距離小於0.1地球半徑，這裡皆省略掉了。

星宿老師

回應來信：

d我也有寫好了，但老師的算法為何沒有計算儒略日???

當初的d沒有換算成儒略日
因為儒略日很龐大（要跨越無限日期和講究經緯度時差）
並且也不知道真的開始著手寫程式了
（而且儒略日使用習慣也不同，如果不要求很完整，有的人會用簡化儒略）
真的要寫完整的程式，的確需要用儒略日
（像是面對上面換曆的問題）


用那個d公式就可以直接轉換成儒略日

2000年一月一日經度0度0時＝2451544.5000（儒略日）

也可以先計算儒略日再轉成d套入公式，而不先算d
或是依據轉換原理去更動公式來配合


不知馮略目前是著手以儒略日為準寫程式嗎？
還是使用d，儒略日僅顯示用？

星宿老師

2000年1月1日0時

日期符號d:::1
黃道傾斜度ecl:::23.4392996437
儒略日:::2451544.5
黃赤交角【以秒為單位】:::84381.406173
黃赤交角【以度為單位】:::23.4392794925

===太陽開頭====
M:::6.23139441146758
w:::4.93824238884639
非正圓近點角E:::6.23051499057306
xv:::0.981904240661308
yv:::-0.0526386176731259
太陽的真近點角v:::-0.053557441082655
太陽的距離r:::0.983314172530426
太陽的真正經度lonsun【弧度】:::4.88468494776374
太陽的真正經度lonsun【角度】:::279.871831757943
地球中心主軸黃道體系座標xs:::0.16858407375364
地球中心主軸黃道體系座標ys:::-0.968754959716761
===太陽結尾====


===月亮開頭====
N:::2.18288061100806
i:::0.0898041713321162
M:::2.24153321662332
w:::5.55412213651163
非正圓近點角E:::2.28307233971804
月亮的xv:::-42.6963622936683
月亮的yv:::45.5455701381712
月亮的v:::2.32391711100731
月亮的r:::62.4289861548564
xh:::-50.0116531182411
yh:::-36.9443348145874
zh:::5.59723054837338
lonecl:::-2.505353456 76472
latecl:::0.0897781146061011
Math.sqrt(xh*xh+yh*yh+zh*zh):::62.4289861548564
lonecl【角度】:::-143.54617926113
latecl【角度】:::5.1439070595714
===月亮結尾====

......我只看了幾條，大致是OK的

小數點位數真多


馮略真的用我給的公式來寫？
還是另外有方法？

我給的程式，並不是最精準的萬年使用法
如果要無限時間延伸，這些公式仍有不足

星宿老師

對了，我換了這個信箱
astrara.0@gmail.com

請以後改用這個信箱聯絡
比較不會有問題

馮略

因為我看了很多程式之後....光是儒略日就有好多種計算

我也找到了天文算法一書....觀看裏面的範例....

只要計算1582年....有些程式就會掛掉....出現儒略日偏差....

占星程式真是辛苦.....重點不是程式....是天文、物理、曆法觀念....

馮略

>程式有人用1582年10月25號12點來判斷

它用2299171儒略日為判斷點，大於此數為格里曆

用天文算法一書的範例，兩者比對，算出的數字都對

這個程式是
http://destiny.xfiles.to/html/astrology/ahelper/ahelper.html

大陸有人改寫成flash
http://www.12sign.cn/htmls/c.htm

但是沒有進行相關歲差修正，所以三王星差很多....

馮略

原文為 星宿一 於 2007-10-6 01:17 發表

>>不知馮略目前是著手以儒略日為準寫程式嗎？
>>還是使用d，儒略日僅顯示用？


是使用d...儒略日僅顯示用

馮略

>>馮略真的用我給的公式來寫？
>>還是另外有方法？

表裏的都用老師的公式算...除了儒略日與黃赤交角

黃赤交角我是列來比較看看...這是網上天文用matlab寫的程式...寫法與天文算法一書一樣

星宿老師

1582年10/25這天是2299171沒錯
往後全都是格里曆，所以不會錯
但是往前則要重新設定了
所以等於沒有涉及、解決換曆問題

所給之程式計算範圍年份是有限的

有的軟體是採取1582該年
9/25=10/5
9/26=10/6
9/27=10/7
9/28=10/8
9/29=10/9
9/30=10/10
10/1=10/11
10/2=10/12
10/3=10/13
10/4=10/14
10/5=10/15
這樣的呈現方式
看起來好像5~14日存在
比較感覺不出中斷
西曆日期有接上
但儒略日和星體度數就得跳了

馮略

曆法真是中西命術的大問題啊....還是塔羅輕鬆多了.....