2020年11月1日日曜日

天保暦の暦法 (8) 月離 (3) 黄道緯度

[江戸頒暦の研究 総目次へ] 

前回までのところで月の黄経の算出が完了した。今回は月の黄緯の算出について。

寛政暦での月の黄緯は、月の昇交点離角 F と、軌道傾斜角 i に基づいて、月の緯度 β は、sinβ=sinisinF を算出するだけだったが、天保暦では、黄緯についても種々の不等項を置いて算出している。


黄道緯度汎数

[新法暦書巻二 推太陰黄道緯度用数]
最大黄白大距五度二十九分三十零秒五十五微
最小黄白大距四度九十九分九十一秒六十七微
黄白大距中数五度一十四分六十一秒一十一微
黄白大距半較一十四分六十九秒四十四微(即最大二均)
[新法暦書巻二 推太陰黄道緯度法]
求黄道緯度汎数「以半径為一率、実月距正交之正弦為二率、黄白大距中数之正弦為三率、求得四率為正弦、検表得黄道緯度汎数。実月距正交初宮至五宮為加、六宮至十一宮為減」
半径を以って一率と為し、実月距正交の正弦、二率と為し、黄白大距中数の正弦、三率と為し、求めて得る四率、正弦と為し、表を検じ黄道緯度汎数を得。実月距正交、初宮より五宮に至るは加と為し、六宮より十一宮に至るは減と為す。
最大黄白大距=5°.293055(=黄白大距中数+黄白大距半較)最小黄白大距=4°.999167(=黄白大距中数黄白大距半較)黄白大距中数=5°.146111黄白大距半較=0°.146944黄道緯度汎数=sin1(sin(黄白大距中数)sin(実月距正交))

月の黄緯を βm、軌道傾斜角(黄白大距)を im、月の黄経 λm の昇交点黄経 Ωm からの離角を F とするとき、
sinβm=sinimsinF
である。これが、黄道緯度汎数であり、寛政暦の月黄緯の計算であればこれで終わりである(ただし、寛政暦では、昇交点真黄経・月の軌道傾斜角に不等を織り込むことにより、sin(F2D) に比例する「緯度の出差」ともいうべき不等が黄緯に自然に織り込まれて計算されてくるが、天保暦では、昇交点真黄経・月の軌道傾斜角にはそのような不等を織り込んでいないので、「緯度の出差」は別途計算することになる(緯度二均にあたる))。

なお、「黄白大距中数」は平均の月の軌道傾斜角であり、軌道傾斜角の不等(緯度二均)を織り込んだ場合の軌道傾斜角の値のレンジが「最小黄白大距」~「最大黄白大距」である。

太陰緯度一均

最大一均一十六秒六十七微
求一均「以半径為一率、三倍実月距正交之正弦為二率、最大一均為三率、求得四率為一均。三倍実月距正交初宮至五宮為減、六宮至十一宮為加」
半径を以って一率と為し、三倍実月距正交の正弦、二率と為し、最大一均、三率と為し、求めて得る四率、為一均。三倍実月距正交、初宮より五宮に至るは減と為し、六宮より十一宮に至るは加と為す。
一均=0°.001667sin(3×実月距正交)

黄道緯度汎数=sin1(sin(黄白大距中数)sin(実月距正交)) は、おおざっぱにいえば、+sinF に比例する項であり、その 3 倍成分が一均である。

が、フーリエ級数展開すると、
黄道緯度汎数=sin1(sin(5°.146111)sinF)+5°.144379sinF0°.001731sin3F+0°.000002sin5F
であり、緯度一均の 0°.001667sin3F は、黄道緯度汎数にすでに織り込まれている 0°.001731sin3F を二重に計算しているものなんじゃないかという疑いがなくもない。

太陰緯度二均

求真月距日「置実月距日、加減十二均及十三均、得真月距日。倍之(満十二宮去之)為倍真月距日」
実月距日を置き、十二均及び十三均を加減し、真月距日を得。これを倍し(満十二宮これを去く)倍真月距日と為す。
求二均引数「置倍真月距日、減実月距正交(不足減者、加十二宮減之)、得二均引数」
倍真月距日を置き、減実月距正交(減に足らざれば、十二宮を加へこれを減ず)、二均引数を得。
求二均「以半径為一率、二均引数之正弦為二率、黄白大距半較為三率、求得四率為二均。二均引数初宮至五宮為加、六宮至十一宮為減」
半径を以って一率と為し、二均引数の正弦、二率と為し、黄白大距半較、三率と為し、求めて得る四率、二均と為す。二均引数、初宮より五宮に至るは加と為し、六宮より十一宮に至るは減と為す。
真月距日=実月距日+太陰経度十二均+太陰経度十三均二均引数=2×真月距日実月距正交二均=+黄白大距半較sin(二均引数)

緯度における出差。緯度の計算での最大の不等項である。寛政暦では黄緯に自然に織り込まれて算出されていた。
二均=+0°.146944sin(2DF)=0°.146944sin(F2D)

太陰緯度三均

最大三均五秒五十六微
求三均引数「置実月距正交、減太陽引数(不足減者、加十二宮減之)、得三均引数」
実月距正交を置き、太陽引数を減じ(減に足らざれば、十二宮を加へこれを減ず)、三均引数を得。
求三均「以半径為一率、三均引数之正弦為二率、最大三均為三率、求得四率為三均。三均引数初宮至五宮為加、六宮至十一宮為減」
半径を以って一率と為し、三均引数の正弦、二率と為し、最大三均、三率と為し、求めて得る四率、三均と為す。三均引数、初宮より五宮に至るは加と為し、六宮より十一宮に至るは減と為す。
三均引数=実月距正交太陽引数三均=+0°.000556sin(三均引数)

緯度三均~六均は、F に、太陽の平均遠点角(太陽引数)l+180° や、月の平均遠点角(真引数) l+180° を加減したものの sin に比例する不等になっている。黄道緯度汎数に対する調整項と言える。
三均=+0°.000556sin(F(l+180°))=0°.000556sin(Fl)

太陰緯度四均

最大四均四十八秒三十三微
求真引数「置白道実行、減最高実行(不足減者、加十二宮減之)、得真引数。倍之(満十二宮去之)為倍真引数。三之(満十二宮去之)為三倍真引数」
白道実行を置き、最高実行を減じ(減に足らざれば、十二宮を加へこれを減ず)、真引数を得。これを倍し(満十二宮これを去く)倍真引数と為す。これを三し(満十二宮去之)三倍真引数と為す。
求四均引数「置実月距正交、減真引数(不足減者、加十二宮減之)、得四均引数」
実月距正交を置き、真引数を減じ(減に足らざれば、十二宮を加へこれを減ず)、四均引数を得。
求四均「以半径為一率、四均引数之正弦為二率、最大四均為三率、求得四率為四均。四均引数初宮至五宮為減、六宮至十一宮為加」
半径を以って一率と為し、四均引数の正弦、二率と為し、最大四均、三率と為し、求めて得る四率、四均と為す。四均引数、初宮より五宮に至るは減と為し、六宮より十一宮に至るは加と為す。
真引数=白道実行最高実行四均引数=実月距正交真引数四均=0°.004833sin(四均引数)

四均=0°.004833sin(F(l+180°))=+0°.004833sin(Fl)

太陰緯度五均

最大五均六十六秒九十四微
求五均引数「置四均引数、減真引数(不足減者、加十二宮減之)、得五均引数」
四均引数を置き、真引数を減じ(減に足らざれば、十二宮を加へこれを減ず)、五均引数を得。
求五均「以半径為一率、五均引数之正弦為二率、最大五均為三率、求得四率為五均。五均引数初宮至五宮為減、六宮至十一宮為加」
半径を以って一率と為し、五均引数の正弦、二率と為し、最大五均、三率と為し、求めて得る四率、五均と為す。五均引数、初宮より五宮に至るは減と為し、六宮より十一宮に至るは加と為す。
五均引数=四均引数真引数(=実月距正交2×真引数)五均=0°.006694sin(五均引数)

五均=0°.006694sin(F2(l+180°))=0°.006694sin(F2l)

太陰緯度六均

最大六均七秒五十微
求六均引数「置五均引数、減真引数(不足減者、加十二宮減之)、得六均引数」
五均引数を置き、真引数を減じ(減に足らざれば、十二宮を加へこれを減ず)、六均引数を得。
求六均「以半径為一率、六均引数之正弦為二率、最大六均為三率、求得四率為六均。六均引数初宮至五宮為加、六宮至十一宮為減」
半径を以って一率と為し、六均引数の正弦、二率と為し、最大六均、三率と為し、求めて得る四率、六均と為す。六均引数、初宮より五宮に至るは加と為し、六宮より十一宮に至るは減と為す。
六均引数=五均引数真引数(=実月距正交3×真引数)六均=+0°.000750sin(六均引数)

六均=+0°.000750sin(F3(l+180°))=0°.000750sin(F3l)

太陰緯度七均

最大七均二十三秒零六微
求七均引数「置二均引数、加太陽引数(満十二宮去之)、得七均引数」
二均引数を置き、太陽引数を加へ(満十二宮これを去く)、七均引数を得。
求七均「以半径為一率、七均引数之正弦為二率、最大七均為三率、求得四率為七均。七均引数初宮至五宮為減、六宮至十一宮為加」
半径を以って一率と為し、七均引数の正弦、二率と為し、最大七均、三率と為し、求めて得る四率、七均と為す。七均引数、初宮より五宮に至るは減と為し、六宮より十一宮に至るは加と為す。
七均引数=二均引数+太陽引数(=2×真月距日実月距正交+太陽引数)七均=0°.002306sin(七均引数)

緯度七均~十一均は、F2D に、太陽の平均遠点角(太陽引数)l+180° や、月の平均遠点角(真引数) l+180° を加減したものの sin に比例する不等になっている。緯度の出差(二均)に対する調整項と言える。
七均=0°.002306sin((2DF)+(l+180°))=0°.002306sin(F2Dl)

太陰緯度八均

最大八均一十零秒二十八微
求八均引数「置二均引数、減太陽引数(不足減者、加十二宮減之)、得八均引数」
二均引数を置き、太陽引数を減じ(減に足らざれば、十二宮を加へこれを減ず)、八均引数を得。
求八均「以半径為一率、八均引数之正弦為二率、最大八均為三率、求得四率為八均。八均引数初宮至五宮為減、六宮至十一宮為加」
半径を以って一率と為し、八均引数の正弦、二率と為し、最大八均、三率と為し、求めて得る四率、八均と為す。八均引数、初宮より五宮に至るは減と為し、六宮より十一宮に至るは加と為す。
八均引数=二均引数太陽引数(=2×真月距日実月距正交太陽引数)八均=0°.001028sin(八均引数)
八均=0°.001028sin((2DF)(l+180°))=0°.001028sin(F2D+l)

太陰緯度九均

最大九均六秒一十一微
求七均引数「置二均引数、加真引数(満十二宮去之)、得九均引数」
二均引数を置き、真引数を加へ(満十二宮これを去く)、九均引数を得。
求九均「以半径為一率、九均引数之正弦為二率、最大九均為三率、求得四率為九均。九均引数初宮至五宮為減、六宮至十一宮為加」
半径を以って一率と為し、九均引数の正弦、二率と為し、最大九均、三率と為し、求めて得る四率、九均と為す。九均引数、初宮より五宮に至るは減と為し、六宮より十一宮に至るは加と為す。
九均引数=二均引数+真引数(=2×真月距日実月距正交+真引数)九均=0°.000611sin(九均引数)
九均=0°.000611sin((2DF)+(l+180°))=0°.000611sin(F2Dl)

太陰緯度十均

最大十均四十一秒六十七微
求十均引数「置二均引数、減真引数(不足減者、加十二宮減之)、得十均引数」
二均引数を置き、真引数を減じ(減に足らざれば、十二宮を加へこれを減ず)、十均引数を得。
求十均「以半径為一率、十均引数之正弦為二率、最大十均為三率、求得四率為十均。十均引数初宮至五宮為加、六宮至十一宮為減」
半径を以って一率と為し、十均引数の正弦、二率と為し、最大十均、三率と為し、求めて得る四率、十均と為す。十均引数、初宮より五宮に至るは加と為し、六宮より十一宮に至るは減と為す。
十均引数=二均引数真引数(=2×真月距日実月距正交真引数)十均=+0°.004167sin(十均引数)

十均=+0°.004167sin((2DF)(l+180°))=+0°.004167sin(F2D+l)

太陰緯度十一均

最大十一均一十六秒六十七微
求十一均引数「置十均引数、減真引数(不足減者、加十二宮減之)、得十一均引数」
十均引数を置き、真引数を減じ(減に足らざれば、十二宮を加へこれを減ず)、十一均引数を得。
求十一均「以半径為一率、十一均引数之正弦為二率、最大十一均為三率、求得四率為十一均。十一均引数初宮至五宮為減、六宮至十一宮為加」
半径を以って一率と為し、十一均引数の正弦、二率と為し、最大十一均、三率と為し、求めて得る四率、十一均と為す。十一均引数、初宮より五宮に至るは減と為し、六宮より十一宮に至るは加と為す。
十一均引数=十均引数真引数(=2×真月距日実月距正交2×真引数)十一均=0°.001667sin(十一均引数)

十一均=0°.001667sin((2DF)2(l+180°))=+0°.001667sin(F2D+2l)

太陰黄道緯度

求太陰黄道緯度「自黄道緯度汎数及一均至十一均、加号相併為加均、減号相併為減均。両均数相減、得太陰黄道緯度。加数大為北、減数大為南」
黄道緯度汎数及び一均より十一均に至る、加号相併せ加均と為し、減号相併せ減均と為す。両均数、相減じ、太陰黄道緯度を得。加数大は北と為し、減数大は南と為す。
太陰黄道緯度=黄道緯度汎数+一均+二均+三均+四均+五均+六均+七均+八均+九均+十均+十一均

黄道緯度汎数に、一均~十一均を加減して、太陰黄道緯度を得る。

水路部式との比較

天保暦の月緯度算出と、水路部式とを比較してみる。

No. 水路部式
[水路部式解釈]
天保均数 天保振幅
1 +5°.1282sin(236°.231+4832°.0202t+B)
[+5.1282sin(F+)]
汎数 +5°.146111
2 +0°.2806sin(215°.147+9604°.0088t)
[+0.2806sin(F+l)]


3 +0°.2777sin(77°.316+60°.0316t)
[0.2777sin(Fl)]
四均 +0°.004833
4 +0°.1732sin(4°.5634073°.322t)
[0.1732sin(F2D)]
二均 0°.146944
5 +0°.0554sin(308°.98+8965°.374t)
[+0.0554sin(F+2Dl)]


6 +0°.0463sin(343°.48+698°.667t)
[0.0463sin(F2D+l)]
十均 +0°.004167
7 +0°.0326sin(287°.90+13737°.362t)
[+0.0326sin(F+2D)]


8 +0°.0172sin(194°.06+14375°.997t)
[+0.0172sin(F+2l)]


9 +0°.0033sin(25°.68845°.31t)
[0.0033sin(F2Dl)]
九均 0°.000611
10 +0°.0088sin(98°.44711°.96t)
[0.0088sin(F2l)]
五均 0°.006944
11 +0°.0082sin(1°.13713°.33t)
[0.0082sin(F2D+l)]
八均 0°.001028
12 +0°.0043sin(322°.4+5470°.66t)
[0.0043sin(F2D+2l)]
十一均 +0°.001667
13 +0°.0042sin(266°.8+18509°.35t)
[+0.0042sin(F+2D+l)]


14 +0°.0034sin(188°.04433°.31t)
[+0.0034sin(F2Dl)]
七均 0°.002306
15 +0°.0025sin(312°.5+8605°.38t)
[+0.0025sin(F+2Dll)]


16 +0°.0022sin(291°.4+13377°.37t)
[+0.0022sin(F+2Dl)]


17 +0°.0021sin(340°.0+1058°.66t)
[0.0021sin(F2D+l+l)]


18 +0°.0019sin(218°.6+9244°.02t)
[+0.0019sin(F+ll)]


19 +0°.0018sin(291°.88206°.68t)
[0.0018sin(F4D+l)]


20 +0°.0018sin(52°.8+5192°.01t)
[0.0018sin(F+l)]


21 +0°.0017sin(168°.7+14496°.06t)
[0.0017sin(3F)]
一均 0°.001667
22 +0°.0016sin(73°.8+420°.02t)
[0.0016sin(Fl+l)]


23 +0°.0015sin(262°.1+9284°.69t)
[0.0015sin(F+D)]


24 +0°.0015sin(31°.7+9964°.00t)
[0.0015sin(F+l+l)]


25 +0°.0014sin(260°.8299°.96t)
[+0.0014sin(Fll)]


26 +0°.0013sin(239°.7+4472°.03t)
[+0.0013sin(Fl)]
三均 0°.000556
27 +0°.0013sin(30°.4+379°.35t)
[+0.0013sin(FD)]


28 +0°.0012sin(304°.9+4812°.68t)
[0.0012sin(λm)]


29 +0°.0012sin(12°.44851°.36t)
[+0.0012sin(λm2F)]


30 +0°.0011sin(173°.0+19147°.99t)
[+0.0011sin(F+3l)]


31 +0°.0010sin(312°.912978°.66t)
[0.0010sin(F4D)]


32 +0°.0008sin(1°.+17870°.7t)
[+0.0008sin(F+4Dl)]


33 +0°.0008sin(190°.+9724°.1t)
[0.0008sin(3Fl)]


34 +0°.0007sin(22°.+13098°.7t)
[+0.0007sin(F+4D2l)]


35 +0°.0006sin(117°.+5590°.7t)
[0.0006sin(3F2D)]


36 +0°.0006sin(47°.13617°.3t)
[0.0006sin(F2D2l)]


37 +0°.0005sin(22°.8485°.3t)
[0.0005sin(F2Dl+l)]


38 +0°.0005sin(150°.+4193°.4t)
[0.0005sin(F+2D2l)]


39 +0°.0004sin(119°.9483°.9t)
[0.0004sin(F3l)]
六均 0°.000750
40 +0°.0004sin(246°.+23281°.3t)
[+0.0004sin(F+2D+2l)]


41 +0°.0004sin(301°.+10242°.6t)
[0.0004sin(F2D+3l)]


42 +0°.0004sin(126°.+9325°.4t)
[0.0004sin(F+2Dl+l)]


43 +0°.0004sin(104°.+14097°.4t)
[0.0004sin(F+2D+l)]


44 +0°.0003sin(340°.+22642°.7t)
[+0.0003sin(F+4D)]


45 +0°.0003sin(270°.+18149°.4t)
[+0.0003sin(F+2D+ll)]


46 +0°.0003sin(358°.3353°.3t)
[0.0003sin(F2D+2l)]


47 +0°.0003sin(148°.+19268°.0t)
[0.0003sin(3F+l)]


正直、あまり合致している感はない。

まあ単純比較はできないのだが。天保暦では、月黄経、太陽黄経は真黄経を用いているのに対し、水路部式では平均黄経を用いており、平均黄経と真黄経とで出る差違は別の不等として切り出している。汎数 sinF について、月黄経に真黄経を用いるとすれば、
汎数sin(F+2esinl)=sinFcos(2esinl)+cosFsin(2esinl)sinF+2ecosFsinl=sinF+esin(F+l)esin(Fl)
となり、No. 2, 3 の sin(F+l),sin(Fl) の項が出てくるわけである。よって天保暦ではこれは汎数に含まれて算出されるわけだから別には切り出されていない。

水路部式には F+2D ベースの不等が含まれているようだ。天保暦では F ベース、F2D ベースのものは含まれているが、F+2D ベースのものはない。

天保一均 0°.001667sin3F は、No. 21 0.0017sin(3F) に対応させてはいるが、これは、汎数、
黄道緯度汎数=sin1(sin(5°.146111)sinF)+5°.144379sinF0°.001731sin3F+0°.000002sin5F
に包含されているように思われるので、やはり二重計上なんじゃないかという気がする。


月の黄経・黄緯が得られたところで、次回は赤経・赤緯への変換について。

また、黄経・黄緯で「左右・上下」の座標が得られたわけだが、三次元的位置を確定するための残りの座標「前後」の位置、つまり、月と地球との間の距離を示す「地半径差」(地平視差 horizontal parallax)の算出について。

[江戸頒暦の研究 総目次へ]

[参考文献]

渋川 景祐; 足立 信頭「新法暦書」 国立公文書館デジタルアーカイブ蔵

渋川 景祐; 足立 信行「新法暦書続編」 国立公文書館デジタルアーカイブ蔵

長沢 工 (1981, 1985)「天体の位置計算 増補版」, 地人書館 ISBN-9784805202258

0 件のコメント:

コメントを投稿