许剑伟先生(福建莆田第十中学) 超好用“寿星”天文万年历.doc

許劍偉先生(福建莆田第十中學) 超好用“壽星”天文萬年曆這是一款真正精確的萬年曆，它提供西元-4712年到西元9999年的西曆及農曆日期查詢功能，其中-721年到1960年的農曆資料已經與張培瑜的三千五百年曆日天象、陳垣的二十史朔閏表、方詩銘的中國史曆日和中西曆日對照表核對；它還含有西元前2000多年以前到今的基本年號；含有二千多個國內城市的經緯度，並且使用者可根據自已的需要擴展經緯度數據。 這是一款採用現代天文演算法製作的天文、曆算程式，可以方便地進行西曆、農曆、回曆三曆之間的轉換；提供了曆譜資料匯出功能；提供精確的日月食等天象的計算功能；含有行星、恒星星曆表計算。 本軟體在最近的幾百年範圍內節氣時刻及日月合朔時刻的平均誤差小於1秒。如果以力學時作為時間尺規，與DE406星曆表比對，太陽座標的最大可能誤差為0.1角秒，月亮座標的最大可能誤差為0.8角秒，平均誤差為最大可能誤差的1/6，這使得精確的日月食計算成為可能，是目前世界上唯一一款整合了精確日月食過程計算的萬年曆。 它還提供日出日沒、月出月沒時間查詢等功能，與國內外著名的天文軟體或天文年曆比較，結果完全一致。此外，天文曆中還即時顯示了“月亮被照亮部分比例”、太陽和月亮的地心黃道座標、地心赤道座標、站心地平座標等。其中日月的地心黃道、赤道座標採用高精度演算法，“月亮被照亮部分比例”的計算採用低精度的演算法，但可以滿足一般需求；地平座標的計算已適當考慮了大氣折射，由於大氣的氣壓及氣溫變化具有許多不確定的因素，所以只考慮平均情況下的大氣折射修正，在本軟體中，地平座標被描述為方位角與高度角，它是以觀測點為中心的座標，所以已經考慮了周日視差修正。螢幕顯示的方位角從正南向西測量。方位角與高度角已轉換到是站心地平座標，並且做了視差修正，同時在地平真緯度大於0時進行大氣折射修正。 本軟體是一款精准的年代跨度大的日曆工具，可作為一般的實用日曆工具，對史學家、考古學家、曆算工作者、天文愛好者均有較大的參考價值。農曆核心演算法類庫（含源碼），源自許劍偉先生的壽星萬年曆(v4.11)，由 Javascript 源碼改寫為 C# 源碼，在 Visual Studio 2008 中測試通過。強烈推薦您下載一份壽星萬年曆的源碼進行對照參考。 真正的“萬年”曆，提供西元 -4712 年到西元 9999年 的日期查詢功能，可讓您告別查表演算法。版權問題 本程式是開源的，你可以使用其中的任意部分代碼，但不得隨意修改“天文演算法(eph.js)”及“農曆演算法(lunar.js)中古曆部分的資料及演算法”。一旦修改可能影響萬年曆的準確性，如果你對天文學不太瞭解而僅憑對曆法的熱情，請不要對此做任何修改，以免弄巧成拙。 如果在你自己開發的軟體中使用了本程式的核心演算法及資料，你可以在你的軟體中申明“資料或演算法來源於本天文曆”，也可以不申明，但不可以申明為它其它來源。如有異義，可與我共內探討。 天文演算法：進行月曆推算時採用了先進的現代天文演算法，它是基於2002年巴黎天文臺西爾特處Jean CHAPRONT 和 Gerard FRANCOU發表的ELP/MPP02月球運動理論，在幾百年範圍內，該理論與當今世界上公認的最精密的LE405/406系列星曆表僅相差3毫角秒，在幾千年範圍相差不超過4角秒。 地標數據：城市經緯度資料庫做了壓縮處理。要擴展經緯資料庫,可使用附件中的【城市經緯壓縮器】生成代碼並替換程式中的JWv陣列。程式中現有的經緯度數據本來只有幾百個，後來網友鄭彬給我一份3000多個城市的資料，比較全面，因此本版使用他的資料。 曆譜說明：1645年農曆八月開始使用時憲曆，七月及之前為授時曆 日食月食：天文曆中提供的圖表，可以精確的確定日月食情況。計算結果與2008年中國天文年曆比對，日食誤差在1秒左右，月食誤差在3秒以內， 日月出沒：受到各地的地形地貌、大氣狀態、海拔高度等的影響，日出日沒、月出月沒時間無法計算得很準確，通常只能精確到1分鐘，因此，軟體中時間顯示精確到到了秒數量級，但這並不表明日出日沒時間的計算達到了這個精度。日月中天時刻受大氣等因素的影響要小得多，所以可以精確到秒。日月的“出、中天、沒”的計算不需要高精度日月位置座標(因為一天中，日月在視野中的位置主要由地球自轉決定)，程式中使用非常低精度的方法計算日月座標仍可以把中天時刻精確到秒。 時區數據：各地時區資料庫做了小量壓縮處理。如要擴展，可使用附件中的【時區合成器】生成代碼並替換程式中的SQv陣列。時區資料庫是“中華農曆網”netghost先生提供。 年號紀年：天文曆中的紀年資料參考萬國鼎中國歷史紀年表，有些朝代存在多組帝王，程式中不一定全部納入。 命理八字：滑鼠經過日曆表，顯示普通的八字資訊。螢幕下方的八字是嚴格的，以立春時刻界定年，以節氣時刻界定月，以本地真太陽時23點界定日(這裡的真太陽時使用低精度演算法，誤差可能達到1秒)，以真太陽時23點、1點、3點界定時辰。注意，這時螢幕上方輸入的時間為當地時間(電腦時間，UTC)。 月球視座標比對 V4.13版本及以後 時間 2008.1.6 00:00:00 TT (JED = 2926.5+2451545) 本程式 視黃經 2565436.31 視黃緯 -45214.12 距離 401817.73千米 DE406 視黃經 2565436.319 視黃緯 -45214.134 距離 401817.6711 中國 視黃經 25654 36.32 視黃緯: -45214.09 2008年中國天文年曆它與JPL不一致 swiss 視黃經: 2565436.319 視黃緯 -45214.090 距離* 401798.6263 本程式 視赤經 17h 00m58.06s 視赤緯 -273833.64 距離 401817.73千米 DE406 視赤經 17h 00m58.061s 視赤緯 -273833.691 距離:401817.6711 中國 視赤經 17h 00m58.061s 視赤緯:-273833.69 2008年中國天文年曆它與JPL一致 swiss 視赤經 17h 00m58.061s 視赤緯 -273833.648 距離*:401798.6263 JPL 視赤經 17h 00m58.0575s視赤緯:-273833.689 距離*:401798.6270 JPL網站查詢(url=/horizons.cgi/horizons.cgi/url) 比較的最後結果： (1)“ELP/mpp02”、“個人利用DE406計算”、“JPL網站查詢”三者結果基本一致。 存在幾個毫角秒的差異主要是歲差參數以及黃赤交角參數的選用值不同造成的。 以上三者座標都在J2000慣性黃道&赤道坐標系中計算並比對。 (2)中國天文年曆的月亮視赤經、視赤緯與JPL的一致。 (3)中國天文年曆的月亮視黃緯與JPL不一致,它採用的黃道與J2000慣性黃道存在0.047角秒的夾角。 (4)2008中國天文年曆的月亮座標參考系存在前後不統一的現象！ (5)swiss的月亮視黃緯、視赤緯與JPL不一致，它採用的黃道與J2000慣性黃道存在0.047角秒的夾角。 (6)swiss與2008中國天文年曆的月亮視黃緯一致，其黃道與J2000慣性黃道存在0.047角秒的夾角。 (7)上面視座標中的“距離*”與“距離”的含義是不相同的。設月亮光行時為T,前者是t-T時刻月球 與t時刻地球的距離,後者為t-T時刻地月距。 (8)此刻本程式誤差為幾十毫角秒 月球黃經計算結果與2008年中國天文年曆等權威資料的比較 V4.12版本及以前 2008年01月01日0h TD +19719 24.43 中國天文年曆 +1971924.91 本程式 2008年01月06日0h TD +25654 36.32 中國天文年曆 +25654 36.11 本程式 2008年01月18日0h TD + 5604 29.83 中國天文年曆 + 5604 29.68 本程式 2100年01月01日0h TD +15724 01.183 瑞士星曆表 +15724 01.96 本程式 2100年01月18日0h TD + 2214 39.400 瑞士星曆表 + 2214 40.47 本程式 2200年01月02日0h TD +10826 45.916 瑞士星曆表+10826 46.12 本程式 月球方位角與高度角 網站座標: L=-11623 =+ 3954 日期時間: 2000年1月1日12h TD(即力學時=0) SkyMap 方位角 34813 16 高度角 -6058 19 本程式 方位角 16813 14 高度角 -6058 19新增內容已將古代真實日曆推算到西元前104年（與實曆相符），西元前104年以前及196