旧暦の計算《計算式編》

本ページは旧暦の計算のサブページです。

1873(明治6)年のグレゴリオ暦への改暦まで日本で使われていた旧暦
本ページでは、朔や二十四節気の時刻を計算するのに必要な諸関数の計算式について記載する。

目次

旧暦の計算

儀鳳暦

月離による補正値

儀鳳暦における月離による補正値 $\DLunar(t)$ は、近点月 $\PLunar = 27\dfrac{743\frac{1}{12}}{1340}\,\day$を周期とする周期関数であり、その具体的な形は下図の折れ線で与えられる。

儀鳳暦における月離による補正値
儀鳳暦における月離による補正値 $\DLunar(t)$

折れ線の頂点の値は下表の通りである。
五紀暦においても儀鳳暦とほぼ同じであるから、併せて掲載する。

儀鳳暦・五紀暦における月離による補正値
近地点からの経過日数
$t$
補正値(儀鳳暦)
$\DLunar(t)$
補正値(五紀暦)
$\DLunar(t)$
近地点 $0\,\day$ $0\,\day$ $0\,\day$
$1\,\day$ $+\dfrac{134}{1340}\,\day$ $+\dfrac{135}{1340}\,\day$
$2\,\day$ $+\dfrac{251}{1340}\,\day$ $+\dfrac{252}{1340}\,\day$
$3\,\day$ $+\dfrac{350}{1340}\,\day$ $+\dfrac{351}{1340}\,\day$
$4\,\day$ $+\dfrac{428}{1340}\,\day$ $+\dfrac{429}{1340}\,\day$
$5\,\day$ $+\dfrac{484}{1340}\,\day$ $+\dfrac{485}{1340}\,\day$
$6\,\day$ $+\dfrac{517}{1340}\,\day$ $+\dfrac{518}{1340}\,\day$
$6\dfrac{1190}{1340}\,\day$ $+\dfrac{526}{1340}\,\day$ $+\dfrac{526}{1340}\,\day$
$7\,\day$ $+\dfrac{526}{1340}\,\day$ $+\dfrac{525}{1340}\,\day$
$8\,\day$ $+\dfrac{512}{1340}\,\day$ $+\dfrac{511}{1340}\,\day$
$9\,\day$ $+\dfrac{474}{1340}\,\day$ $+\dfrac{473}{1340}\,\day$
$10\,\day$ $+\dfrac{412}{1340}\,\day$ $+\dfrac{411}{1340}\,\day$
$11\,\day$ $+\dfrac{327}{1340}\,\day$ $+\dfrac{326}{1340}\,\day$
$12\,\day$ $+\dfrac{223}{1340}\,\day$ $+\dfrac{223}{1340}\,\day$
$13\,\day$ $+\dfrac{102}{1340}\,\day$ $+\dfrac{105}{1340}\,\day$
遠地点 $13\dfrac{1041}{1340}\,\day$ $0\,\day$ $0\,\day$
$14\,\day$ $-\dfrac{29}{1340}\,\day$ $-\dfrac{30}{1340}\,\day$
$15\,\day$ $-\dfrac{157}{1340}\,\day$ $-\dfrac{158}{1340}\,\day$
$16\,\day$ $-\dfrac{272}{1340}\,\day$ $-\dfrac{273}{1340}\,\day$
$17\,\day$ $-\dfrac{367}{1340}\,\day$ $-\dfrac{368}{1340}\,\day$
$18\,\day$ $-\dfrac{441}{1340}\,\day$ $-\dfrac{442}{1340}\,\day$
$19\,\day$ $-\dfrac{493}{1340}\,\day$ $-\dfrac{494}{1340}\,\day$
$20\,\day$ $-\dfrac{521}{1340}\,\day$ $-\dfrac{522}{1340}\,\day$
$20\dfrac{892}{1340}\,\day$ $-\dfrac{525}{1340}\,\day$ $-\dfrac{528}{1340}\,\day$
$21\,\day$ $-\dfrac{525}{1340}\,\day$ $-\dfrac{525}{1340}\,\day$
$22\,\day$ $-\dfrac{505}{1340}\,\day$ $-\dfrac{505}{1340}\,\day$
$23\,\day$ $-\dfrac{461}{1340}\,\day$ $-\dfrac{463}{1340}\,\day$
$24\,\day$ $-\dfrac{393}{1340}\,\day$ $-\dfrac{398}{1340}\,\day$
$25\,\day$ $-\dfrac{304}{1340}\,\day$ $-\dfrac{309}{1340}\,\day$
$26\,\day$ $-\dfrac{196}{1340}\,\day$ $-\dfrac{200}{1340}\,\day$
$27\,\day$ $-\dfrac{71}{1340}\,\day$ $-\dfrac{75}{1340}\,\day$
$27\dfrac{743}{1340}\,\day$ $0\,\day$ $0\,\day$
近地点 $27\dfrac{743\frac{1}{12}}{1340}\,\day$ $0\,\day$ $0\,\day$

近地点/遠地点の通過時刻と、補正値が極大/極小になる時刻を例外として、基本的には1日毎に頂点が打たれている。
実際の計算に当たっては、$\dfrac{1}{1340}\,\day$(儀鳳暦における1分)を単位として丸めて計算すればよい。

日躔による補正値

儀鳳暦における日躔による補正値 $\DSolar(t)$ は、近点年 $\PSolar = 365\dfrac{328}{1340}\,\day$を周期とする周期関数であり、その具体的な形は下図に示す通り24個の2次関数(放物線)を繋ぎ合わせた形である。

儀鳳暦における日躔による補正値
儀鳳暦における日躔による補正値 $\DSolar(t)$

2次関数は、各区間 $[ \tau_s, \tau_e )$ 毎に次式の形で与えられる。

\[ \DSolar(t) = \bigg[ a + b(t - \tau_s) + \dfrac{1}{2}c(t - \tau_s)(t - \tau_s - 1\,\day) \bigg] \cdot \dfrac{1}{1340}\,\day \]

区間の端点 $\tau_s, \tau_e$ 及び係数 $a, b, c$ の値は下表の通りである。

儀鳳暦における日躔による補正値
区間の始点 $\tau_s$
(近日点からの経過日数)
区間の終点 $\tau_e$
(近日点からの経過日数)
初日朓朒ちょうじく
$a$
初日損益
$b$
損益の変化率
$c$
$ 0 \,\day$(定気冬至)$ 14\dfrac{ 911}{1340}\,\day$(定気小寒)$ 0$$+3.9546/\day$$-0.0372/\day^2$
$ 14\dfrac{ 911}{1340}\,\day$(定気小寒)$ 29\dfrac{ 586}{1340}\,\day$(定気大寒)$ +54$$+3.4091/\day$$-0.0372/\day^2$
$ 29\dfrac{ 586}{1340}\,\day$(定気大寒)$ 44\dfrac{ 365}{1340}\,\day$(定気立春)$+100$$+2.8636/\day$$-0.0372/\day^2$
$ 44\dfrac{ 365}{1340}\,\day$(定気立春)$ 59\dfrac{ 144}{1340}\,\day$(定気雨水)$+138$$+2.3181/\day$$+0.0372/\day^2$
$ 59\dfrac{ 144}{1340}\,\day$(定気雨水)$ 73\dfrac{1159}{1340}\,\day$(定気啓蟄)$+176$$+2.8636/\day$$+0.0372/\day^2$
$ 73\dfrac{1159}{1340}\,\day$(定気啓蟄)$ 88\dfrac{ 729}{1340}\,\day$(定気春分)$+222$$+3.4091/\day$$+0.0372/\day^2$
$ 88\dfrac{ 729}{1340}\,\day$(定気春分)$104\dfrac{ 404}{1340}\,\day$(定気清明)$+276$$-3.7220/\day$$+0.0329/\day^2$
$104\dfrac{ 404}{1340}\,\day$(定気清明)$119\dfrac{1315}{1340}\,\day$(定気穀雨)$+222$$-3.2086/\day$$+0.0329/\day^2$
$119\dfrac{1315}{1340}\,\day$(定気穀雨)$135\dfrac{ 782}{1340}\,\day$(定気立夏)$+176$$-2.6952/\day$$+0.0329/\day^2$
$135\dfrac{ 782}{1340}\,\day$(定気立夏)$151\dfrac{ 249}{1340}\,\day$(定気小満)$+138$$-2.1818/\day$$-0.0329/\day^2$
$151\dfrac{ 249}{1340}\,\day$(定気小満)$166\dfrac{1160}{1340}\,\day$(定気芒種)$+100$$-2.6952/\day$$-0.0329/\day^2$
$166\dfrac{1160}{1340}\,\day$(定気芒種)$182\dfrac{ 834}{1340}\,\day$(定気夏至)$ +54$$-3.2086/\day$$-0.0329/\day^2$
$182\dfrac{ 834}{1340}\,\day$(定気夏至)$198\dfrac{ 509}{1340}\,\day$(定気小暑)$ 0$$-3.7220/\day$$+0.0329/\day^2$
$198\dfrac{ 509}{1340}\,\day$(定気小暑)$214\dfrac{ 80}{1340}\,\day$(定気大暑)$ -54$$-3.2086/\day$$+0.0329/\day^2$
$214\dfrac{ 80}{1340}\,\day$(定気大暑)$229\dfrac{ 887}{1340}\,\day$(定気立秋)$-100$$-2.6952/\day$$+0.0329/\day^2$
$229\dfrac{ 887}{1340}\,\day$(定気立秋)$245\dfrac{ 354}{1340}\,\day$(定気処暑)$-138$$-2.1818/\day$$-0.0329/\day^2$
$245\dfrac{ 354}{1340}\,\day$(定気処暑)$260\dfrac{1265}{1340}\,\day$(定気白露)$-176$$-2.6952/\day$$-0.0329/\day^2$
$260\dfrac{1265}{1340}\,\day$(定気白露)$276\dfrac{ 939}{1340}\,\day$(定気秋分)$-222$$-3.2086/\day$$-0.0329/\day^2$
$276\dfrac{ 939}{1340}\,\day$(定気秋分)$291\dfrac{ 510}{1340}\,\day$(定気寒露)$-276$$+3.9546/\day$$-0.0372/\day^2$
$291\dfrac{ 510}{1340}\,\day$(定気寒露)$306\dfrac{ 185}{1340}\,\day$(定気霜降)$-222$$+3.4091/\day$$-0.0372/\day^2$
$306\dfrac{ 185}{1340}\,\day$(定気霜降)$320\dfrac{1304}{1340}\,\day$(定気立冬)$-176$$+2.8636/\day$$-0.0372/\day^2$
$320\dfrac{1304}{1340}\,\day$(定気立冬)$335\dfrac{1083}{1340}\,\day$(定気小雪)$-138$$+2.3181/\day$$+0.0372/\day^2$
$335\dfrac{1083}{1340}\,\day$(定気小雪)$350\dfrac{ 758}{1340}\,\day$(定気大雪)$-100$$+2.8636/\day$$+0.0372/\day^2$
$350\dfrac{ 758}{1340}\,\day$(定気大雪)$365\dfrac{ 328}{1340}\,\day$(定気冬至)$ -54$$+3.4091/\day$$+0.0372/\day^2$

日躔による補正値も、月離の場合と同じく $\dfrac{1}{1340}\,\day$(儀鳳暦における1分)を単位として丸めて計算すればよい。

大衍暦

月離による補正値

大衍暦における月離による補正値 $\DLunar(t)$ は、近点月 $\PLunar = 27\dfrac{1685\frac{79}{80}}{3040}\,\day$を周期とする周期関数であり、その具体的な形は下図の折れ線で与えられる。
※近地点ではなく遠地点を基準に定義されているので、儀鳳暦・五紀暦に対して $\dfrac{\PLunar}{2}$ だけズレている。

大衍暦における月離による補正値
大衍暦における月離による補正値 $\DLunar(t)$

折れ線の頂点の値は下表の通りである。

大衍暦における月離による補正値
遠地点からの経過日数
$t$
補正値(大衍暦)
$\DLunar(t)$
遠地点 $0\,\day$ $0\,\day$
$1\,\day$ $-\dfrac{297}{3040}\,\day$
$2\,\day$ $-\dfrac{556}{3040}\,\day$
$3\,\day$ $-\dfrac{776}{3040}\,\day$
$4\,\day$ $-\dfrac{956}{3040}\,\day$
$5\,\day$ $-\dfrac{1095}{3040}\,\day$
$6\,\day$ $-\dfrac{1192}{3040}\,\day$
$6\dfrac{2701}{3040}\,\day$ $-\dfrac{1240}{3040}\,\day$
$7\,\day$ $-\dfrac{1234}{3040}\,\day$
$8\,\day$ $-\dfrac{1170}{3040}\,\day$
$9\,\day$ $-\dfrac{1064}{3040}\,\day$
$10\,\day$ $-\dfrac{916}{3040}\,\day$
$11\,\day$ $-\dfrac{727}{3040}\,\day$
$12\,\day$ $-\dfrac{498}{3040}\,\day$
$13\,\day$ $-\dfrac{231}{3040}\,\day$
近地点 $13\dfrac{2362}{3040}\,\day$ $0\,\day$
$14\,\day$ $+\dfrac{66}{3040}\,\day$
$15\,\day$ $+\dfrac{355}{3040}\,\day$
$16\,\day$ $+\dfrac{605}{3040}\,\day$
$17\,\day$ $+\dfrac{816}{3040}\,\day$
$18\,\day$ $+\dfrac{987}{3040}\,\day$
$19\,\day$ $+\dfrac{1117}{3040}\,\day$
$20\,\day$ $+\dfrac{1204}{3040}\,\day$
$20\dfrac{2024}{3040}\,\day$ $+\dfrac{1240}{3040}\,\day$
$21\,\day$ $+\dfrac{1222}{3040}\,\day$
$22\,\day$ $+\dfrac{1149}{3040}\,\day$
$23\,\day$ $+\dfrac{1033}{3040}\,\day$
$24\,\day$ $+\dfrac{876}{3040}\,\day$
$25\,\day$ $+\dfrac{678}{3040}\,\day$
$26\,\day$ $+\dfrac{441}{3040}\,\day$
$27\,\day$ $+\dfrac{165}{3040}\,\day$
$27\dfrac{1685}{3040}\,\day$ $0\,\day$
遠地点 $27\dfrac{1685\frac{79}{80}}{3040}\,\day$ $0\,\day$

近地点/遠地点の通過時刻と、補正値が極大/極小になる時刻を例外として、基本的には1日毎に頂点が打たれている。
実際の計算に当たっては、$\dfrac{1}{3040}\,\day$(大衍暦における1分)を単位として丸めて計算すればよい。

日躔による補正値

大衍暦における日躔による補正値 $\DSolar(t)$ は、近点年 $\PSolar = 365\dfrac{743}{3040}\,\day$を周期とする周期関数であり、その具体的な形は下図に示す通り24個の2次関数(放物線)を繋ぎ合わせた形である。
なお、五紀暦においてもこれと同じ関数を用いる。

大衍暦・五紀暦における日躔による補正値
大衍暦・五紀暦における日躔による補正値 $\DSolar(t)$

2次関数は、各区間 $[ \tau_s, \tau_e )$ 毎に次式の形で与えられる。

\[ \DSolar(t) = \bigg[ a + b(t - \tau_s) + \dfrac{1}{2}c(t - \tau_s)(t - \tau_s - 1\,\day) \bigg] \cdot \dfrac{1}{3040}\,\day \]

区間の端点 $\tau_s, \tau_e$ 及び係数 $a, b, c$ の値は下表の通りである。

大衍暦・五紀暦における日躔による補正値
区間の始点 $\tau_s$
(近日点からの経過日数)
区間の終点 $\tau_e$
(近日点からの経過日数)
初日朓朒数
$a$
初日損益
$b$
損益の変化率
$c$
$ 0 \,\day$(定気冬至)$ 14\dfrac{1352}{3040}\,\day$(定気小寒)$ 0$$+13.4524/\day$$-0.1886/\day^2$
$ 14\dfrac{1352}{3040}\,\day$(定気小寒)$ 29\dfrac{ 171}{3040}\,\day$(定気大寒)$+176$$+10.5564/\day$$-0.1634/\day^2$
$ 29\dfrac{ 171}{3040}\,\day$(定気大寒)$ 43\dfrac{2485}{3040}\,\day$(定気立春)$+314$$ +8.0408/\day$$-0.1446/\day^2$
$ 43\dfrac{2485}{3040}\,\day$(定気立春)$ 58\dfrac{2174}{3040}\,\day$(定気雨水)$+418$$ +5.8160/\day$$-0.1318/\day^2$
$ 58\dfrac{2174}{3040}\,\day$(定気雨水)$ 73\dfrac{2250}{3040}\,\day$(定気啓蟄)$+491$$ +3.7987/\day$$-0.1240/\day^2$
$ 73\dfrac{2250}{3040}\,\day$(定気啓蟄)$ 88\dfrac{2700}{3040}\,\day$(定気春分)$+535$$ +1.9265/\day$$-0.1240/\day^2$
$ 88\dfrac{2700}{3040}\,\day$(定気春分)$104\dfrac{ 539}{3040}\,\day$(定気清明)$+551$$ -0.2048/\day$$-0.1178/\day^2$
$104\dfrac{ 539}{3040}\,\day$(定気清明)$119\dfrac{1791}{3040}\,\day$(定気穀雨)$+535$$ -1.9968/\day$$-0.1190/\day^2$
$119\dfrac{1791}{3040}\,\day$(定気穀雨)$135\dfrac{ 391}{3040}\,\day$(定気立夏)$+491$$ -3.7956/\day$$-0.1240/\day^2$
$135\dfrac{ 391}{3040}\,\day$(定気立夏)$150\dfrac{2446}{3040}\,\day$(定気小満)$+418$$ -5.6626/\day$$-0.1324/\day^2$
$150\dfrac{2446}{3040}\,\day$(定気小満)$166\dfrac{1915}{3040}\,\day$(定気芒種)$+314$$ -7.6555/\day$$-0.1436/\day^2$
$166\dfrac{1915}{3040}\,\day$(定気芒種)$182\dfrac{1892}{3040}\,\day$(定気夏至)$+176$$ -9.9405/\day$$-0.1436/\day^2$
$182\dfrac{1892}{3040}\,\day$(定気夏至)$198\dfrac{1869}{3040}\,\day$(定気小暑)$ 0$$-12.0819/\day$$+0.1436/\day^2$
$198\dfrac{1869}{3040}\,\day$(定気小暑)$214\dfrac{1338}{3040}\,\day$(定気大暑)$-176$$ -9.7018/\day$$+0.1324/\day^2$
$214\dfrac{1338}{3040}\,\day$(定気大暑)$230\dfrac{ 353}{3040}\,\day$(定気立秋)$-314$$ -7.5450/\day$$+0.1240/\day^2$
$230\dfrac{ 353}{3040}\,\day$(定気立秋)$245\dfrac{1993}{3040}\,\day$(定気処暑)$-418$$ -5.5634/\day$$+0.1190/\day^2$
$245\dfrac{1993}{3040}\,\day$(定気処暑)$261\dfrac{ 205}{3040}\,\day$(定気白露)$-491$$ -3.7038/\day$$+0.1178/\day^2$
$261\dfrac{ 205}{3040}\,\day$(定気白露)$276\dfrac{1084}{3040}\,\day$(定気秋分)$-535$$ -1.8954/\day$$+0.1178/\day^2$
$276\dfrac{1084}{3040}\,\day$(定気秋分)$291\dfrac{1534}{3040}\,\day$(定気寒露)$-551$$ +0.1783/\day$$+0.1240/\day^2$
$291\dfrac{1534}{3040}\,\day$(定気寒露)$306\dfrac{1610}{3040}\,\day$(定気霜降)$-535$$ +2.0042/\day$$+0.1318/\day^2$
$306\dfrac{1610}{3040}\,\day$(定気霜降)$321\dfrac{1299}{3040}\,\day$(定気立冬)$-491$$ +3.8950/\day$$+0.1446/\day^2$
$321\dfrac{1299}{3040}\,\day$(定気立冬)$336\dfrac{ 573}{3040}\,\day$(定気小雪)$-418$$ +5.9214/\day$$+0.1634/\day^2$
$336\dfrac{ 573}{3040}\,\day$(定気小雪)$350\dfrac{2432}{3040}\,\day$(定気大雪)$-314$$ +8.1610/\day$$+0.1886/\day^2$
$350\dfrac{2432}{3040}\,\day$(定気大雪)$365\dfrac{ 743}{3040}\,\day$(定気冬至)$-176$$+10.9010/\day$$+0.1886/\day^2$

日躔による補正値も、月離の場合と同じく $\dfrac{1}{3040}\,\day$(大衍暦における1分)を単位として丸めて計算すればよい。

五紀暦

月離による補正値

儀鳳暦の項に併せて掲載。

日躔による補正値

大衍暦の項のものと同一。

宣明暦

月離による補正値

宣明暦における月離による補正値 $\DLunar(t)$ は、近点月 $\PLunar = 27\dfrac{4658.19}{8400}\,\day$を周期とする周期関数であり、その具体的な形は下図の折れ線で与えられる。
※基準点を近地点とするか遠地点とするかは半周期毎に切り替わる。

宣明暦における月離による補正値
宣明暦における月離による補正値 $\DLunar(t)$

折れ線の頂点の値は下表の通りである。

宣明暦における月離による補正値
近地点/遠地点からの経過日数
$t$
補正値(宣明暦)
$\DLunar(t)$
近地点 $0\,\day$ $0\,\day$
$1\,\day$ $+\dfrac{830}{8400}\,\day$
$2\,\day$ $+\dfrac{1556}{8400}\,\day$
$3\,\day$ $+\dfrac{2154}{8400}\,\day$
$4\,\day$ $+\dfrac{2618}{8400}\,\day$
$5\,\day$ $+\dfrac{2947}{8400}\,\day$
$6\,\day$ $+\dfrac{3142}{8400}\,\day$
$6\dfrac{7465}{8400}\,\day$ $+\dfrac{3195}{8400}\,\day$
$7\,\day$ $+\dfrac{3188}{8400}\,\day$
$8\,\day$ $+\dfrac{3106}{8400}\,\day$
$9\,\day$ $+\dfrac{2881}{8400}\,\day$
$10\,\day$ $+\dfrac{2515}{8400}\,\day$
$11\,\day$ $+\dfrac{2014}{8400}\,\day$
$12\,\day$ $+\dfrac{1386}{8400}\,\day$
$13\,\day$ $+\dfrac{646}{8400}\,\day$
$13\dfrac{6529}{8400}\,\day$ $0\,\day$
遠地点 $0\,\day$ $0\,\day$
$1\,\day$ $-\dfrac{830}{8400}\,\day$
$2\,\day$ $-\dfrac{1556}{8400}\,\day$
$3\,\day$ $-\dfrac{2162}{8400}\,\day$
$4\,\day$ $-\dfrac{2633}{8400}\,\day$
$5\,\day$ $-\dfrac{2970}{8400}\,\day$
$6\,\day$ $-\dfrac{3172}{8400}\,\day$
$6\dfrac{7465}{8400}\,\day$ $-\dfrac{3225}{8400}\,\day$
$7\,\day$ $-\dfrac{3218}{8400}\,\day$
$8\,\day$ $-\dfrac{3136}{8400}\,\day$
$9\,\day$ $-\dfrac{2912}{8400}\,\day$
$10\,\day$ $-\dfrac{2546}{8400}\,\day$
$11\,\day$ $-\dfrac{2037}{8400}\,\day$
$12\,\day$ $-\dfrac{1394}{8400}\,\day$
$13\,\day$ $-\dfrac{646}{8400}\,\day$
近地点 $13\dfrac{6529}{8400}\,\day$ $0\,\day$

近地点/遠地点の通過時刻と、補正値が極大/極小になる時刻を例外として、基本的には1日毎に頂点が打たれている。
実際の計算に当たっては、$\dfrac{1}{8400}\,\day$(宣明暦における1分)を単位として丸めて計算すればよい。

日躔による補正値

宣明暦における日躔による補正値 $\DSolar(t)$ は、近点年 $\PSolar = 365\dfrac{2055}{8400}\,\day$を周期とする周期関数であり、その具体的な形は下図に示す通り24個の2次関数(放物線)を繋ぎ合わせた形である。

宣明暦における日躔による補正値
宣明暦における日躔による補正値 $\DSolar(t)$

2次関数は、各区間 $[ \tau_s, \tau_e )$ 毎に次式の形で与えられる。

\[ \DSolar(t) = \bigg[ a + b(t - \tau_s) + \dfrac{1}{2}c(t - \tau_s)(t - \tau_s - 1\,\day) \bigg] \cdot \dfrac{1}{8400}\,\day \]

区間の端点 $\tau_s, \tau_e$ 及び係数 $a, b, c$ の値は下表の通りである。

宣明暦における日躔による補正値
区間の始点 $\tau_s$
(近日点からの経過日数)
区間の終点 $\tau_e$
(近日点からの経過日数)
初日朓朒数
$a$
初日損益
$b$
損益の変化率
$c$
$ 0 \,\day$(定気冬至)$ 14\dfrac{4236}{8400}\,\day$(定気小寒)$ 0$$+33.4511/\day$$-0.3695/\day^2$
$ 14\dfrac{4236}{8400}\,\day$(定気小寒)$ 29\dfrac{1072}{8400}\,\day$(定気大寒)$ +449$$+28.0389/\day$$-0.3606/\day^2$
$ 29\dfrac{1072}{8400}\,\day$(定気大寒)$ 43\dfrac{7307}{8400}\,\day$(定気立春)$ +823$$+22.6998/\day$$-0.3519/\day^2$
$ 43\dfrac{7307}{8400}\,\day$(定気立春)$ 58\dfrac{6143}{8400}\,\day$(定気雨水)$+1122$$+17.8923/\day$$-0.4068/\day^2$
$ 58\dfrac{6143}{8400}\,\day$(定気雨水)$ 73\dfrac{6179}{8400}\,\day$(定気啓蟄)$+1346$$+11.7966/\day$$-0.3998/\day^2$
$ 73\dfrac{6179}{8400}\,\day$(定気啓蟄)$ 88\dfrac{7414}{8400}\,\day$(定気春分)$+1481$$ +5.7986/\day$$-0.3998/\day^2$
$ 88\dfrac{7414}{8400}\,\day$(定気春分)$104\dfrac{1450}{8400}\,\day$(定気清明)$+1526$$ -0.2433/\day$$-0.3779/\day^2$
$104\dfrac{1450}{8400}\,\day$(定気清明)$119\dfrac{5085}{8400}\,\day$(定気穀雨)$+1481$$ -6.1254/\day$$-0.3634/\day^2$
$119\dfrac{5085}{8400}\,\day$(定気穀雨)$135\dfrac{1521}{8400}\,\day$(定気立夏)$+1346$$-12.2048/\day$$-0.2987/\day^2$
$135\dfrac{1521}{8400}\,\day$(定気立夏)$150\dfrac{7357}{8400}\,\day$(定気小満)$+1122$$-16.9060/\day$$-0.2919/\day^2$
$150\dfrac{7357}{8400}\,\day$(定気小満)$166\dfrac{5792}{8400}\,\day$(定気芒種)$ +823$$-21.5362/\day$$-0.2854/\day^2$
$166\dfrac{5792}{8400}\,\day$(定気芒種)$182\dfrac{5228}{8400}\,\day$(定気夏至)$ +449$$-26.0498/\day$$-0.2854/\day^2$
$182\dfrac{5228}{8400}\,\day$(定気夏至)$198\dfrac{4664}{8400}\,\day$(定気小暑)$ 0$$-30.3119/\day$$+0.2854/\day^2$
$198\dfrac{4664}{8400}\,\day$(定気小暑)$214\dfrac{3099}{8400}\,\day$(定気大暑)$ -449$$-25.8126/\day$$+0.2919/\day^2$
$214\dfrac{3099}{8400}\,\day$(定気大暑)$230\dfrac{ 535}{8400}\,\day$(定気立秋)$ -823$$-21.2454/\day$$+0.2987/\day^2$
$230\dfrac{ 535}{8400}\,\day$(定気立秋)$245\dfrac{5370}{8400}\,\day$(定気処暑)$-1122$$-17.0296/\day$$+0.3634/\day^2$
$245\dfrac{5370}{8400}\,\day$(定気処暑)$261\dfrac{ 606}{8400}\,\day$(定気白露)$-1346$$-11.4744/\day$$+0.3779/\day^2$
$261\dfrac{ 606}{8400}\,\day$(定気白露)$276\dfrac{3042}{8400}\,\day$(定気秋分)$-1481$$ -5.6429/\day$$+0.3779/\day^2$
$276\dfrac{3042}{8400}\,\day$(定気秋分)$291\dfrac{4277}{8400}\,\day$(定気寒露)$-1526$$ +0.1432/\day$$+0.3998/\day^2$
$291\dfrac{4277}{8400}\,\day$(定気寒露)$306\dfrac{4313}{8400}\,\day$(定気霜降)$-1481$$ +6.1488/\day$$+0.4068/\day^2$
$306\dfrac{4313}{8400}\,\day$(定気霜降)$321\dfrac{3149}{8400}\,\day$(定気立冬)$-1346$$+12.6336/\day$$+0.3519/\day^2$
$321\dfrac{3149}{8400}\,\day$(定気立冬)$336\dfrac{ 984}{8400}\,\day$(定気小雪)$-1122$$+17.8043/\day$$+0.3606/\day^2$
$336\dfrac{ 984}{8400}\,\day$(定気小雪)$350\dfrac{6220}{8400}\,\day$(定気大雪)$ -823$$+23.0590/\day$$+0.3695/\day^2$
$350\dfrac{6220}{8400}\,\day$(定気大雪)$365\dfrac{2055}{8400}\,\day$(定気冬至)$ -449$$+28.4618/\day$$+0.3695/\day^2$

日躔による補正値も、月離の場合と同じく $\dfrac{1}{8400}\,\day$(宣明暦における1分)を単位として丸めて計算すればよい。

貞享暦・宝暦暦・修正宝暦暦

平均黄経に対する月離の進み

貞享暦・宝暦暦・修正宝暦暦における平均黄経に対する月離の進み $\dLunar(t)$ は、近点月 $\PLunar = 27.5546\,\day$を周期とする周期関数であり、その具体的な形は下図に示す通り4個の3次関数を繋ぎ合わせた形である。

平均黄経に対する月離の進み
平均黄経に対する月離の進み $\dLunar(t)$

各区間における3次関数の係数は下表の通りである。

平均黄経に対する月離の進み
区間の始点 $\tau_s$
(近地点からの経過日数)
区間の終点 $\tau_e$
(近地点からの経過日数)
平均黄経に対する月離の進み
$\dLunar(t)$
近地点 $ 0\,\day$$ 6.511958\,\day$$ 1.324 u - 0.052u^2 - 0.005u^3$ 日度
但し $u := t/\day$
$ 6.511958\,\day$遠地点 $\dfrac{\PLunar}{2} = 13.7773\,\day$$ 1.1731u - 0.037u^2 - 0.004u^3$ 日度
但し $u := (\tau_e - t)/\day$
遠地点 $\dfrac{\PLunar}{2} = 13.7773\,\day$$21.042642\,\day$$-1.1731u + 0.037u^2 + 0.004u^3$ 日度
但し $u := (t - \tau_s)/\day$
$21.042642\,\day$近地点 $\PLunar = 27.5546\,\day$$-1.324 u + 0.052u^2 + 0.005u^3$ 日度
但し $u := (\tau_e - t)/\day$

平均黄経に対する日躔の進み

貞享暦・宝暦暦・修正宝暦暦における平均黄経に対する日躔の進み $\dSolar(t)$ は、近点年 $\PSolar$を周期とする周期関数であり、その具体的な形は下図に示す通り4個の3次関数を繋ぎ合わせた形である。
(貞享暦では $\PSolar = 365.256696\,\day$、宝暦暦では $\PSolar = 365.256556\,\day$、修正宝暦暦では $\PSolar = 365.256626\,\day$)

平均黄経に対する日躔の進み
平均黄経に対する日躔の進み $\dLunar(t)$

各区間における3次関数の係数は下表の通りである。
(暦法に依り $\PSolar$ の値は異なるが、これに因る $\dLunar(t)$ の値への影響は $10^{-6}\,\day$程度である)

平均黄経に対する日躔の進み
区間の始点 $\tau_s$
(近日点からの経過日数)
区間の終点 $\tau_e$
(近日点からの経過日数)
平均黄経に対する日躔の進み
$\dSolar(t)$
近日点 $ 0 \,\day$$ 89.253920\,\day$$ 0.043600u - 0.0002000u^2 - 0.00000034u^3$ 日度
但し $u := t/\day$
$ 89.253920\,\day$遠日点 $\dfrac{\PSolar}{2}$$ 0.041198u - 0.0001764u^2 - 0.00000031u^3$ 日度
但し $u := (\tau_e - t)/\day$
遠日点 $\dfrac{\PSolar}{2}$$\PSolar - 89.253920\,\day$$-0.041198u + 0.0001764u^2 + 0.00000031u^3$ 日度
但し $u := (t - \tau_s)/\day$
$\PSolar - 89.253920\,\day$近日点 $\PSolar $$-0.043600u + 0.0002000u^2 + 0.00000034u^3$ 日度
但し $u := (\tau_e - t)/\day$

寛政暦・天保暦

寛政暦・天保暦では、平均黄経を介さずに日時から直接黄経を求める。
寛政暦・天保暦の暦法書に規定された黄経の計算は非常に複雑なのであるが、おまつのブログ(おまつさん)というブログにて詳説されているので、正確な計算方法はそちらを参照されたい。

実のところ、寛政暦・天保暦はそれまでの暦法と比べて格段に精度が高いので、時刻系の違い(平均太陽時ではなく真太陽時を用いること)にさえ注意すれば、次節(明治改暦以降の旧暦)に示す現代天文学による略算式を用いて計算しても然程大きな誤差は生じない。
現代天文学による略算式を用いて計算した場合に朔や節気の日付が史実の寛政暦・天保暦と異なるのは、下表に示す6つのみである。

略算式と寛政暦・天保暦との差異
差異箇所 次節の略算式による日付
(グレゴリオ暦)
史実の暦における日付
(グレゴリオ暦)
寛政暦1802年7月朔(第22283朔)
1819年5月朔(第22492朔)
1824年5月朔(第22553朔)
1828年10月朔(第22608朔)
天保暦1848年立春(第44331節気)
1852年清明(第44431節気)

計算を簡便に済ませたければ、現代天文学による略算式を用いて計算した上で、上表の朔・節気の日付のみ補正すれば良いだろう。
(節気の日付は補正せずとも旧暦の日付の決定には影響しない)

明治改暦以降の旧暦

明治改暦以降の旧暦では、何らかの暦法書に規定された計算式ではなく、実際の月や太陽の黄経が作暦の基準である。

実際に月や太陽の黄経を求めるに当たっては、何らかの略算式を用いることになる。
よく知られた略算式としては以下のようなものがある。

本稿では、前述のおまつさんの作成された水路部式の改訂版という略算式を紹介する。
「水路部式の改訂版」とは言うが、VSOP87 及び ELP2000-82B を基にした高精度なものである。

太陽の黄経

おまつさんの略算式によると、太陽の黄経(真黄経)$\lSolar(t)$ は次式で計算できる。
ここで、次式右辺の $u$ は、時刻 $t$ に対応するユリウス年数 $\JY(t)$ である。

\[ \begin{align*} \lSolar(t) &= \dummypl 0.0000025\degree u \sin(240.97024\degree + 719.9874571\degree u) \\ &\dummyeq + 0.0001181\degree u \sin(243.44584\degree + 359.9937286\degree u) \\ &\dummyeq + 0.00006\degree \sin(129.01901\degree +    8.9049329\degree u) \\ &\dummyeq + 0.00006\degree \sin( 85.80080\degree +  122.9662205\degree u) \\ &\dummyeq - 0.00006\degree \sin(109.92923\degree +   38.6827252\degree u) \\ &\dummyeq + 0.00006\degree \sin(334.46952\degree +    0.4075762\degree u) \\ &\dummyeq + 0.00006\degree \sin(145.88880\degree +  900.7377128\degree u) \\ &\dummyeq + 0.00006\degree \sin(126.43578\degree +  268.9465583\degree u) \\ &\dummyeq - 0.00006\degree \sin( 76.40919\degree + 9625.3576239\degree u) \\ &\dummyeq + 0.00007\degree \sin(126.98151\degree +    0.0561683\degree u) \\ &\dummyeq + 0.00007\degree \sin(152.05282\degree +    1.1903425\degree u) \\ &\dummyeq + 0.00008\degree \sin(285.44371\degree +  168.5907355\degree u) \\ &\dummyeq + 0.00009\degree \sin(137.73067\degree +   12.2211379\degree u) \\ &\dummyeq + 0.00012\degree \sin(230.80908\degree +  347.7725906\degree u) \\ &\dummyeq + 0.00012\degree \sin(197.10816\degree +  319.3175611\degree u) \\ &\dummyeq + 0.00012\degree \sin(  5.38791\degree +  145.7784780\degree u) \\ &\dummyeq + 0.00014\degree \sin(109.75127\degree +  314.3692135\degree u) \\ &\dummyeq + 0.00016\degree \sin(108.03170\degree +  628.9402869\degree u) \\ &\dummyeq + 0.00016\degree \sin(198.78622\degree +   45.6245150\degree u) \\ &\dummyeq + 0.00018\degree \sin( 65.11345\degree +  675.5532846\degree u) \\ &\dummyeq + 0.00020\degree \sin(257.27739\degree +    0.0038566\degree u) \\ &\dummyeq + 0.00028\degree \sin(209.06815\degree -   44.4341725\degree u) \\ &\dummyeq + 0.00029\degree \sin(352.56522\degree + 1079.9811857\degree u) \\ &\dummyeq - 0.00037\degree \sin(200.99619\degree +  720.0153950\degree u) \\ &\dummyeq + 0.00043\degree \sin(235.14760\degree +  315.5595560\degree u) \\ &\dummyeq + 0.00045\degree \sin(157.53705\degree +  299.2956151\degree u) \\ &\dummyeq + 0.00049\degree \sin(248.97860\degree -   22.8122575\degree u) \\ &\dummyeq + 0.00052\degree \sin(332.82704\degree -    1.5067827\degree u) \\ &\dummyeq + 0.00057\degree \sin( 29.80522\degree +  337.1814711\degree u) \\ &\dummyeq + 0.00069\degree \sin(153.57565\degree +   90.3751278\degree u) \\ &\dummyeq + 0.00073\degree \sin(333.28344\degree -   30.3490567\degree u) \\ &\dummyeq + 0.00076\degree \sin(132.52960\degree +  659.2893436\degree u) \\ &\dummyeq + 0.00134\degree \sin( 81.51535\degree +  225.1844282\degree u) \\ &\dummyeq + 0.00153\degree \sin(343.13042\degree +  450.3688564\degree u) \\ &\dummyeq + 0.00180\degree \sin(242.22020\degree - 4452.6711152\degree u) \\ &\dummyeq + 0.00196\degree \sin(287.91793\degree -    0.2018598\degree u) \\ &\dummyeq + 0.00200\degree \sin(247.22117\degree +  329.6446718\degree u) \\ &\dummyeq - 0.00478\degree \sin(125.03373\degree -   19.3413626\degree u) \\ &\dummyeq + 0.01999\degree \sin(355.04476\degree +  719.9874571\degree u) \\ &\dummyeq + 1.91463\degree \sin(357.52586\degree +  359.9937286\degree u) \\ &\dummyeq + 0.000000030\degree u^2 + 360.0076974\degree u + 719650.46075\degree \end{align*} \]

上式では、加減算する項の順番を原文とは逆順にしている。
これは、コンピュータの浮動小数点数演算で計算する際の情報落ち誤差[注1]を軽減するためである。
(より正確には、各項の値をそれぞれ求めた上で、絶対値の小さな項から順に加減算するべきである)

また、原文では $280.46075\degree$ であった定数項を上式では $719650.46075\degree$ としている。
これは、第0節気(西暦0年冬至)において $\lSolar(t) = 0\degree$ となるようにするための補正項である。
これにより、各整数 $m$ に対し、次が成り立つ。

$\lSolar(t) = 15m\degree$ なる時刻 $t$ が第 $m$ 節気

太陽の平均黄経

太陽の平均黄経 $\lbarSolar(t)$ は、前項の真黄経の式のうち周期項を除いたものである。

\[ \begin{align*} \lbarSolar(t) &= 0.000000030\degree u^2 + 360.0076974\degree u + 719650.46075\degree \end{align*} \]

月の黄経

おまつさんの略算式によると、月の黄経 $\lLunar(t)$ は次式で計算できる。
ここで、次式右辺の $u$ は、時刻 $t$ に対応するユリウス年数 $\JY(t)$ である。

\[ \begin{align*} \lLunar(t) &= \dummypl 0.0000047\degree u \sin(357.52909\degree + 359.9905029\degree u) \\ &\dummyeq + 0.00029\degree \sin(280.58970\degree + 23221.3082593\degree u) \\ &\dummyeq + 0.00030\degree \sin( 19.58410\degree -    40.6729418\degree u) \\ &\dummyeq + 0.00032\degree \sin( 53.08650\degree +  4792.6428976\degree u) \\ &\dummyeq + 0.00032\degree \sin(310.02141\degree +  5491.9696821\degree u) \\ &\dummyeq + 0.00033\degree \sin(148.09739\degree + 18089.3290801\degree u) \\ &\dummyeq + 0.00033\degree \sin(113.48956\degree -  3494.7184317\degree u) \\ &\dummyeq + 0.00034\degree \sin( 96.37637\degree -  3814.0359929\degree u) \\ &\dummyeq + 0.00035\degree \sin( 70.34233\degree +  9584.6502944\degree u) \\ &\dummyeq + 0.00037\degree \sin( 21.00755\degree +   720.0153950\degree u) \\ &\dummyeq + 0.00038\degree \sin(263.38263\degree -   120.0629979\degree u) \\ &\dummyeq + 0.00040\degree \sin(226.68534\degree -   398.7076173\degree u) \\ &\dummyeq + 0.00049\degree \sin(332.07641\degree +  5091.3062375\degree u) \\ &\dummyeq + 0.00052\degree \sin( 66.12900\degree - 17450.6939578\degree u) \\ &\dummyeq + 0.00054\degree \sin(179.85287\degree + 19087.9547053\degree u) \\ &\dummyeq + 0.00055\degree \sin(246.36331\degree + 22582.6731370\degree u) \\ &\dummyeq + 0.00060\degree \sin(128.37285\degree +  1118.6886231\degree u) \\ &\dummyeq + 0.00069\degree \sin(216.69729\degree +   278.6446194\degree u) \\ &\dummyeq + 0.00070\degree \sin(264.20498\degree -  4853.3345598\degree u) \\ &\dummyeq + 0.00071\degree \sin(139.90503\degree +  4052.0076705\degree u) \\ &\dummyeq + 0.00076\degree \sin(336.05544\degree -  7906.7166051\degree u) \\ &\dummyeq + 0.00081\degree \sin(188.19236\degree + 14037.3214096\degree u) \\ &\dummyeq + 0.00089\degree \sin(321.41315\degree -  8586.0246692\degree u) \\ &\dummyeq + 0.00111\degree \sin(276.47024\degree + 19208.0177032\degree u) \\ &\dummyeq + 0.00122\degree \sin(201.09222\degree - 12678.7052814\degree u) \\ &\dummyeq + 0.00159\degree \sin(242.24385\degree + 18569.3825809\degree u) \\ &\dummyeq + 0.00177\degree \sin(  4.11946\degree +  4013.2905561\degree u) \\ &\dummyeq + 0.00196\degree \sin(125.04550\degree -    19.3413618\degree u) \\ &\dummyeq + 0.00205\degree \sin( 74.32136\degree -  3413.3725482\degree u) \\ &\dummyeq + 0.00207\degree \sin(175.05819\degree +   719.9810058\degree u) \\ &\dummyeq + 0.00212\degree \sin( 87.45553\degree +  9903.9678555\degree u) \\ &\dummyeq + 0.00224\degree \sin(299.35814\degree -  8185.3612245\degree u) \\ &\dummyeq + 0.00235\degree \sin(252.81324\degree +  9224.6597915\degree u) \\ &\dummyeq + 0.00239\degree \sin(211.75548\degree +   998.6256252\degree u) \\ &\dummyeq + 0.00260\degree \sin( 72.71937\degree - 13797.3939046\degree u) \\ &\dummyeq + 0.00269\degree \sin(272.39734\degree +  9183.9868497\degree u) \\ &\dummyeq + 0.00367\degree \sin(349.18961\degree +  5410.6237986\degree u) \\ &\dummyeq + 0.00386\degree \sin(111.40009\degree + 17810.6844607\degree u) \\ &\dummyeq + 0.00396\degree \sin( 60.24759\degree -     1.3184887\degree u) \\ &\dummyeq + 0.00399\degree \sin(145.62648\degree + 18449.3195830\degree u) \\ &\dummyeq + 0.00404\degree \sin( 13.13417\degree + 13317.3404037\degree u) \\ &\dummyeq + 0.00478\degree \sin(305.03343\degree -    19.3413626\degree u) \\ &\dummyeq + 0.00499\degree \sin(295.37912\degree +  4812.6616181\degree u) \\ &\dummyeq + 0.00516\degree \sin(197.11319\degree +   319.3175611\degree u) \\ &\dummyeq + 0.00677\degree \sin( 53.22914\degree +  9265.3327332\degree u) \\ &\dummyeq + 0.00789\degree \sin(261.73408\degree -  4493.3440569\degree u) \\ &\dummyeq + 0.00855\degree \sin(338.52634\degree -  8266.7071080\degree u) \\ &\dummyeq + 0.01003\degree \sin( 44.88965\degree + 14315.9660289\degree u) \\ &\dummyeq + 0.01067\degree \sin(203.56313\degree - 13038.6957844\degree u) \\ &\dummyeq + 0.01098\degree \sin(308.41941\degree -  4892.0516742\degree u) \\ &\dummyeq + 0.01253\degree \sin(141.50702\degree + 14436.0290269\degree u) \\ &\dummyeq + 0.01533\degree \sin(130.84376\degree +   758.6981202\degree u) \\ &\dummyeq + 0.03038\degree \sin(312.49231\degree +  5131.9791792\degree u) \\ &\dummyeq + 0.03472\degree \sin(117.85002\degree +  4452.6711152\degree u) \\ &\dummyeq + 0.04092\degree \sin(137.43412\degree +  4411.9981734\degree u) \\ &\dummyeq + 0.04576\degree \sin(301.82905\degree -  8545.3517274\degree u) \\ &\dummyeq + 0.05332\degree \sin( 10.66326\degree + 13677.3309066\degree u) \\ &\dummyeq + 0.05707\degree \sin( 76.79227\degree -  3773.3630511\degree u) \\ &\dummyeq + 0.05879\degree \sin(214.22639\degree +   638.6351223\degree u) \\ &\dummyeq + 0.11433\degree \sin(  6.54381\degree +  9664.0403505\degree u) \\ &\dummyeq + 0.18512\degree \sin(177.52909\degree +   359.9905029\degree u) \\ &\dummyeq + 0.21362\degree \sin(269.92643\degree +  9543.9773526\degree u) \\ &\dummyeq + 0.65831\degree \sin(235.70005\degree +  8905.3422303\degree u) \\ &\dummyeq + 1.27401\degree \sin( 79.26317\degree -  4133.3535540\degree u) \\ &\dummyeq + 6.28877\degree \sin(134.96312\degree +  4771.9886763\degree u) \\ &\dummyeq - 0.000000133\degree u^2 + 4812.6788118\degree u + 9620588.31645\degree \end{align*} \]

項の順番については太陽の黄経のときと同様である。

また、原文では $218.31645\degree$ であった定数項を上式では $9620588.31645\degree$ としている。
これは、第0朔(0年11月朔)において $\lLunar(t) = \lSolar(t)$ となるようにするための補正項である。
これにより、各整数 $m$ に対し、次が成り立つ。

$\lLunar(t) - \lSolar(t) = 360m\degree$ なる時刻 $t$ が第 $m$ 朔

出典・参考資料

脚注