黃赤交角（黃赤大距）算法

[丁易]

黃赤交角精確的算法，與太陽到山，日影測向等有著密切的關係，下列程式算法給需要的人做參考。
黃赤交角的常數是23度26分21.448秒（黃赤大距）若欲求得精確之數，可用下列公式計算：
ε黃赤交角＝23°27′8〞.26－46〞.845Ｔ－0〞.0059Ｔ2＋0〞.181Ｔ3
Ｔ為1900年1月0.5日後的儒略周日紀，如1991年12月22日的儒略周日數是2448613。
2448613－2415020＝33593（日）
33593÷36525＝0.919835紀（Ｔ）
則⊿ε（黃赤交角）＝46.845－0.059Ｔ2＋0.0018Ｔ3＝43.093366弧秒
43.093366弧秒÷3600＝0.011970（度）
23°27′8〞.26－0.011970（度）＝23°26′25〞.1
黃赤交角的數值每年都在變化。例如2002 年是23.4390（23°26′20.4〞） 度，到了2010 年是23.4380（23°26′16.8〞） 度。總的來說，其變化是週期性的。約37158 年為一週期，變化幅度約為1.612 度。2008 年到達中間點。

[yfsum]

:emm12: :emm22:

[bbq2005]

Thanks for sharing.

[talang]

[ 這篇文章最後由talang在 2005/07/21 08:25pm 第 1 次編輯 ]

谢谢丁易兄提供计算方法，请教儒略周日纪的资料何处能找到或购买，请兄长赐教，十分感谢！
另贴中有几处不明，请指教，谢谢！
1、33593÷36525＝0.919835紀（Ｔ）按计算应等于：0.91972621492 所得数值丁易兄有点偏差。
2、⊿ε（黃赤交角）＝46.845－0.059Ｔ2＋0.0018Ｔ3＝43.093366弧秒 所得数值有差。
按丁易兄提供公式应为：
ε黃赤交角＝23°27′8〞.26－46〞.845Ｔ－0〞.0059Ｔ2＋0〞.181Ｔ3
⊿ε（黃赤交角）＝（46.845×0.919835）-（0.059×0.919835×2）+（0.0018×0.919835×3）=42.986097154弧秒 所得数值与丁易兄的也稍有不同。
还请丁易兄指教，是为谢！

[mars15]

[ 這篇文章最後由mars15在 2005/09/14 07:47am 第 1 次編輯 ]

計算「儒略周日」的網站
http://aa.usno.navy.mil/data/docs/JulianDate.html
http://www.physics.sfasu.edu/astro/javascript/julianday.html
http://www.starlightccd.com/walter/bytebook/javascript/jdcalc.htm
http://www.csgnetwork.com/juliandaydate.html
請參考

[ Last edited by yfsum on 2005-10-25 at 08:45 PM ]

[1154248456]

谢谢分享！