塔羅異次元

 找回密碼
 註冊
搜索
查看: 5758|回復: 11

[占星學] 老師請收一下email,另外格里曆與儒略曆問題

[複製鏈接]
發表於 2007-10-3 22:08:43 | 顯示全部樓層 |閱讀模式
老師請收一下email,有占星問題想問,另外格里曆與儒略曆問題

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

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

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

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

為何是如此啊?

謝謝

很想上老師的歷法課....不過要看運氣了....看有沒有放假
發表於 2007-10-6 00:51:38 | 顯示全部樓層
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年(公元一年之前一年,就是公元前一年)
發表於 2007-10-6 01:04:36 | 顯示全部樓層
回應來信:

四、修正月球的攝動

首先計算:(參考軌道要件)
    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地球半徑,這裡皆省略掉了。




發表於 2007-10-6 01:17:29 | 顯示全部樓層
回應來信:

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


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


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

2000年一月一日經度0度0時=2451544.5000(儒略日)

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


不知馮略目前是著手以儒略日為準寫程式嗎?
還是使用d,儒略日僅顯示用?
發表於 2007-10-6 01:30:54 | 顯示全部樓層
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的

小數點位數真多


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

我給的程式,並不是最精準的萬年使用法
如果要無限時間延伸,這些公式仍有不足
發表於 2007-10-6 02:33:58 | 顯示全部樓層
對了,我換了這個信箱
astrara.0@gmail.com

請以後改用這個信箱聯絡
比較不會有問題
 樓主| 發表於 2007-10-6 12:45:58 | 顯示全部樓層

回覆 #4 星宿一 的文章

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

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

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

占星程式真是辛苦.....重點不是程式....是天文、物理、曆法觀念....
 樓主| 發表於 2007-10-6 12:51:25 | 顯示全部樓層
>程式有人用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 12:53:11 | 顯示全部樓層
原文為 星宿一 於 2007-10-6 01:17 發表



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


是使用d...儒略日僅顯示用
 樓主| 發表於 2007-10-6 12:55:42 | 顯示全部樓層
>>馮略真的用我給的公式來寫?
>>還是另外有方法?

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

黃赤交角我是列來比較看看...這是網上天文用matlab寫的程式...寫法與天文算法一書一樣
發表於 2007-10-6 17:15:21 | 顯示全部樓層
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日存在
比較感覺不出中斷
西曆日期有接上
但儒略日和星體度數就得跳了
 樓主| 發表於 2007-10-6 22:15:35 | 顯示全部樓層
曆法真是中西命術的大問題啊....還是塔羅輕鬆多了.....
您需要登錄後才可以回帖 登錄 | 註冊

本版積分規則

小黑屋|神祕網域

GMT+8, 2019-6-26 20:25 , Processed in 0.019887 second(s), 17 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回復 返回頂部 返回列表